The discretization approximation method commonly used to simulate the dynamics of quantum systems coupled to the environment in continuum often suffers from the periodically partial recovery of the initial state because of the effect of finite dimension, dubbed the recurrence. To address this issue, we propose a generalization of the discretization approximation method into the complex frequency space based on complex Gauss quadratures. An effective Hamiltonian can be established by this way, which is non-Hermitian and demonstrates the complex energy modes with a negative imaginary part, describing the dissipation of the system. This method is applied to examine the dynamics in two exactly solvable models, the dephasing model and the single-excitation dissipative dynamics in the Aubry–André–Harper model. By comparison with the exact numerics and analytical results, it is found that our approach not only significantly reduces the effect of recurrence and improves the effectiveness of calculation but also provides a unique perspective into the dynamics of open systems from the point of complex energy levels. Furthermore, we establish a simple relationship between the parameters in computation and the effectiveness of simulation by analyzing the computational error.
The discretization approximation method commonly used to simulate the dynamics of quantum systems coupled to the environment in continuum often suffers from the periodically partial recovery of the initial state because of the effect of finite dimension, dubbed the recurrence. To address this issue, we propose a generalization of the discretization approximation method into the complex frequency space based on complex Gauss quadratures. An effective Hamiltonian can be established by this way, which is non-Hermitian and demonstrates the complex energy modes with a negative imaginary part, describing the dissipation of the system. This method is applied to examine the dynamics in two exactly solvable models, the dephasing model and the single-excitation dissipative dynamics in the Aubry–André–Harper model. By comparison with the exact numerics and analytical results, it is found that our approach not only significantly reduces the effect of recurrence and improves the effectiveness of calculation but also provides a unique perspective into the dynamics of open systems from the point of complex energy levels. Furthermore, we establish a simple relationship between the parameters in computation and the effectiveness of simulation by analyzing the computational error.
The study of phase transitions in dissipative quantum systems based on the Liouvillian is often hindered by the difficulty of constructing a time-local master equation when the system-environment coupling is strong. To address this issue, the complex discretization approximation for the environment is proposed to study the quantum phase transition in the extended Jaynes-Cumming model with an infinite number of boson modes. This approach yields a non-Hermitian effective Hamiltonian that can be used to simulate the dynamics of the spin. It is found that the ground state of this effective Hamiltonian determines the spin dynamics in the single-excitation subspace. Depending on the opening of the energy gap and the maximum population of excitations on the spin degree of freedom, three distinct phases can be identified: fast decaying, localized, and stretched dynamics of the spin. This approach can be extended to multiple excitations, and similar dynamics were found in the double-excitation subspace, indicating the robustness of the single-excitation phase.
Understanding the dynamics of dissipative quantum systems, particularly beyond the weak coupling approximation, is central to various quantum applications. While numerically exact methods provide accurate solutions, they often lack the analytical insight provided by theoretical approaches. In this study, we employ the recently developed method dubbed the effective Hamiltonian theory to understand the dynamics of system-bath configurations without resorting to a perturbative description of the system-bath coupling energy. Through a combination of mapping steps and truncation, the effective Hamiltonian theory offers both analytical insights into signatures of strong couplings in open quantum systems and a straightforward path for numerical simulations. To validate the accuracy of the method, we apply it to two canonical models: a single spin immersed in a bosonic bath and two noninteracting spins in a common bath. In both cases, we study the transient regime and the steady state limit at nonzero temperature and spanning system-bath interactions from the weak to the strong regime. By comparing the results of the effective Hamiltonian theory with numerically exact simulations, we show that although the former overlooks non-Markovian features in the transient equilibration dynamics, it correctly captures non-perturbative bath-generated couplings between otherwise non-interacting spins, as observed in their synchronization dynamics and correlations. Altogether, the effective Hamiltonian theory offers a powerful approach for understanding strong coupling dynamics and thermodynamics, capturing the signatures of such interactions in both relaxation dynamics and in the steady state limit.
Non-Markovian effects are important in modeling the behavior of open quantum systems arising in solid-state physics, quantum optics as well as in study of biological and chemical systems. A common approach to the analysis of such systems is to approximate the non-Markovian environment by discrete bosonic modes thus mapping it to a Lindbladian or Hamiltonian simulation problem. While systematic constructions of such modes have been proposed in previous works [D. Tamascelli et al, PRL (2012), A. W. Chin et al, J. of Math. Phys (2010)], the resulting approximation lacks rigorous convergence guarantees. In this paper, we initiate a rigorous study of the convergence properties of these methods. We show that under some physically motivated assumptions on the system-environment interaction, the finite-time dynamics of the non-Markovian open quantum system computed with a sufficiently large number of modes is guaranteed to converge to the true result. Furthermore, we show that, for most physically interesting models of non-Markovian environments, the approximation error falls off polynomially with the number of modes. Our results lend rigor to numerical methods used for approximating non-Markovian quantum dynamics and allow for a quantitative assessment of classical as well as quantum algorithms in simulating non-Markovian quantum systems.
Abstract In the conventional quantum mechanics of conserved systems, Hamiltonian is assumed to be a Hermitian operator. However, when it comes to quantum systems in presence of dissipation and/or noise, including open quantum optical systems, the strict hermiticity requirement is nor longer necessary. In fact, it can be substantially relaxed: the non-Hermitian part of a Hamiltonian is allowed, in order to account for effects of dissipative environment, whereas its Hermitian part would be describing subsystem’s energy. Within the framework of the standard approach to dissipative phenomena based on a master equation for the reduced density operator, we propose a replacement of the hermiticity condition by a more general condition of commutativity between Hermitian and anti-Hermitian parts of a Hamiltonian. As an example, we consider a dissipative two-mode quantum system coupled to a single-mode electromagnetic wave, where we demonstrate that the adjoint-commutativity condition does simplify the parametric space of the model.
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.
We propose a numerical technique based on a combination of short-iterative Lanczos and exact diagonalization methods, suitable for simulating the time evolution of the reduced density matrix of a single qubit interacting with an environment. By choosing a mode discretization method and a flexible bath states truncation scheme, we are able to include in the physical description multiple-excitation processes, beyond weak coupling and Markov approximations. We apply our technique to the simulation of three different model Hamiltonians, which are relevant in the field of adiabatic quantum computation. We compare our results with those obtained on the basis of the widely used Lindblad master equation, as well as with well-known exact and approximated approaches. We show that our method is able to recover the thermodynamic behavior of the qubit-bath system, beyond the Born-Markov approximation. Finally, we show that even in the case of the adiabatic quantum annealing of a single qubit the bath can be beneficial in reaching the reduced system ground state.