Abstract
We have studied the decoherence mechanism in a fermion and scalar quantum field theory with the Yukawa interaction in the Minkowski spacetime, using the non-equilibrium effective field theory formalism appropriate for open systems. The scalar field is treated as the system whereas the fermions as the environment. As the simplest realistic scenario, we assume that an observer measures only the Gaussian 2-point correlator for the scalar field. The cause of decoherence and the subsequent entropy generation is the ignorance of information stored in higher-order correlators, Gaussian and non-Gaussian, of the system and the surrounding. Using the 2-loop 2-particle irreducible effective action, we construct the renormalised Kadanoff–Baym equation, i.e., the equation of motion satisfied by the 2-point correlators in the Schwinger–Keldysh formalism. These equations contain the non-local self-energy corrections. We then compute the statistical propagator in terms of the 2-point functions. Using the relationship of the statistical propagator with the phase space area, we next compute the von Neumann entropy, as a measure of the decoherence or effective loss of information for the system. We have obtained the variation of the entropy with respect to various relevant parameters. We also discuss the qualitative similarities and differences of our results with the scenario when both the system and the environment are scalar fields.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
A fundamental mark of quantum mechanics is the superposition principle, which leads naturally to the phenomenon of coherence and interference, as appropriate for isolated or ideal quantum systems. In contrast, realistic quantum systems are never completely isolated from their environment. For such open quantum systems, there is always a source of decoherence due to the interaction or entanglement of the system with the environment or the surrounding [1,2,3]. Such entanglement may influence what we observe locally upon measurement of the system, even from a classical point of view. Quantum decoherence leads to quantum to classical transition and ensures consistency between quantum and classical predictions for the observed system [4].
Quantum decoherence is rapidly gaining interest in the research community, chiefly in the context of interacting quantum field theories [5, 6]. Decoherence should be closely related to the loss or ignorance of information of an open quantum system and as we have mentioned above, is usually formulated in the setup of system plus environment. It can be characterised via various correlations, like the mutual information, discord and the entanglement entropy. Decoherence can be studied in many different scenarios. A general model of decoherence for a non-relativistic quantum particle interacting with a weak stochastic gravitational perturbation was studied in [7, 8]. Decoherence generation from an accelerated time-delay source for an inertial observer is studied in [9]. The decoherence via the gravitational interaction of the dark matter with its environment, consisting of ordinary matter is analysed in [10]. Perhaps one of the most interesting outcome of the decoherence mechanism could be the classicalisation of the primordial inflationary quantum field theoretic perturbations, leading to the large scale structures in the sky we observe today. We refer our reader to [11,12,13,14,15,16,17,18,19,20,21] and references therein for some such discussions.
As of the issue of quantifying decoherence, different strategies exist in the non-equilibrium quantum field theory, e.g. [5, 22,23,24,25,26] and also the references therein. In the conventional approach to compute decoherence, one traces out the inaccessible environmental degrees of freedom to obtain a reduced density matrix of the system, which in general is a mixed one. Using this reduced density matrix one quantifies the decoherence generated in the system in terms of the von Neumann entropy, e.g. [25, 27,28,29,30,31,32,33,34,35,36,37,38]. Another way to tackle the decoherence problem is to consider the master equation approach [39,40,41,42]. In this work, we shall instead be interested in using the correlator approach proposed in [43,44,45], in order to compute the decoherence in terms of the von Neumann entropy due to observer’s lack of ability to know all correlation functions of the system and the environment, eventually characterising the effective loss of information for the system.
Precisely, there can be n-point correlators which can be generated using the n-particle irreducible effective action [5, 46, 47], capturing the information of interaction between the system and its surrounding. In a realistic scenario however, it is certainly impossible for an observer to quantify the correlators of all orders exactly. Thus one needs to consider a practical scenario allowing us to compute only some finite order correlators such as the 2-point or 4-point Gaussian ones. The ignorance of the higher-order Gaussian and non-Gaussian correlators of the system and surrounding leads to the lack of information, generation of quantum decoherence and hence a non-zero entropy defined in some suitable manner. However, we also note that due to progresses in cold atom experiments, simulating quantum field theory models and measuring higher-order correlations becomes an experimentally relevant problem. For example in [48], multipoint correlation functions in both in and out of equilibrium quantum field theories have been experimentally studied to show the deviation of Gaussianity due to the presence of interaction.
In this work we have considered an open quantum field theory at zero temperature, consisting of bosons and fermions, in order to study the decoherence and subsequent entropy generation via the approach proposed in [43,44,45]. Specifically, we have considered a scalar field coupled to fermions via the Yukawa interaction, treating the scalar as the system and the fermions as the surrounding. We shall work in the flat spacetime, and the present work is to be understood as a warm up exercise before we attempt this problem in the context of the early inflationary universe paradigm, where the late time non-perturbative secular effects may be present. Some earlier analysis on open quantum systems with scalars and fermions can be found in [49,50,51,52].
The rest of the paper is organised as follows. In Sect. 2, we have described the model we are considering and the assumptions we are making along with their justifications. Further, following [43,44,45] and references therein, we have defined the von Neumann entropy of the system in terms of the phase space area. In Sect. 2.1, we have outlined the in–in or the Schwinger–Keldysh or the closed time path formalism which will be required to compute the relationship between the statistical propagator with that of the Wightman functions, i.e. the two point functions relevant for our computations. In Sect. 3 we have computed the 2-loop 2-particle irreducible (2PI) effective action using the Schwinger–Keldysh formalism and have found out the Kadanoff–Baym equations for our model. These are basically the equations of motion for the loop corrected two point functions in an interacting quantum field theory, containing self energy corrections. In Sect. 3.1, we have renormalised the Kadanoff–Baym equations by renormalising the self-energy of the scalar field. Finally in Sect. 4, we have computed the phase space area, the statistical propagator and the entropy of the system, which can be thought of as a quantifier of the decoherence generated on the system due to the interaction with the surrounding and the subsequent ignorance of the higher order correlators. The variation of the entropy and phase space area with respect to the rest mass of the system as well as the Yukawa coupling strength is obtained. Finally, in Sect. 5 we conclude our work. The fermions will be taken to be massless, for the sake of simplicity of computation. We shall use the technical formalism of [44, 45], developed in the context when both system and environment are scalar fields.
We shall work with the mostly positive signature of the metric in \(d=(4-\epsilon )\)-dimensional \((\epsilon =0^+)\) Minkowski spacetime and will set \(c=\hbar =1\) throughout.
2 The basic setup
We consider a hermitian massive scalar field \(\phi (x)\) coupled to fermions by the Yukawa interaction,
Since we are working with the mostly positive signature of the metric, the \(\gamma ^{\mu }\)’s satisfy the anti-commutation,
As we have mentioned earlier, the scalar \(\phi (x)\) will play the role of the system, interacting with the environment \(\psi (x)\), \({\overline{\psi }}(x)\). We assume that the environment is at zero temperature and is in its vacuum state.
We assume that our observer (observing the system) can measure two point correlation functions only. Note that a tree level two point function is essentially Gaussian. A quantum corrected such function may or may not be Gaussian. Higher correlations containing the effect of interaction, such as the three point correlators, are non-Gaussian. We have to work on perturbative corrections to the self energy and will restrict ourselves to one loop, \({\mathcal {O}}(g^2)\). Note also that the two loop \({\mathcal {O}}(g^4)\) diagrams of Fig. 1 contains backreaction of the system on the environment, which will be ignored. Restricting ourselves only to one loop seems justified at least when the Yukawa coupling constant g is not too strong and second, in particular when there is no secular effect at late times. Note also that since the fermions are taken to be massless, there can be decay of the system (i.e., the scalar) into fermion-anti-fermion pairs starting at \({\mathcal {O}}(g^2)\). However, this decay involves three point correlators between the system and the surrounding, with on shell process also going on within the surrounding which is unobserved by the observer. Hence we shall ignore such decay in our computations. As we have stated in the preceding section, the lack of observer’s ability to measure higher order correlators due to the system-environment interaction essentially leads to lack of information for the system and it decoheres. This lack of information or decoherence will be quantified by the von Neumann entropy, as follows.
A quantum system can be represented by the density matrix that holds all information about the system. One can calculate various correlators from the density matrix. In particular, different two point correlators are given as [53]
where the first one is the Wightman function, \(\pi (x)={\dot{\phi }}(x)\) is the momentum conjugate to the scalar field and the quantity \(F_{\phi }(x,x')\) is called the statistical propagator defined as
Thus the statistical propagator is basically the Wightman function symmetrised in x and \(x'\). For a given density operator \(\rho (t)\) (pure or mixed), the statistical propagator tells us about how the states are occupied. One can also relate this to the average particle density. We shall focus only on the above two point correlators and their quantum corrections, as a practical scenario.
In the spatial momentum space, the statistical propagator reads,
One also defines the phase space area for each spatial momentum mode as the Fourier transform of the quantity
given by
where we have used the equal time limits of Eqs. (2) and (4). In order to see the analogy of this construction with that of ordinary quantum mechanics, let us recall the generalised uncertainty relation (with \(\hbar =1\)),
Combining the two above, we write
Thus the quantity \(\Delta \ge 1\), interpreted as the phase space area, can be thought of as a measure of the impurity of the quantum state. One plausible way to understand the increase in \(\Delta \) is the transfer of momentum between the system and the environment, thereby increasing the momentum uncertainty. Putting these all in together, the loss or ignorance of information due to the inaccessibility of all the correlations in the system, environment and between them, is characterised via the von Neumann entropy for our field theoretic continuum system [53],
One can also relate the phase space area with the statistical particle number density per mode as
We shall see below that \(\Delta _{|{\vec {k}}|}(t)\) becomes identity in the absence of the Yukawa interaction.
For our purpose, we need to compute the various two-point functions in the in–in or the Schwinger–Keldysh formalism, which we outline below.
2.1 Propagators in the in–in formalism
The Schwinger–Keldysh or the in–in formalism is useful for studying the quantum dynamics of a system in a non-equilibrium scenario [54, 55]. Using this, one can meaningfully compute the causal expectation value of an operator with respect to some suitable initial state. Precisely, one follows the time evolution of an operator from some initial state at \(t = t_0\), without having the knowledge of appropriate late time states. Evolution of an operator from \(t_0\) to t requires a time-ordered evolution and then it is brought back to the initial time following an anti-time ordered evolution t to \(t_0\), as depicted in Fig. 2.
Thus in the in–in formalism, the expectation value of an operator O(t) with respect to some initial density operator \(\rho (t_{0})\) (defined in the Heisenberg picture) is given by,
where \({\overline{T}}\) stands for the anti-time ordering, and H(t) is the Hamiltonian. The generating functional for the corresponding path integral for Eq. (1), subject to Fig. 2 is given by
where for notational convenience we have written \([{\mathcal {D}}\psi ]={\mathcal {D}}\bar{\psi }{\mathcal {D}}\psi \), \([J_{+}^{\psi } \psi ^{+}]= J_{+}^{\psi } \psi ^{+}+\bar{\psi }^{+}\bar{J_+}^{{\bar{\psi }}^+} \) etc. The fields \(\phi ,\, \psi ,\,\bar{\psi } \) have their sources \(J^{\phi }\), \(J^{\psi }\) \({\bar{J}}^{{\bar{\psi }}}\) respectively. Note that there are two species of fields and sources for each category. The \(+\) sign denotes forward evolution in time whereas the − sign denotes backward evolution. The \(\delta \)-function ensures the field configurations are the same on the final hypersurface at \(t=t_f\). The expectation values of n-point functions can be found by taking functional differentiation of the generating functional Eq. (11) with respect to the sources J’s, for instance
With the help of the above, we now define the following propagators for the scalar field,
where we have taken \(t>t'\) above. The time ordered and the anti-time ordered propagators are respectively the Feynman and anti-Feynman propagators and the rest are the two Wightman functions. One can also write
These propagators and Wightman functions also satisfy the following properties,
The retarded and the advanced propagators, useful for our later purpose, are respectively defined as
One also defines the spectral two point function and the statistical propagator as (e.g. [44]),
The spectral function can tell us about the states of the system and the spectrum, but it does not hold any information on how these states are occupied. The statistical propagator on the other hand, gives information about how the states are populated. Hence the latter is much relevant to study decoherence and entropy [22] (also references therein).
The free Feynman, anti-Feynman and the Wightman functions satisfy
The propagators for the fermionic field can be defined in a likewise manner, keeping in mind the anti-commutations satisfied by them and the replacement of the d’Alembertian \(\partial ^2\) by . We shall consider a zero temperature field theory, and hence various averages such as in Eq. (13) are to be understood simply as the vacuum expectation values below.
With these equipments, we are now ready to compute the loop correction to the two-point functions and to find out the Kadanoff–Baym equation, i.e. the equation of motion satisfied by them.
3 Derivation of the Kadanoff–Baym equations
The Kadanoff–Baym equations are integro-differential equations satisfied by the two-point functions in an interacting quantum field theory [56]. It is well known that the one particle irreducible effective action, when varied with respect to some background field, gives rise to the quantum corrected field equation. Likewise, the 2PI effective action, when varied with respect to the two point functions, yields equation of motion satisfied by them known as the Kadanoff–Baym equations. Since the 2PI effective action contains quantum corrections, we obtain extensions of free theory equations like Eq. (19), essentially containing the effect of loops. One can obtain the 2PI effective action as a double Legendre transform from the generating functional for connected Green functions with respect to the linear source J and also another quadratic source [22, 47, 57,58,59]. These equations contain the effects of the non-local self energy.
We expand the effective action corresponding to Eq. (1) up to two loop order as
where \(s,s'=\pm \), and \(\imath {}S^{ss^{\prime }}_{\psi }\) are the fermion propagators. \(\Gamma ^{(2)}\) denotes the 2PI contribution to the effective action at two loop, as shown in Fig. 3. Note however that since we are considering massless fermions, the tadpoles vanish. Even with a massive fermion, the tadpoles can be completely renormalised away in the flat spacetime. Thus we need to consider only the sunset like diagram of Fig. 3. We now explicitly write down the parts of Eq. (20) containing quantum corrections as
where the last line corresponds to the first of Fig. 3.
We now vary the effective action with respect to the bosonic and the fermionic propagators, \(i\Delta _{\phi }\) and \(iS_{\psi }\), to obtain the equations of motion
explicitly giving respectively
We next multiply Eqs. (23a) and (23b) respectively by \(s \imath \Delta ^{s^{\prime }s^{\prime \prime }}_{\phi }(x^{\prime },x^{\prime \prime })\) and \(s \imath {}S^{s^{\prime }s^{\prime \prime }}_{\psi }(x^{\prime },x^{\prime \prime })\) (from the right), and then integrate over \(x^{\prime }\) and sum over \(s^{\prime }=\pm \). Recalling that the integration over the product of the propagator and its inverse gives a \(\delta \)-function, we obtain the one-loop Kadanoff–Baym equations
where we have abbreviated the one loop the self-energies as
Since we are considering the scalar to be our system which is observed, we shall consider only the scalar self energy, Eq. (25a). The corresponding Feynman diagram is given by the first of Fig. 1. Accordingly, by expanding the summations, we rewrite Eq. (24a) as
Setting \(g=0\) above makes the self energies vanishing, thereby reproducing the free theory results of Eq. (19). Note also that even though the self energies appearing in the above equations are of one loop order, we may integrate these equations to find out the propagators, eventually non-perturbative in the coupling constant. In other words, the Kadanoff–Baym equations gives a framework to resum the self energies. This seems to be in particular useful in the context of the primordial cosmic inflation, where late time secular effects may be present, necessitating resummation, e.g. [60,61,62].
We shall solve Eq. (26) by going to the momentum space. We define the Fourier transform
in terms of which Eq. (26) become
where \(k^2=k_{\mu }k^{\mu }=-k_{0}^{2}+|{\vec {k}}|^{2}\). Solving these coupled algebraic equations, we can find out various two point functions. For example, on subtracting Eq. (28b) from Eq. (28c), we have the momentum space expression for the advanced propagator Eq. (17b),
where the advanced self-energy \(\imath M^{\textrm{a}}_{\phi }(k)\) is given by
Likewise we find the Wightman functions
where the retarded self-energy \(\imath M^{\textrm{r}}_{\phi }(k)\) is given by
Substituting Eq. (29) into Eqs. (30a) and (30b), we finally obtain the momentum space expression of the statistical propagator Eq. (18),
Recall that the statistical propagator will yield the expression of the phase space area and entropy, Eqs. (5) and (8). Thus in order to compute the statistical propagator, we need to determine various self-energies, as appearing in Eq. (31). However, note that the above expressions are not renormalised. Hence we shall use the renormalised self energies in Eq. (31), in order to compute the entropy. Also in particular, note that since the self energies are \({\mathcal {O}}(g^2)\), the only coupling constant dependence of the above expression comes in the denominator of the terms within the parenthesis. Thus the expression for the statistical propagator is actually non-perturbative and contains the resummed self-energy corresponding to the series of one-loop diagrams (i.e., the first of Fig. 1), owing to the Kadanoff–Baym equations.
3.1 The retarded self-energy and its renormalisation
In this section, we compute the self-energy \(\imath M_{\phi }^{ss^{\prime }}(x,x')\) and find out the renormalised retarded self-energy, \(\imath M_{\phi ,{\textrm{ren}}}^{\textrm{r}}(x,x')= \imath M_{\phi ,{\textrm{ren}}}^{++}(x,x') - \imath M_{\phi }^{+-}(x,x')\), to be useful for our future purpose. Due to the subtraction of the two propagators, it is easily done in coordinate space. Necessary techniques in order to deal with such coordinate space computations can be seen in, e.g. [63, 64].
The fermion propagator \(\imath {}S^{ss^{\prime }}_{\psi }(x,x')\) is obtained by acting on the massless scalar field propagator,
The Poincaré invariant biscalar distance functions with appropriate \(i\epsilon \) prescription, \(\Delta x_{ss^{\prime }}^{2}(x,x')\), necessary for the in–in formalism are defined as
The one loop scalar self-energy \(\imath M_{\phi }^{++}(x,x')\) is readily found from Eqs. (25a) and (32)
Similarly, we can find out the other self-energies \(\imath M_{\phi }^{--}(x,x')\), \(\imath M_{\phi }^{+-}(x,x')\) and \(\imath M_{\phi }^{-+}(x,x')\) using the suitable \(i\epsilon \) prescriptions as given in Eq. (33). We will now identify the divergence of Eq. (34). Note first that for an arbitrary exponent \(\alpha \), we have
Furthermore, we can write
For the other distance functions of Eq. (33), we have
We now rewrite Eq. (34) using Eqs. (35) and (36a) as
where \(\mu \) is an arbitrary mass scale. We now Taylor expand the terms inside the curly brackets around \(d=4\) to obtain,
The first term of the above expression contains an ultraviolet divergence around \(d=4\) which we have separated and the second term contains a non-local contribution to the self-energy. Since the divergence contains a \(\partial ^2\), we have to add a scalar field strength renormalisation counterterm \((\phi \rightarrow (1+\delta Z)^{1/2}\phi ).\) This yields in the action an additional kinetic term
where we have ignored a total divergence. The amputated version of the Feynman diagram corresponding to the above term involves two functional differentiations with respect to the two scalar field operators, yielding two \(\delta \)-functions. Integrating either of them, and choosing
we remove the divergence of Eq. (38). Note that such renormalisation could also be performed at the level of the Kadanoff–Baym equations, Eq. (26). In that case, the counterterm contribution is achieved by replacing the self energies \(\imath M_{\phi }^{++}(x,x')\) and \(\imath M_{\phi }^{--}(x,x')\), which contain divergent contributions, by the amputated counterterm contribution corresponding to \(\delta Z\) given above.
\(\imath M_{\phi }^{--}(x,x')\) is given by just the complex conjugation of Eq. (38)
The divergence appearing in the above expression can be tackled as above, with the counterterm Eq. (39). Thus the renormalised expressions for the above two self-energies are given by
Also, choosing the suitable pole prescription from Eq. (33), \(\imath M_{\phi }^{+-}(x,x')\) and \(\imath M_{\phi }^{-+}(x,x')\) are determined by
Note that \(\imath M_{\phi }^{+-}(x,x')\) and \(\imath M_{\phi }^{-+}(x,x')\) do not need any renormalisation as they do not contain any divergence around \(d = 4\), as can be readily verified from Eq. (36c). From Eqs. (41) and (42), we now have the renormalised retarded self-energy
Using now
we can write the nonlocal terms of Eq. (43) as,
Breaking now the logarithms into real and complex parts, we have
Putting these all in together the renormalised retarded self-energy, Eq. (43), takes the form
The step functions appearing above ensure that the retarded self-energy is non-vanishing only if \(\Delta t>0\) and \(\Delta t^2-\Delta x^2>0\). This ensures the expected causal characteristics of the retarded self-energy.
As we have stated earlier, we wish to make a momentum space computation of the entropy and hence the statistical propagator, Eq. (31). The relevant renormalised expressions, including that of various self energies are found in Appendices A and B. In particular, the Fourier transform of Eq. (47) is done in Appendix B. Using these ingredients, we wish to compute the von Neumann entropy in the following section.
4 Phase space area and entropy
Renormalised expression for the statistical propagator is given by the renormalised version of Eq. (31)
The renormalisation of the various self-energies has been performed in the preceding section and Appendices A and B. Substituting now Eqs. (66), (69) and (71) into the above equation, we find after some algebra
We have plotted the above statistical propagator with respect to the dimensionless variable \(k^0/\mu \) in Figs. 4 and 5. Note that as the Yukawa coupling g gets smaller, \(F_{\phi }(k)\) approaches a \(\delta \)-function dispersion, as expected from Eq. (49). However as g increases, the \(\delta \)-function peak becomes broadened to a quasi-particle peak like that of the Breit–Wigner kind. Further broadening of the peak with increasing g implies that the resonance becomes broadened and we can no longer sensibly talk about a quasi-particle. Also, these plots show that with the increase in \(m/\mu \), the peak gets shifted towards the higher values of \(k^0/\mu \) with a decrease in the value of the statistical propagator, which means that the corresponding state becomes less populated.
Now in order to compute the entropy, we would use the definition of the same in terms of the phase space area and the three momentum, Eqs. (5) and (8). We shall assume the mass of the scalar field to be time independent. However, one can also consider time-dependent mass as a signature of a non-equilibrium system, as has been considered in [44]. A time dependent mass function will break the time translation invariance, inducing an explicit proper time dependence on the statistical propagator. We also note that in standard quantum field theories, one generally talks about the early or late times when the system is in equilibrium and it is found in some eigenstates of the free Hamiltonian. For a system which is out-of-equilibrium however, an exact distinction between energy states is unclear. A standard approach in such scenario is the adiabatic approximation, in which one specifies a reference set of approximate states under the assumption of a slowly varying dynamical background. Using the projection of system’s evolution onto these approximate states, one may hope to study its dynamics at intermediate times. There are schemes to truncate the adiabatic expansion of the system to form the aforementioned approximate basis set at intermediate times. For example in [44, 45], the Bogoliubov transformations are used to achieve the same.
For turning on the perturbation non-adiabatically which we have not considered here, we consider replacing the coupling constant g with \(g\rightarrow g\,\theta (t-t_0)\). The step function changes the limit of integration Eq. (26) from \(-\infty \) to \(t_0\) and \(t_0\) to some final time. The self energies vanish in the absence of interaction and would only contribute for the second time interval. Thus it is clear that the Schwinger–Keldysh contours presented in Fig. 2 are not equivalent to this abrupt switching case. Note that turning on the interaction abruptly may certainly lead to particle creation, which may significantly modify the statistical propagator. Some relevant discussion on this can be seen in e.g. [44].
Using now the Fourier transforms, we have from Eq. (49)
Substituting now Eq. (49) into the above integrals, we have evaluated them numerically. For example, for \(|{\vec {k}}|/\mu =1\), \(m/\mu =2\) and \(g=0.5\), we find the numerical value of the phase space area Eq. (5), to be
which is indeed greater than unity and hence indicating a non-vanishing von Neumann entropy, as dictated by Eq. (8). We have further analysed the variation of the phase space area with respect to the dimensionless system mass \(m/\mu \), by fixing all the other parameters in Fig. 6. Thus \(\Delta \) decreases monotonically with increasing \(m/\mu \) and asymptotically reaches unity, indicating very small or almost vanishing entropy. This corresponds to the fact that with the increasing mass, the system becomes more stable, i.e., it becomes difficult for the surrounding to disturb a heavier system.
From Eq. (8), we obtain the non-vanishing von Neumann entropy corresponding to the phase space area of Eq. (51) (\(|{\vec {k}}|/\mu =1\), \(m/\mu =2\) and \(g=0.5\))
We have also analysed the variation of the entropy with respect to the system mass by fixing all the other parameters in Fig. 7. Due to the aforementioned reason as that of the phase space area, the entropy is decreasing with the increasing \(m/\mu \).
As we have emphasised earlier, such non-vanishing entropy is due to the coupling between the system and the environment along with the ignorance of all kinds of correlations between them for an observer located in the practical world. The non-trivial statistical propagator gives up phase space for the system field that previously was inaccessible to it. Such increase in the accessible phase space for the system in turn, implies that less information about the system field is accessible to the observer, and hence we observe a non-vanishing von Neumann entropy. Due to this reason, as our results show, the increase in the coupling strength increases the von Neumann entropy.
5 Conclusion
This work studies the decoherence and subsequent entropy generation using the non-equilibrium effective field theory for the Yukawa interaction in the Minkowski spacetime, using the correlator approach of [43,44,45]. The scalar is treated as the system, whereas the fermions as the surrounding and it is assumed that the surrounding or the environment is in its vacuum state. We have assumed as the simplest realistic scenario that the observer measures only the leading and next to the leading order, two point correlators for the system. The ignorance about all the other correlators yields the generation of entropy.
We have constructed the Kadanoff–Baym equations from the 2-loop 2PI effective action, i.e. the quantum corrected equations satisfied by the two point correlators in the in–in formalism in Sect. 3. In Sect. 3.1, we renormalise our results and finally compute the phase space area, the statistical propagator and the entropy in Sect. 4. Note that all the results found in the preceding section are explicitly independent of time, owing to the fact that we have not taken any time dependent background field or time dependent mass. Discussion on the latter scenario can be seen in [44].
For a fixed value of all the relevant parameters such as mass, coupling, momentum and energy, we observe a higher numerical value of phase space area and entropy in (cf., Sect. 4), compared to the scenario when both the system and the surrounding are scalars [44]. This is due to the fact that the Yukawa coupling is dimensionless whereas the cubic coupling considered in [44] is dimensionfull, leading to different momentum dependence in the statistical propagator Eq. (49). Although we note that the qualitative aspects of the variation of the entropy (e.g. Fig. 7) for both the cases are similar.
Finally, we note that the entire above analysis might seem to rest upon the implicit assumption that the observer cannot measure the correlations corresponding to the fermions. For the case when the fermions represent a thermal bath and the scalar is initially at zero temperature, such distinction between the system and the environment seems obvious. For the zero temperature case we have considered presently, however, such distinction might seem albeit arbitrary. For example, we do not have any explicit hierarchy of scales of physical quantities here which decides such partitioning. Moreover, we could also have computed the correlators corresponding to the fermions instead. Thus for the zero temperature field theory, the above analysis is based upon the fact that the observer only measures the correlations for the scalar field, and the effect of fermions comes only as virtual particles inside a loop. Such implicit assumption is justified only when the system is much ‘small’ compared to the surrounding, i.e., when it can be considered as a ‘bath’ at zero temperature. This seems to have some qualitative similarity with the standard formalism of tracing out the fermionic degrees of freedom and look into the effective scalar field dynamics. Now, since we are ignoring a part of the theory, it would lead to an entropy, as we have obtained above. An obvious extension of this would be to go to the two loop self energies of Fig. 1, which contains some effect of system’s backreaction onto the surrounding. Perhaps more importantly, one should also consider the three-point scalar-fermion-anti-fermion correlators containing the effect of the decay of the scalar. This would involve constructing an effective action generating those correlators. We reserve these works as future tasks.
As we have emphasised earlier, understanding such decoherence mechanism can be very important in the early inflationary universe scenario, which is in fact our chief motivation. Even though this also must be a zero-temperature case, the short and the super-Hubble wavelength parts of the scalar field might give us some clue to naturally identify the system and the surrounding. Some related works, using formalism different from this paper can be seen in [65,66,67]. Most importantly, for a massless minimal scalar in an inflationary background, we may expect the appearance of late time, non-perturbative secular effects. Resumming them non-perturbatively using the Kadanoff–Baym equations could be a challenging task. We hope to return to these issues in our future publications.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: We have not used any data for this work.]
References
M. Schlosshauer, Quantum decoherence. Phys. Rep. 831, 1 (2019). arXiv:1911.06282 [quant-ph]
C. Kiefer, E. Joos, Decoherence: concepts and examples. Lect. Notes Phys. 517, 105–128 (1999). arXiv:quant-ph/9803052
W.H. Zurek, Environment-induced superselection rules. Phys. Rev. D 26, 1862 (1982)
M. Schlosshauer, The quantum-to-classical transition and decoherence. arXiv:1404.2635 [quant-ph]
E. Calzetta, B.L. Hu, Nonequilibrium Quantum Field Theory (Cambridge University Press, Cambridge, 2008)
E. Calzetta, B.L. Hu , Correlations, decoherence, dissipation, and noise in quantum field theory. arXiv:hep-th/9501040
L. Asprea, G. Gasbarri, A. Bassi, Gravitational decoherence: a general non relativistic model. Phys. Rev. D 103, 104041 (2021). arXiv:1905.01121 [quant-ph]
L. Asprea, A. Bassi, H. Ulbricht, G. Gasbarri, On the decoherence effect of a stochastic gravitational perturbation on scalar matter and the possibility of its interferometric detection. Phys. Rev. Lett. 126, 200403 (2021). arXiv:1912.12732 [gr-qc]
S. Onoe, D. Su, T. Ralph, Particle production and apparent decoherence due to an accelerated time-delay. Phys. Rev. D 98, 036011 (2018). arXiv:1806.00929 [quant-ph]
I. Allali, M.P. Hertzberg, Gravitational decoherence of dark matter. JCAP 07, 056 (2020). arXiv:2005.12287 [gr-qc]
A. Ashtekar, A. Corichi, A. Kesavan, Emerge Universe 2018, 4(6), 71nce of classical behavior in the early universe. Phys. Rev. D 102, 023512 (2020). arXiv:2004.10684 [gr-qc]
M. Rotondo, Y. Nambu, The decoherence and interference of cosmological arrows of time for a de Sitter universe with quantum fluctuations. Universe 4(6), 71 (2018). arXiv:1805.02346 [gr-qc]
D. Campo, Decoherence and entropy of primordial fluctuations. I. Formalism and interpretation. Phys. Rev. D 78, 065044 (2008). arXiv:0805.0548 [hep-th]
D. Campo, R. Parentani, Decoherence and entropy of primordial fluctuations II. The entropy budget. Phys. Rev. D 78, 065045 (2008). arXiv:0805.0424 [hep-th]
E. Calzetta, B.L. Hu, Quantum fluctuations, decoherence of the mean field, and structure formation in the early Universe. Phys. Rev. D 52, 6770 (1995)
S. Brahma, O. Alaryani, R. Brandenberger, Entanglement entropy of cosmological perturbations. Phys. Rev. D 102, 043529 (2020). Published 31 August 2020
T. Janssen, T. Prokopec, A graviton propagator for inflation. Class. Quantum Gravity 25, 055007 (2008). arXiv:0707.3919 [gr-qc]
P. Friedrich, T. Prokopec, Entropy production in inflation from spectator loops. Phys. Rev. D 100, 083505 (2019). arXiv:1907.13564 [astro-ph.co]
T. Markkanen, Decoherence can relax cosmic acceleration: an example. JCAP 09, 022 (2017). https://doi.org/10.1088/1475-7516/2017/09/022. arXiv:1610.06637 [gr-qc]
B.L. Hu, Quantum statistical processes in the early universe. Vistas Astron. 37, 391 (1993). https://doi.org/10.1016/0083-6656(93)90066-S. arXiv:gr-qc/9302031
B.L. Hu, Statistical mechanics and quantum cosmology. Elsevier, 233–252 (1991). https://doi.org/10.1016/B978-0-444-88903-4.50025-1. arXiv:grqc/9511079
J. Berges, Introduction to nonequilibrium quantum field theory. AIP Conf. Proc. 739, 3–62 (2005). arXiv:hep-ph/0409233
D. Glavan, T. Prokopec, A pedestrian introduction to non-equilibrium QFT
J. Luczka, Non-Markovian stochastic processes: colored noise. Chaos 15, 026107 (2005)
D. Boyanovsky, Information loss in effective field theory: entanglement and thermal entropies. Phys. Rev. D 97, 065008 (2018). arXiv:1801.06840 [hep-th]
F.C. Lombardo, Influence functional approach to decoherence during inflation. Braz. J. Phys. 35, 391–396 (2005)
D. Boyanovsky, Effective field theory during inflation. II. Stochastic dynamics and power spectrum suppression. Phys. Rev. D 93, 043501 (2016). arXiv:1511.06649 [astro-ph.CO]
S. Bhattacharya, S. Chakrabortty, H. Hoshino, S. Kaushal, Background magnetic field and quantum correlations in the Schwinger effect. Phys. Lett. B 811, 135875 (2020). arXiv:2005.12866 [hep-th]
M.S. Ali, S. Bhattacharya, S. Chakrabortty, S. Kaushal, Fermionic Bell violation in the presence of background electromagnetic fields in the cosmological de Sitter spacetime. Phys. Rev. D 104, 125012 (2021). arXiv:2102.11745v1 [hep-th]
S. Bhattacharya, S. Chakrabortty, S. Goyal, Dirac fermion, cosmological event horizons and quantum entanglement. Phys. Rev. D 101, 085016 (2020). arXiv:1912.12272 [hep-th]
S. Bhattacharya, H. Gaur, N. Joshi, Some measures for fermionic entanglement in the cosmological de Sitter spacetime. Phys. Rev. D 102, 045017 (2020). arXiv:2006.14212 [hep-th]
S. Bhattacharya, N. Joshi, Entanglement degradation in multi-event horizon spacetimes. Phys. Rev. D 105, 065007 (2022). arXiv:2105.02026 [hep-th]
S. Kaushal, Schwinger effect and a uniformly accelerated observer. Eur. Phys. J. C 82, 872 (2022). arXiv:2201.03906 [hep-th]
F. Marquardt, A. Patmann, Introduction to dissipation and decoherence in quantum systems. arXiv:0809.4403 [quant-ph]
S. Deffner, Quantum entropy production in phase space. EPL 103, 30001 (2013). arXiv:1307.3183 [cond-mat.stat-mech]
A. Giraud, J. Serreau, Decoherence and thermalization of a pure quantum state in quantum field theory. Phys. Rev. Lett. 104, 230405 (2010)
J. Berges, S. Borsanyi, J. Serreau, Thermalization of fermionic quantum fields. Nucl. Phys. B 660, 00261 (2003). arXiv:0212404 [hep-ph]
T. Prokopec, M.G. Schmidt, J. Weenink, The Gaussian entropy of fermionic systems. Ann. Phys. 327, 3138 (2012). arXiv:1204.4124 [hep-th]
R.Z. Shaisultanov, Back reaction in scalar QED, Langevin equation and decoherence functional. arXiv:hep-th/9509154
T.J. Hollowood, J.I. McDonald, Decoherence, discord and the quantum master equation for cosmological perturbations. Phys. Rev. D 95, 103521 (2017)
C. Anastopoulos, B.L. Hu, A master equation for gravitational decoherence: probing the textures of spacetime. Class. Quantum Gravity 30, 165007 (2013)
C.H. Chou, T. Yu, B.L. Hu, Exact master equation and quantum decoherence of two coupled harmonic oscillators in a general environment. Phys. Rev. E 77, 011112 (2008). https://doi.org/10.1103/PhysRevE.77.011112. arXiv:quant-ph/0703088
J.F. Koksma, T. Prokopec, M.G. Schmidt, Entropy and correlators in quantum field theory. Ann. Phys. 325, 1277 (2010)
J.F. Koksma, T. Prokopec, M.G. Schmidt, Decoherence in an interacting quantum field theory: the vacuum case. Phys. Rev. D 81, 065030 (2010). arXiv:0910.5733 [hep-th]
J.F. Koksma, T. Prokopec, M.G. Schmidt, Decoherence in an interacting quantum field theory: thermal case. Phys. Rev. D 83, 085011 (2011). arXiv:1102.4713 [hep-th]
B. Fayzullaev, Singular structure of the QED effective action. Int. J. Mod. Phys. Conf. Ser. 49, 1960006 (2019). arXiv:1804.02313 [hep-th]
E. Calzetta, B.L. Hu, Nonequilibrium quantum fields: closed time path effective action, Wigner function and Boltzmann equation. Phys. Rev. D 37, 2878 (1988)
I. Kukuljan, S. Sotiriadis, G. Takacs, Correlation functions of the quantum sine-Gordon model in and out of equilibrium. Phys. Rev. Lett. 121, 110402 (2018)
A. Nüßeler, I. Dhand, S.F. Huelga, M.B. Plenio, Efficient simulation of open quantum systems coupled to a fermionic bath. Phys. Rev. B 101(15), 155134 (2020). arXiv:1909.09589 [quant-ph]
K. Enqvist, J. Hogdahl, Scalar condensate decay in a fermionic heat bath in the early universe. JCAP 09, 013 (2004). arXiv:hep-ph/0405299
A. Karmakar, G. Gangopadhyay, Decoherence without dissipation due to fermionic bath. Phys. Scr. 85, 045008 (2012)
J. Lankinen, J. Malmi, I. Vilja, Fermionic decay of a massive scalar in the early Universe. Eur. Phys. J. C 80(6), 502 (2020). arXiv:1904.05084 [gr-qc]
J.F. Koksma, T. Prokopec, M.G. Schmidt, Ann. Phys. 325, 1277–1303 (2010). arXiv:1002.0749 [hep-th]
J.S. Schwinger, Brownian motion of a quantum oscillator. J. Math. Phys. 2, 407 (1961)
L.V. Keldysh, Diagram technique for nonequilibrium processes. Zh. Eksp. Teor. Fiz. 47, 1515 (1964) [Sov. Phys. JETP 20, 1018 (1965)]
L.P. Kadanoff, G. Baym, Quantum Statistical Mechanics (Benjamin Press, New York, 1962)
J.M. Cornwall, R. Jackiw, E. Tomboulis, Effective action for composite operators. Phys. Rev. D 10, 2428 (1974)
R. Jackiw, Functional evaluation of the effective potential. Phys. Rev. D 9, 1686 (1974)
A. Lesniewski, Effective action for the Yukawa(2) quantum field theory. Commun. Math. Phys. 108, 437–467 (1987)
J.A. Cabrer, D. Espriu, Secular effects on inflation from one-loop quantum gravity. Phys. Lett. B 663, 361–366 (2008). arXiv:0710.0855 [gr-qc]
N.C. Tsamis, R.P. Woodard, Stochastic quantum gravitational inflation. Nucl. Phys. B 724, 295–328 (2005). arXiv:gr-qc/0505115
T. Brunier, V.K. Onemli, R.P. Woodard, Two loop scalar self-mass during inflation. Class. Quantum Gravity 22, 59–84 (2005). arXiv:gr-qc/0408080
L.D. Duffy, R.P. Woodard, Yukawa scalar self-mass on a conformally flat background. Phys. Rev. D 72, 024023 (2005). arXiv:hep-ph/0505156
S.P. Miao, R.P. Woodard, Leading log solution for inflationary Yukawa. Phys. Rev. D 74, 044019 (2006). arXiv:gr-qc/0602110
T.J. Hollowood, J.I. McDonald, Decoherence, discord and the quantum master equation for cosmological perturbations. Phys. Rev. D 95(10), 103521 (2017). arXiv:1701.02235 [gr-qc]
D. Boyanovsky, Imprint of entanglement entropy in the power spectrum of inflationary fluctuations. Phys. Rev. D 98(2), 023515 (2018). arXiv:1804.07967 [astro-ph.CO]
P. Friedrich, T. Prokopec, Entropy production in inflation from spectator loops. Phys. Rev. D 100(8), 083505 (2019). arXiv:1907.13564 [astro-ph.CO]
M.E. Peskin, D.V. Schroeder, An Introduction to Quantum Field Theory (Westview Press, Boulder, 1995)
M. Abramowitz, I. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (National Bureau of Standards, Gaithersburg, 1964)
Acknowledgements
The authors would like to thank anonymous referee for a careful critical reading of the manuscript and for making various useful comments.
Author information
Authors and Affiliations
Corresponding author
Additional information
Sourav Bhattacharya: On lien from IIT Ropar, Punjab, India.
Appendices
Appendix A: Calculation of renormalised \(\imath M_{\phi ,\textrm{ren}}^{++}(k)\) and \(\imath M_{\phi ,\textrm{ren}}^{--}(k)\)
In this appendix, we compute the renormalised self-energies \(\imath M_{\phi }^{++}\) and \(\imath M_{\phi }^{--}\) in Fourier space, to be useful for our future purpose. From Eq. (25a), we have in momentum space
where we have used Feynman’s trick (e.g. [68]). Here p and k are the internal and external momentum corresponding to fermion and scalar propagators respectively. We have also used \(q=p+xk\), as the new integration variable. Performing now the usual Euclideanisation of q in the last integral of the above equation, we have
The integral can now straightforwardly be performed, yielding
As of the coordinate space expression Eq. (40), the divergence can be absorbed in the scalar field strength renormalisation. Accordingly we have
\(\imath M_{\phi }^{--}(k)\) on the other hand, will be the negative of the complex conjugation of the above, follows directly from Eq. (41),
We note that the above expression is similar to that of the case when the environment is also a scalar [44], however in our case we have one extra factor of \(k^2\) multiplied with the logarithm.
Appendix B: Computations for \(\imath M_{\phi ,{\textrm{ren}}}^{\textrm{r}}(k)\), \(\imath M_{\phi ,{\textrm{ren}}}^{\textrm{a}}(k)\), \(\imath M^{+-}_{\phi }(k)\) and \(\imath M^{-+}_{\phi }(k)\)
In this appendix, we wish to find out the momentum space renormalised expressions of the retarded and advanced self energies as well as the self energies corresponding to the Wightman functions \(\imath M^{+-}_{\phi }(k)\) and \(\imath M^{-+}_{\phi }(k)\), as dictated by Eq. (31).
The renormalised expression for the retarded self-energy in coordinate space was found in Eq. (47). We shall now take the Fourier transform of it. We do it in two steps, for the sake of convenience of calculations. We first take the Fourier transformation of Eq. (47) with respect to its spatial part only
The last integral is given by the special function [69],
where \({\text {si}}(x)\) and \({\text {ci}}(x)\) are the sine and cosine integral functions, given by [69]
Using Eq. (58) into Eq. (57), we arrive at
Note that the \(\theta (\Delta t)\) appearing above commutes with one of the \((\partial ^{2}_{t}+|{\vec {k}}|^{2})\) operators because the term in square brackets is proportional to \((\Delta t)^2\) when \(\Delta t \rightarrow 0\). Accordingly, by taking one \((\partial ^{2}_{t}+|{\vec {k}}|^{2})\), we have
We next take the Fourier transform of the above equation with respect to time as well
where we have introduced \(\epsilon =0^+\) at necessary places, in order to regularise the integral. In order to evaluate the above equation further, we make use of the following equations [69]
Using these, Eq. (62) evaluates to
To evaluate the remaining integrals, we use, for real \(\alpha , \beta \), [69]
Putting these all in together, we obtain
The additive constant \(2\gamma _E\) appearing in the above equation can be further absorbed in a scalar field strength renormalisation counterterm as of Appendix A. Thus we have the final expression
We next wish to identify the self-energies \(\imath M^{+-}_{\phi }(k)\) and \(\imath M^{-+}_{\phi }(k)\), corresponding to the Wightman functions. The simplest way to achieve this is to use the retarded self-energy Eq. (67), and \(\imath M^{++}_{\phi ,{\textrm{ren}}}(k)\) derived in Appendix A. The relationship between \(\imath M^{++}_{\phi ,{\textrm{ren}}}(k)\), \(\imath M^{--}_{\phi ,{\textrm{ren}}}(k)\), \(\imath M^{\textrm{r}}_{\phi ,{\textrm{ren}}}(k)\), \(\imath M^{+-}_{\phi }(k)\) and \(\imath M^{-+}_{\phi }(k)\) is given by the renormalised version of the equation appearing below Eq. (30b),
Using now Eqs. (55), (56), (66) and (68), we find out the Wightman self-energies
Likewise, we can find out the renormalised advanced self-energy by using the renormalised version of the equation appearing below Eq. (29),
Using Eqs. (55), (56) and (69), we have
Being equipped with all these, we have computed the renormalised statistical propagator in momentum space as quoted in the main text, Eq. (48).
Finally, as a check of consistency, from Eqs. (29), (30), (67) and Eq. (71), we compute the Wightman functions,
Using
it is easy to see that in the limit \(g \rightarrow 0\), the above expressions agree with the free Wightman functions
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3. SCOAP3 supports the goals of the International Year of Basic Sciences for Sustainable Development.
About this article
Cite this article
Bhattacharya, S., Joshi, N. & Kaushal, S. Decoherence and entropy generation in an open quantum scalar-fermion system with Yukawa interaction. Eur. Phys. J. C 83, 208 (2023). https://doi.org/10.1140/epjc/s10052-023-11357-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-023-11357-6