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 $$\rho (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.
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.
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.
Based on the references [P1–P6], we investigate in this cumulative PhD thesis the dynamics and transport in closed and open interacting spin chain models. We focus on (i) the connection between the open and the closed system dynamics [P1, P2], (ii) the influence of various system parameters on the transport behavior in closed and open systems [P3, P5], and in addtion (iii) the dynamics and transport in closed models [P4, P6]. As for (i), we explore the connection of the dynamics in open systems, modeled by the Lindblad equation, and the correlation functions in the closed system [P1]. Building on this connection for magnetization transport, we study whether the time evolution of the open quantum system can be described in terms of classical correlation functions generated by the Hamiltonian equations of motion for real vectors [P2]. Second, we study in (ii) the influence of different parameters on the diffusion constants in closed and open systems [P3, P5]. In Ref. [P3], we study the impact of integrability-breaking perturbations in both closed and open spin chain systems. While we find a continuous change of the diffusion constant with perturbation strength in the closed system, in the open system, in the weak-perturbation limit, we indicate the possibility of a diverging diffusion constant. Also, we investigate the influence of the system-bath coupling on the transport behavior in boundary-driven open systems [P5]. We observe that the diffusion constants depend on the system-bath coupling and that in the nonequilibrium steady state, they inevitably retain finite-size effects, limiting the computation of long-time transport coefficients. Additional for (iii), in Ref. [P4], we first examine the impact of perturbations on the dynamics in closed spin chain models. We employ projector-operator techniques to understand, on perturbative grounds, the nontrivial damping of the standard correlation function. Second, we also discuss a dynamical quantum typicality approach for studying the energy transport in bipartite systems at low temperatures [P6].
In this Thesis we discussed several issues related to the dynamics of isolated quantum many-body systems. The brief overview of Chapter 1 highlighted a series of open fundamental theoretical questions which are crucial for the investigation of nonequilibrium phenomena. Among them, we mainly focused on the emergence of thermal behaviors and on the characterization of dynamical phase transitions and of the associated universal dynamics. To this end, we analyzed some paradigmatic models, such as the Fermi-Hubbard model, the quantum Ising chain, and the O(N) model. -Chapter 2: Absence of thermalization in a Fermi liquid. In this Chapter we report two results that question the common expectation of thermalization in a Landau-Fermi liquid after an interaction quench. We first consider the perturbative expansion in the interaction strength of the long-wavelength structure factor S(q) and show that it does not satisfy the hypothesis that steady-state averages correspond to thermal ones. In particular, S(q) has an analytical component ∼ const. + O(q^2) for q → 0, compatible with thermalization, but it retains also a nonanalytic term ∼ |q| characteristic of a Fermi liquid at zero-temperature. In real space, this nonanalyticity turns in a power-law decay of the density-density correlations, in contrast with the exponential decay associated to thermalization. We next analyze the low-density case, where one can obtain results valid at any order in interaction but at leading in the density, and find that in the steady state the momentum distribution jump at the Fermi surface is strictly finite, though smaller than at equilibrium. -Chapter 3: Nonadiabatic stationary behavior in a driven Ising chain. In this Chapter we discuss the emergence of a nonadiabatic behavior in the dynamics of the order parameter in a low-dimensional system subject to a linear ramp of one of its parameters within a gapped phase, which should be the most favorable situation for an adiabatic evolution. We study in details the dynamics of the spontaneous magnetization m_x(t) in a quantum Ising chain after a linear variation in time of the transverse field within the ordered phase. In particular, focusing on the value of the order parameter at the end of the ramp m_x(τ), we show that the smaller the switching rate of the transverse field is the closer m_x(τ) gets to its ground state value. Nonetheless, this apparently small correction to adiabaticity eventually leads to the disruption of the order exponentially fast in the subsequent time evolution, no matter how slowly the ramp is performed. -Chapter 4: Aging and coarsening in isolated quantum systems after a quench. This Chapter analyzes the nonequilibrium dynamics of an isolated quantum system after a sudden quench to the dynamical critical point, where the emergence of scaling and universal exponents is expected, due to the absence of time and length scales. We explore these features for a bosonic interacting scalar field theory with O(N) symmetry in the large-N limit, where the model is exactly solvable and the exponents and scaling forms of the relevant two-time correlation functions can be analytically calculated. Moreover, we provide evidence of the emergence of a dynamical scaling behavior also for quenches below the dynamical critical point, associated with coarsening. We find that the latter case is characterized by the same scaling functions as those describing the critical case, yet with different exponents. -Chapter 5: Dynamical transitions and statistics of excitations. In this Chapter we study the dynamics of the O(N) model (introduced in the previous Chapter) resulting from a different protocol: a linear ramp of one of its parameters. We find that the presence of a dynamical phase transition, as well as its critical properties, are robust against the change of the protocol. We show that a characterization based on the critical dimensions and exponents would suggest that the dynamical phase transition is analogous to the equilibrium thermal one. However, its nonequilibrium nature becomes evident analyzing the statistics of excitations produced in the ramp process. In particular, the critical properties are encoded in the fluctuations in the number of excitations, which display qualitatively different behaviors depending on the ramp being performed above, at, or below the dynamical critical point. These behaviors bear no dependence on the duration of the protocol.
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.
Often quantum systems are not isolated and interactions with their environments must be taken into account. In such open quantum systems these environmental interactions can lead to decoherence and dissipation, which have a marked influence on the properties of the quantum system. In many instances the environment is well-approximated by classical mechanics, so that one is led to consider the dynamics of open quantum-classical systems. Since a full quantum dynamical description of large many-body systems is not currently feasible, mixed quantum-classical methods can provide accurate and computationally tractable ways to follow the dynamics of both the system and its environment. This review focuses on quantum-classical Liouville dynamics, one of several quantum-classical descriptions, and discusses the problems that arise when one attempts to combine quantum and classical mechanics, coherence and decoherence in quantum-classical systems, nonadiabatic dynamics, surface-hopping and mean-field theories and their relation to quantum-classical Liouville dynamics, as well as methods for simulating the dynamics.
Often quantum systems are not isolated and interactions with their environments must be taken into account. In such open quantum systems these environmental interactions can lead to decoherence and dissipation, which have a marked influence on the properties of the quantum system. In many instances the environment is well-approximated by classical mechanics, so that one is led to consider the dynamics of open quantum-classical systems. Since a full quantum dynamical description of large many-body systems is not currently feasible, mixed quantum-classical methods can provide accurate and computationally tractable ways to follow the dynamics of both the system and its environment. This review focuses on quantum-classical Liouville dynamics, one of several quantum-classical descriptions, and discusses the problems that arise when one attempts to combine quantum and classical mechanics, coherence and decoherence in quantum-classical systems, nonadiabatic dynamics, surface-hopping and mean-field theories and their relation to quantum-classical Liouville dynamics, as well as methods for simulating the dynamics.