Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2412.13890

Lindblad dynamics of open multi-mode bosonic systems: Algebra of bilinear superoperators, exceptional points and speed of evolution

Andrei Gaidash, Alexei D. Kiselev, Козубов Антон Владимирович, George P. Miroshnichenko

Published Dec 18, 2024• 0 citations
arXivHTMLPDF

Abstract

We develop the algebraic method based on the Lie algebra of quadratic combinations of left and right superoperators associated with matrices to study the Lindblad dynamics of multimode bosonic systems coupled a thermal bath and described by the Liouvillian superoperator that takes into account both dynamical (coherent) and environment mediated (incoherent) interactions between the modes. Our algebraic technique is applied to transform the Liouvillian into the diagonalized form by eliminating jump superoperators and solve the spectral problem. The temperature independent effective non-Hermitian Hamiltonian, $\hat{H}_{eff}$, is found to govern both the diagonalized Liouvillian and the spectral properties. It is shown that the Liouvillian exceptional points are represented by the points in the parameter space where the matrix, $H$, associated with $\hat{H}_{eff}$ is non-diagonalizable. We use our method to derive the low-temperature approximation for the superpropagator and to study the special case of a two mode system representing the photonic polarization modes. For this system, we describe the geometry of exceptional points in the space of frequency and relaxation vectors parameterizing the intermode couplings and, for a single-photon state, evaluate the time dependence of the speed of evolution as a function of the angles characterizing the couplings and the initial state.

DOI
10.48550/arxiv.2412.13890

Related papers

Algebraic approach for investigation of a multi-mode quantum system dynamics

Andrei Gaidash, Козубов Антон Владимирович, Alexei D. Kiselev +1 more

We introduce algebraic approach for superoperators that might be useful tool for investigation of quantum (bosonic) multi-mode systems and its dynamics. In order to demonstrate potential of proposed method we consider multi-mode Liouvillian superoperator that describes relaxation dynamics of a quantum system (including thermalization and intermode coupling). Considered algebraic structure of superoperators that form Liouvillian and their algebraic properties allows us to diagonilize multi-mode Liouvillian to find its spectrum. Also it allows to derive linear by mean number of thermal (environmental) photons approximation for time-evolution superoperator that keeps amount of considered dimensions in Fock space finite (assuming initial amount of dimensions finite) that might be helpful regarding entanglement dynamics problems. Conjugate Liouvillian is considered as well in order to perform analysis in Heisenberg picture, it can be implemented for multi-time correlation functions derivation.

arXiv:2207.01383•Jul 4, 2022•0 citations

Speed of Evolution and Correlations in Multi-Mode Bosonic Systems

Alexei D. Kiselev, Ranim Ali, Andrei V. Rybin

We employ an exact solution of the thermal bath Lindblad master equation with the Liouvillian superoperator that takes into account both dynamic and environment-induced intermode couplings to study the speed of evolution and quantum speed limit (QSL) times of a open multi-mode bosonic system. The time-dependent QSL times are defined from quantum speed limits, giving upper bounds on the rate of change of two different measures of distinguishability: the fidelity of evolution and the Hilbert-Schmidt distance. For Gaussian states, we derive explicit expressions for the evolution speed and the QSL times. General analytical results are applied to the special case of a two-mode system where the intermode couplings can be characterized by two intermode coupling vectors: the frequency vector and the relaxation rate vector. For the system initially prepared in a two-mode squeezed state, dynamical regimes are generally determined by the intermode coupling vectors, the squeezing parameter and temperature. When the vectors are parallel, different regimes may be associated with the disentanglement time, which is found to be an increasing (a decreasing) function of the length of the relaxation vector when the squeezing parameter is below (above) its temperature-dependent critical value. Alternatively, we study dynamical regimes related to the long-time asymptotic behavior of the QSL times, which is characterized by linear time dependence with the proportionality coefficients defined as the long-time asymptotic ratios. These coefficients are evaluated as a function of the squeezing parameter at varying temperatures and relaxation vector lengths. We also discuss how the magnitude and orientation of the intermode coupling vectors influence the maximum speed of evolution and dynamics of the entropy and the mutual information.

Entropy•Dec 5, 2022•2 citations

Liouvillian exceptional points of any order in dissipative linear bosonic systems: Coherence functions and switching between PT and anti- PT symmetries

Ievgen I. Arkhipov, Adam Miranowicz, Fabrizio Minganti +1 more

Usually, when investigating exceptional points (EPs) of an open Markovian bosonic system, one deals with spectral degeneracies of a non-Hermitian Hamiltonian (NHH), which can correctly describe the system dynamics only in the semiclassical regime. A recently proposed quantum Liouvillian framework [Minganti et al., Phys. Rev. A 100, 062131 (2019)] enables the complete determination of the dynamical properties of such systems and their EPs (referred to as Liouvillian EPs, or LEPs) in the quantum regime by taking into account the effects of quantum jumps, which are ignored in the NHH formalism. Moreover, the symmetry and eigenfrequency spectrum of the NHH become a part of much larger Liouvillian eigenspace. As such, the EPs of an NHH form a subspace of the LEPs. Here we show that once an NHH of a dissipative linear bosonic system exhibits an EP of a certain finite order $n$, it immediately implies that the corresponding LEP can become of any higher order $m\ensuremath{\ge}n$ defined in the infinite Hilbert space. These higher-order LEPs can be identified by the coherence and spectral functions at the steady state. The coherence functions can offer a convenient tool to probe extreme system sensitivity to external perturbations in the vicinity of higher-order LEPs. As an example, we study a linear bosonic system of a bimodal cavity with incoherent mode coupling to reveal its higher-order LEPs; particularly, of second and third order via first- and second-order coherence functions, respectively. Accordingly, these LEPs can be additionally revealed by squared and cubic Lorentzian spectral lineshapes in the power and intensity-fluctuation spectra. Moreover, we demonstrate that these EPs can also be associated with spontaneous parity-time ($\mathcal{PT}$) and anti-$\mathcal{PT}$-symmetry breaking in the system studied. These symmetries can be switched in the output fields (the so-called supermodes) of an additional linear coupler with a properly chosen coupling strength. Thus, we show that the initial loss-loss dynamics for the supermodes can be equivalent to the balanced gain-loss evolution.

arXiv:2006.03557•Physical Review A•Sep 11, 2020•73 citations

Third quantization for bosons: symplectic diagonalization, non-Hermitian Hamiltonian, and symmetries

Steven Kim, Fabian Hassler

Abstract Open quantum systems that interact with a Markovian environment can be described by a Lindblad master equation. The generator of time-translation is given by a Liouvillian superoperator acting on the density matrix of the system. As the Fock space for a single bosonic mode is already infinite-dimensional, the diagonalization of the Liouvillian has to be done on the creation- and annihilation-superoperators, a process called ‘third quantization’. We propose a method to solve the Liouvillian for quadratic systems using a single symplectic transformation. We show that the non-Hermitian effective Hamiltonian of the system, next to incorporating the dynamics of the system, is a tool to analyze its symmetries. As an example, we use the effective Hamiltonian to formulate a -‘symmetry’ of an open system. We describe how the inclusion of source terms allows us to obtain the cumulant generating function for observables such as the photon current.

arXiv:2304.02367•Journal of Physics A Mathematical and Theoretical•Aug 17, 2023•7 citations

Dynamics of interacting oscillators in the presence of environments

Juan Nicolas Moreno Tarquino

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

•Oct 23, 2024•0 citations

Lindblad dynamics of open multimode bosonic systems: Algebra of quadratic superoperators, exceptional points, and speed of evolution

Andrei Gaidash, Alexei D. Kiselev, Козубов Антон Владимирович +1 more

Physical Review A•Jun 2, 2025•4 citations
We develop the algebraic method based on the Lie algebra of quadratic combinations of left and right superoperators associated with matrices to study the Lindblad dynamics of multimode bosonic systems coupled a thermal bath and described by the Liouvillian superoperator that takes into account both dynamical (coherent) and environment mediated (incoherent) interactions between the modes. Our algebraic technique is applied to transform the Liouvillian into the diagonalized form by eliminating jump superoperators and solve the spectral problem. The temperature independent effective non-Hermitian Hamiltonian, H^eff\hat{H}_{eff}H^eff​, is found to govern both the diagonalized Liouvillian and the spectral properties. It is shown that the Liouvillian exceptional points are represented by the points in the parameter space where the matrix, HHH, associated with H^eff\hat{H}_{eff}H^eff​ is non-diagonalizable. We use our method to derive the low-temperature approximation for the superpropagator and to study the special case of a two mode system representing the photonic polarization modes. For this system, we describe the geometry of exceptional points in the space of frequency and relaxation vectors parameterizing the intermode couplings and, for a single-photon state, evaluate the time dependence of the speed of evolution as a function of the angles characterizing the couplings and the initial state.
We introduce algebraic approach for superoperators that might be useful tool for investigation of quantum (bosonic) multi-mode systems and its dynamics. In order to demonstrate potential of proposed method we consider multi-mode Liouvillian superoperator that describes relaxation dynamics of a quantum system (including thermalization and intermode coupling). Considered algebraic structure of superoperators that form Liouvillian and their algebraic properties allows us to diagonilize multi-mode Liouvillian to find its spectrum. Also it allows to derive linear by mean number of thermal (environmental) photons approximation for time-evolution superoperator that keeps amount of considered dimensions in Fock space finite (assuming initial amount of dimensions finite) that might be helpful regarding entanglement dynamics problems. Conjugate Liouvillian is considered as well in order to perform analysis in Heisenberg picture, it can be implemented for multi-time correlation functions derivation.
We employ an exact solution of the thermal bath Lindblad master equation with the Liouvillian superoperator that takes into account both dynamic and environment-induced intermode couplings to study the speed of evolution and quantum speed limit (QSL) times of a open multi-mode bosonic system. The time-dependent QSL times are defined from quantum speed limits, giving upper bounds on the rate of change of two different measures of distinguishability: the fidelity of evolution and the Hilbert-Schmidt distance. For Gaussian states, we derive explicit expressions for the evolution speed and the QSL times. General analytical results are applied to the special case of a two-mode system where the intermode couplings can be characterized by two intermode coupling vectors: the frequency vector and the relaxation rate vector. For the system initially prepared in a two-mode squeezed state, dynamical regimes are generally determined by the intermode coupling vectors, the squeezing parameter and temperature. When the vectors are parallel, different regimes may be associated with the disentanglement time, which is found to be an increasing (a decreasing) function of the length of the relaxation vector when the squeezing parameter is below (above) its temperature-dependent critical value. Alternatively, we study dynamical regimes related to the long-time asymptotic behavior of the QSL times, which is characterized by linear time dependence with the proportionality coefficients defined as the long-time asymptotic ratios. These coefficients are evaluated as a function of the squeezing parameter at varying temperatures and relaxation vector lengths. We also discuss how the magnitude and orientation of the intermode coupling vectors influence the maximum speed of evolution and dynamics of the entropy and the mutual information.
Usually, when investigating exceptional points (EPs) of an open Markovian bosonic system, one deals with spectral degeneracies of a non-Hermitian Hamiltonian (NHH), which can correctly describe the system dynamics only in the semiclassical regime. A recently proposed quantum Liouvillian framework [Minganti et al., Phys. Rev. A 100, 062131 (2019)] enables the complete determination of the dynamical properties of such systems and their EPs (referred to as Liouvillian EPs, or LEPs) in the quantum regime by taking into account the effects of quantum jumps, which are ignored in the NHH formalism. Moreover, the symmetry and eigenfrequency spectrum of the NHH become a part of much larger Liouvillian eigenspace. As such, the EPs of an NHH form a subspace of the LEPs. Here we show that once an NHH of a dissipative linear bosonic system exhibits an EP of a certain finite order nnn, it immediately implies that the corresponding LEP can become of any higher order $m\ensuremath{\ge}n$ defined in the infinite Hilbert space. These higher-order LEPs can be identified by the coherence and spectral functions at the steady state. The coherence functions can offer a convenient tool to probe extreme system sensitivity to external perturbations in the vicinity of higher-order LEPs. As an example, we study a linear bosonic system of a bimodal cavity with incoherent mode coupling to reveal its higher-order LEPs; particularly, of second and third order via first- and second-order coherence functions, respectively. Accordingly, these LEPs can be additionally revealed by squared and cubic Lorentzian spectral lineshapes in the power and intensity-fluctuation spectra. Moreover, we demonstrate that these EPs can also be associated with spontaneous parity-time (PT\mathcal{PT}PT) and anti-PT\mathcal{PT}PT-symmetry breaking in the system studied. These symmetries can be switched in the output fields (the so-called supermodes) of an additional linear coupler with a properly chosen coupling strength. Thus, we show that the initial loss-loss dynamics for the supermodes can be equivalent to the balanced gain-loss evolution.
Abstract Open quantum systems that interact with a Markovian environment can be described by a Lindblad master equation. The generator of time-translation is given by a Liouvillian superoperator acting on the density matrix of the system. As the Fock space for a single bosonic mode is already infinite-dimensional, the diagonalization of the Liouvillian has to be done on the creation- and annihilation-superoperators, a process called ‘third quantization’. We propose a method to solve the Liouvillian for quadratic systems using a single symplectic transformation. We show that the non-Hermitian effective Hamiltonian of the system, next to incorporating the dynamics of the system, is a tool to analyze its symmetries. As an example, we use the effective Hamiltonian to formulate a -‘symmetry’ of an open system. We describe how the inclusion of source terms allows us to obtain the cumulant generating function for observables such as the photon current.
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