[
Abstract
We study the bath dynamics in the dephasing model of a twostate quantum system (qubit) coupled to an environment of harmonic oscillators. This model was shown [Morozov et al., Phys. Rev. A, 2012 85, 022101] to admit the analytic solution for the qubit and environment dynamics. Using this solution, we derive the exact expression for the bath reduced density matrix in the presence of initial qubitenvironment correlations. We obtain the nonequilibrium phonon distribution function and discuss in detail the time behavior of the bath energy. It is shown that only the inclusion of dynamic correlations between the qubit and the bath ensures the proper time behavior of the quantity which may be interpreted as the ‘‘environment energy’’.
\keywordsopen quantum systems, qubitenvironment correlations \pacs03.65.Yz
Abstract
Проведено дослiдження динамiки термостату у випадку моделi з розфазуванням, що описує дворiвневу квантову систему (кубiт), яка взаємодiє з гармонiчними осциляторами зi свого оточення. Ця модель має аналiтичний розв’язок [Morozov et al., Phys. Rev. A, 2012 85, 022101] як для спiнових змiнних, так i для змiнних термостату. Використовуючи цей розв’язок, отримано аналiтичний вираз для приведеної матрицi густини термостату при наявностi початкових кореляцiй. Отримано нерiвноважну функцiю розподiлу фононiв та детально дослiджено часову еволюцiю енергiї термостату. Показано, що лише належне врахування динамiчних кореляцiй мiж кубiтом та його оточенням забезпечує збереження величини, яку слiд вважати ‘‘енергiєю оточення’’.
\keywordsквантовi вiдкритi системи, кореляцiї, кубiт
201316334001
\doinumber10.5488/CMP.16.34001
\articletypeRapid Communication
Bath dynamics]Bath dynamics in an exactly solvable qubit model with initial
qubitenvironment correlations
V.V. Ignatyuk, V.G. Morozov]V.V. Ignatyuk\refaddrlabel1,
V.G. Morozov\refaddrlabel2
\addresses
\addrlabel1 Institute for Condensed Matter Physics
of the National Academy of Sciences of Ukraine,
1 Svientsitskii St., 79011 Lviv, Ukraine\addrlabel2 Moscow State Technical University of Radioengineering, Electronics, and
Automation,
78 Vernadsky Prospect, 119454 Moscow, Russia
\authorcopyrightV.V. Ignatyuk, V.G. Morozov, 2013
1 Introduction
The dynamics of open quantum systems has attracted a great deal of interest over the last few decades. Its importance is dictated by the prospects of applications in quantum optics, quantum computing, quantum measurements and control [1, 2], as well as by the necessity of a deeper understanding of the theory itself [3, 4, 5, 6, 7, 8]. The dynamics of open quantum systems was studied in several aspects: (i) the effect of initial correlations between an open system and its environment has been investigated in [9, 10]; (ii) a new viewpoint concerning the nature and the measure of nonMarkovianity has been presented in [11, 12]; (iii) the effect of nonequilibrium environment on quantum coherence and the level populations has been considered in [13, 14].
The latter problem is of a particular interest. Usually, when constructing a master equation for the reduced density matrix of an open system, one considers the bath to be at thermal equilibrium [1], even though there is a buildup of dynamical correlations [4, 9] caused by the entanglement of quantum states. Moreover, even if the effect of nonequilibrium environment on the system behavior is not neglected [13], the intrinsic bath dynamics is beyond consideration, since the corresponding bath variables are always integrated out from the equations of motion. Although this approach seems to be quite natural as long as one studies solely the open system dynamics, the investigation of the bath evolution itself can undoubtedly be an interesting problem, yielding some useful hints about how to deal with more realistic systems, especially with those which possess slow relaxation to equilibrium and do not admit exact solutions.
In this brief report, we show the relevance of proper (or intrinsic) bath dynamics. The paper is structured as follows. In section 2 we derive an exact expression for the bath reduced density matrix in the socalled dephasing model, describing a twostate system (qubit) coupled to a bosonic bath [5, 6, 7]. In section 3 the nonequilibrium distribution function for the bath modes (phonons) is calculated. Special attention is paid to the description of time evolution of the phonon energy and the correlation energy in the ‘‘qubitbath’’ system. In the last section we discuss the results and draw final conclusions.
2 Bath density matrix in the dephasing model
We consider a simple version of the spinboson model describing a twostate system (qubit) () coupled to the bath () of harmonic oscillators [5, 6, 7, 8]. In the ‘‘spin’’ representation for the qubit, the total Hamiltonian of the model is written as follows (in our units )
(1) 
where is the energy difference between the excited and the ground states of the qubit, and is one of the Pauli matrices , , . Note that the operator can also be written in the basis of the ground and the excited states of the twolevel system as . Bosonic creation and annihilation operators and correspond to the th bath mode with frequency , and are the coupling constants.
The distinctive feature of the dephasing model (1) is that the average populations of the qubit states do not depend on time, and hence there is no relaxation to the complete equilibrium between the qubit and the environment. In other words, the model is nonergodic^{1}^{1}1Another special feature of the model (1) is that the qubit Hamiltonian and the sum are integrals of motion.. Note, however, that the Heisenberg picture operators evolve in time, leading to a nontrivial decay of the coherences . Thus, we have a unique situation, where the system relaxation may be interpreted physically as ‘‘pure’’ decoherence and the entropy exchange [8] rather than the energy dissipation.
It was shown in [7] that equations of motion for the Heisenberg picture operators , , and can be solved exactly with the results
(2) 
(3) 
where
(4) 
We shall use the exact expressions (2)–(3) to evaluate the reduced density matrix of the bath. We start with the obvious relation
(5) 
for the nonequilibrium density matrix of the composite system () and assume that the initial density matrix has the form [7]
(6) 
where the projector can be expressed in terms of the Bloch vector [7, 8]:
(7) 
where , , and denotes the identity matrix. In equation (6) with is the state vector of the qubit. The constrained initial density matrix of the bath, , is given by ^{2}^{2}2 The equation (8) can be derived by using obvious relations , and the following properties [7] of the Hamiltonian (1): , .
(8) 
where . The bath Hamiltonians and the corresponding partition functions are defined as
(9) 
Physically, the density matrix (6) corresponds to a situation where at times the open system is in thermal equilibrium with its environment , and at time zero one makes a perfect (selective) measurement on the system only. As a result [1], the system is prepared in some pure state .
Making use of equations (6)–(7) and noting that is an integral of motion, we can recast equation (5) into the form
(10) 
where we have introduced the notation
(11) 
The bath density matrix is obtained from (2) by taking the trace over the qubit states,
(12) 
Since does not contain the spin operators , the terms with do not contribute to (12). Expressing the probabilities in terms of in a usual way, , , it is straightforward to manipulate the timedependent bath density matrix (12) to
(13)  
where
(14) 
denotes a phonon part of in the interaction picture (with the replacement , see equation (5)), whereas the quantity
(15) 
is the nonequilibrium correction to the phonons energy (see also the next section for discussion).
The expression for the bath density matrix becomes much simpler if one neglects the initial correlations in the system by taking a direct product
(16) 
instead of (6). Proceeding in a similar way, after some algebra one obtains
(17) 
where denotes the nonequilibrium contribution to the correlation energy in the interaction picture. It is seen from equation (17) that even in this case there is a dynamical buildup of correlations in the system. We touch upon this point in the next section, when analyzing the relevance of nonequilibrium correlations.
3 Phonon nonequilibrium distribution function and energy
To gain some insight into the time behavior of the bath modes (phonons), let us first calculate the phonon distribution function . This can be done in two equivalent ways: either using the exact expression (13) for the bath density matrix, or averaging the bosonic Heisenberg picture operators (3) over the initial state of the composite system. Here, we shall follow the latter procedure which is simpler. Taking the initial density matrix in the form (6) and then applying the unitary transformation technique [7] (or the method of the displaced harmonic oscillator modes [15]), after some straightforward algebra one obtains
(18) 
where
(19) 
is the initial phonon distribution function, and the function
(20) 
represents the contribution of initial correlations.
With equations (18) and (19), it is easy to calculate the time evolution of the nonequilibrium phonon energy . The final result is conveniently written in terms of the bath spectral density which is introduced by the wellknown rule [1, 2, 7]
(21) 
After simple manipulations we arrive at
(22) 
where the initial phonon energy is given by
(23) 
Here the last term occurs due to initial correlations in the system.
Usually [1, 2, 7], the spectral density function is chosen in the form
(24) 
where and is a dimensionless coupling constant. This formula ensures both a proper lowfrequency behavior of and a cutoff at high frequencies (). The case is usually called the ‘‘Ohmic’’ case, the case ‘‘superOhmic’’, and the case ‘‘subOhmic’’. Using expression (24), it is possible to analyze the time behavior of for different , but in this brief report we would like to discuss only one physically interesting point related to the result (22) for the phonon energy.
At first glance, the fact that the phonon energy (22) depends on time may appear as an apparent paradox. Indeed, on the one hand, one may conclude that there is an energy exchange between the qubit and the bath. On the other hand, the qubit Hamiltonian commutes with the total Hamiltonian (1), and hence the qubit energy does not depend on time. To explain this paradox, let us calculate the nonequilibrium correlation energy , where is the interaction term in equation (1) (in the Heisenberg picture) and the average is taken over the initial state (6). After some algebra, which we omit, we obtain
(25) 
Combining this expression with equation (22), it is easy to check that the sum is a timeindependent quantity (see also a footnote on page 2). We see that the dynamics of the correlation energy exactly compensates the time dependence of the nonequilibrium phonon energy , ensuring the energy conservation law. Physically, the sum is precisely the quantity which should be interpreted as the environment energy.
One more remark is to the point. It can be seen from equations (20), (22), and (25) that both the nonequilibrium phonon energy and the correlation energy do not depend on time under conditions . Note in this connection that the correlational contribution to the generalized decoherence function [7] vanishes for the same values of the mean inversion population of the levels, manifesting a close relationship between the essentially nonequilibrium behavior of the correlation energy and the onset of the additional channel of decoherence in the system.
4 Conclusions
Here, we present a summary of the results and discuss their relation to some problems in the dynamics of open quantum systems.
We have derived exact expressions (13) and (17) for the bath density matrix in the model (1) which describes the dephasing mechanism of decoherence in a qubit interacting with a bosonic environment. To the best of our knowledge, a derivation of a bath density matrix has never been performed in the theory of open quantum systems. The explicit form of could be essential, for instance, when constructing master equations (especially nonMarkovian) and taking into account the intrinsic dynamics of the environment along with the equation of motion for . Such an approach would modify the wellknown ZwanzigNakajima projection technique [1, 3, 9] where the bath degrees of freedom are ‘‘eliminated’’. We believe that this modification is quite natural in the case of the finite size of the bath, when all the environmental modes are involved in the composite system dynamics, and a backflow of energy (information) from the bath to the open system is essential. Thus, the exact solutions (13) and (17) can give a valuable insight into general properties of the dynamics of decoherence and can serve as a step toward consistent derivation of master equations ensuring regular behavior of composite systems on all timescales and for strong coupling regimes.
Our analysis of the phonon energy in section 3 illustrates the special role of dynamic correlations between an open system and its environment. We have seen that a ‘‘naive’’ picture with the ‘‘energy exchange between the qubit and the bath’’ is inadequate (even in the case of weak coupling), and only the proper inclusion of nonequilibrium correlations ensures the conservation of the total energy.
References
 [1] Breuer H.P., Petruccione F., The Theory of Open Quantum Systems, Oxford University, Oxford, 2007.
 [2] Weiss U., Quantum Dissipative Systems, World Scientific, Singapore, 2008.
 [3] Vacchini B., Breuer H.P., Phys. Rev. A, 2010, 81, 042103; \doi10.1103/PhysRevA.81.042103.
 [4] Semin V., Sinayskiy I., Petruccione F., Phys. Rev. A, 2012, 86, 062114; \doi10.1103/PhysRevA.86.062114.
 [5] Łuczka J., Physica A, 1990, 167, 919; \doi/10.1016/03784371(90)902998.
 [6] Unruh W.G., Phys. Rev. A, 1995, 51, 992; \doi10.1103/PhysRevA.51.992.
 [7] Morozov V.G., Mathey S., Röpke G., Phys. Rev. A, 2012, 85, 022101; \doi10.1103/PhysRevA.85.022101.
 [8] Morozov V.G., Röpke G., Condens. Matter Phys., 2012, 15, 43004; \doi10.5488/CMP.15.43004.
 [9] Fischer J., Breuer H.P., Phys. Rev. A, 2007, 76, 052119; \doi10.1103/PhysRevA.76.052119.
 [10] Smirne A., Breuer H.P., Piilo J., Vacchini B., Phys. Rev. A, 2010, 82, 062114; \doi10.1103/PhysRevA.82.062114.
 [11] Clos G., Breuer H.P., Phys. Rev. A, 2012, 86, 012115; \doi10.1103/PhysRevA.86.012115.
 [12] Vacchini B., Smirne A., Laine E.M., Piilo J., Breuer H.P., New J. Phys., 2011, 13, 093004; \doi10.1088/13672630/13/9/093004.

[13]
Shresta S., Anastopoulos C., Dragulescu A., Hu B.L., Phys.
Rev. A, 2005, 71, 022109;
\doi10.1103/PhysRevA.71.022109.  [14] Martens C., J. Phys. B: At. Mol. Opt. Phys., 2012, 45, 154008; \doi10.1088/09534075/45/15/154008.
 [15] Chaudhry A.Z, Gong J., Phys. Rev. A, 2013, 87, 012129; \doi10.1103/PhysRevA.87.012129.
Динамiка термостату для точної моделi кубiту при наявностi початкових кореляцiй з його оточенням
В.В. Iгнатюк\refaddrlabel1, В.Г. Морозов\refaddrlabel2
\addresses
\addrlabel1 Iнститут фiзики конденсованих систем НАН України,
вул. Свєнцiцького, 1, 79011 Львiв, Україна
\addrlabel2 Московський державний технiчний унiверситет
радiоелектронiки та автоматики,
просп. Вернадського, 78, 119454 Москва, Росiя