We study the bath dynamics in the dephasing model of a two-state 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 qubit-environment correlations. We obtain the non-equilibrium 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’’.\keywords
open quantum systems, qubit-environment correlations \pacs03.65.Yz
Проведено досл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т
Bath dynamics]Bath dynamics in an exactly solvable qubit model with initial
V.V. Ignatyuk, V.G. Morozov]V.V. Ignatyuk\refaddrlabel1,
\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
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 non-Markovianity has been presented in [11, 12]; (iii) the effect of non-equilibrium 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 , even though there is a build-up of dynamical correlations [4, 9] caused by the entanglement of quantum states. Moreover, even if the effect of non-equilibrium environment on the system behavior is not neglected , 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 so-called dephasing model, describing a two-state system (qubit) coupled to a bosonic bath [5, 6, 7]. In section 3 the non-equilibrium 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 ‘‘qubit-bath’’ 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 spin-boson model describing a two-state 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 )
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 two-level 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 non-ergodic111Another 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 non-trivial 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  rather than the energy dissipation.
It was shown in  that equations of motion for the Heisenberg picture operators , , and can be solved exactly with the results
for the non-equilibrium density matrix of the composite system () and assume that the initial density matrix has the form 
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 222 The equation (8) can be derived by using obvious relations , and the following properties  of the Hamiltonian (1): , .
where . The bath Hamiltonians and the corresponding partition functions are defined as
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 , the system is prepared in some pure state .
where we have introduced the notation
The bath density matrix is obtained from (2) by taking the trace over the qubit states,
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 time-dependent bath density matrix (12) to
denotes a phonon part of in the interaction picture (with the replacement , see equation (5)), whereas the quantity
is the non-equilibrium 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
instead of (6). Proceeding in a similar way, after some algebra one obtains
where denotes the non-equilibrium contribution to the correlation energy in the interaction picture. It is seen from equation (17) that even in this case there is a dynamical build-up of correlations in the system. We touch upon this point in the next section, when analyzing the relevance of non-equilibrium correlations.
3 Phonon non-equilibrium 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  (or the method of the displaced harmonic oscillator modes ), after some straightforward algebra one obtains
is the initial phonon distribution function, and the function
represents the contribution of initial correlations.
With equations (18) and (19), it is easy to calculate the time evolution of the non-equilibrium phonon energy . The final result is conveniently written in terms of the bath spectral density which is introduced by the well-known rule [1, 2, 7]
After simple manipulations we arrive at
where the initial phonon energy is given by
Here the last term occurs due to initial correlations in the system.
where and is a dimensionless coupling constant. This formula ensures both a proper low-frequency behavior of and a cut-off at high frequencies (). The case is usually called the ‘‘Ohmic’’ case, the case ‘‘super-Ohmic’’, and the case ‘‘sub-Ohmic’’. 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 non-equilibrium 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
Combining this expression with equation (22), it is easy to check that the sum is a time-independent quantity (see also a footnote on page 2). We see that the dynamics of the correlation energy exactly compensates the time dependence of the non-equilibrium 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 non-equilibrium 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  vanishes for the same values of the mean inversion population of the levels, manifesting a close relationship between the essentially non-equilibrium behavior of the correlation energy and the onset of the additional channel of decoherence in the system.
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 non-Markovian) and taking into account the intrinsic dynamics of the environment along with the equation of motion for . Such an approach would modify the well-known Zwanzig-Nakajima 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 back-flow 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 non-equilibrium correlations ensures the conservation of the total energy.
-  Breuer H.-P., Petruccione F., The Theory of Open Quantum Systems, Oxford University, Oxford, 2007.
-  Weiss U., Quantum Dissipative Systems, World Scientific, Singapore, 2008.
-  Vacchini B., Breuer H.-P., Phys. Rev. A, 2010, 81, 042103; \doi10.1103/PhysRevA.81.042103.
-  Semin V., Sinayskiy I., Petruccione F., Phys. Rev. A, 2012, 86, 062114; \doi10.1103/PhysRevA.86.062114.
-  Łuczka J., Physica A, 1990, 167, 919; \doi/10.1016/0378-4371(90)90299-8.
-  Unruh W.G., Phys. Rev. A, 1995, 51, 992; \doi10.1103/PhysRevA.51.992.
-  Morozov V.G., Mathey S., Röpke G., Phys. Rev. A, 2012, 85, 022101; \doi10.1103/PhysRevA.85.022101.
-  Morozov V.G., Röpke G., Condens. Matter Phys., 2012, 15, 43004; \doi10.5488/CMP.15.43004.
-  Fischer J., Breuer H.-P., Phys. Rev. A, 2007, 76, 052119; \doi10.1103/PhysRevA.76.052119.
-  Smirne A., Breuer H.-P., Piilo J., Vacchini B., Phys. Rev. A, 2010, 82, 062114; \doi10.1103/PhysRevA.82.062114.
-  Clos G., Breuer H.-P., Phys. Rev. A, 2012, 86, 012115; \doi10.1103/PhysRevA.86.012115.
-  Vacchini B., Smirne A., Laine E.-M., Piilo J., Breuer H.-P., New J. Phys., 2011, 13, 093004; \doi10.1088/1367-2630/13/9/093004.
Shresta S., Anastopoulos C., Dragulescu A., Hu B.L., Phys.
Rev. A, 2005, 71, 022109;
-  Martens C., J. Phys. B: At. Mol. Opt. Phys., 2012, 45, 154008; \doi10.1088/0953-4075/45/15/154008.
-  Chaudhry A.Z, Gong J., Phys. Rev. A, 2013, 87, 012129; \doi10.1103/PhysRevA.87.012129.
Динамiка термостату для точної моделi кубiту при наявностi початкових кореляцiй з його оточенням
В.В. Iгнатюк\refaddrlabel1, В.Г. Морозов\refaddrlabel2
\addrlabel1 Iнститут фiзики конденсованих систем НАН України,
вул. Свєнцiцького, 1, 79011 Львiв, Україна
\addrlabel2 Московський державний технiчний унiверситет
радiоелектронiки та автоматики,
просп. Вернадського, 78, 119454 Москва, Росiя