We numerically study the relaxation dynamics of impurity-host systems, focusing on the presence of long-lived metastable states in the nonequilibrium dynamics after an initial excitation of the impurities. In generic systems, an excited impurity coupled to a large bath at zero temperature is expected to relax and approach its ground state over time. However, certain exceptional cases exhibit metastability, where the system remains in an excited state on timescales largely exceeding the typical relaxation time. We study this phenomenon for three prototypical impurity models: a tight-binding quantum model of independent spinless fermions on a lattice with two stub impurities, a classical-spin Heisenberg model with two weakly coupled classical impurity spins, and a tight-binding quantum model of independent electrons with two classical impurity spins. Through numerical integration of the fundamental equations of motion, we find that all three models exhibit similar qualitative behavior: complete relaxation for nearest-neighbor impurities and incomplete or strongly delayed relaxation for next-nearest-neighbor impurities. The underlying mechanisms leading to this behavior differ between models and include impurity-induced bound states, emergent approximately conserved local observables, and exact cancellation of local and nonlocal dissipation effects. Published by the American Physical Society 2024
We numerically study the relaxation dynamics of impurity-host systems, focusing on the presence of long-lived metastable states in the nonequilibrium dynamics after an initial excitation of the impurities. In generic systems, an excited impurity coupled to a large bath at zero temperature is expected to relax and approach its ground state over time. However, certain exceptional cases exhibit metastability, where the system remains in an excited state on timescales largely exceeding the typical relaxation time. We study this phenomenon for three prototypical impurity models: a tight-binding quantum model of independent spinless fermions on a lattice with two stub impurities, a classical-spin Heisenberg model with two weakly coupled classical impurity spins, and a tight-binding quantum model of independent electrons with two classical impurity spins. Through numerical integration of the fundamental equations of motion, we find that all three models exhibit similar qualitative behavior: complete relaxation for nearest-neighbor impurities and incomplete or strongly delayed relaxation for next-nearest-neighbor impurities. The underlying mechanisms leading to this behavior differ between models and include impurity-induced bound states, emergent approximately conserved local observables, and exact cancellation of local and nonlocal dissipation effects. Published by the American Physical Society 2024
Impurities coupled to superconductors offer a controlled platform to understand the interplay between superconductivity, many-body interactions, and nonequilibrium physics. In the equilibrium situation, local interactions at the impurity induce a transition from the spin-singlet to the spin-doublet ground state, resulting in a supercurrent sign reversal ($0\text{\ensuremath{-}}\ensuremath{\pi}$ transition). In this work, we apply the exact time-dependent density matrix renormalization group method to simulate the transient dynamics of such superconducting systems. We also use a perturbative approximation to analyze their properties at longer times. These two methods agree for a wide range of parameters. In a phase-biased situation, the system gets trapped in a metastable state characterized by a lower supercurrent compared to the equilibrium case. We show that local Coulomb interactions do not provide an effective relaxation mechanism for the initially trapped quasiparticles. In contrast, other relaxation mechanisms, such as coupling to a third normal lead, make the impurity spin relax for parameter values corresponding to the equilibrium 0 phase. For parameters corresponding to the equilibrium $\ensuremath{\pi}$ phase the impurity converges to a spin-polarized stationary state. Similar qualitative behavior is found for a voltage-biased junction, which provides an effective relaxation mechanism for the trapped quasiparticles in the junction.
An important step towards a comprehensive understanding of far-from-equilibrium dynamics of quantum many-body systems is the identification of unifying features that are independent of microscopic details of the system. We experimentally observe such robust features in the magnetization relaxation dynamics of disordered Heisenberg XX, XXZ, and Ising Hamiltonians. We realize these Heisenberg spin models with tunable anisotropy parameter and power-law interactions in an ensemble of Rydberg atoms by encoding the spin in suitable Rydberg state combinations. We consistently observe stretched-exponential relaxation of magnetization for all considered spin models, collapsing onto a single curve after appropriate rescaling of time. This robust short-time relaxation behavior is explained by a perturbative treatment that exploits the strong disorder in pairwise couplings, which leads to a description in terms of approximately independent pairs of spins. In numerical simulations of small systems, we show that these pairs of spins constitute approximate local integrals of motion, which remain at least partially conserved on a timescale exceeding the duration of the relaxation dynamics of the magnetization. Published by the American Physical Society 2024
Open quantum systems with nearly degenerate energy levels have been shown to exhibit long-lived metastable states in the approach to equilibrium, even when modelled with certain Lindblad-form quantum master equations. This is a result of dramatic separation of timescales due to differences between Liouvillian eigenvalues. These metastable states often have nonzero coherences that die off only in the long-time limit once the system reaches thermal equilibrium. We examine two distinct situations that give rise to this effect: one in which dissipative dynamics couple together states only within a nearly degenerate subspace, and one in which they give rise to jumps over finite-energy splittings, between nearly degenerate subspaces. We find, in each case, that a change of basis can often lead to a representation that more naturally captures the impact of the system-bath interaction than does the energy eigenbasis, revealing that separate timescales are associated with separate processes (e.g., decoherence into a nonenergy eigenbasis, decay of population correlations to the initial state). This approach is paired with the inspection of quantum trajectories, which further provide intuition as to how open system evolution is characterized when coherent oscillations, thermal relaxation, and decoherence all occur simultaneously. Published by the American Physical Society 2024
The nonequilibrium dynamics of an impurity immersed with a finite velocity into a one-dimensional system of weakly interacting bosons is studied within the framework of the time-dependent Gross-Pitaevskii equation. We uncover and characterize different regimes of relaxation dynamics. We find that the final impurity velocity remains constant in a large interval of sufficiently big and realistic initial velocities. The underlying physical mechanism is the emission of the dispersive density shock waves that carry away the excess of the initial impurity momentum, while locally the system remains in the same stationary state. In contrast, a heavy impurity with the same coupling constant relaxes differently and the regime of constant final velocity disappears. Furthermore, a fast heavy impurity exhibits damped velocity oscillations in time before reaching a stationary state. This process is accompanied by the oscillations of the local depletion of the boson density around the impurity, until their positions coincide and they continue the motion together. Decreasing the impurity-boson coupling or increasing the strength of repulsion between bosons, the oscillations get amplified. In the case of a heavy impurity with the mass bigger than the critical one, the ground-state energy as a function of momentum exhibits cusps and metastable branches. We show that they manifest themselves by a soliton emission, a considerable slowing down of the relaxation, and a change of the impurity direction of motion with respect to the initial one.
We present a general theory of classical metastability in open quantum systems. Metastability is a consequence of a large separation in timescales in the dynamics, leading to the existence of a regime when states of the system appear stationary, before eventual relaxation toward a true stationary state at much larger times. In this work, we focus on the emergence of classical metastability, i.e., when metastable states of an open quantum system with separation of timescales can be approximated as probabilistic mixtures of a finite number of states. We find that a number of classical features follow from this approximation, for the manifold of metastable states, long-time dynamics between them, and symmetries of the dynamics. Namely, those states are approximately disjoint and thus play the role of metastable phases, the relaxation toward the stationary state is approximated by a classical stochastic dynamics between them, and weak symmetries correspond to their permutations. Importantly, the classical dynamics is observed not only on average but also at the level of individual quantum trajectories: We show that time coarse-grained continuous measurement records can be viewed as noisy classical trajectories, while their statistics can be approximated by that of the classical dynamics. Among others, this explains how first-order dynamical phase transitions arise from metastability. Finally, to verify the presence of classical metastability in a given open quantum system, we develop an efficient numerical approach that delivers the set of metastable phases together with the effective classical dynamics. Since the proximity to a first-order dissipative phase transition manifests as metastability, the theory and tools introduced in this work can be used to investigate such transitions---which occur in the large size limit---through the metastable behavior of many-body systems of moderate sizes accessible to numerics.
Theoretical studies of electron and energy transport in nanostructures usually involve nonequilibrium quantum impurity models, which feature a locally interacting quantum system (impurity) coupled to noninteracting bosonic and/or electronic baths each at its own equilibrium. In these models, the interplay of strong correlation, nonequilibrium effects and dissipation leads to rich and complex phenomena. As a nonequilibrium strongly correlated many-body problem, these models, except at some rare cases, have no exact analytical solution, and hence their treatment heavily relies upon numerical techniques. The goal of this thesis is to develop numerically accurate and relatively low-cost techniques for treatment of strongly correlated open nonequilibirum quantum systems. Such impurity solvers are in high demand both as standalone methods for simulation of response in nanoscale junctions and as part of divide-and-conquer schemes such as dynamical mean field theory (DMFT).For nonequilibirum open quantum systems, two characteristic energy scales governing the physics are intra-system interactions and strength of the system-baths couplings. Thus, systems strongly coupled to their baths with weak intra-system interactions or weakly coupled systems with strong local interactions can be treated within developed perturbation theory approaches. My research focuses on development of nonequilibirum dual techniques, which are capable to accurately treat situations where perturbative techniques are inapplicable due to absence of a small parameter in the system and which reduce to standard perturbative expansions in the limits of weak coupling or weak intra-system interaction. \tDual techniques rely on ability to solve exactly a simplified reference system. The latter includes exactly all intra-system interactions and accounts approximately for system-bath couplings. Dual techniques (described below) allow to account for deviation from the correct hybridization function in a organized and controlled manner of superperturbative expansion. Here, I use Markovian Lindbald quantum master equation (QME) as a solver for the reference system. In the thesis I discuss possibility of mapping between true non-Markov dynamics and simpler Markov Lindblad QME (mapping between physical and reference systems), numerical methods for the mapping and for solution of the resulting Lindblad QME, and application of the mapping in novel dual formulations: auxiliary-QME dual fermion and dual boson approaches. \tFirst I analytically prove that the continuous fermionic baths that the impurity is coupled to can be mapped onto a set of dissipative auxiliary modes under Lindblad QME. Based on this mapping, complicated non-Markovian dynamics of the open fermionic system can be treated by a much simpler Markovian Lindblad QME. However, nice mapping generally requires many auxiliary modes, and numerically only a small number of modes can be exactly diagonalized due to exponential increase of the Hilbert space. To perturbatively treat the difference between the true bath and the small-size auxiliary one, we introduce dual-fermion method. It can directly target the steady state and avoid long-time propagation. Besides current-voltage relation, it can also provide spectral functions, indicating itself as a potential impurity solver for DMFT. Furthermore, we generalize dual-fermion method to consider additionally the effects of bosonic environment; hence it is called dual-boson method. It can be used for studying electron-phonon interaction and light-matter interaction in nanojunctions. By comparing with numerically exact results we argue that both methods are very accurate and relatively cheap. \tFinally, matrix product state representation of the Lindblad QME enables us to treat a large number of auxiliary modes, which provide an exponentially improved mapping. Adding counting field to the auxiliary QME, we can compute the generating function of time-dependent full-counting statistics of electron transport, which yields not only electronic current but also higher-order cumulants (e.g., zero-frequency noise).