Abstract The Lindblad equation describes the dissipative time evolution of a density matrix that characterizes an open quantum system in contact with its environment. The widespread ensemble interpretation of a density matrix requires its time evolution to be linear. However, when the dynamics of the density matrix is of a quantum system results not only from the interaction with an external environment, but also with other quantum systems of the same type, the ensemble interpretation is inappropriate and nonlinear dynamics arise naturally. We therefore study the dynamical behavior of nonlinear Lindblad equations using the example of a two-level system. By using techniques developed for classical dynamical systems we show that various types of bifurcations and even chaotic dynamics can occur. As specific examples that display the various types of dynamical behavior, we suggest explicit models based on systems of interacting spins at finite temperature and exposed to a magnetic field that can change in dependence of the magnetization. Due to the interaction between spins, which is treated at mean-field level, the Hamiltonian as well as the transition rates of the Lindblad equation become dependent on the density matrix.
Abstract The Lindblad equation describes the dissipative time evolution of a density matrix that characterizes an open quantum system in contact with its environment. The widespread ensemble interpretation of a density matrix requires its time evolution to be linear. However, when the dynamics of the density matrix is of a quantum system results not only from the interaction with an external environment, but also with other quantum systems of the same type, the ensemble interpretation is inappropriate and nonlinear dynamics arise naturally. We therefore study the dynamical behavior of nonlinear Lindblad equations using the example of a two-level system. By using techniques developed for classical dynamical systems we show that various types of bifurcations and even chaotic dynamics can occur. As specific examples that display the various types of dynamical behavior, we suggest explicit models based on systems of interacting spins at finite temperature and exposed to a magnetic field that can change in dependence of the magnetization. Due to the interaction between spins, which is treated at mean-field level, the Hamiltonian as well as the transition rates of the Lindblad equation become dependent on the density matrix.
The Lindblad equation describes the dissipative time evolution of a density matrix that characterizes an open quantum system in contact with its environment. The widespread ensemble interpretation of a density matrix requires its time evolution to be linear. However, when the density matrix is obtained from a mean field theory of interacting quantum systems or from a top-down control by a changing classical environment, the ensemble interpretation is inappropriate and nonlinear dynamics arise naturally. We therefore study the dynamical behavior of nonlinear Lindblad equations using the example of a two-level system. By using techniques developed for classical dynamical systems we show that various types of bifurcations and even chaotic dynamics can occur. We also discuss experimental situations for which our results could be relevant.
The dynamics of the one-dimensional spin model (s = 1/2) XXZ in the quantum “quenching” mode is investigated on a classical computer using the Qiskit and LindbladMPO software. The quantitative description of the time behavior of local magnetizations and spin correlation functions are proposed to be considered as a standard for comparative analysis with similar results previously obtained on existing noisy intermedium-scale quantum (NISQ) computers. Quantum quenching of the spin system is created at the initial time, t=0 , by an abrupt (nonadiabatic) change of the spin Hamiltonian parameter. The subsequent ( t>0 ) time evolution of the reduced density matrix ρ(t) of the system is traced from the solutions of the quantum Liouville equation in the case of an isolated system and its generalization in the form of the Lindblad equation for an open system interacting with the environment. In the case of an open quantum system, dissipation processes are caused by the influence of fluctuating environmental fields on the spin system.
We study a periodically driven macrospin system with anisotropic long-range interactions and collective dissipation, described by a Lindblad master equation. In the thermodynamic limit (N→∞), a mean-field treatment yields classical equations of motion, whose dynamics are characterized via the maximal Lyapunov exponent (MLE). Focusing on the thermodynamic limit, we map out chaotic, quasiperiodic, and periodic phases via bifurcation diagrams, MLEs, and Fourier spectra of evolved observables, identifying classic period-doubling bifurcations and fractal boundaries in the regions of attractors. Finite-size quantum simulations in the Dicke basis reveal that while both quantum and classical systems exhibit diverse dynamical phases, finite-size effects suppress some behaviors present in the thermodynamic limit. The sign of λmax serves as a key indicator of convergence between quantum and classical dynamics, which agree over timescales up to the Lyapunov time. Analysis of the density matrix shows that convergence occurs only when its nonzero elements are sharply localized. However, the nonconvergence does not imply a fundamental difference between quantum and classical dynamics: in chaotic regimes, although the evolution orbits of quantum and classical systems show significant differences, quantum evolution becomes mixed and diffusively explores the Hilbert space, signaling quantum chaos, which can be confirmed by the delocalized nature of the density matrix.
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.
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..
Abstract This chapter introduces the density operator and the Liouville equation and outlines how reduced dynamical description of open quantum systems naturally leads to the appearance of relaxation processes. (1) The density operator and the quantum Liouville equation (the density matrix for a pure system, statistical mixtures, representations, coherences, thermodynamic equilibrium). (2) An example: The time evolution of a two-level system in the density matrix formalism. (3) Reduced descriptions. (4) Time evolution equations for reduced density operators: The quantum master equation (projection operators, The Nakajima–Zwanzig equation, Derivation of the quantum master equation using the thermal projector, The quantum master equation in the interaction representation, The Markovian limit—the Redfield equation, The secular approximation and the Lindblad equation). (5) The Lindblad equation in the ‘Heisenberg representation’. (6) Relaxation of a two-level system and a harmonic oscillator revisited. (7) The optically driven two-level system in a thermal environment—the Bloch equations. (8) The oscillating dipole model of a two-level molecule. (9) Analogy of a coupled two-level system to a spin ½ system in a magnetic field.