Next Article in Journal
Starlikeness of New General Differential Operators Associated with q-Bessel Functions
Previous Article in Journal
UAV-Assisted Three-Dimensional Spectrum Mapping Driven by Spectrum Data and Channel Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Dynamics of Characteristic and One-Point Correlation Functions of Multi-Mode Bosonic Systems: Exactly Solvable Model

1
Laboratory of Quantum Processes and Measurements, ITMO University, Kadetskaya Line 3b, 199034 Saint Petersburg, Russia
2
Faculty of Photonics, ITMO University, Kadetskaya Line 3b, 199034 Saint Petersburg, Russia
3
Center of Information Optical Technology, ITMO University, Birzhevaya Line 14a, 199034 Saint Petersburg, Russia
*
Authors to whom correspondence should be addressed.
Symmetry 2021, 13(12), 2309; https://doi.org/10.3390/sym13122309
Submission received: 14 October 2021 / Revised: 23 November 2021 / Accepted: 1 December 2021 / Published: 3 December 2021

Abstract

:
In this communication we study dynamics of the open quantum bosonic system governed by the generalized Lindblad equation with both dynamical and environment induced intermode couplings taken into account. By using the method of characteristics we deduce the analytical expression for the normally ordered characteristic function. Analytical results for one-point correlation functions describing temporal evolution of the covariance matrix are obtained.

1. Introduction

Temporal evolution of open quantum systems is generally governed by completely positive and trace-preserving dynamical linear maps known as the quantum channels parameterized by time [1]. Under certain conditions, these maps can be described by master equations for the reduced density matrix derived using different assumptions and approximations [2,3,4,5,6]. It is well known that, within the Markov approximation, the master equations can be written in the so-called Lindblad form [7,8,9]. Such equation is also sometimes referred to as the Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) equation and its form preserves complete positivity of the dynamics.
Though the Lindblad-type master equations has been extensively discussed [10,11,12,13], the bulk of the mathematical methods developed for analysis of bosonic systems are applicable to the special case of single-mode Lindblad Equations [14,15,16,17,18,19]. These techniques cannot be employed to treat its multi-mode generalizations. In this paper, the most general form of the multi-mode Lindblad equation which is described, e.g., in [20] will be our primary concern.
The multi-mode bosonic systems play a leading role in quantum optics and quantum information with continuous variables [21,22]. In particular, there is a family of models that represents the most studied special case of two coupled oscillators [23,24,25,26]. Typically, in such models, it is assumed that relaxation does not involve interaction between the center-of-mass and relative-distance modes, so that the relaxation part of Lindbladian is a sum of two the commuting relaxation superoperators. In this communication we shall relax this assumption and study the generalized multi-mode Lindblad equation that, in particular, takes into account nontrivial effects induced by intermode couplings arising from the interaction between the system and the environment.
So we deal with a generalization of the model of a linear chain of coupled harmonic oscillators experiencing Markovian losses recently studied using Lie algebras [27]. Description of quantum dynamics of such quantum systems is important for understanding the quantum effects that occur in non-Hermitian systems which might pave the way to explore and apply quantum parity-time ( PT ) symmetry [28] (for recent reviews see, e.g., [29]) and quantum exceptional points [30,31].
The structure of this communication is as follows. In Section 2, after introducing generalized GKSL master equation for multi-mode bosonic systems interacting with the thermal bath we solve the dynamical equation for the normally ordered characteristic function χ N using the method of characteristics. The formula for χ N is then utilized to deduce the expressions for the one-point correlators describing dynamics of the second order moments that enter the covariance matrix. In Section 3 we discuss the results and make some concluding remarks.

2. Results

Our starting point is the Lindblad equation for the density matrix describing the Markovian quantum dynamics of the N-mode open bosonic systemis of the general form [20]:
ρ ^ t = L ρ ^ = n , m = 1 N i Ω n m C a ^ n a ^ m ρ ^ + Δ n m * C a ^ n a ^ m ρ ^ Δ n m C a ^ m a ^ n ρ ^ + n , m = 1 N K n m D a ^ m a ^ n ρ ^ + L n m D a ^ n a ^ m ρ ^ + M n m * D a ^ n a ^ m ρ ^ + M n m D a ^ m a ^ n ρ ^
written in terms of two superoperators given by
C A ^ B ^ : ρ ^ C A ^ B ^ ρ ^ = [ A ^ B ^ , ρ ^ ] ,
D A ^ B ^ : ρ ^ D A ^ B ^ ρ ^ = 2 A ^ ρ ^ B ^ B ^ A ^ ρ ^ ρ ^ B ^ A ^ = [ A ^ , ρ ^ B ^ ] + [ A ^ ρ ^ , B ^ ] ,
where the dagger denotes Hermitian conjugation, ρ ^ is the density matrix representing the quantum state; a ^ n ( a ^ n ) is the creation (annihilation) operator of the nth mode; [ A ^ , B ^ ] = A ^ B ^ B ^ A ^ stands for the commutator; Ω n m , K n m and L n m are the elements of the frequency and relaxation matrices, respectively. All these matrices are assumed to be Hermitian: Ω = Ω , K = K and L = L .
From Equation (1) combined with the algebraic identities
Tr ( S ^ C A ^ B ^ ρ ^ ) = [ S ^ , A ^ ] B ^ + A ^ [ S ^ , B ^ ] ,
Tr ( S ^ D A ^ B ^ ρ ^ ) = B ^ [ S ^ , A ^ ] [ S ^ , B ^ ] A ^
we have
χ N ( α ˜ , t ) t = L ^ χ N ( α ˜ , t ) ,
where the operator L ^ on the right hand side of Equation (6)
L ^ = n , m = 1 N i Ω n m D ^ m n ( ) + n , m = 1 N Δ n m ( C ^ m n ( + ) α n * α m * ) + Δ n m * ( [ C ^ m n ( + ) ] * α n α m ) n , m = 1 N ( K n m L n m ) D ^ m n ( + ) + 2 L n m α m α n * ) n , m = 1 N M n m ( C ^ m n ( ) + α n * α m * ) + M n m * ( [ C ^ m n ( ) ] * + α n α m )
is expressed in terms of the differential operators given by
D ^ m n ( ± ) = α m α n ± α n * α m * ,
C ^ m n ( ± ) = α m * α n ± α n * α m .
Temporal evolution of the characteristic function χ N is governed by the dynamical Equation (6) supplemented with the initial condition
χ N ( α ˜ , 0 ) χ ini ( α ˜ ) = d μ ( β ) e [ ( α * , β ) ( α , β * ) ] P ( β ˜ , 0 ) ,
d μ ( β ) = n = 1 N d μ ( β i ) , d μ ( β i ) = d 2 β i π ,
where
( α * , β ) i = 1 N α i * β i
and χ N is related to the Glauber–Sudarshan P function (quasidistribution), P ( β ˜ , t ) , that enters the P-representation of the density matrix:
ρ ^ ( t ) = d μ ( β ) P ( β ˜ , t ) β β .
We can now employ the method of characteristics [32] to solve the above initial value problem. According to this method, we begin with the system of characteristic equations
α n t = m = 1 N Q n m ( 11 ) α m m = 1 N Q n m ( 12 ) α m * ,
that can be conveniently put into the matrix form
α ˜ t t α α * = Q α α * , Q = Q 11 Q 12 Q 12 * Q 11 * ,
where the block structure of Q is defined by the submatrices given by
Q 11 = i Ω Γ , Q 12 = Δ s M a ,
Γ = K L , Δ s , a = Δ ± Δ T , M s , a = M ± M T .
Solution of the system (15) written in the matrix form as follows
α ˜ α α * = U ( t ) α 0 α 0 * ,
U ( t ) = e Q t = U 11 ( t ) U 12 ( t ) U 12 * ( t ) U 11 * ( t ) ,
where α 0 = α ( 0 ) . It is not difficult to obtain the solution along the characteristic curves (18) given by
χ N ( α ˜ 0 , t ) = exp 0 t ( α ˜ 0 * , G ( τ ) α ˜ 0 ) d τ χ ini ( α ˜ 0 ) ,
G ( t ) = U ( t ) G U ( t ) , G = L 1 2 ( Δ s M s ) 1 2 ( Δ s * M s * ) L * .
We can now express α ˜ 0 in terms of α ˜ with the help of Equation (18): α ˜ 0 = U ( t ) α ˜ to transform Formula (20) into the final expression for the characteristic function:
χ N ( α ˜ , t ) = e ( α ˜ * , B ( t ) α ˜ ) χ ini ( U ( t ) α ˜ ) ,
B ( t ) = B 11 ( t ) B 12 ( t ) B 12 * ( t ) B 11 * ( t ) = 0 t G ( t τ ) d τ .
We can now use formula for the characteristic function (22) to compute time dependence of averages that can be regarded as one-point correlation functions. Derivatives of this function can be easily evaluated giving the expressions for mean values of normally ordered operators. In particular, analytical expressions for second order moments can be readily derived in the following general form:
a ^ i a ^ j ( t ) = 2 χ N ( α , t ) α i α j * | α = 0 = 2 B j i ( 11 ) ( t ) n , m U n i ( 11 ) ( t ) U m j ( 12 ) ( t ) a ^ n a ^ m ( 0 ) n , m [ U n i ( 12 ) ( t ) U m j ( 11 ) ( t ) ] * a ^ n a ^ m ( 0 )
+ n , m U n i ( 11 ) ( t ) [ U m j ( 11 ) ] * ( t ) a ^ n a ^ m ( 0 ) + n , m [ U n i ( 12 ) ( t ) ] * U m j ( 12 ) ( t ) a ^ m a ^ n ( 0 ) , a ^ i a ^ j ( t ) = 2 χ N ( α , t ) α i * α j * | α = 0 = 2 B j i ( 12 ) ( t ) + n , m U n i ( 12 ) ( t ) U m j ( 12 ) ( t ) a ^ n a ^ m ( 0 ) + n , m [ U n i ( 11 ) ( t ) U m j ( 11 ) ( t ) ] * a ^ n a ^ m ( 0 )
n , m U n i ( 12 ) ( t ) [ U m j ( 11 ) ] * ( t ) a ^ n a ^ m ( 0 ) n , m [ U n i ( 11 ) ( t ) ] * U m j ( 12 ) ( t ) a ^ m a ^ n ( 0 ) .
Relations (24) and (25) govern time dependence of the covariance matrix for our bosonic system [22].

3. Conclusions

In this communication we have obtained the expression for the normally ordered characteristic function of the open multi-mode bosonic system governed by the Lindblad equation taken in the general form. This equation accounts for the effects of dynamical and environment induced intermode couplings and squeezing. The analytical results are applied to obtain formulas describing time dependence of the averages representing the one-point correlation functions that enter the covariance matrix.
Our technique being directly applicable to quantum dynamics of continuous variables [22] might also provide a useful method to investigate Lindblad dynamics of the systems of coupled oscillators with PT symmetry previously studied in both classical [33,34] and quantum [35,36] regimes. The analytical approach can be readily extended to study a number of problems such as quantum dynamics of mixed polarization states transmitted in lossy quantum communication channels [37,38,39] and the problems related to controlled quantum dynamics in a realistic setup involving open environments.
Our concluding remark concerns quantum navigation which is an important class of controlled quantum dynamics whereby the objective is to transport one quantum state into another, or to generate quantum gates, in the shortest possible time under the influence of an uncontrollable external field. Problems of this kind can be thought of as representing the quantum counterpart of the classical Zermelo navigation problem of finding the time-optimal control that takes a ship from one location to another, under the influence of external wind or currents [40,41]. In a forthcoming publication we will apply our results to the Zermelo navigation problem for open multimode bosonic systems.

Author Contributions

Research: A.D.K., R.A., A.V.R.; writing—original draft preparation: A.D.K., R.A.; writing—review and editing: A.D.K., A.V.R. All authors have read and agreed to the published version of the manuscript.

Funding

The work of A.D.K. and A.V.R. was supported by the Russian Science Foundation under grant No. 20-11-20226.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Caruso, F.; Giovannetti, V.; Lupo, C.; Mancini, S. Quantum channels and memory effects. Rev. Mod. Phys. 2014, 86, 1203. [Google Scholar] [CrossRef]
  2. Carmichael, H. An Open Systems Approach to Quantum Optics; Springer: Berlin/Heidelberg, Germany, 1993; p. 179. [Google Scholar]
  3. Breuer, H.-P.; Petruccione, F. The Theory of Open Quantum Systems; Oxford University Press: Oxford, UK, 2002; p. 625. [Google Scholar]
  4. Rivas, A.; Huelga, S.F. Open Quantum Systems: An Introduction; SpringerBriefs in Physics; Springer: Berlin/Heidelberg, Germany, 2012; p. 97. [Google Scholar]
  5. Rivas, Á.; Huelga, S.F.; Plenio, M.B. Quantum non-Markovianity: Characterization, quantification and detection. Rep. Prog. Phys. 2014, 77, 094001. [Google Scholar] [CrossRef] [PubMed]
  6. de Vega, I.; Alonso, D. Dynamics of non-Markovian open quantum systems. Rev. Mod. Phys. 2017, 89, 015001. [Google Scholar] [CrossRef] [Green Version]
  7. Kossakowski, A. On quantum statistical mechanics of non-Hamiltonian systems. Rep. Math. Phys. 1972, 3, 247–274. [Google Scholar] [CrossRef]
  8. Lindblad, G. On the generators of quantum dynamical semigroups. Commun. Math. Phys. 1976, 48, 119–130. [Google Scholar] [CrossRef]
  9. Gorini, V.; Kossakowski, A.; Sudarshan, E.C.G. Completely positive dynamical semigroups of N-level systems. J. Math. Phys. 1976, 17, 821. [Google Scholar] [CrossRef]
  10. Pearle, P. Simple derivation of the Lindblad equation. Eur. J. Phys. 2012, 33, 805–822. [Google Scholar] [CrossRef]
  11. Albash, T.; Boixo, S.; Lidar, D.A.; Zanardi, P. Quantum adiabatic Markovian master equations. New J. Phys. 2012, 14, 123016. [Google Scholar] [CrossRef]
  12. McCauley, G.; Cruikshank, B.; Bondar, D.I.; Jacobs, K. Accurate Lindblad-form master equation for weakly damped quantum systems across all regimes. npj Quantum Inf. 2020, 6, 74. [Google Scholar] [CrossRef]
  13. Manzano, D. A short introduction to the Lindblad master equation. AIP Adv. 2020, 10, 025106. [Google Scholar] [CrossRef]
  14. Arévalo-Aguilar, L.M.; Moya-Cessa, H. Solution to the master equation for a quantized cavity mode. Quantum Semiclassical Opt. J. Eur. Opt. Soc. Part B 1998, 10, 671–674. [Google Scholar] [CrossRef]
  15. Klimov, A.B.; Romero, J.L. An algebraic solution of Lindblad-type master equations. J. Opt. B Quantum Semiclassical Opt. 2003, 5, S316–S321. [Google Scholar] [CrossRef]
  16. Lu, H.-X.; Yang, J.; Zhang, Y.-D.; Chen, Z.-B. Algebraic approach to master equations with superoperator generators of su(1, 1) and su(2) Lie algebras. Phys. Rev. A 2003, 67, 024101. [Google Scholar] [CrossRef]
  17. Tay, B.A.; Petrosky, T. Biorthonormal eigenbasis of a Markovian master equation for the quantum Brownian motion. J. Math. Phys. 2008, 49, 113301. [Google Scholar] [CrossRef] [Green Version]
  18. Honda, D.; Nakazato, H.; Yoshida, M. Spectral resolution of the Liouvillian of the Lindblad master equation for a harmonic oscillator. J. Math. Phys. 2010, 51, 072107. [Google Scholar] [CrossRef] [Green Version]
  19. Tay, B. Eigenvalues of the Liouvillians of quantum master equation for a harmonic oscillator. Phys. A Stat. Mech. Its Appl. 2020, 556, 124768. [Google Scholar] [CrossRef]
  20. Benatti, F.; Floreanini, R. Entangling oscillators through environment noise. J. Phys. A Math. Gen. 2006, 39, 2689–2699. [Google Scholar] [CrossRef]
  21. Braunstein, S.L.; van Loock, P. Quantum information with continuous variables. Rev. Mod. Phys. 2005, 77, 513. [Google Scholar] [CrossRef] [Green Version]
  22. Serafini, A. Quantum Continuous Variables: A Primer of Theoretical Methods; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2017; p. 350. [Google Scholar]
  23. Hiroshima, T. Decoherence and entanglement in two-mode squeezed vacuum states. Phys. Rev. A 2001, 63, 022305. [Google Scholar] [CrossRef] [Green Version]
  24. Paz, J.P.; Roncaglia, A.J. Dynamics of the Entanglement between Two Oscillators in the Same Environment. Phys. Rev. Lett. 2008, 100, 220401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Linowski, T.; Gneiting, C.; Rudnicki, Ł. Stabilizing entanglement in two-mode Gaussian states. Phys. Rev. A 2020, 102, 042405. [Google Scholar] [CrossRef]
  26. Vendromin, C.; Dignam, M.M. Continuous-variable entanglement in a two-mode lossy cavity: An analytic solution. Phys. Rev. A 2021, 103, 022418. [Google Scholar] [CrossRef]
  27. Teuber, L.; Scheel, S. Solving the quantum master equation of coupled harmonic oscillators with Lie-algebra methods. Phys. Rev. A 2020, 101, 042124. [Google Scholar] [CrossRef]
  28. Bender, C.M.; Boettcher, S. Real spectra in non-hermitian hamiltonians having symmetry. Phys. Rev. Lett. 1998, 80, 5243. [Google Scholar] [CrossRef] [Green Version]
  29. Christodoulides, D.; Yang, J. (Eds.) Parity-Time Symmetry and Its Applications, Springer Tracts in Modern Physics; Springer: Singapore, 2018; Volume 280, p. 578. [Google Scholar]
  30. Heiss, W.D. The physics of exceptional points. J. Phys. A Math. Theor. 2012, 45, 444016. [Google Scholar] [CrossRef]
  31. Özdemir, Ş.K.; Rotter, S.; Nori, F.; Yang, L. Parity-time symmetry and exceptional points in photonics. Nat. Mater. 2019, 18, 783–798. [Google Scholar] [CrossRef]
  32. Shearer, M.; Levy, R. Partial Differential Equations: An Introduction to Theory and Applications; Princeton University Press: New Jersey, NJ, USA, 2015; p. 300. [Google Scholar]
  33. Bender, C.M.; Gianfreda, M.; Klevansky, S.P. Systems of coupledPT-symmetric oscillators. Phys. Rev. A 2014, 90, 022114. [Google Scholar] [CrossRef] [Green Version]
  34. Tsoy, E.N. Coupled oscillators with parity-time symmetry. Phys. Lett. A 2017, 381, 462–466. [Google Scholar] [CrossRef]
  35. Xu, X.-W.; Liu, Y.-X.; Sun, C.-P.; Li, Y. Mechanical symmetry in coupled optomechanical systems. Phys. Rev. A 2015, 92, 013852. [Google Scholar] [CrossRef] [Green Version]
  36. Li, W.; Li, C.; Song, H. Theoretical realization and application of parity-time-symmetric oscillators in a quantum regime. Phys. Rev. A 2017, 95, 023827. [Google Scholar] [CrossRef] [Green Version]
  37. Gaidash, A.; Kozubov, A.; Miroshnichenko, G. Dissipative dynamics of quantum states in the fiber channel. Phys. Rev. A 2020, 102, 023711. [Google Scholar] [CrossRef]
  38. Kozubov, A.V.; Gaidash, A.A.; Miroshnichenko, G.P.; Kiselev, A.D. Quantum dynamics of mixed polarization states: Effects of environment-mediated intermode coupling. J. Opt. Soc. Am. B 2021, 38, 2603. [Google Scholar]
  39. Kiselev, A.D.; Ali, R.; Rybin, A.V. Lindblad Dynamics and Disentanglement in Multi-Mode Bosonic Systems. Entropy 2021, 23, 1409. [Google Scholar] [CrossRef] [PubMed]
  40. Brody, D.C.; Gibbons, G.W.; Meier, D.M. Time-optimal navigation through quantum wind. New J. Phys. 2015, 17, 033048. [Google Scholar] [CrossRef]
  41. Brody, D.C.; Longstaff, B. Evolution speed of open quantum dynamics. Phys. Rev. Res. 2019, 1, 033127. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kiselev, A.D.; Ali, R.; Rybin, A.V. Dynamics of Characteristic and One-Point Correlation Functions of Multi-Mode Bosonic Systems: Exactly Solvable Model. Symmetry 2021, 13, 2309. https://doi.org/10.3390/sym13122309

AMA Style

Kiselev AD, Ali R, Rybin AV. Dynamics of Characteristic and One-Point Correlation Functions of Multi-Mode Bosonic Systems: Exactly Solvable Model. Symmetry. 2021; 13(12):2309. https://doi.org/10.3390/sym13122309

Chicago/Turabian Style

Kiselev, Alexei D., Ranim Ali, and Andrei V. Rybin. 2021. "Dynamics of Characteristic and One-Point Correlation Functions of Multi-Mode Bosonic Systems: Exactly Solvable Model" Symmetry 13, no. 12: 2309. https://doi.org/10.3390/sym13122309

APA Style

Kiselev, A. D., Ali, R., & Rybin, A. V. (2021). Dynamics of Characteristic and One-Point Correlation Functions of Multi-Mode Bosonic Systems: Exactly Solvable Model. Symmetry, 13(12), 2309. https://doi.org/10.3390/sym13122309

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop