Time-convolutionless cosmological master equations:
Late-time resummations and decoherence for non-local kernels
Suddhasattwa Brahma1***e-mail address: [email protected], Jaime Calderón-Figueroa2†††e-mail address: [email protected] and Xiancong Luo1‡‡‡e-mail address: [email protected]
1Higgs Centre for Theoretical Physics, School of Physics and Astronomy,
University of Edinburgh, Edinburgh EH9 3FD, UK
2Astronomy Centre, University of Sussex, Falmer, Brighton, BN1 9QH, UK
Abstract
We revisit a simple toy model of two scalar fields in de Sitter space, playing the roles of “system” and “environment” degrees of freedom, which interact with each other. We show that there are secular divergences in physically relevant observables which arise solely from the non-Markovian part of the memory kernel, contrary to popular belief that secular growth typically comes from local terms in evolution equations. Nevertheless, we show that these terms can still be non-perturbatively resummed, using the time-convolutionless master equation formalism, which improves upon previous approximations. At the same time, there are other physical quantities in the same model that are dominated by local terms in the memory kernel. Therefore, we conclude that, for cosmological backgrounds, either the dissipation or the noise kernel can end up being dominated by non-local terms depending on the nature of the system-environment coupling.
1 Introduction
Inflation – the most widely accepted paradigm for the early-universe – predicts that the large scale structure of the cosmos has its origins in the quantum vacuum fluctuations of the scalar field responsible for the background accelerated expansion [starobinsky1980new, fang1980entropy, Guth:1980zm, Sato:1981qmu, Linde:1981mu, Albrecht:1982wi]. However, a deeper dive into this formalism shows that one needs to break up the same inflaton field into the homogeneous background, sourcing the evolution of the universe, and its quantum fluctuations, responsible for structure formation [Mukhanov:1982nu]. A natural question that follows is what is the backreaction [brandenberger2002back, PerreaultLevasseur:2014ziv] of the sub-horizon quantum fluctuations on the effective field theory (EFT) of the ‘classicalized’ super-Hubble modes [Cohen:2020php, burgess2015eft]? Another way to phrase this problem is as follows: The equal-time in-in correlations of a massless scalar field in quasi-de Sitter (dS) space give us the spectra of the adiabatic perturbations which can be verified in the temperature fluctuations in the CMB [collaboration2020planck]. Nevertheless, issues of infrared (IR) singularities and late-time secular divergences [burgess2015eft, Burgess:2015ajz, Cespedes:2023aal] are know to afflict such computations when considering loop corrections [senatore2010loops, Assassi:2012et].
These questions are typically tackled using Starobinsky’s stochastic formulation of inflation[Starobinsky:1979ty, Morikawa:1989xz, nambu1988stochastic, nambu1989stochastic, kandrup1989stochastic, nakao1988stochastic, mollerach1991stochastic, linde1994big, starobinsky1994equilibrium]. The general idea behind it is that gradient terms are expensive once the wavelength of a perturbation becomes super-Hubble, and thus its evolves essentially independently outside the horizon. This allows one to write down a Fokker-Planck equation for the probability distribution of the scalar field, coarse-grained over each Hubble volume. Its evolution is a balance between the classical drift due to the quasi-dS scalar potential and ‘quantum kicks’, modeled by a Gaussian white noise, due to the sub-horizon fluctuations. Since it offers a gradient expansion for the super-Hubble density perturbation, instead of the one in terms of its amplitude, it can make predictions regarding the tail of its probability distribution where linear perturbation theory is bound to fail. And in this way, it can make nontrivial predictions for, say, the production of primordial black-holes during inflation (see, for instance, [Vennin:2020kng, vennin2024quantum, Panagopoulos:2019ail]).
Stochastic inflation is supposed to also be precisely the non-perturbative formalism designed to cure the IR issues of massless fields during inflation. Since it is expected to provide the leading order effect in the resummation of late-time IR divergences in cosmology111For partial results in this direction, see [Seery:2007we, Enqvist:2008kt, seery2009parton, Burgess:2009bs, Seery:2010kh, Gautier:2013aoa, Guilleux:2015pma, Gautier:2015pca, Hardwick:2017fjo, Markkanen:2017rvi, LopezNacir:2019ord, Gorbenko:2019rza, Mirbabayi:2019qtx, Adshead:2020ijf, Moreau:2020gib]., its assumptions need to be carefully examined. For instance, one of the primary approximations for stochastic inflation is made in the computation of the quantum noise due to the short wavelength fluctuations. Generally, it is typically assumed to be white [Finelli:2008zg, Onemli:2015pma, Calzetta:1995ys, Garbrecht:2013coa, Vennin:2015hra, burgess2015eft, Baumgart:2019clc, Baumgart:2020oby], although ‘colored’ corrections are expected to come from non-Gaussian contributions [Cohen:2021fzf, Cohen:2021jbo]. While it has been argued that both the choice of the ‘sharp’ window function [Mahbub:2022osb, Andersen:2021lii] (to demarcate the short and long wavelength scales) as well as the ‘test scalar’ approximation [Boyanovsky:1994me] can invalidate it, nevertheless this ‘Markovian’ assumption222At the very outset, we want to emphasize that we do not use Markovian and time-locality interchangeably in this draft, as is sometimes erroneously done in the cosmology literature. Also, we loosely say that a white noise implies a Markovian system in the classical sense even though we will define what we mean by Markovianity in a quantum setup as having a semi-group property later on [Prudhoe:2022pte]. of having a white noise for the integrated-out quantum degrees of freedom (dofs) remains a ubiquitous one in cosmology (e.g., [Mirbabayi:2020vyt, Pattison:2019hef, Gong:2019yyz]) with a few notable exceptions [shandera2018open, Pinol:2020cdp, Zarei:2021dpb, Figueroa:2021zah].
What is even more, this assumption in stochastic inflation has often been used as an argument for inferring about the origin of secular terms that are to be expected when computing loop corrections during inflation. There is a rich literature on computing one-graviton loop corrections to, say, a conformally coupled scalar field propagator in dS space [Glavan1, Glavan2, Boran:2014xpa] in which it has been explicitly shown that the secular terms are the ones which come from time-local terms in the evolution equation whereas the non-local terms decay and do not contribute at late-times. Similar results have also been obtained for the one-loop graviton correction to the vacuum polarization for dynamical photons in a dS background [wang2015excitation]. Such secular enhancement is generally found to independent of the nonlocal terms in a variety of models [Glavan3, Glavan4, Glavan5, Degueldre:2013hba, miao2006gravitons, kahya2007quantum, Prokopec:2002uw], and is hence free of any memory effects. These results lend support to the general expectation that stochastic inflation captures all the leading IR secular terms via a local equation, and thus the Markovian approximation is a sufficient one for cosmology.
In this paper, we will revisit an old model [Boyanovsky:2015tba] to show that indeed secular divergences can just as easily come from the non-local part of the memory kernel – an object that captures correlations of the integrated-out fields. More concretely, we will show explicitly that secularly divergent terms come from non-local terms in the (dissipation) part of the kernel and, consequently, that it is not necessary that late-time behaviour is always dominated by local terms. More importantly, we shall demonstrate that such secular divergences, arising from non-Markovian terms, can nevertheless be exactly resummed, without invoking any additional, ad hoc approximations, using an appropriate master equation formalism. In our example, we shall give the physical interpretation of different co-efficients of the master equation, corresponding to different parts of the kernel, which will be dominated by local and non-local terms at late-times. The main purpose of our work is thus to show that:
-
1.
Secular divergences can just as easily stem from non-Markovian terms and, more importantly, such terms can still be resummed at late times following a precise algorithm that does not involve any arbitrary approximations.
-
2.
The memory kernel, corresponding to the integrated-out (or coarse-grained) fields, in the same model can affect different physical quantities differently. More specifically, local and non-local parts of the kernel can become dominant for different physical observables.
In the next section, we briefly discuss master equations in cosmology before setting up the one relevant for us in Sec-3. We discuss physical interpretations of coefficients of the master equations as they pertain to dissipative corrections to cosmological observables and diffusion terms leading to decoherence of the quantum fields. In Sec-4, we compute corrections to the power spectrum and the purity of the (system) density matrix and show how they are dominated by the non-local and local parts of the kernel, respectively. Finally, we conclude by summarizing our findings and outlining their applications to future work.
2 Open quantum systems in cosmology
To achieve the above objectives, we shall adopt an open EFT approach [calzetta2009nonequilibrium, breuer2002theory] to the problem, an idea that has been gaining considerable importance over the last decade in cosmology [Burgess:2022rdo, burgess2016open, Colas:2022hlq, Colas:2022kfu, Colas:2023wxa, Colas:2024lse, Colas:2024xjy, Colas:2024ysu, Salcedo:2024smn, Boyanovsky:2015xoa, Boyanovsky:2015jen, Boyanovsky:2015tba, hollowood2017decoherence, shandera2018open, Choudhury:2018ppd, Brahma:2020zpk, Brahma:2021mng, Brahma:2022yxu, Burgess:2022nwu, banerjee2023thermalization, Kaplanek:2022xrr, Cao:2022kjn, Prudhoe:2022pte, Kading:2022hhc, Kading:2022jjl, Kading:2023mdk, Alicki:2023rfv, Alicki:2023tfz, Creminelli:2023aly, Keefe:2024cia, Bowen:2024emo, crossley2017effective, Colas:2022hlq, Burgess:2015ajz]. The main realization here is that open quantum systems are particularly suited to spacetimes with horizons, such as the one for accelerated expansion. There is a sector of the Hilbert space that is hidden from the observer, thus being dubbed “environment” , which can still exchange energy and information with the “system” sub-sector. The usefulness of using an open quantum system is that it is able to coarse-grain the effects of the high energy or sub-Hubble dofs on the dynamics of the long-wavelength modes when there is no exclusion principle (such as the conservation of energy) available.
Hence, the master equation aims to capture the evolution of the system dofs by integrating (or tracing) out the environment modes. One of the main advantages of this formalism is that it can capture reliable late-time dynamics of the system by overcoming the secular growth associated with time-dependent backgrounds in cosmology. The underlying reason for secular growth at late times is that gravitational interactions are universal, and once turned on, it is not possible to isolate the system from the ever-present medium, namely, gravity. This implies that even in the weakly-interacting regime, for long time-intervals small effects can accumulate to lead to the failure of perturbation theory. Resummation of such late-time divergences typically implies deriving the master equation at some order in perturbations theory, say , and then treating it as the bona fide dynamical evolution equation that needs to be solved non-perturbatively without performing any expansion [Colas:2022hlq]. This means that the master equation would contain all term at and some terms at .
The details of much of this discussion can be found in, say [breuer2002theory, Colas:2022hlq], and we only very briefly touch upon the time-convolutionless master equation formalism that we will employ for our model. The density matrix, corresponding to the entire system undergoes unitary evolution governed by the Louville-von Neumann equation:
(1) |
where the interaction Hamiltonian has operators corresponding to both and . For the sub-sector, this equation can be put in an exact, closed form which determines the evolution of the reduced density matrix corresponding to , known as the Nakajima-Zwanzig master equation. However, it is an integro-differential equation that has exactly the same level of difficulty as the Louville-von Neumann equation, and thus needs to be solved using some approximation scheme.
The first of these is the so-called Born approximation, which assumes a weak coupling between and and expands in this interaction paarameter . One typically also assumes that the initial density matrix is in a factored state: . Although this is not an essential assumption, it is a natural one for cosmology since we will start off both the system and environment modes in the Bunch-Davies vacuum state. At , this results in a master equation of the form:
(2) |
Since the RHS of the above equation has depending on , rather than on alone, this equation is not amenable to late-time resummations. This is because assuming that the master equation generates dynamics non-perturbatively, even though it is derived at , works only because it is not dependent on some specific value of time. The next popular approximation is the Markovian one where the above master equation reduces to the standard Lindblad form. However, in this work, we shall not invoke this restrictive limit which basically assumes that the master equation generates dynamics that fulfills a semi-group property, i.e., where .
Instead, we will work with the time-convolutionless approximation, based on a cumulant expansion (see [breuer2002theory] for details). In this case, one can rewrite the master equation, at leading order (hereafter the TCL2 master equation), as
(3) |
Thus, at this order, the only difference is that in the TCL2 master equation is evaluated at instead of . It can be shown that the terms dropped in this equation, as compared to (2), is at or higher. More importantly, the TCL formalism allows one to capture non-Markovian effects even though we end up working with a time-local equation. And this is exactly one of the technical issues that we wish to clarify in this work, and correct some associated errors in the literature on this topic.
Note that the reader might think that the substitution of in the above master equation is nothing but an application of the Markovian approximation to the system. Indeed, this line of thinking has been mistakenly employed before [Boyanovsky:2015tba, hollowood2017decoherence]. One of the goals of this work is to underscore that despite the (time-)local nature of the TCL2 equation, it remains capable of adequately describing a non-Markovian system. However, a subtlety arises in the form of a residual artifact of the “memory” of the initial state remains in the lower limit of the integral on the RHS. In other words, when we compute the coefficients of the TCL2 master equation, in terms of expressions quantifying diffusion and dissipation, we shall have to take particular care to drop such terms with explicit dependence on such initial time . Such terms have been called “spurious” earlier in [Colas:2022hlq], and it will be shown that the TCL2 cannot be resummed unless such terms are dropped by hand. Going beyond these technical aspects, the main message is that the system described by the TCL2 master equation is not a Markovian one and we will show how previous results were only approximate since it did not involve a proper analyses of such spurious terms, and how we improve upon those results in our work.
3 A toy-model for sub- and super-Hubble mode-coupling
Our model consists of a system comprising a minimally coupled scalar field () and an environment composed of a conformally coupled scalar field (), with an interaction term which has previously been studied in [Boyanovsky:2015tba, Boyanovsky:2015jen, hollowood2017decoherence]. The reason behind studying this model is two-fold. Firstly, the associated memory kernel of the environment field has the advantage of cleanly splitting into two different contributions, one that is clearly time-local while the other non-local, as we will show later on. Secondly, the conformally-coupled “environment” field is a good proxy for the sub-Hubble modes since the latter are least influenced by the background curvature. On the other hand, the minimally coupled “system” scalar would also be assumed to be massless, thereby playing the role of the adiabatic degree of freedom. Therefore, this setup can play the role of a toy-model for stochastic inflation. Nonetheless, it is important to keep in mind that this is still a model of two interacting test scalars in a dS background: where is the scalar factor, the Hubble parameter and denotes the conformal time.
The action, written in terms of (canonical) Mukhanov-Sasaki fields, is given by
(4) |
where primes denote derivatives with respect to conformal time, . From the above action, we can obtain the system and environment Hamiltonians, together with the interaction potential, as:
(5) | ||||
(6) |
where is the momentum conjugate of , and is that of , respectively. From here on, operators in the Schrödinger picture will be denoted with hats, while interaction picture operators will be denoted with tildes. Following [Boyanovsky:2015tba], the interaction term has been normal-ordered as , where is the initial (product) density matrix. Given the free Hamiltonians, the Heisenberg equations of motion can be evaluated as
(7) |
with solutions, assuming Bunch-Davies initial conditions, given by
(8) |
and similary for their conjugate momenta:
(9) |
Consequently, the quantum operators for the system field (and its conjugate momenta) read
(10) | ||||
(11) |
where and are the annihilation and creation operators at , respectively.
Finally, let us express these fields at arbitrary times in a time-local way, as required in the TCL framework. For this, we ‘solve’ (10) and (11) for and at , and subsequently substitute these solutions back into the corresponding equations at , yielding
(12) | ||||
(13) |
Plugging the mode functions above, one obtains
(14) | ||||
(15) |
3.1 Building the TCL2 master equation
Following the recipe outlined in Sec-2, some of the details of which can be found in [Colas:2022hlq, us:QSR], we can write down the TCL2 master equation as333We have used and interchangably in this paper.
(16) |
where the memory kernel is given by444To write the solution of the integral in the quoted form, one invokes the Sokhotski–Plemelj theorem, which states where denotes the principal part which shall be expressed more explicitly later on.
(17) |
Given the system field appears at most at quadratic order in the the full Lagrangian, the evolution equation for the density matrix can be written as a sum over independent momentum modes without any mode-coupling. This is another technical advantage of choosing the type of coupling that is being considered in this model (as opposed to, say, a coupling term).
It will be more compact to express the field and its conjugate momentum as coordinates of phase space, denoted as , . Using (12),(13), the resulting master equation can be written as
(18) |
where all the operators are evaluated at , and the functions are given by the following integrals:
(19) | ||||
(20) | ||||
(21) | ||||
(22) |
We leave the solutions of these integrals for the next section.
At this point, let us consider the lower limit of the (formal) antiderivative of the above functions. Following the convention adopted in [Colas:2022hlq], we can rewrite the RHS of each of the above equations as , the abstract index denoting the different functions for the four different quantities above. The second term, coming from the lower limit of the integral, has the memory of the initial state and is the one that has been termed ‘spurious’ in [Colas:2022hlq]. In the Appendix, we compute these terms explicitly and show that they dominate at late-times near the conformal boundary. What this implies is that keeping such terms would not allow for the TCL2 master equation to be resummed. Moreover, as has also been noticed in [Colas:2023wxa], such terms are absent in a perturbative treatment. In other words, if one were to solve the master equation by expanding the density matrix to a fixed order in , then these terms get cancelled amongst themselves. However, when solving the TCL2 equation non-perturbatively by considering it as the authentic dynamical generator, we only keep some of the terms at orders higher than , and therefore such spurious terms need to be removed by hand. However, getting rid of these spurious terms is not the same as imposing a Markovian approximation. As will be shown later on, some of the master equation coefficients in (19)-(22) have contributions solely from the non-local part of the memory kernel and hence, by definition, has a non-Markovian origin. It is important to point out that previous papers where this model has been studied lacks any discussion of these spurious terms and we shall explain later on in the next section why their results are, therefore, necessarily approximate when considering non-perturbative resummations.
With these expressions at hand, the master equation is written in a rather economical and suggestive way,
(23) |
where is the antisymmetric matrix
(24) |
and and are symmetric matrices, whose entries are given by (19)-(22). Noticeably, there is no term present, a consequence of the specific form of the interaction involving only configuration field-field coupling. On the other hand, in interactions involving time derivatives of the fields, one anticipates the emergence of such terms (as that would involve a field-momentum coupling [Colas:2022kfu]). The master equation, presented in this manner, bears resemblance to stochastic equation describing (quantum) Brownian motion, such as the one obtained from the Caldeira-Leggett model [Colas:2022hlq, Bhattacharyya:2024duw].
For practical purposes, it is more useful to work in the Schrödinger picture where the master equation takes the form:
(25) |
with the effective quadratic Hamiltonian given by
(26) |
Finally, we can rewrite the master equation in a Lindblad–like form as
(27) |
where is known as the dissipator matrix.
3.2 The dissipator matrix: Dissipation and Diffusion
Let us focus on the coefficients of the master equation, some of which appear in the effective Hamiltonian and the others as entries of the dissipator matrix. Before presenting the explicit form of these expressions, let us delve into the physical implications of these terms. The real and imaginary part of the memory kernel correspond to the noise and the dissipation kernel, respectively555For thermal states, the real and imaginary parts are related to each other via the fluctuation-dissipation theorem. Such dynamics takes place in scenarios like warm inflation [Berera:1995wh, Berera:1995ie].. In terms of the RHS of the dissipator matrix written above, corresponds to the noise kernel while the dissipation kernel, respectively. Both these terms signify non-unitary contributions to the reduced density matrix due to interactions with the environment field. On the other hand, in (26) corresponds to a unitary term renormalizing the energy spectrum of the system Hamiltonian. Since the term is responsible for diffusion, it leads to decoherence of the quantum modes associated with the system dofs. On the other hand, the is the dissipative or friction term that controls the amount of energy transfer between the system and the environment and would be reflected in the equal-time correlation functions, e.g., the configuration field power-spectrum. ( plays the dual role of renormalizing the comoving Hubble parameter.) The utility of our model is that it would be clear which of these quantities are dominated by the local and the non-local parts of the kernel, and thereby which physical observables would reflect the non-Markovianity of the model.
Keeping in mind the physical significance of , and , we wish to find their explicit form using the integrals in (19)-(22). Although intricate, these integrals can indeed be evaluated analytically for this model. Particularly delicate are the terms arising from the principal value in the kernel (see (17)), i.e., from the time non-local term. One strategy to address this, and to isolate the associated UV-divergence, is to express the principal value as:
(28) |
and subsequently expand around . We present the result of these integrals derived through this method, while relegating the expressions of the spurious functions, originating from the lower limits of the integrals, to the Appendix. In the following, we have highlighted in blue those terms arising from the local part of the kernel while the rest all originate from the non-local bit. Furthermore, the UV-divergent piece (also coming from the non-local part of the kernel) has been marked in red. The resulting expressions are:
(29) | ||||
(30) | ||||
(31) | ||||
(32) |
We show the evolution of these functions in Figs. 1 and 2. The depicted plots exclusively come from the finite terms where, as promised, we have simply dropped the spurious contributions for now. Note that we have appropriately rescaled these quantities with the correct factors of to make them dimensionless. For generating the plots and performing subsequent numerical analyses, we have selected values of such that the system evolves for 20 foldings before horizon crossing. We verified that extending this time does not alter our results. denotes the value of the scale-factor when the mode crosses the horizon; thus, the zero of the horizontal axis denotes horizon-crossing in the following plots.
![Refer to caption](x1.png)
![Refer to caption](x2.png)
![Refer to caption](x3.png)
![Refer to caption](x4.png)
Let us first take a close look at the expression for as given in (29). The term coming from the local part of the kernel apparently looks subdominant () to the nonlocal part (with a piece going as ). However, on tracking the evolution of these terms separately, it is clearly seen that the time-local contribution consistently dominates, becoming more pronounced over timescales of a few foldings after horizon exit (see Fig. 3). One can show that the Taylor expansion of around is:
(33) |
What this indicates is that one of the diffusion terms () is dominated by the contribution coming from the local part of the memory kernel (at all times and certainly after horizon-crossing). This will be important for our discussion about decoherence later on. On the other hand, only has contributions from the non-local part of the kernel since the term coming from the local bit turns out to be zero for our model. Next let us look at the expression of and from (30) and (32), respectively. These have contributions coming solely from the non-local part of the memory kernel which shows that the dissipation in this model is necessarily non-Markovian. However, does contain a UV-divergent term that needs to be taken care of by invoking a mass counter-term. Note that this implies that this mass term has to counteract long-time correlations such that they get screened. Or, conversely, we might say that the non-local origin of this UV-divergent term means that the time non-local history dependence of the correlations in the environment induce a mass term for the system. does not have any such UV-divergent piece. Finally, note that both and also display secular divergence even though these now have a non-local origin.
![Refer to caption](x5.png)
4 Cosmological Observables
The purpose of tracking the time evolution of the (reduced) density matrix through the master equation is to compute observables, particularly equal-time -point correlation functions. In this work, we primarily focus on 2-point correlators, where every relevant combination is encapsulated by the following matrix:
(34) |
where each term on the RHS is to be understood as the -mode of their respective field. Observables are identified as the expectation values of these operators, represented by . This is known as the covariance matrix, and the differential equations governing the evolution of its entries are dubbed transport equations.
4.1 The transport equations
To derive the transport equations, let us begin by considering the expectation value of an arbitrary operator , which in the Schrödinger picture is given by
(35) |
Then, the evolution equation of the expectation value can be found through
(36) |
Using the TCL2 master equation in any of its forms, say (25), and leveraging the cyclic property of the trace, we can cast (36) as a combination of expectation values (underscoring the importance of a time convolutionless master equation). In doing so, and repeatedly using
(37) |
where , we obtain
(38) |
The first two terms on the RHS come from the unitary evolution, while the latter two are associated with diffusion and dissipation, respectively, as explained in the previous subsection. For parity conserving Hamiltonians, one can safely assume , yielding the simple set of transport equations
(39) |
or, in matrix form,
(40) |
where it is understood that all of the terms are computed at a specific momentum mode .
4.1.1 A simple consistency check
Before diving into solving the transport equations – something that can only be achieved through numerical simulations – let us verify that the equations yield sensible expressions for the free theory, where analytical exact results are readily available. For such a scenario, we know the elements of the covariance matrix, which are merely combinations of the mode functions described in (8) and (9), and their complex conjugates. Upon calculation, we find:
(41) |
which perfectly matches with the free-theory transport equations:
(42) |
On a cursory glance, the transport equation for , namely the configuration field () power spectrum, looks remarkably similar in structure to the analogous equation in (40). However, this resemblance does not imply that interactions have no impact on the power spectrum. Quite the contrary, the influence of interactions is encapsulated in the coupling of to other equations, where the dissipative and diffusion coefficients appear explicitly.
4.2 Numerical Solutions for the covariance matrix
As previously mentioned, the system of transport equations (40) requires numerical techniques to obtain a solution. Alternative approaches, like the one implemented in [Boyanovsky:2015tba, Brahma:2021mng], involve approximations in the super-Hubble regime by ignoring the contribution of the so-called decaying modes to the system. As a by-product, we shall show that these approximations render results in excellent agreement with the exact (numerical) results obtained here.
Figs. 4-6 depict the solutions to the transport equations for a range of couplings. In our code, we simulated 20 folds of evolution in the sub-horizon regime, ensuring that extending this duration yields consistent final results once the mode has long exited the horizon and is in the super-Hubble regime. As before, we have rescaled the quantities – in this case, the covariance matrix elements – to make them dimensionless. Additionally, the value of was chosen based on the same reasoning discussed earlier, and is the instant of horizon-crossing.
![Refer to caption](x6.png)
![Refer to caption](x7.png)
![Refer to caption](x8.png)
4.2.1 Non-perturbative resummation and spurious terms
One aspect of our solutions that might not be immediately apparent is the inherent non-perturbative resummation of any potential secular divergences. This feature has long been recognized as one of the advantages of using master equation techniques for open quantum systems [Chaykov:2022zro, burgess2008decoherence, Burgess:2018sou, Kaplanek:2019dqu, Brahma:2021mng, Boyanovsky:2015tba, Boyanovsky:2015jen, Burgess:2015ajz, Colas:2022hlq, Colas:2024xjy, Burgess:2024eng]. An analytical scheme which has been used to find resummations for master equations follows a super-horizon approximation which, in effect, calls for dropping all terms proportional to the decaying mode outside the horizon [Boyanovsky:2015tba, Brahma:2021mng]. Specifically, on applying this to this model resulted in an amplitude for the spectrum of the configuration field (, for a given ,) in the approximate form [Boyanovsky:2015tba, Boyanovsky:2015jen]:
(43) |
where
(44) |
and represents the correlation between the free–theory growing modes, and is given by at horizon exit. Further, represents a renormalized mass, whose origin (in our framework) is the counterterm that subtracts the divergent term in (see (30)). The double exponential ensures a finite result since the second exponential –– having a quadratic dependence on the number of folds –– will dominate at late times.
To compare our results with the above analytical expression, must be fully determined. The form presented in (44) arises from equating the renormalization scale to the initial time at the beginning of inflation when the density matrix gets prepared (). The reason behind this choice shall be discussed shortly. Our treatment of the divergence in differs in that we chose a counterterm that precisely cancels out the infinity, leaving only the finite parts shown in (30). Consequently, no term analogous to contributes to our result, leading us to drop it and consider for comparison purposes666In the final expression in [Boyanovsky:2015tba], to provide numerical estimates, is chosen such that the mode would have exited the horizon 50 e-folds before the end of inflation, i.e., . On the other hand, the renormalization scheme there had eliminated the constant term in the exponent.. Although our implementation of the renormalization procedure varies slightly from previous results [Boyanovsky:2015jen] – mainly due to the fact that a numerical, yet exact, computation that we have executed here would have been significantly difficult otherwise – note that this small differences appear in the term linear in and should not matter for the actual resummation which follows from the double exponential.
With this caveat in mind, we quantify the consistency of our results with those of [Boyanovsky:2015tba] by defining the relative deviation as:
(45) |
where is given by (43), and is the corresponding result obtained with our approach. This function is plotted in Fig. 7, showing a remarkable agreement for small , with expected discrepancies for , the limit in which both the TCL2 approximation and the perturtabative approach in [Boyanovsky:2015tba] break down. To provide further evidence with respect to the consistency of both approaches, we have fitted our results to a double exponential function as described in (43), which ultimately forms an exponential of a quadratic function of . The outcomes of this fitting process are summarized in Table 1, focussing on the fit for the quadratic function as it is the most important one for resummation purposes.
Several insights can be drawn from these fits. Most importantly, the quadratic coefficient in aligns closely with the prediction in [Boyanovsky:2015tba, Boyanovsky:2015jen], even for the strongest couplings. As mentioned above, this coefficient is crucial for resummation, indicating that both approaches are comparable in this regard. The reason behind the discrepancy for the linear term, due to the slightly different ways of implementing the renormalization schemes, has been explained above. Also, it is worth noting that the analytical and numerical approaches predict almost identical coefficients for the term independent of any foldings (). For the former, this value corresponds to the logarithm of the last term in (43), , which aligns well with the fitted functions.
![Refer to caption](x9.png)
Having shown that the super-horizon approximation that leads to an analytical non-perturbative resummation works quite well when compared to our exact numerical solution, we are now in a position to explain why that is the case and also why does one need an additional approximation in the former case at all. Recall that there was no proper analysis of the spurious terms done in [Boyanovsky:2015tba, Brahma:2021mng], and yet the results of such a resummation procedure compares quite well to our exact numerical solutions. The reason is the judicious choice of the renormalization scale which was identified with in order to minimize the effects from the lower limits of the time integrals involved in the calculations. Although this process is reminiscent to removing spurious terms, it does not quite completely eliminate the contribution of the lower limit of the integral to the master equation. To do so, one needs to implement a further approximation, that of dropping all the decaying modes in the final result. In this way, all contributions of the spurious terms were effectively eliminated in that approach. Although this works fine if one is concerned with computing the value of an observable at the final conformal boundary, the approximation becomes progressively worse at earlier times, when evaluating the evolution of the observable over time as one goes away from the limit. On the other hand, in the TCL2 approach, one simply identifies terms which would cancel in a perturbative approach and contain remnants of the initial state, and drops them from the very beginning. In this way, no further approximations are necessary to derive dynamics of cosmological observables. Moreover, apart from not having to impose additional assumptions, the TCL resummation has been shown to better match with exact results for solvable models [Colas:2022hlq] than the super-horizon one.
4.3 Purity as the determinant of the covariance matrix
Beside computing the in-in correlation functions, open EFTs are key to studying non-unitary effects such as decoherence (see [Colas:2024xjy, brandenberger1990classical, Calzetta:1995ys, barvinsky1999decoherence, lombardo2005decoherence, lombardo2005influence, martineau2007decoherence, prokopec2007decoherence, burgess2008decoherence, sharman2007decoherence, campo2008decoherence, anastopoulos2013master, Nelson:2016kjm, martin2018non, Oppenheim:2022xjr, DaddiHammou:2022itk, Sharifian:2023jem, Ning:2023ybc, joos2014decoherence, Burgess:2024eng, Bhattacharyya:2024duw, Colas:2024lse, Burgess:2022nwu, Micheli:2023qnc, Ning:2023ybc, Kaplanek:2019dqu, DaddiHammou:2022itk, hollowood2017decoherence, burgess2008decoherence, Chandran:2018wwc, Bhattacharya:2022wpe, Bhattacharya:2023twz, Bhattacharya:2023xvd] for a partial list). In particular, in this paper, we will consider the behaviour of a specific quantum informatic measure, namely purity, for each mode777As mentioned earlier, for our model different system modes are decoupled from each other at the level of the TCL2 master equation, which greatly simplifies the task at hand., which signals the amount of quantum entanglement between the system and the environment. By analysing the purity of the reduced , we can study the decoherence (and possibly more exotic phenomena of recoherence or purity-freezing, as a consequence of non-Markovianity [Colas:2022kfu, Colas:2024xjy, Colas:2024ysu]) of the system quantum field.
The purity of the system, for a particular momentum mode, is defined as:
(46) |
where is the reduced density matrix, corresponding to the system dofs, as before. When , the system is in a pure state and there is no entanglement between system and environment. When , the system density matrix is maximally mixed. As a side-note, recall that the purity is indeed independent under field-redifinitions as the trace indicates. From hereon, we will drop the subscript to avoid clutter since the evaluation can be done mode by mode.
One can simply solve the master equation (3.1) to solve for the reduced density matrix and use that to compute the evolution of the purity. However, an equivalent step would be to use the covariance matrix to do so. More specifically, it is known that for Gaussian states, we can use the covariance matrix of the system sector 888We introduce the superscript here for the covariance matrix to emphasize that this corresponds to the system field alone and is not the full covariance matrix. to calculate the purity of the system [Ferraro:2005hen]:
(47) |
Solving from the master equation by truncating at the TCL2 order is the same as using the covariance matrix to compute the purtiy, while ignoring contributions from higher order correlators, and thus our approximation is self-consistent.
Although we can solve for each of the components of the covariance matrix individually and plug them into the above equation, in practice, it is simpler to derive the equation for the determinant of covariance matrix from (39) as
(48) |
We choose initial conditions such that the field is in the Bunch-Davies vacuum. Explicitly, these are given by , which implies that as .
With these initial conditions, we solve for the above differential equation for the determinant of and use that to solve for . For a mode that spends e-folds before horizon exit, we solve the above equation and compare the dynamical behavior of purity for different values of as a ratio over . As shown in Fig. 8, the system evolves rapidly to a mixed state due to its coupling with the environment. We conclude that decoherence proceeds very efficiently for this model supporting previous findings [hollowood2017decoherence] derived in a different manner.
![Refer to caption](x10.png)
However, this raises an important question. In [Colas:2022kfu, Colas:2024ysu], it was shown that the purity decreases from 1 and then bounces back up to some value, when considering the (Gaussian) coupling of an adiabatic mode with an entropic one. This was rightly pointed out to be an effect of the non-Markovian nature of the system-environment coupling for an accelerating background. In the previous sub-sections, as was also done in [hollowood2017decoherence] using different approximations, we have demonstrated that the corrections to comes from the non-local part of the memory kernel, and thus it is established that this model is also non-Markovian. Why do we then not see such recoherence or purity-freezing in this case?
This can be understood if we focus on different terms appearing in (48). The first term on the RHS of that equation is characterized by diffusion, and is the dominant contribution to this term due to the initial conditions mentioned above. (Recall that is zero for our model and , which does have a non-local origin, multiplies the off-diagonal terms in .) It can be easily seen that the dynamics is primarily driven by this first term on the RHS which dominates over the second piece featuring the dissipative term (at least, until very late-times once the system classicalizes). Furthermore, the dissipative term forms the homogeneous part of the equation and since is negative for this model, the overall sign becomes positive for this term. Consequently, this yields an exponentially growing solution for the determinant, indicative of a decline in purity, i.e., decoherence. However, the main reason why we see no recoherence or purity-freezing is that the diffusion term , which dominates the evolution of the determinant of , is, in turn, dominated by the local part of the memory kernel as shown in Fig. 3 and in (33). Consequently, it also indicates that if we were to consider a model in which the coupling for a time-dependent background such that the diffusion terms are primarily derived from the non-local part of the kernel, then we might be able to see such exotic phenomena such as recoherence [Colas:2024xjy].
To drive home this message even more forcefully, let us solve (48) separately with only local and non-local components, where we note that while the non-local effects leave an oscillating flow of information between the system and environment, decoherence is purely driven by the local component of . We need to zoom in the plot (Fig. 9) for the evolution of the purity to see this. We plot sub-horizon evolution of the purity for this since, on artificially only keeping the non-local contributions to the transport equations, the plot for diverges after horizon-crossing. Nevertheless, when all the terms are kept (plot corresponding to ), the purity is perfectly well-behaved and shows efficient decoherence, and closely follows the plot for purity if we had (again, artificially) dropped all the non-local terms in the transport equations.
![Refer to caption](x11.png)
For comparison, we also evaluate the purity of the reduced density matrix using perturbation theory (PT). Since has been shown to contain all the terms appearing in and some other higher order terms [Colas:2023wxa], we get the equation by terminating the equation at order. Accordingly, the perturbative transport equation, at leading order, is given by
(49) |
The calculation of the perturbative purity shows that solution is almost identical to solution, which is a signature of validity for method. It also points to the real advantage of using the TCL formalism lies in its ability to resum secular divergences which appear for some physical observables, such as as shown in the previous sub-section when computed perturbatively, and not so much for purity in this model.
![Refer to caption](x12.png)
Finally, following the expressions of the spurious term , as computed in the Appendix, we also consider the solution for purity on including these terms. Since the spurious terms dominate at late-times, the purity naturally behaves radically differently in this case. In general, if we do not drop these terms, the system will go through incomplete decoherence and stabilize itself around some non-zero value of purity, as shown in Fig. 11, which is another physical motivation to drop such terms.
![Refer to caption](x13.png)
5 Conclusions
What have we learnt from revisiting this toy model? Firstly, we have shown that secular late-time divergences can appear from the non-local part of the memory kernel and show up in the quantum corrections to the power spectrum of the system field. This is made abundantly clear by the fact that even if we drop the local term in the kernel, the secular divergences in the dissipative terms persist in our model. This is in stark contrast to previous computations for quantum corrections to mode functions of scalar or photon fields due to graviton loops, where the local terms in the evolution equation always dominate at late-times and were responsible for secular divergences. More explicitly, for our model, the solution for the power spectrum is numerical. The leading order expansion of this quantity in the coupling constant is in general what one would have obtained in perturbation theory (in the absence of any resummations). Looking up Table 1, the leading order expansion can be written (for, say, ), as , which clearly denotes secularly divergent terms, the superscript ‘PT’ reminding us that this would be the result in perturbation theory999This, of course, can also be demonstrated from the results of [Boyanovsky:2015tba] (50) . Thus, quantum corrections to a physical observable can indeed have secularly growing contributions from non-Markovian terms and, what is more, the time-convolutionless form of the master equation is nevertheless able to dynamically resum such terms and render finite results. Although some of these observations have been made in the past, we correct some misconceptions regarding the Markovian approximation, namely that having a time-local master equation does not necessarily imply that the system must be Markovian.
This gives us a segue into how we have also shed light on an important aspect of such resummations which have been ignored in the past – spurious terms which arise from the lower (time) limit of the kernel integration as remnants of the initial state. This can be very cleanly handled using the TCL formalism and we have shown why such terms must be dropped by hand in order to have a meaningful resummation of physical observables. It makes sense to drop such terms since they are cancelled, between themselves, in a perturbative treatment to any order. Furthermore, we have resolved the interesting paradox of why previous computations of the power spectrum [Boyanovsky:2015jen, Boyanovsky:2015tba] yielded accurate results even when such a treatment of spurious terms were lacking. This was due to a judicious choice of the renormalization scale coupled with an additional ad hoc approximation of dropping the decaying mode contributions in the final answer. Although this approximation is good one for computing the power spectrum at the conformal boundary, it falls short in tracking the evolution of the elements of the covariance matrix when going to earlier times. The TCL formalism does not demand such approximations and render them redundant by applying a consistent analysis of the spurious terms. Another reason to have more faith in the TCL2 resummation over the super-horizon one employed in [Brahma:2021mng, Boyanovsky:2015tba] is that for exactly solvable Gaussian models, the former has been shown to have a closer fit with exact results [Colas:2022hlq].
Secondly, what our model shows is that although the memory kernel has a non-local part, the diffusion term which controls decoherence is dominated by the local term. In this context, we clarify that unlike previous claims, there is no emergent Markovian behaviour after horizon-crossing. First, this is evidenced by the fact that the corrections to the power spectrum, coming from non-local part of the kernel, survive in the limit. Just because the system density matrix attains “positivity” at late-times [hollowood2017decoherence], which is a quantitative measure of classicality of the system, does not mean that system has become Markovian. In fact, as has been shown in [hollowood2017decoherence] itself, one of the eigenvalues of the dissipator matrix remains negative at all times, thus clearly satisfying the definition of a non-Markovian system. Rather, what happens is that is dominated throughout by the contributions coming to it from local terms, which get rather pronounced at after horizon-crossing and at late times. Although plays the crucial role for decoherence, its contribution to the quantum correction to the power spectrum is negligible and thus different physical observables can depend on the local and non-local pieces of the memory kernel. One of the main reasons for the belief that late-time secular growth in inflationary models typically come from local terms is stochastic inflation since the latter invokes a white noise for modelling the quantum diffusion due to the sub-horizon modes. Our work shows that there might be a case for considering the non-local corrections to the noise term for stochastic inflation as well. Indeed, in this model, it is the dissipation terms that appear solely from the non-local terms and hence it will be interesting to consider a model in the future where the diffusion (or the noise) terms depend exclusively on the non-Markovian piece of the memory kernel. We leave constructing such a realistic model for future work [US:non-local_D22]. Our results also find that purity is not a sufficiently useful tool to probe the Markovianity of the model. In [Colas:2022kfu, Colas:2024ysu], it was shown that a non-Markovian kernel can lead to recoherence. What we find is that this is, however, not a smoking gun signal of non-Markovianity.
To the best of our knowledge, this is the first application the TCL formalism to a cosmological model involving a non-linear interaction generalizing beyond strictly Gaussian models. Our results, although sometimes technical, would serve well as benchmarks against applying this particular master equation formalism to future realistic models for inflation.
Appendix: Spurious Terms
As mentioned in Sec-3, the lower limits of integrals (19)-(22) give rise to so-called spurious terms. These terms, initially studied in the context of the cosmological Caldeira-Leggett model [caldeira1981influence, banerjee2023thermalization], were found to dominate at late times for all coefficients, which is unexpected for time-local (or quasi-time-local) master equations. In the same context, these terms were shown to be absent for perturbative solutions to any given order, as explained earlier. If these spurious terms are considered, they could lead to a breakdown of the resummation of the power spectrum, contradicting the exact results available for such a model.
Clearly, it is crucial to check the behaviour of these terms within the model we are working with. This not only completes the TCL2 description, but also adds this is necessary when considering non-Gaussian models such as this one, for which exact results are not available. An expectation that was articulated in [Colas:2023wxa] was that for more realistic interactions, going beyond the Gaussian case, it might be true that such spurious terms (coming from the lower limit of the integrals) would get suppressed. We want to show that this is certainly not the case for our model, and this has consequences for the non-perturbative resummation as has been explained in Sec-4.
For the reader’s convenience, we show again the explicit form of the master equation coefficients (keeping our color code for the local and divergent terms):
(51) |
(52) |
(53) |
(54) |
where are the expressions coming from the lower limit of the integrals. As for the upper limits, Analytical expressions are available for these functions, namely
(55) | ||||
(56) | ||||
(57) | ||||
(58) |
where we have isolated the nonzero parts in the limit . Notice how even after doing this, there is a surviving contribution at late times, which we should now check to see if they are dominant over the others. Expanding these functions in the limit , we get:
(59) |
Compare this to the non-spurious late-time behaviour of the master equation coefficients, given by
(60) | ||||
(61) | ||||
(62) | ||||
(63) |
Notice how the spurious terms are at least comparable to these expressions, with only the spurious part of being subdominant. Fig. 12 presents the solution of the transport equation for , comparing cases with and without the inclusion of spurious terms. It is immediately noticeable that when spurious terms are excluded, the solutions increase at a lower rate, a signal of resummation. Further, for the time range considered, the solutions obtained by including spurious terms look insensitive to the choice of , as illustrated in Fig. 13.
![Refer to caption](x14.png)
![Refer to caption](x15.png)
Acknowledgements
SB thanks Drazen Glavan for interesting discussions regarding secular divergences from non-local terms that spurred this work. The authors also thank Thomas Colas for numerous discussions about open quantum systems in cosmology.
SB is supported in part by the Higgs Fellowship and by the UK Science and Technology Facilities Council (STFC) Consolidated Grant “Particle Physics at the Higgs Centre”. JCF is supported by the STFC under grant ST/X001040/1. XL is supported in part by the Program of China Scholarship Council (Grant No. 202208170014).
For the purpose of open access, the author has applied a Creative Commons Attribution (CC BY) licence to any Author Accepted Manuscript version arising from this submission.