Abstract We formulate and discuss the explicit computation of dynamic correlation functions in open quadratic fermionic systems which are driven and dissipated by Lindblad jump processes that are linear in canonical fermionic operators. Dynamic correlators are interpreted in terms of a local quantum quench where the pre-quench state is a non-equilibrium steady state, i.e. a fixed point of the Liouvillian. As an example, we study the XY spin 1/2 chain and quadratic Majorana chains with boundary Lindblad driving, whose dynamics exhibit asymmetric (skewed) light cone behaviour. We also numerically treat the two-dimensional XY model and the XY spin chain with additional Dzyaloshinskii–Moriya interactions. The latter exhibits a new non-thermal phase transition which can be understood in terms of bifurcations of the quasi-particle dispersion relation. Finally, considering in some detail the periodic quadratic fermionic ring with dissipation at a single (arbitrary) site, we present analytical expressions for first order corrections (in the strength of dissipation) to the spectrum and non-equilibrium steady state (NESS) correlation functions.
Abstract We formulate and discuss the explicit computation of dynamic correlation functions in open quadratic fermionic systems which are driven and dissipated by Lindblad jump processes that are linear in canonical fermionic operators. Dynamic correlators are interpreted in terms of a local quantum quench where the pre-quench state is a non-equilibrium steady state, i.e. a fixed point of the Liouvillian. As an example, we study the XY spin 1/2 chain and quadratic Majorana chains with boundary Lindblad driving, whose dynamics exhibit asymmetric (skewed) light cone behaviour. We also numerically treat the two-dimensional XY model and the XY spin chain with additional Dzyaloshinskii–Moriya interactions. The latter exhibits a new non-thermal phase transition which can be understood in terms of bifurcations of the quasi-particle dispersion relation. Finally, considering in some detail the periodic quadratic fermionic ring with dissipation at a single (arbitrary) site, we present analytical expressions for first order corrections (in the strength of dissipation) to the spectrum and non-equilibrium steady state (NESS) correlation functions.
Abstract The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k -point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.
The dynamical correlations of nonlocal operators in general quadratic open fermion systems is still a challenging problem. Here we tackle this problem by developing a new formulation of open fermion many-body systems, namely, the characteristic function approach. Illustrating the technique, we analyze a finite Kitaev chain with boundary dissipation and consider anyon-type nonlocal excitations. We give explicit formula for the Green's functions, demonstrating an asymmetric light cone induced by the anyon statistical parameter and an increasing relaxation rate with this parameter. We also analyze some other types of nonlocal operator correlations such as the full counting statistics of the charge number and the Loschmidt echo in a quench from the vacuum state. The former shows clear signature of a nonequilibrium quantum phase transition, while the later exhibits cusps at some critical times and hence demonstrates dynamical quantum phase transitions.
We obtain exact formulas for the time-dependence of a few physical observables for the open XX spin chain with Lindbladian dynamics. Our analysis is based on the fact that the Lindblad equation for an arbitrary open quadratic system of N N fermions is explicitly solved in terms of diagonalization of a 4N\times4N 4 N × 4 N matrix called structure matrix by following the scheme of the third quantization. We mainly focus on the time-dependence of magnetization and spin current. As a short-time behavior at a given site, we observe the plateau regime except near the center of the chain. Basic features of this are explained by the light-cone structure created by propagations of boundary effects from the initial time, but we can explain their more detailed properties analytically using our exact formulas. On the other hand, after the plateau regime, the magnetization and spin current exhibit a slow decay to the steady state values described by the Liouvillian gap. We analytically establish its O(N^{-3}) O ( N − 3 ) scaling and also determine its coefficient.
We study spectral and steady-state properties of generic Markovian dissipative systems described by quadratic fermionic Liouvillian operators of the Lindblad form. The Hamiltonian dynamics is modeled by a generic random quadratic operator, i.e., as a featureless superconductor of class D, whereas the Markovian dissipation is described by M M random linear jump operators. By varying the dissipation strength and the ratio of dissipative channels per fermion, m=M/(2N_F) m = M / ( 2 N F ) , we find two distinct phases where the support of the single-particle spectrum has one or two connected components. In the strongly dissipative regime, this transition occurs for m=1/2 m = 1 / 2 and is concomitant with a qualitative change in both the steady-state and the spectral gap that rules the large-time dynamics. Above this threshold, the spectral gap and the steady-state purity qualitatively agree with the fully generic (i.e., non-quadratic) case studied recently. Below m=1/2 m = 1 / 2 , the spectral gap closes in the thermodynamic limit and the steady-state decouples into an ergodic and a nonergodic sector yielding a non-monotonic steady-state purity as a function of the dissipation strength. Our results show that some of the universal features previously observed for fully random Liouvillians are generic for a sufficiently large number of jump operators. On the other hand, if the number of dissipation channels is decreased the system can exhibit nonergodic features, rendering it possible to suppress dissipation in protected subspaces even in the presence of strong system-environment coupling.
We study the quench dynamics of noninteracting fermionic quantum many-body systems that are subjected to Markovian drive and dissipation and are described by a quadratic Liouvillian which has parity-time (PT) symmetry. In recent work, we have shown that such systems relax locally to a maximum entropy ensemble that we have dubbed the PT-symmetric generalized Gibbs ensemble (PTGGE), in analogy to the generalized Gibbs ensemble that describes the steady state of isolated integrable quantum many-body systems after a quench. Here, using driven-dissipative versions of the Su-Schrieffer-Heeger (SSH) model and the Kitaev chain as paradigmatic model systems, we corroborate and substantially expand upon our previous results. In particular, we confirm the validity of a dissipative quasiparticle picture at finite dissipation by demonstrating light cone spreading of correlations and the linear growth and saturation to the PTGGE prediction of the quasiparticle-pair contribution to the subsystem entropy in the PT-symmetric phase. Further, we introduce the concept of directional pumping phases, which is related to the non-Hermitian topology of the Liouvillian and based upon qualitatively different dynamics of the dual string order parameter and the subsystem fermion parity in the SSH model and the Kitaev chain, respectively: depending on the postquench parameters, there can be pumping of string order and fermion parity through both ends of a subsystem corresponding to a finite segment of the one-dimensional lattice, through only one end, or there can be no pumping at all. We show that transitions between dynamical pumping phases give rise to a new and independent type of dynamical critical behavior of the rates of directional pumping, which are determined by the soft modes of the PTGGE. Published by the American Physical Society 2024
Quadratic Lindbladians encompass a rich class of dissipative electronic and bosonic quantum systems, which have been predicted to host new and exotic physics. In this study, we develop a Lindblad-Keldysh spectroscopic response formalism for open quantum systems that elucidates their steady-state response properties and dissipative phase transitions via finite-frequency linear and non-linear probes. As illustrative examples, we utilize this formalism to calculate the (1) density and dynamic spin susceptibilities of a boundary driven XY model at and near criticality, (2) linear and non-linear optical responses in Bernal bilayer graphene coupled to dissipative leads, and (3) steady state susceptibilities in a bosonic optical lattice. We find that the XY model spin density wavelength diverges with critical exponent 1/2, and there are gapless dispersive modes in the dynamic spin response that originate from the underlying spin density wave order; additionally the dispersing modes of the weak and ultra-strong dissipation limits exhibit a striking correspondence since the boundary dissipators couple only weakly to the bulk in both cases. In the optical response of the Bernal bilayer, we find that the diamagnetic response can decrease with increasing occupation, as opposed to in closed systems where the response increases monotonically with occupation; we study the effect of second harmonic generation and shift current and find that these responses, forbidden in centrosymmetric closed systems, can manifest in these open systems as a result of dissipation. We compare this formalism to its equilibrium counterpart and draw analogies between these non-interacting open systems and strongly interacting closed systems.