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:2610.05454

Game of Nons: Non-Gaussian Nonunitary Environments with Nonlinear Coupling

Artem Karasev, Alexander Teretenkov

quant-phSubmitted Oct 4, 2026• Announced Oct 6, 2026• 0 citations
arXivHTMLPDF

Abstract

We develop a perturbative framework for systems coupled quadratically to bosonic environments with quadratic Hamiltonians and linear Lindblad operators. The method applies to both single-mode and multimode environments whose initial states may be non-Gaussian and nonstationary. By combining linear evolution in a doubled operator space with a characteristic-function representation, we express ordered multitime environmental correlation functions in terms of derivatives of the initial characteristic function. These correlation functions enter a recursive construction of the time-convolutionless generator at arbitrary finite order. We illustrate the method with a qubit coupled to a damped mode via a two-photon interaction. Initial Fock and thermal states with identical first and second field moments produce different second-order generators through their fourth moments.

Related papers

Non-Perturbative Treatment of Open-System Multi-Time Expectation Values in Gaussian Bosonic Environments

Andrea Smirne, Dario Tamascelli, James Lim +2 more

We determine the conditions for the equivalence between the multi-time expectation values of a general finite-dimensional open quantum system when interacting with, respectively, an environment undergoing a free unitary evolution or a discrete environment under a free evolution fixed by a proper Gorini-Kossakowski-Lindblad-Sudarshan generator. We prove that the equivalence holds if both environments are bosonic and Gaussian and if the one- and two-time correlation functions of the corresponding interaction operators are the same at all times. This result leads to a non-perturbative evaluation of the multi-time expectation values of operators and maps of open quantum systems interacting with a continuous set of bosonic modes by means of a limited number of damped modes, thus setting the ground for the investigation of open-system multi-time quantities in fully general regimes.

arXiv:2209.00293•Open Systems & Information Dynamics•Dec 1, 2022•10 citations

Recursive perturbation approach to time-convolutionless master equations: Explicit construction of generalized Lindblad generators for arbitrary open systems

Alessandra Colla, Heinz‐Peter Breuer, Giulio Gasbarri

We develop a recursive perturbative expansion for the time-convolutionless generator of an open quantum system in a generalized Lindblad form. This formulation provides a systematic approach to derive the generator at arbitrary order while preserving a Lindblad-like structure, without imposing assumptions on the system or environment beyond an initially uncorrelated state. The generator is written, at all orders, in a canonical form, which also corresponds to the minimal dissipation condition, which uniquely specifies the decomposition of the generator into Hamiltonian and dissipative contributions. To validate the method and show its effectiveness in addressing non-Markovian dynamics and strong-coupling effects, we compute the generator explicitly up to fourth order.

arXiv:2506.04095•Physical Review A•Nov 24, 2025•3 citations

On Non-Gaussian Linear Quantum Optics

Alexander Evgen’evich Teretenkov

Abstract This paper studies a class of single-mode bosonic models of open quantum systems in which the dynamics of the creation and annihilation operators remains linear despite the strongly non-Gaussian nature of the induced dissipation. It is shown that, in contrast to isolated bosonic systems—where linear dynamics uniquely corresponds to quadratic Hamiltonians and Gaussian states—linearity of the operator dynamics in open systems does not guarantee Gaussianity of the state evolution. We construct quantum channels and Gorini–Kossakowski–Sudarshan–Lindblad generators for which the dynamics of the creation and annihilation operators is linear and closed, while Gaussian states are driven beyond the convex hull of Gaussian states during the evolution. In particular, the emergence of negativity of the Wigner function at a finite evolution time is demonstrated.

Bulletin of the Russian Academy of Sciences Physics•Jul 1, 2026•1 citation

Nonperturbative Treatment of non-Markovian Dynamics of Open Quantum Systems

Dario Tamascelli, Andrea Smirne, S. F. Huelga +1 more

We identify the conditions that guarantee equivalence of the reduced dynamics of an open quantum system (OQS) for two different types of environments-one a continuous bosonic environment leading to a unitary system-environment evolution and the other a discrete-mode bosonic environment resulting in a system-mode (nonunitary) Lindbladian evolution. Assuming initial Gaussian states for the environments, we prove that the two OQS dynamics are equivalent if both the expectation values and two-time correlation functions of the environmental interaction operators are the same at all times for the two configurations. Since the numerical and analytical description of a discrete-mode environment undergoing a Lindbladian evolution is significantly more efficient than that of a continuous bosonic environment in a unitary evolution, our result represents a powerful, nonperturbative tool to describe complex and possibly highly non-Markovian dynamics. As a special application, we recover and generalize the well-known pseudomodes approach to open-system dynamics.

arXiv:1709.03509•Physical Review Letters•Jan 19, 2018•243 citations

Non-Markovian qubit dynamics in nonequilibrium environments

Xiang-Ji Cai

We theoretically study the non-Markovian dynamics of qubit systems coupled to nonequilibrium environments with nonstationary and non-Markovian statistical properties. The reduced density matrix of the single qubit system satisfies a closed third-order differential equation with all the higher-order environmental correlations taken into account and the reduced density matrix of the two qubit system can be expressed as the Kraus representation in terms of the tensor products of the single qubit Kraus operators. We derive the relation between the entanglement and nonlocality of the two qubit system which are both closely associated with the decoherence function. We identify the threshold values of the decoherence function to ensure the existences of the concurrence and nonlocal quantum correlations for a given evolution time when the two qubit system is initially prepared in the composite Bell states and the extended Werner states, respectively. It is shown that the environmental nonstationary feature can suppress the decoherence and disentanglement dynamics and can reduce the coherence and entanglement revivals in non-Markovian dynamics regime. In addition, it is shown that the environmental non-Markovian feature can enhance the coherence revivals in the single decoherence dynamics and the entanglement revivals in the two qubit disentanglement dynamics, respectively. Furthermore, the environmental nonstationary and non-Markovian features can enhance the nonlocality of the two qubit system.

arXiv:2008.00640•Aug 3, 2020•0 citations

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

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

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.

arXiv:2412.13890•Dec 18, 2024•0 citations
We develop a perturbative framework for systems coupled quadratically to bosonic environments with quadratic Hamiltonians and linear Lindblad operators. The method applies to both single-mode and multimode environments whose initial states may be non-Gaussian and nonstationary. By combining linear evolution in a doubled operator space with a characteristic-function representation, we express ordered multitime environmental correlation functions in terms of derivatives of the initial characteristic function. These correlation functions enter a recursive construction of the time-convolutionless generator at arbitrary finite order. We illustrate the method with a qubit coupled to a damped mode via a two-photon interaction. Initial Fock and thermal states with identical first and second field moments produce different second-order generators through their fourth moments.
We determine the conditions for the equivalence between the multi-time expectation values of a general finite-dimensional open quantum system when interacting with, respectively, an environment undergoing a free unitary evolution or a discrete environment under a free evolution fixed by a proper Gorini-Kossakowski-Lindblad-Sudarshan generator. We prove that the equivalence holds if both environments are bosonic and Gaussian and if the one- and two-time correlation functions of the corresponding interaction operators are the same at all times. This result leads to a non-perturbative evaluation of the multi-time expectation values of operators and maps of open quantum systems interacting with a continuous set of bosonic modes by means of a limited number of damped modes, thus setting the ground for the investigation of open-system multi-time quantities in fully general regimes.
We develop a recursive perturbative expansion for the time-convolutionless generator of an open quantum system in a generalized Lindblad form. This formulation provides a systematic approach to derive the generator at arbitrary order while preserving a Lindblad-like structure, without imposing assumptions on the system or environment beyond an initially uncorrelated state. The generator is written, at all orders, in a canonical form, which also corresponds to the minimal dissipation condition, which uniquely specifies the decomposition of the generator into Hamiltonian and dissipative contributions. To validate the method and show its effectiveness in addressing non-Markovian dynamics and strong-coupling effects, we compute the generator explicitly up to fourth order.
Abstract This paper studies a class of single-mode bosonic models of open quantum systems in which the dynamics of the creation and annihilation operators remains linear despite the strongly non-Gaussian nature of the induced dissipation. It is shown that, in contrast to isolated bosonic systems—where linear dynamics uniquely corresponds to quadratic Hamiltonians and Gaussian states—linearity of the operator dynamics in open systems does not guarantee Gaussianity of the state evolution. We construct quantum channels and Gorini–Kossakowski–Sudarshan–Lindblad generators for which the dynamics of the creation and annihilation operators is linear and closed, while Gaussian states are driven beyond the convex hull of Gaussian states during the evolution. In particular, the emergence of negativity of the Wigner function at a finite evolution time is demonstrated.
We identify the conditions that guarantee equivalence of the reduced dynamics of an open quantum system (OQS) for two different types of environments-one a continuous bosonic environment leading to a unitary system-environment evolution and the other a discrete-mode bosonic environment resulting in a system-mode (nonunitary) Lindbladian evolution. Assuming initial Gaussian states for the environments, we prove that the two OQS dynamics are equivalent if both the expectation values and two-time correlation functions of the environmental interaction operators are the same at all times for the two configurations. Since the numerical and analytical description of a discrete-mode environment undergoing a Lindbladian evolution is significantly more efficient than that of a continuous bosonic environment in a unitary evolution, our result represents a powerful, nonperturbative tool to describe complex and possibly highly non-Markovian dynamics. As a special application, we recover and generalize the well-known pseudomodes approach to open-system dynamics.
We theoretically study the non-Markovian dynamics of qubit systems coupled to nonequilibrium environments with nonstationary and non-Markovian statistical properties. The reduced density matrix of the single qubit system satisfies a closed third-order differential equation with all the higher-order environmental correlations taken into account and the reduced density matrix of the two qubit system can be expressed as the Kraus representation in terms of the tensor products of the single qubit Kraus operators. We derive the relation between the entanglement and nonlocality of the two qubit system which are both closely associated with the decoherence function. We identify the threshold values of the decoherence function to ensure the existences of the concurrence and nonlocal quantum correlations for a given evolution time when the two qubit system is initially prepared in the composite Bell states and the extended Werner states, respectively. It is shown that the environmental nonstationary feature can suppress the decoherence and disentanglement dynamics and can reduce the coherence and entanglement revivals in non-Markovian dynamics regime. In addition, it is shown that the environmental non-Markovian feature can enhance the coherence revivals in the single decoherence dynamics and the entanglement revivals in the two qubit disentanglement dynamics, respectively. Furthermore, the environmental nonstationary and non-Markovian features can enhance the nonlocality of the two qubit system.
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.