MASTER EQUATION FOR CORRELATORS OF NORMAL-ORDERED FIELD MODE OPERATORS Текст научной статьи по специальности «Физика»

NANOSYSTEMS: Vatutin A.D., et al. Nanosystems:

PHYSICS, CHEMISTRY, MATHEMATICS Phys. Chem. Math., 2022,13 (6), 628-631.

http://nanojournal.ifmo.ru

Original article DOI 10.17586/2220-8054-2022-13-6-628-631

Master equation for correlators of normal-ordered field mode operators

Alexander D. Vatutin", George P. Miroshnichenko6, Alexander I. Trifanovc ITMO University, Kronverkskiy, 49, St. Petersburg, 197101, Russia

° advatutin@yandex.ru, bgpmirosh@gmail.com, calextrifanov@gmail.com Corresponding author: Alexander D. Vatutin, advatutin@yandex.ru

PACS 03.65.Yz, 03.65.Ca

Abstract We study the master equation for open quantum systems in the alternative form, preserving the

normal form of the averaged normal-ordered operators. We give an example of using this equation for the

correlators of normal-ordered field mode operators. We explore the properties of the system of linear equations

for the higher-order field operators based on the example of a two-mode bosonic system.

Keywords master equation, normal-ordered correlators, two-mode bosonic system

For citation Vatutin A.D., Miroshnichenko G.P., Trifanov A.I. Master equation for correlators of normal-

ordered field mode operators. Nanosystems: Phys. Chem. Math., 2022,13 (6), 628-631.

1. Introduction

Quantum information processing is transfer, storage, and conversion of quantum information. It is essential for quantum protocols like quantum teleportation [1], quantum computing [2], quantum key distribution [3], dense coding, and quantum memory [4]. There are plenty of approaches to various problems for these protocols, but there is no general model that describes all kinds of them.

Models of open quantum systems are intended for solving problems of transport and storage of quantum information. This approach can be described in terms of completely positive trace-preserving (CPTP), or trace-non-increasing linear mappings known as quantum channels. The quantum master equation is the most effective way to determine the dynamics of open systems. More precisely, we use the Lindblad-type master equation [5, 6]. This equation was studied using physical methods [7-9] and mathematical techniques for the case of the single-mode Lindblad equation [6,10-17]. In the case of multi-mode bosonic systems, the mathematical approach was used to construct the Fock-like eigenstates of the Lindblad superoperators using the Lie algebras [18].

This approach simplifies the determination of the dynamics of the averaged moments of the Stokes operator and the polarization of light [19]. These parameters are sufficient for a fairly accurate description of many physical processes. However, these averaged moments of the Stokes operator are inadequate in cases when other statistical parameters than the wave average intensity are recorded [20]. Hence, the investigation of the higher-order moment correlators is required for a more detailed description of these quantum processes. We offer a more convenient representation of the Lindblad equation for the calculation of higher-order moments. The equation saves the normal form of the operators: if the original averaged operator was represented in the normal form, then each term in the equation has the form of the average of the operator in the normal form.

2. Model

The starting point is the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) form of the master equation for polarized light in the two-mode bosonic system in optical fibers:

dp = dt

- i ^ p] - ^ ^rn,m(^(ut + l)(a\lamp + pa,lam - 2ampa)n) + nT(ama)np + pama)n - 2alpam)j ,

n,m n,m

(1)

where indexes n,m € {1, 2}, dagger denotes the Hermitian conjuration, p is the density matrix of a quantum state, a,n and an are the creation and the annihilation operators of the n-th mode, Q and r are the frequency and the relaxation matrices [21], and nT is the mean number of thermal photons:

1

ut ,

' HQ o \ i

kBT 1 - 1

© Vatutin A.D., Miroshnichenko G.P., Trifanov A.I., 2022

where Qo is the free-space frequency, h is the reduced Planck constant, T is the environment temperature. kB is the Boltzmann constant. The frequency (Q) and the relaxation (r) matrices are Hermitian: Q = Q^, r = T. The positivity condition for the relaxation matrix has the following (conventional) form: z*Tz > 0,z e C2.

3. Alternative form for the Lindblad equation

Let us apply some operator A to the both parts of equation (1) and calculate the trace of the obtained operators. Note that the average A of the operator A is given by the following formula A = (a!) = Tr(Ap). Keeping in mind that the trace is invariant under cyclic permutations, we obtain the following relation:

^ = Tr ( - 2 2 Qn,m(Âalamp - a)namAp)-

- 2^2 rn,m(nT + V)(Aa)namp + a,[amAp - 2a!nAamp)-

n,m

1

- 2 2 rn,mnT (Aamalp + amalAp - 7amAalp)^ . (2)

Our goal is to bring the equation to such a form that if the operator A is in the normal form then the equation contains operators in the normal form only. Consider each of the sums in (2), starting with the first one:

7

n,m

Tr( - 2 2 ün,m(Aalamp - alamAp)^ =

Tr( - 2 2 ^n,m(Aa)nam - a)nAam + alAam - alamA)f^ =

[A, al]am - iü

n,man [am, A])p). (3)

2

n,m

As for the second sum, one subtracts and adds a)namAp to obtain commutators: 1

- 2 2 rn,m(nT + 1)(Âanamp + anamAp - 7alJiam^ =

n,m

- 2 2(nT + 1)rn,^([A,al]am + al[am,A(4)

Let us consider the third sum in more details following the procedure used above: 1

Tr( - 2 2 rn,mnT(Aamanp + amanÂp - 7amAanp)j =

n,m

Tr( - 2 2 rn,mnT(Aama,n - ama,nA + 2[am, [an, A]] + , A]am)^j =

n,m

Tr( - 2 rn,mnT([am, [a^A^p - 2 2 rn,mnT(Aa^m - altamA + A]am)^j =

n,m n,m

Tr( rn,mnT[[A, ai],am]p

)+Tr{2J2 rn nT ([Â,an\am + an[am,A])pj. (5)

Substituting (3), (4) and (5) in (2), we obtain the following relation for averaging of any operator A:

dA ~3t

2

(([A,an]am^ (-iün,m - rn,m) + ^ an[a,m,A]^ (iün,m - rn,m))--nt 2 rn,m{[ [[A, an], am]^ . (6)

n,m

One can note that this form of the master equation saves the normal form of the operators. An example showing this property is presented below.

1

630

A. D. Vatutin, G. P. Miroshnichenko, A. I. Trifanov

4. Second normal-ordered moments

apaq

Let us assume the operator A = apaq has the form of a field operator. Hence, the commutators in the expression (6)

can be rewritten as follows:

[apaq, 4] = Sqnap, [[apaq' aniL am] 3qn^pmi

[amo apaq] &pmaq

and equation (6) takes the form:

^ dt ^ = 2 ^ {a\am) ( iQq,m - rq,m) + ^ ^ (a^q) (iQn,p - ^,p) + U^qp. (7)

m n

Note that all operations in (6) are very simple. All transformations are reduced to counting of three simple commutators. The commutators in the first and the second sums preserve the order of the operators, and the double commutator in the third sum lowers the order by 2. This property takes place for any field operator. Furthermore, it is valid for any operator, if we expand it in respect to the basis of the field operators.

Let the set of operators Ai(Bi) consist of the normal-ordered operators. Let any operator Ai be a product of a operators of annihilation and fi operators of creations in any possible variants of modes; any operator Bi is a product of a - 1 operators of annihilation and fi - 1 operators of creations in any possible variants of modes. Then the system of differential equations for the set Jii derived from equation (6) gives one the expression containing two sets of operators A-i and Bi. This means that the equations of higher orders depend on the lower ones recursively. Example (7) shows this property in the simplest case (a = fi =1).

We can mention the amount of non-recurring expressions in the system. Only an amount of possible variations of operators of the creation (annihilation) in the first mode affects the quantity of all variants, because of [ap, a^q] = [ap, aq] = 0. Hence, the amount of equations is equal to (a + 1)(fi +1).

5. Conclusion

In this paper, we presented the equation for the normal-ordered moments and its property of saving the normal order. We started with the master equation in the GKSL form, converted it to the convenient form for the normal-ordered moments (6), and gave the example for the two-mode bosonic system. This is a logical continuation for the description of open quantum systems. We will try to continue the idea in terms of hidden polarization of light.

References

[1] Pirandola S., Eisert J., Weedbrook C., Furusawa A., Braunstein S.L. Advances in quantum teleportation. Nat. Photon., 2015, 9, P. 641-652.

[2] Gyongyosi L., Imre S. A survey on quantum computing technology. Comput. Sci. Rev., 2019, 31, P. 51-71 .

[3] Pirandola S., Andersen U., Banchi L., Berta M., Bunandar D., Colbeck R., Englund D., Gehring T., Lupo C., Ottaviani C., Pereira J. Advances in quantum cryptography. Adv. Opt. Photon., 2020,12, P. 1012-1236.

[4] Lvovsky A.I., Sanders B.C., Tittel W. Optical quantum memory. Nat. Photon., 2009, 3, P. 706-714.

[5] Gaidash A., Kozubov A., Miroshnichenko G. Dissipative dynamics of quantum states in the fiber channel. Phys. Rev. A, 2020, 102, 023711.

[6] Bonetti J., Hernandez S.M., Grosz D.F. Master equation approach to propagation in nonlinear fibers. Opt. Lett., 2021, 46, P. 665-668.

[7] Pearle P. Simple derivation of the Lindblad equation. Eur. J. Phys., 2012, 33, P. 805-822.

[8] Albash T., Boixo S., Lidar D.A., Zanardi P. Quantum adiabatic Markovian master equations. New J. Phys., 2012,14, 123016.

[9] 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.

[10] Miroshnichenko G. Decoherence of a one-photon packet in an imperfect optical fiber. Bull. Russ. Acad. Sci., 2018, 82, P. 1550-1555.

[11] Kozubov A., Gaidash A., Miroshnichenko G. Quantum model of decoherence in the polarization domain for the fiber channel. Phys. Rev. A, 2019, 99, 053842.

[12] Miroshnichenko G.P. Hamiltonian of photons in a single-mode optical fiber for quantum communications protocols. Opt. Spectrosc., 2012, 112, P. 777-786.

[13] 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.

[14] Tay B.A., Petrosky T. Biorthonormal eigenbasis of a Markovian master equation for the quantum Brownian motion. J. Math. Phys., 2008, 49, 113301.

[15] 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.

[16] Tay B.A. Eigenvalues of the Liouvillians of quantum master equation for a harmonic oscillator. Physica A, 2020, 556, 124768.

[17] Shishkov V.Y., Andrianov E.S., Pukhov A.A., Vinogradov A.P., Lisyansky A.A. Perturbation theory for Lindblad superoperators for interacting open quantum systems. Phys. Rev. A, 2020, 102, 032207.

[18] Teuber L., Scheel S. Solving the quantum master equation of coupled harmonic oscillators with Lie-algebra methods. Phys. Rev. A, 2020, 101, 042124.

[19] Gaidash A., Kozubov A., Miroshnichenko G., Kiselev A.D. Quantum dynamics of mixed polarization states: effects of environment-mediated intermode coupling. JOSA B, 2021, 38 (9), 425226.

[20] Klyshko D.M. Polarization of light: fourth-order effects and polarization-squeezed states. Zh. Eksp. Teor. Fiz., 1997,111, P. 1955-1983.

[21] Miroshnichenko G.P. Parameterization of an interaction operator of optical modes in a single-mode optical fiber. Nanosystems: Physics, Chemistry, Mathematics, 2015, 6 (6), P. 857-865.

Submitted 24 October 2022; revised 31 October 2022; accepted 13 November 2022

Information about the authors:

Alexander D. Vatutin - ITMO University, Kronverkskiy, 49, St. Petersburg, 197101, Russia; ORCID 0000-0003-1853-2171; advatutin@yandex.ru,

George P. Miroshnichenko - ITMO University, Kronverkskiy, 49, St. Petersburg, 197101, Russia; ORCID 0000-0002-4265-8818; gpmirosh@gmail.com

Alexander I. Trifanov - ITMO University, Kronverkskiy, 49, St. Petersburg, 197101, Russia; ORCID 0000-0002-5109-2086; alextrifanov@gmail.com

Conflict of interest: the authors declare no conflict of interest.