A system of N spin-1/2 particles interacting with a thermal reservoir is used as a pedagogical example for advanced undergraduate and graduate students. We introduce and illustrate some methods, approximations, and phenomena related to dissipation and nonlinearity in many-particle physics. We start our analysis from the dynamical Sudarshan-Lindblad quantum master equation for the density operator of a system S interacting with a thermal reservoir R. We derive the quantum version of the so-called Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) equations such that the master equation can be decomposed in a hierarchical set of N - 1 equations (N > 1). The hierarchy is broken by introducing the mean-field approximation and reducing the problem to a nonlinear single particle system. In this scenario, the Hamiltonian is nonlinear (i.e., it depends on the state of <S), although the superoperator responsible for the dissipation and decoherence of S remains unaffected. To provide a useful tool to students: (1) we discuss the physical approximations involved, (2) we derive the analytical solution to the mean values equations of motion resulting from the Hamiltonian, (3) we solve analytically the master equation in the stationary regime, (4) we obtain and discuss the solution of the nonlinear master equation, numerically, and finally, (5) we discuss the master equation beyond the mean-field approximation and show how to introduce higher order quantum correlations that have been previously neglected.
A system of N spin-1/2 particles interacting with a thermal reservoir is used as a pedagogical example for advanced undergraduate and graduate students. We introduce and illustrate some methods, approximations, and phenomena related to dissipation and nonlinearity in many-particle physics. We start our analysis from the dynamical Sudarshan-Lindblad quantum master equation for the density operator of a system S interacting with a thermal reservoir R. We derive the quantum version of the so-called Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) equations such that the master equation can be decomposed in a hierarchical set of N - 1 equations (N > 1). The hierarchy is broken by introducing the mean-field approximation and reducing the problem to a nonlinear single particle system. In this scenario, the Hamiltonian is nonlinear (i.e., it depends on the state of <S), although the superoperator responsible for the dissipation and decoherence of S remains unaffected. To provide a useful tool to students: (1) we discuss the physical approximations involved, (2) we derive the analytical solution to the mean values equations of motion resulting from the Hamiltonian, (3) we solve analytically the master equation in the stationary regime, (4) we obtain and discuss the solution of the nonlinear master equation, numerically, and finally, (5) we discuss the master equation beyond the mean-field approximation and show how to introduce higher order quantum correlations that have been previously neglected.
A system of N spin-1/2 particles interacting with a thermal reservoir is used as a pedagogical example for advanced undergraduate and graduate students. We introduce and illustrate some methods, approximations, and phenomena related to dissipation and nonlinearity in many-particle physics. We start our analysis from the dynamical Sudarshan-Lindblad quantum master equation for the density operator of a system S interacting with a thermal reservoir R. We derive the quantum version of the so-called Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) equations such that the master equation can be decomposed in a hierarchical set of N−1 equations (N>1). The hierarchy is broken by introducing the mean-field approximation and reducing the problem to a nonlinear single particle system. In this scenario, the Hamiltonian is nonlinear (i.e., it depends on the state of S), although the superoperator responsible for the dissipation and decoherence of S remains unaffected. To provide a useful tool to students: (1) we discuss the physical approximations involved, (2) we derive the analytical solution to the mean values equations of motion resulting from the Hamiltonian, (3) we solve analytically the master equation in the stationary regime, (4) we obtain and discuss the solution of the nonlinear master equation, numerically, and finally, (5) we discuss the master equation beyond the mean-field approximation and show how to introduce higher order quantum correlations that have been previously neglected.
A recently proposed Markov approach provides Lindblad-type scattering superoperators, which ensure the physical (positive-definite) character of the many-body density matrix. We apply the mean-field approximation to such a many-body equation, in the presence of one- and two-body scattering mechanisms, and we derive a closed equation of motion for the electronic single-particle density matrix, which turns out to be nonlinear as well as non-Lindblad. We prove that, in spite of its nonlinear and non-Lindblad structure, the mean-field approximation does preserve the positive-definite character of the single-particle density matrix, an essential prerequisite of any reliable kinetic treatment of semiconductor quantum devices. This result is in striking contrast with conventional (non-Lindblad) Markov approaches, where the single-particle mean-field equations can lead to positivity violations and thus to unphysical results. Furthermore, the proposed single-particle formulation is extended to the case of quantum systems with spatial open boundaries, providing a formal derivation of a recently proposed density-matrix treatment based on a Lindblad-like system-reservoir scattering superoperator.
Reservoir engineering has proven to be a practical approach to control open quantum systems, preserving quantum coherence by appropriately manipulating the reservoir and system-reservoir interactions. In this context, for systems comprised of different parts, it is common to describe the dynamics of a subsystem of interest by making an adiabatic elimination of the remaining components of the system. This procedure often leads to an effective master equation for the subsystem that is not in the well-known Gorini-Kossakowski-Lindblad-Sudarshan form (here called standard Lindblad form). Instead, it has a more general structure (here called generalized Lindblad form), which explicitly reveals the dissipative coupling between the various components of the subsystem. In this work, we present a set of dynamical equations for the first and second moments of the canonical variables for linear systems, bosonic and fermionic, described by master equations in a generalized Lindblad form. Our method is efficient and allows one to obtain analytical solutions for the steady state. Further, we include as a review some covariance matrix methods for which our results are particularly relevant, paying special attention to those related to the measurement of entanglement. Finally, we prove that the Duan criterion for entanglement is also applicable to fermionic systems.
An open quantum system is a quantum system that interacts with some environment whose degrees of freedom have been coarse grained away. This model describes non-equilibrium processes more general than scattering-matrix formulations. Furthermore, the microscopically-derived environment provides a model of noise, dissipation and decoherence far more general than Markovian (white noise) models. The latter are fully characterized by Lindblad equations and can be motivated phenomenologically. Non-Markovian processes consistently account for backreaction with the environment and can incorporate effects such as finite temperature and spatial correlations. We consider linear systems with bilinear coupling to the environment, or quantum Brownian motion, and nonlinear systems with weak coupling to the environment. For linear systems we provide exact solutions with analytical results for a variety of spectral densities. Furthermore, we point out an important mathematical subtlety which led to incorrect master-equation coefficients in earlier derivations, given nonlocal dissipation. For nonlinear systems we provide perturbative solutions by translating the formalism of canonical perturbation theory into the context of master equations. It is shown that unavoidable degeneracy causes an unfortunate reduction in accuracy between perturbative master equations and their solutions. We also extend the famous theorem of Lindblad, Gorini, Kossakowski and Sudarshan on completely positivity to non-Markovian master equations. Our application is primarily to model atoms interacting via a common electromagnetic field. The electromagnetic field contains correlations in both space and time, which are related to its relativistic (photon-mediated) nature. As such, atoms residing in the same field experience different environmental effects depending upon their relative position and orientation. Our more accurate solutions were necessary to assess sudden death of entanglement at zero temperature. In contrast to previous claims, we found that all initial states of two-level atoms undergo finite-time disentanglement. We were also able to access regimes which cannot be described by Lindblad equations and other simpler methods, such as near resonance. Finally we revisit the infamous Abraham-Lorentz force, wherein a single particle in motion experiences backreaction from the electromagnetic field. This leads to a number of well-known problems including pre-acceleration and runaway solutions. We found a more a more-suitable open-system treatment of the nonrelativistic particle to be perfectly causal and dissipative without any extraneous requirements for finite size of the particle, weak coupling to the field, etc..
In the limit of large particle numbers and low densities systems of cold atoms can be effectively described as macroscopic single-particle systems in a mean-field approximation. In the case of a Bose–Hubbard system, modelling bosons on a discrete lattice with on-site interactions, this yields a discrete nonlinear Schrödinger equation of Gross–Pitaevskii type. It has been recently shown that the correspondence between the Gross–Pitaevskii equation and the Bose–Hubbard system breaks down for complex extensions. In particular, for a Bose–Hubbard dimer with complex on-site energy the mean-field approximation yields a generalized complex nonlinear Schrödinger equation. Conversely, a Gross–Pitaevskii equation with complex on-site energies arises as the mean-field approximation of many-particle Lindblad dynamics rather than a complex extension of the Bose–Hubbard system. Here we address the question of how the mean-field description is modified in the presence of a complex-valued particle interaction term for a Bose–Hubbard dimer. We derive the mean-field equations of motion leading to nonlinear dissipative Bloch dynamics, related to a nontrivial complex generalization of the nonlinear Schrödinger equation. The resulting dynamics are analysed in detail. It is shown that depending on the parameter values there can be up to six stationary states, and for small values of the interaction strength there are limit cycles. Furthermore, we show how a Gross–Pitaevskii equation with a complex interaction term can be derived as the mean-field approximation of a Bose–Hubbard dimer with an additional Lindblad term modelling two-particle losses.
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.