Decoherence and entropy generation in an open quantum scalar-fermion system with Yukawa interaction

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.

B Computations for ıM r φ,ren (k), ıM a φ,ren (k), ıM +− φ (k) and ıM −+ φ (k) 17 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 Section 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 Section 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 Section 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 Section 3.1, we have renormalised the Kadanoff-Baym equations by renormalising the self-energy of the scalar field.Finally in Section 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 Section 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 − )-dimensional ( = 0 + ) Minkowski spacetime and will set c = = 1 throughout.

The basic setup
We consider a hermitian massive scalar field φ(x) coupled to fermions by the Yukawa interaction, Since we are working with the mostly positive signature of the metric, the γ µ 's satisfy the anti-commutation, As we have mentioned earlier, the scalar φ(x) will play the role of the system, interacting with the environment ψ(x), ψ(x).We assume that the environment is at zero temperature and is in its vacuum state.
Figure 1: Self energy diagrams for the scalar field at one and two loop order for the Yukawa coupling.The broaken lines denote fermions whereas the solid lines denote the scalar.We shall restrict our computations to O(g 2 ) only.See main text for discussion.
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, O(g 2 ).Note also that the two loop 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 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, π(x) = φ(x) is the momentum conjugate to the scalar field and the quantity F φ (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 ρ(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 where we have used the equal time limits of Eq. ( 2) and Eq. ( 4).In order to see the analogy of this construction with that of ordinary quantum mechanics, let us recall the generalised uncertainty relation (with = 1), (pure state) Combining the two above, we write Thus the quantity ∆ ≥ 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 ∆ 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 ∆ | 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.

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 nonequilibrium 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 ρ(t 0 ) (defined in the Heisenberg picture) is given by, where 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 The fields φ, ψ, ψ have their sources J φ , J ψ J ψ 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 δ-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 ∂ 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.

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 = ±, and ıS ss ψ are the fermion propagators.Γ (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 Γ (2) [ı∆ ss φ , ıS ss ψ ] = 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∆ φ and iS ψ , to obtain the equations of motion We next multiply Eq. (23a) and Eq.(23b) respectively by sı∆ s s φ (x , x ) and sıS s s ψ (x , x ) (from the right), and then integrate over x and sum over s = ±.Recalling that the integration over the product of the propagator and its inverse gives a δ-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 Solving these coupled algebraic equations, we can find out various two point functions.For example, on subtracting Eq. (28a) from Eq. (28c), we have the momentum space expression for the advanced propagator Eq. (17b), where the advanced self-energy ıM a φ (k) is given by where the retarded self-energy ıM r φ (k) is given by 29) into Eq.(30a) and Eq.(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, Eq. ( 5), Eq. ( 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 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.

The retarded self-energy and its renormalisation
In this section, we compute the self-energy ıM ss φ (x, x ) and find out the renormalised retarded self-energy, ıM r φ,ren (x, x ) = ıM ++ φ,ren (x, x ) − ıM +− φ (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 ıS ss ψ (x, x ) is obtained by acting i / ∂ on the massless scalar field propagator, The Poincarré invariant biscalar distance functions with appropriate i prescription, ∆x 2 ss (x, x ), necessary for the in-in formalism are defined as The one loop scalar self-energy ıM ++ φ (x, x ) is readily found from Eq. (25a) and Eq. ( 32) Similarly, we can find out the other self-energies ıM −− φ (x, x ), ıM +− φ (x, x ) and ıM −+ φ (x, x ) using the suitable i prescriptions as given in Eq. ( 33).We will now identify the divergence of Eq. (34).Note first that for an arbitrary exponent α, we have Furthermore, we can write For the other distance functions of Eq. ( 33), we have We now rewrite Eq. ( 34) using Eq. ( 35) and Eq.(36a) as where µ 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 ∂ 2 , we have to add a scalar field strength renormalisation counterterm (φ → (1 + δZ) 1/2 φ).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 δ-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 ıM ++ φ (x, x ) and ıM −− φ (x, x ), which contain divergent contributions, by the amputated counterterm contribution corresponding to δZ given above.
ıM −− φ (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), ıM +− φ (x, x ) and ıM −+ φ (x, x ) are determined by Note that ıM +− φ (x, x ) and ıM −+ φ (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 Eq. ( 41) and Eq. ( 42), we now have the renormalised retarded self-energy Using now ln ∆x 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 ∆t > 0 and ∆t 2 − ∆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 Appendix A and Appendix 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.

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 Appendix A, Appendix B. Substituting now Eq. ( 67), Eq. ( 69) and Eq. ( 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 /µ in Fig. 4 and Fig. 5.Note that as the Yukawa coupling g gets smaller, F φ (k) approaches a δ-function dispersion, as expected from Eq. ( 49).However as g increases, the δ-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/µ, the peak gets shifted towards the higher values of k 0 /µ 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, Eq. ( 5), Eq. ( 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 → g θ(t − t 0 ).The step function changes the limit of integration Eq. ( 26) from −∞ to t 0 Eq. ( 49), with respect to the dimensionless variable k 0 /µ, where we have taken | k|/µ = 1, m/µ = 1 and the Yukawa coupling strengths g = 0.5 (black curve), g = 1.3 (red curve).Eq. ( 49), with respect to the dimensionless variable k 0 /µ, where we have taken | k|/µ = 1, m/µ = 1.3 and the Yukawa coupling strengths g = 0.5 (green curve), g = 1.3 (blue curve).Comparison with Fig. 4 shows lesser numerical value of the statistical propagator with increasing m/µ.See main text for discussion.
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 eg.[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 | k|/µ = 1, m/µ = 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/µ, by fixing all the other parameters in Fig. 6.Thus ∆ decreases monotonically with increasing m/µ 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) (| k|/µ = 1, m/µ = 2 and g = 0.5) S ≈ 1.27 (52) We have also analysed the variation of the entropy with respect to the system mass by fixing all the other parameters have taken | k|/µ = 1 and, g = 0.5 (blue curve), g = 1 (red curve) and g = 1.5 (green curve).As of the phase space area, Fig. 6, the entropy also decreases monotonically with increasing m/µ.See main text for discussion.
in Fig. 7. Due to the aforementioned reason as that of the phase space area, the entropy is decreasing with the increasing m/µ.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.

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 Section 3. In 3.1, we renormalise our results and finally compute the phase space area, the statistical propagtor and the entropy in Section 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., Section 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 correllators 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.
A Calculation of renormalised ıM ++ φ,ren (k) and ıM −− φ,ren (k) In this appendix, we compute the renormalised self-energies ıM ++ φ and ıM −− φ 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 ıM −− φ (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.
B Computations for ıM r φ,ren (k), ıM a φ,ren (k), ıM +− φ (k) and ıM −+ φ (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 ıM +− φ (k) and ıM −+ φ (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 ıM r φ,ren (| k|, t, t ) = The last integral is given by the special function [69], where si(x) and ci(x) are the sine and cosine integral functions, given by [69] si Using Eq. ( 58) into Eq.( 57), we arrive at where we have introduced = 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] Putting these all in together, we obtain The additive constant 2γ 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 ıM r φ,ren (k) = We next wish to identify the self-energies ıM +− φ (k) and ıM −+ φ (k), corresponding to the Wightman functions.The simplest way to achieve this is to use the retarded self-energy Eq. ( 67), and ıM ++ φ,ren (k) derived in Appendix A. The relationship between ıM ++ φ,ren (k), ıM −− φ,ren (k), ıM r φ,ren (k), ıM +− φ (k) and ıM −+ φ (k) is given by the renormalised version of the equation appearing below Eq. (30b), Using now Eq. ( 55), Eq. ( 56), Eq. ( 67) and Eq. ( 68), we find out the Wightman self-energies Likewise, we can find out the the renormalised advanced self-energy by using the renormalised version of the equation appearing below Eq. ( 29 Being equipped with all these, we have computed the renormalised statistical propagator in momentum space as quoted in the main text, Eq. (48).

Figure 2 :
Figure 2: Schwinger-Keldysh contours with finite and infinite initial and final times.

Figure 3 :
Figure 3: Contributions to the 2PI effective action up to two loop order.The solid lines denote φ-propagators, whereas the dashed lines correspond to ψ-propagators.The tadpoles will not make any contribution.See main text for discussion.

Figure 4 :
Figure 4: Variation of the statistical propagator,

Figure 5 :
Figure 5: Variation of the statistical propagator,

4 ΔFigure 6 :Figure 7 :
Figure6: Variation of the dimensionless phase space area ∆, with respect to the dimensionless mass parameter m/µ of the system, where we have taken | k|/µ = 1 and, for g = 0.5 (blue curve), g = 1 (red curve) and g = 1.5 (green curve).See main text for discussion.