In this report, we study the time evolution of a two-site Hubbard model when immersed in different environments by applying the Lindblad Master Equation. The equation was derived from microscopic dynamics and analysed in terms of symmetries and eigenvalues. Three environments were modelled; one which enabled the loss and gain of particles, one modelling continuous observation of the particle count and the final environment preserving particle-hole symmetry. The spectral gap of the Liouvillian, which determines the asymptotic decay rate, was numerically calculated for each environment with different parameter values (J, U, μ) for the Hubbard Hamiltonian. The ensuing time evolution was also simulated using a fourth-order Runge-Kutta scheme. For the third set of operators, a closure of the spectral gap was observed at U=2μ, as well as a non-decoherent steady state. This was attributed to the particle-hole symmetry of the Hubbard model. To test this hypothesis, Lindblad operators highlighting this symmetry were also simulated, showing expected results. An outlook for further study was formulated to address remaining questions concerning the interpretation of these results.
In this report, we study the time evolution of a two-site Hubbard model when immersed in different environments by applying the Lindblad Master Equation. The equation was derived from microscopic dynamics and analysed in terms of symmetries and eigenvalues. Three environments were modelled; one which enabled the loss and gain of particles, one modelling continuous observation of the particle count and the final environment preserving particle-hole symmetry. The spectral gap of the Liouvillian, which determines the asymptotic decay rate, was numerically calculated for each environment with different parameter values (J, U, μ) for the Hubbard Hamiltonian. The ensuing time evolution was also simulated using a fourth-order Runge-Kutta scheme. For the third set of operators, a closure of the spectral gap was observed at U=2μ, as well as a non-decoherent steady state. This was attributed to the particle-hole symmetry of the Hubbard model. To test this hypothesis, Lindblad operators highlighting this symmetry were also simulated, showing expected results. An outlook for further study was formulated to address remaining questions concerning the interpretation of these results.
We investigate the dynamical properties of one-dimensional dissipative Fermi-Hubbard models, which are described by the Lindblad master equations with site-dependent jump operators. The corresponding non-Hermitian effective Hamiltonians with pure loss terms possess parity-time (PT) symmetry if we compensate the system additionally with an overall gain term. By solving exactly the two-site Lindblad equation with fixed dissipation, we find that the dynamics of the rescaled density matrix shows an instability as the interaction increases over a threshold, which can be equivalently described in the scheme of non-Hermitian effective Hamiltonians. This instability is also observed in multisite systems and closely related to the PT-symmetry breaking accompanied by appearance of complex eigenvalues of the effective Hamiltonian. Moreover, we unveil that the dynamical instability of the antiferromagnetic Mott phase comes from the PT-symmetry breaking in highly excited bands, although the low-energy effective model of the non-Hermitian Hubbard model in the strongly interacting regime is always Hermitian. We also provide a quantitative estimation of the time for the observation of dynamical PT-symmetry breaking which could be probed in experiments.
Abstract Relaxation processes in quantum systems coupled to external environments represent one of the most fundamental nonequilibrium phenomena in condensed matter physics. The Lindblad master equation provides a powerful framework for characterizing such open quantum dynamics. In this work, we systematically investigate how different types of quantum jump operators and system geometries influence the Liouvillian gap and the properties of the nonequilibrium steady state (NESS) in finite-size systems. We demonstrate that, due to the intricate structure of the Liouvillian superoperator, multiple NESSs with unphysical characteristics can emerge. The physically meaningful steady state must instead be understood as a superposition of these NESSs that collectively satisfy the required physical constraints. Furthermore, we find that the Liouvillian gap does not necessarily increase monotonically with the system–environment coupling strength. Instead, it can exhibit a nontrivial peak structure, corresponding to a minimum in the relaxation time. The magnitude of this peak is closely related to the symmetry properties of the system. Our results provide a deeper understanding of nonequilibrium behavior in finite quantum systems and offer new insights into the design and control of open quantum dynamics.
Markovian open quantum systems display complicated relaxation dynamics. The spectral gap of the Liouvillian characterizes the asymptotic decay rate toward the steady state, but it does not necessarily give a correct estimate of the relaxation time because the crossover time to the asymptotic regime may be too long. We here give a rigorous upper bound on the transient decay of autocorrelation functions in the steady state by introducing the symmetrized Liouvillian gap. The standard Liouvillian gap and the symmetrized one are identical in an equilibrium situation but differ from each other in the absence of the detailed balance condition. It is numerically shown that the symmetrized Liouvillian gap always gives a correct upper bound on the decay of the autocorrelation function, but the standard Liouvillian gap does not.
This thesis explores the dynamics of linearly coupled quantum and classical oscillators interacting with environments. In the quantum regime, such dynamics is described in terms of the Liouvillian superoperator, which leads to non-unitary time evolution of the system, specifically of the density matrix. In particular, when the correlation time of the environment is significantly shorter than the characteristic timescale of the system, the dynamics are classified as Markovian. Therefore, the equation of motion for the density matrix that fulfils this property is known as the Lindblad master equation. Systems of coupled cavities and waveguides in the presence of gain and loss can be described by this equation and, more importantly, by non-Hermitian Hamiltonians. In particular, when two coupled cavities exhibit balance between gain and loss, the system is PT -symmetric. However, for certain values in the parameter space, this symmetry is broken and the non-Hermitian Hamiltonian is no longer diagonalizable. Such values are denoted as exceptional points. In the case of large arrays of coupled quantum oscillators, the dynamics description becomes more challenging since the dimensions of the density matrix grow exponentially. To address this, we study the dynamics of these arrays via the quantum state diffusion approach. In such approach, the Lindblad master equation is unraveled by a stochastic Schrödinger equation. Then, we decouple this equation by diagonalizing the non-Hermitian Hamiltonian that describes the system. Such diagonalization leads to a transformation of a system of coupled oscillators into a decoupled one, offering a great reduction in the basis size. For systems of classical oscillators with dissipative nearest-neighbor coupling, the dynamics leads to synchronized motion. We study how an array of such oscillators reaches the synchronized state. This behavior is found when the system is dissipation free, which implies that all the eigenvectors of the dynamical matrix are damped, except one. To measure the degree of synchronization, we introduce the inverse of participation ratio, whose maximum value is associated with a delocalized state and, therefore, fully synchronized motion. For larger arrays of oscillators, we examine the case where the frequency fluctuations of each oscillator are obtained from a Gaussian distribution with zero mean and a given standard deviation. Then, we derive a scaling law for the largest value of the mentioned standard deviation that still allows the oscillators fully synchronize their motion.:1 Introduction 2 Open quantum systems 3 Dynamics of two coupled cavities in the presence of losses 4 Unraveling of Lindblad equation for system of oscillators 5 Synchronization of dissipatively coupled harmonic oscillators 6 Conclusions and outlook A Derivation of equations of motion for the amplitudes B Quantum state diffusion unraveling for N decoupled oscillators C Model for disspative coupling D Perturbation theory for disorder in the frequencies References
In open quantum systems, reduced dynamics is commonly described by a master equation, whose Liouvillian gap closing (LGC) typically signals the emergence of decoherence-free subspace. By contrast, the dynamics of the full system-environment compound is governed by the underlying Hamiltonian spectrum, where bound states in the continuum (BICs) can protect long-lived quantum resources. Despite these parallel perspectives, the relation between LGC and BIC formation has remained largely unexplored. Here we bridge this gap in a paradigmatic giant-atom waveguide platform and show that the occurrence of LGC necessarily benchmarks the presence of a BIC in the full Hamiltonian description. By engineering the giant-atom geometry, we further demonstrate rich dynamical regimes-including Rabi oscillations, fractional decay, and complete exponential relaxation-depending on the number of supported BICs, which can be tuned from three to zero. Remarkably, when two BICs become frequency-degenerate, the long-time dynamics approaches a steady state rather than exhibiting persistent oscillations. Our results establish a direct spectral-dynamical connection between effective Markovian and underlying non-Markovian descriptions, and provide a route toward flexible control of open-system dynamics.
We investigate theoretically the time- and frequency-domain two-particle correlations of a driven-dissipative Bose-Hubbard model at and near a dissipative phase transition (DPT). We compute the Hanbury Brown--Twiss (HBT)-type two-particle temporal correlation function ${g}^{2}(\ensuremath{\tau})$ which, as a function of time delay $\ensuremath{\tau}$, exhibits oscillations with frequencies determined by the imaginary part of the Liouvillian gap. As the gap closes near a transition point, the oscillations at that point die down. For parameters slightly away from the transition point, the HBT correlations show oscillations from superbunching to antibunching regimes. We show that the Fourier transform of HBT correlations into the frequency domain provide information about DPT and Liouvillian dynamics. We numerically solve the many-body Lindblad master equation and calculate the Wigner distribution of the system in the steady state to ascertain the DPT. Below a certain drive strength, the Fourier transform shows a two-peak structure, while above that strength it exhibits either a Lorentzian-like single-peak structure or a structure with two dips. The width of the single-peak structure is minimal at the phase-transition point and the peak of this structure always lies at zero frequency. The positions of the two symmetrical peaks in case of a two-peak structure are given by the imaginary parts of the Liouvillian gap while their half width at half maximum is given by the real part of the gap. The positions and widths of the two dips are also related to low-lying eigenvalues of the Liouvillian operator. We discuss quantum statistical properties of the model in terms of the HBT correlation function and its Fourier transform.