Decoherence and Landauer’s principle in qubit-cavity quantum-field-theory interaction

We consider quantum decoherence and Landauer’s principle in qubit-cavity quantum field theory (QFT) interaction, treating the qubit as the system and cavity QFT as the environment. In particular, we investigate the changes that occur in the system with a pure initial state and environment during the decoherence process, with or without energy dissipation, and compare the results with the case in which the initial state of the system is a mixed state and thus decoherence is absent. When we choose an interaction Hamiltonian such that the energy and coherence of the system change simultaneously, the population change of the system and the energy change are the same when the initial state is mixed. However, the decoherence terms increase the von Neumann entropy of the system. In this case the energy change and decoherence of the system are not independent physical processes. The decoherence process maintains unitarity. On the other hand, if the interaction Hamiltonian does not change the energy of the system, there is only the decoherence effect. The environment will be a distribution in the basis of the displaced number state and always increases the energy. Landauer’s principle is satisfied in both cases.

The question of how a seemingly classical world can emerge from quantum mechanics is one of the greatest foundational problems in modern physics.In quantum mechanics, particles are described by a wave function satisfying the deterministic, linear Schrödinger equation therefore in principle we can compute the time evolution of the state given the initial state of the system and its Hamiltonian.However, the principle of superposition, which is the basic tenet of the quantum mechanics reflected in the linearity of the Schrödinger equation, predicts a probabilistic interpretation of our universe in macroscopic world, thus bringing a long and fascinating history of the measurement problem [1].After all, If the universe is governed by the fundamental principles of quantum mechanics at the microscopic level, why does the macro-world appear to be classical?One of the reasons for this problem is our understanding was built on a tacit but incorrect assumption: the information about a system can be acquired without changing its state.This has led to the idealization of closed systems in physics.However, macroscopic systems are never isolated from their environments.In fact, a truly ideal closed system means that we cannot establish interactions with it in any way, and therefore cannot obtain any information about it.As emphasized by Zeh and Zurek, the Schrödinger equation applies only to a closed system, and the contribution from the environment must be included [2][3][4][5].The relevant theoretical framework for our study is therefore the theory of open quantum systems [6].
This leads us to one of the keywords in the title of our work: decoherence.Decoherence, or more specifically, environment-induced decoherence, represents the loss of quantum coherence, a measure of the definite phase relation between different states of the system.In open quantum systems, the coupling to the environment defines and determines the physical properties of the system.Since almost every system is loosely coupled with other systems, decoherence can be viewed as the loss of information from one system to another.When we consider one of the systems in isolation, the evolution of the system is non-unitary, although the whole system is undergoing a unitary process.The process of decoherence is also the key component in explaining how the classical world emerges from the quantum regime.In the dynamical description of quantum-to-classical transition, the notion of "classicality" is understood as an emergent concept [7].
In many situations, we are used to treating the environment in quantum evolution as something whose details are unimportant and out of our control, whose quantum nature is often neglected.This is despite -in principlethe environment, regardless of its actual content, should also be treated quantum mechanically.There are some justifications for this approach.For example, the detailed dynamics of decoherence can be very complicated, and since interactions are almost ubiquitous, we do not even need to deliberately introduce the environment -even cosmic rays and the microwave background radiation could lead to rapid decoherence of macroscopic bodies.Thus earlier studies have focused on aspects that are under our controls.However, since decoherence has become relevant, perhaps as the key obstacle, to practical applications such as quantum cryptography and quantum computation, it is necessary to further study the issue of the control of decoherence, including the backreaction effects brought to the environment.
In the present work we consider decoherence and Landauer's principle in qubit-cavity quantum field theory (QFT) interaction.The qubit is arguably the simplest quantum system and the basis for quantum computing and quantum cryptography, whereas QFT is the description of Nature that contains all the fundamental interactions except for gravitation (whose quantization is still beyond our understanding).In particular, we will consider during the decoherence process of the system (qubit), with or without dissipation, the changes that occur in the environment (QFT).
Landauer's principle is one of the bridges between thermodynamics and quantum information.It relates the entropy change of a system and the energy consumption of the environment, and it also provides a theoretical limit throughout the evolution [8][9][10].Landauer's principle is based on the following four assumptions: (i) both the system S and environment E are described by Hilbert spaces, (ii) the environment is initially in a thermal state ρ E = e −β ĤE /Tr(e −β ĤE ), where ĤE is the Hamiltonian of the environment and β is the inverse temperature, (iii) system and environment are uncorrelated initially ρ SE = ρ S ⊗ ρ E , (iv) the process proceeds by unitary evolution ρ SE = U ρ SE U † .If all the four assumptions are satisfied, Landauer's principle can be expressed in a form reminiscent of the first law of thermodynamics: ∆Q T E ∆S. ( The quantity ∆Q := Tr ĤE (ρ E − ρ E ) is the heat transferred to the environment E, where ρ E and ρ E denote the final and initial state of E respectively.Here ∆S := S(ρ S )−S(ρ S ) is the difference in the von Neumann entropy between the initial state ρ S and the final state ρ S of the system S.The quantity ρ S/E := Tr E/S [ρ SE ] is the reduced density matrix.Since the process of decoherence is inevitably accompanied by the change in von Neumann entropy, it is natural to ask what changes in the environment will occur accordingly.In addition, the process of decoherence may also be accompanied (or not accompanied) by energy changes in the system S, so we are interested in how this might affect the evolution of the whole system.In the present work we will answer these questions by analyzing the different forms of interactions between qubit and cavity QFT.Henceforth we shall adopt the natural unit system, setting c = = k B = 1 in all the analytical calculations and numerical analyses.

II. ANALYSES OF THE MODELS
In the qubit-cavity QFT interaction the total Hamiltonian Ĥtotal describing our system consists of three terms: Ĥtotal = ĤS + ĤE + Ĥint .The first term ĤS is the free Hamiltonian of the system S and in our case it is just a qubit so we can choose ĤS = Ω 2 σ z , where σ z is the Pauli matrix.The second term ĤE = ∞ j=1 ω j a † j a j is the free Hamiltonian of the cavity QFT, and Ĥint = λη(τ )mφ[x(τ )] is the interaction Hamiltonian, in which λ is a weak coupling constant, and τ denotes the proper time.Here η(τ ) is the so-called "switching function" that controls the interaction, m is the monopole moment of the qubit, and φ[x(τ )] is the field operator at the position of the qubit in the cavity.
Before continuing the analysis we give a short discussion of the interaction Hamiltonian Ĥint to gain some intuition.In interaction picture the field operator φ[x(τ )] reads where the expression of u j [x(τ )] depends on the boundary conditions of the cavity.Thus we can conclude that ĤE and Ĥint do not commute with each other, regardless of the other quantities in Ĥint .The energy expectation value of the cavity QFT is not a conserved quantity in the interaction.
On the other hand, the monopole moment m should be the linear combination of the Pauli matrices σ i (i = x, y, z).For example, if we choose m = σ x (the case σ y is similar), then in the interaction picture m(τ , where |0 , |1 are the ground and excited states of the qubit, respectively.This means that the populations of the qubit density matrix can be exchanged between the ground and excited states, thus causing the change in energy.We can also obtain this by noticing that [ ĤS , Ĥint ] = 0.If we have chosen m = σ z instead, we would have [ ĤS , Ĥint ] = 0, and consequently the qubit does not undergo energy changes.In the present work we will discuss these two cases separately.Since Landauer's principle requires the initial state of the environment to be thermal, the density matrix for E takes the form [11] where for each integral value of j, n j ∈ [0, ∞), and nj := 1/ e ω j T E − 1 is the average photon number.For the qubit, we choose its initial state to be |ψ 0 = √ 1 − p|0 + √ p|1 , in which both √ 1 − p and √ p are chosen to be real numbers.Thus the initial density matrix of the qubit is given by Note that this density matrix differs from (1 − p)|0 0| + p|1 1|.Although both density matrices correspond to the same probability distribution, for the case of (1−p)|0 0|+ p|1 1|, the population may be interpreted as classical probabilities, and the density matrix is not coherent.The main problem studied in this paper is the effect of the changes of the off-diagonal terms of the qubit on the whole system, while the diagonal terms of the qubit change or remain unchanged in the interaction.
A. The case of m = σx In this case we assume that the coupling constant λ is small, thus allowing the use of perturbation methods.The time evolution operator from time τ = 0 to τ = T is given by the Dyson series [12] Û (T, 0 , so the density matrix at a time τ = T is (7) We can expand ρ T order by order as where Let us first consider the ρ T term, in which Û (1) can be written as where I ±,j is defined as Here we have already set η(τ ) = 1 for 0 τ T .Since our initial state ( 4) is thermal and contains only the diagonal terms, the operators a j and a † j from Û (1) acting on ρ 0 would only give rise to |n j − 1 n j | and |n j + 1 n j | (ρ 0 Û (1) † is also similar).Since they are off-diagonal terms, when we subsequently perform the partial trace over the basis of cavity QFT to obtain the reduced density matrix of the qubit, we would necessarily obtain zero.This means that qubit-cavity QFT interaction has no effect at the λ order.In the language of QFT, it is just the one point function φ(x) = 0 for the thermal state.We can also easily see that the contributions from all oddorder λ 2n−1 (n = 1, 2, 3...) vanish.However, we need to emphasize that this is not a general result; it depends on the initial state of the QFT.If the initial state already contains some off-diagonal terms, such as in the case of coherent state, the odd-order λ 2n−1 interaction may create some diagonal terms and the λ order would play the leading role.In the present work we are considering weak coupling and initial state of E to be diagonal, so we shall restrict to at most the λ 2 terms and omit higher order ones.
Next we consider the λ 2 order term ρ T .For the term Û (1) ρ 0 Û (1) † , we already know from ( 12) there are four terms in Û (1) , so Û (1) ρ 0 Û (1) † would have sixteen terms in total.All we need to do is write out all the terms and pick the ones that are not zero after carrying out the partial trace.For example, we can have |n j + 1 n j + 1| and |n j − 1 n j − 1| in the field part, thus performing the partial trace would yield a non-vanishing result, so we can obtain the correction of the qubit density matrix.On the other hand, each term also contains |I −,j | 2 , |I +,j | 2 , etc.These quantities play the role of two-point function of the cavity QFT.As we previously mentioned, u j [x(τ )] depends on the boundary condition.In this work we choose Dirichlet boundary condition and set u j [x(τ )] ∼ sin[k n x(τ )].If the qubit is at rest x(τ ) = constant, then t = τ and u j [x(τ )] gives a constant value.For I −,j , as Ω = ω j , the integrand will just be the value of u j [x(τ )] at x(τ ) = constant, thus I −,j is proportional to the time T .In other words, the qubit and the cavity QFT resonate at the frequency Ω (we refer to this as the "resonance effect").However, for Ω = ω j or the case of I +,j , the integration gives 1−e i(±Ω+ω j )T ±Ω+ωj u j (x), and the contributions from (13) becomes smaller and quickly decays for larger values of ±Ω + ω j .An analogous phenomenon in classical mechanics was reported in [13].In summary, we will select only those terms in Û (1) ρ 0 Û (1) † which are non-vanishing after partial tracing and furthermore must contain |I −,j | 2 .
Similarly, for both Û (2) ρ 0 and ρ 0 Û (2) † , which each contain sixteen terms; we will pick the terms in the same manner.After some tedious calculation, by taking partial trace over the basis of cavity QFT we have the reduced density matrix of the qubit as where and Similarly, for the cavity QFT, partial tracing over the basis of the qubit yields where we have ignored the off-diagonal terms in E, because the energy of the free cavity QFT is only related to the diagonal terms.
If we compare the results here with the case of mixed state (1 − p)|0 0| + p|1 1| being the initial state in [14], where the qubit evolution takes the form of we find that the change of the cavity QFT and δp of the qubit are the same in both cases, but the the off-diagonal term δd in ( 14) is new.This may be confusing for the coherence in qubit does not affect the ∆Q.The reason here lies in the form we have chosen for the interaction Hamiltonian.In order to obtain ∆Q, we need to know the reduced density matrix of the E by taking partial trace over the basis of qubit.Taking Û (1) ρ 0 Û (1) † for example, the matrix associated with the qubit can be written roughly as m(τ )ρ 0 m(τ ).By some algebraic calculations we will find that the off-diagonal terms in the initial state of the qubit does not affect the trace of the matrix.Similarly we can also examine Û (2) ρ 0 and ρ 0 Û (2) † , and we can also find the off-diagonal terms does not affect the trace.However, the off-diagonal terms of the initial state and their variations do affect the von Neumann entropy of the qubit and are thus relevant for Landauer's principle.
The first thing we can observe is that δd is non-negative, which means that quantum coherence can only decrease during the interaction.According to the analysis of |I −,j | 2 , we know that the decoherence process will always continue, dominated by the resonance effect.Furthermore, since the initial state of the qubit is pure, the decoherence process must increase the von Neumann entropy.Diagonalizing ρ S we have the von Neumann entropy of the final state where is the eigenvalues of the ρ S .Omitting the higher-order terms O(δp 2 , δd 2 ), we get Although δp could be negative, in turns out that p − is always non-negative.To see this, we simply insert the expressions ( 15) and ( 16) into the above formula of p − , we have We can easily verify that this is indeed non-negative (the terms in the square bracket is a quadratic in p 2 , with negative discriminant).This also means that in this case the change of energy and decoherence are not independent physical processes.The energy change must be accompanied by the decoherence process to maintain unitarity.
Next we analyze the change in the von Neumann entropy.In Sec.II(B) of [14], Landauer's principle ∆Q T E ∆S has been verified for the qubit with the mixed initial state (1 − p)|0 0| + p|1 1|, and since we have the same ∆Q and δp, if we can prove the entropy change in our case is always smaller than the case of ( 19), Landauer's principle is naturally satisfied.We will find that this is indeed the case.For a non-vanishing p, the entropy change of process ( 19) is On the other hand, if p = 0, the von Neumann entropy will be This is exactly our case, in which p − plays the role of δp.The derivative of the function −(1 − p) ln (1 − p) − p ln p increases when p gets smaller in the range 0 < p < 0.5, and it becomes divergent in the limit p → 0. Thus in our The cases of T E = 1 T = 100 in single-mode approximation.In both figures we set L = 1.234,Ω = ω20, p = 0.2, x = 0.52345 in natural units.In each figure from top to bottom the curves correspond to ∆Q/T E , ∆S in process (19), and ∆S in our case with decoherence, respectively.
case the eigenvalues p ± always give a larger von Neumann entropy and thus yield a smaller and negative ∆S.When p = 0, the two cases coincide.
In Fig. 1 we present explicit examples of the quantities ∆Q/T E and ∆S for both cases to illustrate this result.The two figures corresponds to T E = 1 and T E = 100 respectively.Here we apply the single-mode approximation, focusing on the resonance effect.In each figure from top to bottom the curves correspond to ∆Q/T E , ∆S in process (19), and ∆S in our case with decoherence, respectively.For the case of T E = 1, both ∆Q/T E and ∆S in process (19) are positive, while in the decoherence case ∆S is negative.For the case of T E = 100, ∆Q/T E and ∆S are both negative, but ∆S in the decoherence case is the smallest, which verifies the validity of Landauer's principle.
The case of m = σ y is similar to m = σ x .It can also produce the variation in populations and can trigger decoherence.The calculation and conclusions are the same as m = σ x , so we will not discuss it here.

B. The case of m = σz
The case of m = σ z is different from the m = σ x case, because with [ ĤS , Ĥint ] = 0, the qubit does not undergo energy changes, and since σ z remains unchanged in the interaction picture, we also do not have the e ±iΩτ term.This is why there is no resonance effect in this case.The commutator of Ĥint at different times is not an operator, but just a c-number [15] [ Ĥint (τ ), Ĥint (τ )] = 2i j |u j | 2 sin(ω j (τ − τ )), (26) and the time evolution operator differs only by a timedependent phase from the one obtained by casting the operators in their natural order where the phase ) describes the "additional" motion due to the time ordering.Performing the integral we have where It is not difficult to see that the above formula has the same form as I ±,j .Since σ z in the interaction picture does not contain e ±iΩτ term and we do not have the resonance effect, we need to consider the composition of all modes.The density matrix of the whole system can be directly obtained as We can conclude that the populations of the qubit are unaffected, and the off-diagonal terms are and where is the displacement operator for the j-th field mode, while Tr E (D j (±α j )ρ E ) is the Wigner characteristic function.
Since the thermal state ρ E is Gaussian, we can easily note that it represents a Gaussian function, which immediately leads to the expression We can define the above formula as the suppression factor χ(τ ).On the other hand, the reduced density matrix of the environment is For each j, the quantity ρ E can be described as a distribution in the basis of number state |n j .Thus thus we can obtain the energy of the environment as and the energy change of the environment as which is always non-negative.On the other hand, since the initial state of the qubit is pure, the difference ∆S := S(ρ S ) − S(ρ S ) is always non-positive, thus Landauer's principle is again naturally satisfied.
In FIG. 2 we present an example of the suppression factor χ(τ ).Since we are considering discrete models and |α j | 2 has periodicity, we can observe the phenomenon of Poincaré recurrences.Similarly from FIG. 3 we can also observe that ∆Q/T E is always non-negative while ∆S is non-positive, and Poincaré recurrences are also present.

III. CONCLUSION
In the present work we extend the analysis of Landauer's principle in qubit-cavity QFT interaction [14] to the case that includes decoherence.We analyzed the decoherence of the qubit in the evolutionary process in two cases: m = σ x and m = σ z , corresponding to with and without energy dissipation respectively.The initial state of the qubit is choosen to be |ψ while the cavity QFT is thermal.The initial density matrix is different from the mixed density matrix (1 − p)|0 0| + p|1 1| since it contains coherence (the offdiagonal terms), while the latter may be interpreted as classical probabilities.
In the m = σ x case (similarly for the m = σ y case), we find the energy change of the cavity QFT and population change of the qubit is identical to the case of (1−p)|0 0|+ p|1 1|.The decoherence process corresponds to a new variation of the off-diagonal terms.The energy change must be accompanied by the decoherence to maintain unitarity.Furthermore, since the initial state of the qubit is pure, the interaction can increase the von Neumann entropy more than in the (1−p)|0 0|+p|1 1| case, which has a smaller ∆S.Since Landauer's principle has been proved in the mixed case, it also holds naturally in our case.
In the m = σ z case, we observe that [ ĤS , Ĥint ] = 0, so the qubit does not undergo energy changes.Since the commutator of Ĥint at different times is just a cnumber, the time evolution operator differs only by a timedependent phase from the one obtained by casting the operators in their natural order.We can directly calculate the evolution of the whole system and perform the partial trace.We verify that the populations of the qubit does remain unchanged, while its off-diagonal terms have a suppression factor.The cavity QFT, on the other hand, becomes the distribution in the basis of the displaced number state.The difference ∆S is always non-positive, since the initial state is pure, while ∆Q is always nonnegative, thus Landauer's principle is naturally satisfied.The phenomenon of Poincaré recurrences can also be observed.
In [14] the qubit in accelerated motion (Unruh effect) was also studied.Since the contribution of the acceleration comes from the |I −,j | 2 term, this does not affect any of our conclusions.Recently decoherence was proposed as a mean to detect the Unruh effect [17].We believe that the analysis of qubit-QFT interaction and decoherence can be applied to more contexts, both theoretical and experimental, to deepen our understanding of quantum phenomena.If we use harmonic oscillators instead of qubit to be the system, the interaction can also be solved non-perturbatively in the framework of Gaussian quantum mechanics [18][19][20].
We previously mentioned that QFT formulation can describe all fundamental interactions except gravity.However, gravity can also induce decoherence either via classical or quantum gravitational fluctuations, although a clear understanding of gravitational decoherence is still lacking.See, e.g., [21,22] for reviews and [23] for recent development in this area.Future space-based experiments could help to further understand gravitational decoherence and possibly checking the predictions of different quantum theories [24].It is worth remarking that, just as classical thermodynamics is intimately related to classical gravity (e.g., Jacobson derived Einstein field equations assuming the Bekenstein-Hawking entropy and the first law of thermodynamics [25]), quantum information is also expected to play important roles in quantum gravity [26].In fact, recently it was discovered that black holes would eventually decohere any quantum superpositions [27].Recent studies also included the storage of quantum information at finite temperature in holographic conformal field theories [28].As Landauer's principle is one of the bridges between thermodynamics and quantum information, understanding it better in more realistic quantum interactions, such as those involving decoherence, is a first tiny step towards the grand goal of quantum gravity.We hope to pursue these gravitational connections in our future research.
Finally, let us mention that in the present work we have only considered the original bound of Landauer's principle: ∆Q T E ∆S.If we add more conditions to the previous four assumptions mentioned in the Introduction, we may get a tighter bound [29][30][31][32].The motivation of the these works, including the original paper of Reeb and Wolf [10], are to give some estimates of the evolutionary process of a certain class of systems, i.e. the corresponding bound.These are more general and applicable to a wider range of discussions, allowing us to have a general idea of system evolution without knowing all detailed information about the whole system.In the present work, however, the system is uniquely determined, so that in principle the corresponding physical quantities can be determined directly.We can directly calculate ∆Q and ∆S and thus verify Landauer's principle.