A system of spins coupled to a bath is a traditional setup in open quantum systems. Through Heisenberg's equation, the spin dynamics can be modeled by a set of first-order differential equations. Interpreting the terms as colored noise and non-Markovian damping, one can write them as quantummechanical Heisenberg-Langevin (HL) equations. These are notoriously difficult to solve because of the high dimensionality of the Hilbert space. Classical generalized Langevin equations, involving non-Markovian damping and colored noise, are well understood and can be treated numerically with relative ease. Thus, a classical ansatz can be made by substituting quantum expectation values with classical functions. This allows the application of standard methods developed for classical stochastic dynamical systems to tackle spin dynamics. However, this approach is uncontrolled and should be benchmarked against known quantum dynamics. In this investigation, a Hamiltonian for spin dynamics is modified to obtain a setup analogous to the Weisskopf-Wigner (WW) theory of spontaneous emission, enabling a comparison of the results. This will be compared for T = 0 and with a slight adaptation in the high-temperature limit.
A system of spins coupled to a bath is a traditional setup in open quantum systems. Through Heisenberg's equation, the spin dynamics can be modeled by a set of first-order differential equations. Interpreting the terms as colored noise and non-Markovian damping, one can write them as quantummechanical Heisenberg-Langevin (HL) equations. These are notoriously difficult to solve because of the high dimensionality of the Hilbert space. Classical generalized Langevin equations, involving non-Markovian damping and colored noise, are well understood and can be treated numerically with relative ease. Thus, a classical ansatz can be made by substituting quantum expectation values with classical functions. This allows the application of standard methods developed for classical stochastic dynamical systems to tackle spin dynamics. However, this approach is uncontrolled and should be benchmarked against known quantum dynamics. In this investigation, a Hamiltonian for spin dynamics is modified to obtain a setup analogous to the Weisskopf-Wigner (WW) theory of spontaneous emission, enabling a comparison of the results. This will be compared for T = 0 and with a slight adaptation in the high-temperature limit.
Abstract We present a general numerical simulation method to solve non-Markovian dynamics of an open quantum system influenced by quantum Brownian motion. Based on the determined memory kernel function, this method enables the resolution of non-Markovian dynamics for a wide range of system Hamiltonians and spectral densities. The system dynamics are described by exact integro-differential operator equations without any common approximations and they are simulated in this work by definite-number equations with stochastic initial conditions. This approach ensures the applicability of mature numerical methods and maintains computational complexity that remains largely invariant, even when dealing with more complex models. The high accuracy of our simulation is evident from a comparison with the results obtained from corresponding exact master equations, underscoring the reliability and precision of our method.
Phase-space methods are powerful tools for the treatment of Markovian open quantum systems: they map the reduced dynamics of a system S, in interaction with an environment E, exactly onto Langevin equations for c-number stochastic variables, as opposed to Heisenberg-Langevin equations for operators. Langevin equations provide analytical insight in key regimes and excel at handling strong nonlinearities and couplings, where other methods often falter. Extending phase-space methods to non-Markovian dynamics, however, has remained a long-standing challenge. Here we address this gap by applying phase-space representations to the full S+E system; integrating out the environmental degrees of freedom then yields a general Langevin framework for S that incorporates both deterministic and stochastic contributions from E. Normally ordered representations, such as the Glauber-Sudarshan P representation and its positive variant due to Drummond and Gardiner, lead to Langevin equations in which (i) non-Markovian effects emerge exclusively in the deterministic terms, via a memory kernel, and (ii) noise contributions vanish when E is initially in the vacuum state. To demonstrate the power of this framework, we address the paradigmatic problem of delayed coherent feedback, in which the system is driven by its own past state, and study its impact on the laser linewidth: we recover the narrowing observed well above threshold and predict an enhanced narrowing just above it. Crucially, the number of stochastic variables scales linearly with the system size, making the framework suitable for problems ranging from a few degrees of freedom to genuinely many-body systems. This opens the way to the systematic study of non-Markovian driven-dissipative quantum systems using the same analytical and numerical tools that have long made phase-space methods so successful in the Markovian regime.
Previous years researchers began to simulate open quantum system, taking into account the interaction between system and the environment. One approach to deal with this problem is to use the density matrix within the Liouville-von-Neumann formalism or the Markovian variant the Lindblad equations. Another way is to use a stochastic approach where a random force is added to the system. The benefit of the stochastic approach is to solve the dynamics of the system with less time and memory than the density matrix approaches. In this project we want to develop a stochastic approach that can deal with the stochastic wave functions approach. We did this on a 2-level system and found that it works well when comparing to a density matrix approach. Next, we tested a quantum particle connect to a bath of harmonic oscillators using the stochastic approach. We found that a friction term is necessary and applied it. Like in the classical Langevin equations the friction constant and the random force fluctuations are related by the fluctuation-dissipation constant. We showed that with friction the dynamics decays to an ensemble with energy of Egs+kBT. However, we also found here are problems. The system seems to absorb energy indefinitely if the temperature is higher than the zero point energy or if the system is a Morse oscillator. Thus more research is required to make this method work.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.
The stochastic Liouville-von Neumann (SLN) equation describes the dynamics of an open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. The interaction with the environment is represented by complex colored noises which drive the system, and whose correlation functions are set by the properties of the environment. We present a number of schemes capable of generating colored noises of this kind that are built on a noise amplitude reduction procedure [Imai et al., Chem. Phys. 446, 134 (2015)CMPHC20301-010410.1016/j.chemphys.2014.11.014], including two analytically optimized schemes. In doing so, we pay close attention to the properties of the correlation functions in Fourier space, which we derive in full. For some schemes the method of Wiener filtering for deconvolutions leads to the realization that weakening causality in one of the noise correlation functions improves numerical convergence considerably, allowing us to introduce a well-controlled method for doing so. We compare the ability of these schemes, along with an alternative optimized scheme [Schmitz and Stockburger, Eur. Phys. J.: Spec. Top. 227, 1929 (2019)1951-635510.1140/epjst/e2018-800094-y], to reduce the growth in the mean and variance of the trace of the reduced density matrix, and their ability to extend the region in which the dynamics is stable and well converged for a range of temperatures. By numerically optimizing an additional noise scaling freedom, we identify the scheme which performs best for the parameters used, improving convergence by orders of magnitude and increasing the time accessible by simulation.
In traditional open quantum systems, the baths are usually traced out so that only the system information is left in the equations of motion. However, recent studies reveal that using only the system degrees of freedom can be insufficient. In this work, we develop a stochastic c-number Langevin equation method which can conveniently access the bath information. In our approach, the studied quantities are the expectation values of operators which can contain both system operators and bath operators. The dynamics of the operators of interest is formally divided into separate system and bath parts, with auxiliary stochastic fields. After solving the independent stochastic dynamics of the system part and the bath part, we can recombine them by taking the average over these stochastic fields to obtain the desired quantities. Several applications of the theory are highlighted, including the pure dephasing model, the spin-boson model, and an optically excited quantum dot coupled to a bath of phonons.