# Quantum Brownian Representation for the Quantum Field Modes

Authors: Daniel Arteaga
Journal: Advances in High Energy Physics
 (2009)

---

## Abstract

When analyzing the particle-like excitations in quantum field theory it is natural to regard
the field mode corresponding to the particle momentum as an open quantum system, together
with the opposite momentum mode. Provided that the state of the field is stationary,
homogeneous, and isotropic, this scalar two-mode system can be equivalently represented in
terms of a pair of quantum Brownian oscillators under a Gaussian approximation. In other
words, the two-mode system behaves as if it were interacting linearly with some effective
environment. In this paper we build the details of the effective linear coupling and the
effective environment, and argue that this quantum Brownian representation provides a
simple, universal, and nonperturbative characterization of any single particle-like excitation.
As immediate applications of the equivalence, we reanalyze the interpretation of the self-energy
in terms of decay rates in a general background state and present the master equation
for the field mode corresponding to the particle momentum.


---

## Body

## 1. Introduction and Motivation

Quantum fields can be regarded from the viewpoint of open quantum systems [1–4]. The degrees of freedom of the field which are relevant for the physical problem in question constitute the reduced subsystem and the rest form the environment. If there are several fields in interaction, and the object of interest is a particular field, it is natural to trace over the other environment fields. For instance, in electrodynamics one can study the so-called Euler-Heisenberg effective action for the photons [5], considering the electrons as the environment, or the complementary case, in which the electrons are taken as the system of interest and the photons are integrated out [6–9]. Similarly, in stochastic gravity [10–21] the system of interest is the gravitational field and the matter fields are integrated out. In many other circumstances it is natural to consider as the reduced system the modes of the quantum field which are below some ultraviolet cutoff, with the ultraviolet modes constituting the environment. This approach has been used, for instance, in studying bubble nucleation [22, 23], analyzing decoherence in field theory [24–26] and in inflationary cosmology [27–30]. As these references illustrate, the open quantum system point of view has often provided new tools and insights to different field theory problems, specially when dealing with states different than the Minkowski vacuum.The propagation of a particle-like excitation over a given field background is another situation in which there is a natural system-environment separation: the field mode corresponding to the particle momentum is the object of interest, and the rest of modes of the field, along with any other field in interaction, form the environment. In the Minkowski vacuum, the Källén-Lehmann representation of the propagator [5, 31–33] provides a complete description of single particle excitations from the field theory perspective, but there was no equivalent analysis in thermal backgrounds [34–38], or in curved spacetimes [39–42]. The open quantum system approach proves a useful tool when reanalysing the propagation of particle-like excitations in nonvacuum or nonMinkowski backgrounds [43, 44].Therefore, in this paper we elaborate on the open quantum description of the field mode corresponding to the particle momentum, with the goal of developing useful techniques for the analysis of particle and quasiparticle excitations. We concentrate on the basic methods and tools, and give example applications; the general analysis of the particle-like excitations in nonvacuum or nonflat backgrounds is left for separate publications [43, 44].We would like to consider the field mode as the system and the remaining modes as the environment. However, it proves difficult to implement this system-environment division directly, since the mode-decomposed fieldϕp is not real, but is a complex quantity obeying the contraint ϕp=ϕ-p* and acting on both the Hilbert space sector with momentum p and the sector with momentum  -p. Instead of focusing on a single mode, given that the field naturally links modes with opposite momentum, we will choose as the system of interest any two modes with a given opposite momentum, and as the environment the remaining modes of the field, as well as the modes of any other field in interaction. Namely, given a single scalar field ϕ, the system degrees of freedom are the two modes ϕp and ϕ-p, and the other modes ϕq, with q≠±p, form the environment. A similar system-environment division has been used in [45–47] in an inflationary context.The main goal of this paper is to show that this two-mode system behaves, under certain assumptions, as if it was an open quantum system interacting linearly with some effective environment, even if the original field interaction is nonlinear. The paradigm of linear open quantum system is the quantum Brownian motion (QBM) model, whose system of interest is a nonrelativistic particle interacting linearly with an infinite bath of harmonic oscillators. This model has had many applications in different contexts, among which one may mention the quantum to classical transition [48, 49], the escape from a potential well [50–53], the Unruh effect [54, 55] or quantum optics [3, 56]. In an influential paper Caldeira and Leggett [57] applied the influence functional model of Feynman and Vernon [58, 59] to the QBM model. The QBM model can be generalised to encompass the general class of linearly interacting open quantum systems [60, 61].The results of this paper can bee seen as an explicit implementation of the statement by Hu and Matacz [55] that the motion of a Brownian particle can be used to depict the behaviour of a single quantum field mode.We will allow for stationary and isotropic, but otherwise arbitrary, states for the field (a Gaussian approximation will be also assumed later on). Field theory with arbitrary field states can be studied within the closed time path (CTP), orin-in, method, which was originally proposed by Schwinger [62] and Keldysh [63]. The most characteristic feature of the CTP method, in contrast to the conventional  in-out  method, is the doubling of the number of degrees of freedom.The paper is organised as follows. In Section2, the relevant system-environment separation is discussed. In Section 3, the central section of this paper, we present and analyse the QBM representation for the field modes, given the structure of the generating functional and the two-point propagators. In Section 4, we illustrate the utility of the QBM representation with two-example applications. First, we rederive the interpretation of the imaginary part of the self-energy. Second, we build and analyze the master equation for the field modes corresponding to the particle momentum. The main body of the paper ends with Section 5, where we summarise and discuss the main results. In order not to break the continuity of the exposition, some details of the analysis are left for the appendices, which also provide background reference material which fixes the notation and makes the paper relatively self-contained. In Appendix A, we briefly present the CTP approach to field theory and apply it to the analysis of the structure of the two-point propagators, and in Appendix B, we introduce the theory of linear open quantum systems, focusing also in the structure of the two-point propagators. The analysis of the propagators done in the appendices is important for the discussion in Section 3.Throughout the paper we work with a system of natural units withℏ=c=1, denote quantum mechanical operators with a hat, and use a volume-dependent normalisation in the definition of the field modes (see (2.1) below). The same symbol will be used for a quantity and its Fourier transform whenever there is no danger of confusion.


## 2. Field Modes Regarded as Open Quantum Systems

Let us now present the relevant system-environment separation. For concreteness, we consider a self-interacting field theory model consisting of a single scalar fieldϕ, although results can be straightforwardly extended to any number of fields. The field ϕ can be decomposed in modes according to
(2.1)ϕp=1V∫d3xe-ip·xϕ(x),
where V is the volume of the space, a formally infinite quantity. The factor V-1/2 in the definition of ϕp is chosen so that the propagators verify
(2.2)G+(t,t′;p)=∫d3xe-ip·x〈0|ϕ̂(t,x)ϕ̂(t′,0)|0〉=〈0|ϕ̂p(t)ϕ̂-p(t′)|0〉.As stated in the introduction, given a particular momentump≠0, the system is composed by the two modes ϕp and ϕ-p, and the environment is composed by the other modes of the field, ϕq, with q≠±p. Should there be other fields in interaction of any arbitrary spin, the modes of these additional fields would also form part of the environment.The Hilbert space can be decomposed asℋ=ℋsys⊗ℋenv, where in turn ℋsys=ℋk⊗ℋ-k. Notice that this separation does not correspond to the Fock space decomposition. The entire system is in a state ρ̂; the state of the reduced system is ρ̂s=Trenvρ̂, and the state of the environment is  ρ̂e=Trsysρ̂. Generally speaking, the state for the entire system is not a factorised product state (i.e., ρ̂≠ρ̂s⊗ρ̂e).The action can be decomposed asS=Ssys+Scount+Senv+Sint, where Ssys is the renormalized system action,
(2.3a)Ssys=∫dt(ϕ̇pϕ̇-p-Ep2ϕpϕ-p),Scount is the appropriate counterterm action,
(2.3b)Scount=∫dt{(𝒵p-1)ϕ̇pϕ̇-p-[𝒵p(p2+m02)-Ep2]ϕpϕ-p},
(with m0 being the bare mass of the field), and Senv and Sint are the environment and interaction actions, respectively, which depend on the particular field theory model. Notice that we have allowed for an arbitrary rescaling of the field ϕ→ϕ/𝒵p1/2 and for an arbitrary frequency of the two-mode system Ep (which needs not be necessarily of the form (p2+m2)1/2).Let us draw our attention on the field rescaling and the frequency renormalisation. Since it is always possible to freely move finite terms from the system to the counterterm action andvice versa, both the field rescaling and the frequency renormalisation should be taken into account even if no infinities appeared in the perturbative calculations. A physical criterion needs to be chosen in order to fix the values of these two parameters. In the vacuum, such criterion is provided by the on-shell renormalisation scheme; in nonvacuum situations, it is investigated in [43], and will be briefly discussed in Section 4. Notice that the form of 𝒵p and Ep is not necessarily dictated by the Lorentz symmetry: even if the countertems which remove the infinities from the vacuum theory also remove the infinities in general field states, there can be finite Lorentz-breaking contributions. Anyway for the purposes of most of this paper the values of 𝒵p and Ep are not relevant and will be left unspecified (some comments will be made in Section 4.2 though).The system variablesϕp and ϕ-p are complex quantities verifying ϕ-p*=ϕp. We can construct real degrees of freedom by introducing the following change of variables:
(2.4)ϕΣ=12(ϕp+ϕ-p),ϕΔ=-i2(ϕp-ϕ-p).
In terms of the real variables, the system action can be re-expressed as
(2.5)Ssys=12∫dt(ϕ̇Σ2-Ep2ϕΣ2+ϕ̇Δ2-Ep2ϕΔ).
We could alternatively have directly obtained these two real degrees of freedom by working with the sine and cosine Fourier transform [55]. However we prefer to work with the exponential Fourier transform to make manifest the momentum conservation properties.Most information on the reduced quantum system can be extracted from the set of correlation functions, or equivalently from the CTP generating functional for the reduced system, which can be written as(2.6)Z[ja,α]=exp[-12!∫dtdt′jαa(t)jβb(t)Gabαβ(t,t′)-14!∫dtdt′dt′′dt′′′jαa(t)jβb(t′)jγc(t′′)jδd(t′′′)Gabcd(C)αβγδ(t,t′,t′′,t′′′)+⋯].
This expression is somewhat cumbersome and needs some clarification. Latin indices a,b,c,… are CTP indices and take the values 1 and 2, indicating, respectively, the forward and backwards time branches characteristic of the CTP formalism (see Appendix A). Greek indices α,β,γ,… take the two values +p and -p and make reference to the two field modes ϕp and ϕ-p. An Einstein summation convention is used both for Latin and Greek indices. The propagator Gabαβ(t,t′) is the 2-point propagator connecting CTP indices a and b, whose external legs correspond to particles with momenta α and β. When the state is translation-invariant and isotropic, momentum conservation imposes:
(2.7)Gab(+p)(-p)(t,t′)=Gab(-p)(+p)(t,t′)=Gab(t,t′;p),Gab(+p)(+p)(t,t′)=Gab(-p)(-p)(t,t′)=0.
In turn, Gabcd(C)αβγδ(t,t′,t′′,t′′′) is the connected part of the four-point correlation function having external legs with momenta α, β, γ, and δ. For translation-invariant states momentum conservation implies that only when momentum is balanced (i.e., two incoming and two outgoing external legs) the correlation function is nonvanishing. Terms with a higher number of external legs behave similarly.The open quantum system is nonlinear, and a systematic treatment of the generating functional can be done by using the tools of nonlinear open quantum systems. We follow a different path in the next section.


## 3. Quantum Brownian Motion Analogy

Let us reconsider the generating functional (2.6). It depends on the n-point correlation functions, with n being arbitrarily large. However, in many situations one is interested in properties which only depend on the two-point correlation functions. In other cases one is doing a perturbative expansion of the generating functional, and connected higher-order correlation functions are usually also of higher-order in the expansion parameter. Finally there are situations in which one only has access to the two-point correlation functions and expects (or simply hopes) that connected higher-order correlation functions are subdominant. In any of this situations one can be tempted to approximate the generating functional by the following Gaussian expression:
(3.1)Z[ja,α]≈exp[-12!∫dtdt′jαa(t)jβb(t)Gabαβ(t,t′)].
This Gaussian approximation can be controlled in the framework or the large-N expansion (where N is the number of scalar fields) [64–66]. We should emphasise that it does not necessarily imply any perturbative expansion in the coupling parameter, nor any free field approximation.The Gaussian approximation of the generating functional can be expanded to(3.2)Z[ja,p,ja,-p]=exp{-12∫dtdt′[jpa(t)j-pb(t′)Gab(t,t′)+jpa(t′)j-pb(t)Gab(t,t′)]},
or equivalently, in terms of the real variables ϕΣ and ϕΔ and their corresponding classical sources,
(3.3)Z[ja,Σ,ja,Δ]=exp{-12∫dtdt′[jΣa(t)jΣb(t′)Gab(t,t′)+jΔa(t)jΔb(t′)Gab(t,t′)]},
or
(3.4)Z[ja,Σ,ja,Δ]=Z[ja,Σ]Z[ja,Δ],Z[ja]=exp[-12∫dtdt′ja(t)jb(t′)Gab(t,t′)].
Equation (3.4) shows that for translation-invariant states within the Gaussian approximation the reduced two-mode state effectively behaves as two decoupled quantum mechanical degrees of freedom.The Gaussian approximation implies that the expression of the generating functional of a reduced two-mode in terms of the two-point propagators coincides with that of a QBM model, provided that the system is isotropic and translation invariant: compare (B.27) with (3.4). Moreover, by comparing (B.26) with (A.20), we also realise that the structure of the two-point propagators is identical in both cases. Notice that this latter fact is independent of the Gaussian approximation.Therefore, we conclude that, assuming homogeneity, isotropy, and stationarity,there is an equivalent QBM system for every scalar two-mode pair treated under the Gaussian approximation. In other words, within the Gaussian approximation, every two-mode of a given quantum field theory can be described in terms of a pair of quantum Brownian particles interacting linearly with some effective environment. We must stress that, similarly as the linear interaction does not coincide with the real coupling, the effective environment does not coincide with the real environment. The precise details of the equivalence are summarised in Table 1 for the particular case of the λϕ4 theory. We will see below that both the strength of the linear coupling and the state of the equivalent environment depend on the details of the original environment.Table 1
Detail of the equivalent linear QBM system for aλϕ4 quantum field theory. The λϕ4 model has been chosen for concreteness, but the correspondence would be analogous for any other field theory model. The symbol (=) indicates that the original and equivalent quantities are indeed identical despite the name change.
Original systemEquivalent linear QBMSystemTwo field modesTwo identical oscillatorsSystem d.o.f.ϕp, ϕ-p2 copies ofqEnvironmentOther modes 3-d field2 1-d fieldEnv. d.o.f.ϕq, q≠±pφpFrequencyEpΩ (=)Coupling(λ/3!)∑qq'ϕ±pϕqϕq'ϕ∓p-q-q'g2∑pℐ(p)q̇φp2-point functionGab(t,t′;p)Gab(t,t′)(=)Let us investigate on this correspondence. On the one hand, as  is shown in Appendix  B, the effect of the environment in the QBM system is fully encoded in two kernels, the dissipation kernelD(t,t′)—or its closely related counterpart H(t,t′); see (B.13b)—and the noise kernel N(t,t′). The dynamics of the quantum Brownian particle can be determined once the frequency and the noise and dissipation kernels are known. On the other hand, as shown in Appendix A, in a quantum field theory the two-point correlation functions are fully characterised by the frequency Ep and the retarded and Hadamard self-energies, ΣR(t,t′) and Σ(1)(t,t′), respectively. By comparing again (B.25) and (B.26) with (A.20) we realise that the precise analogy goes as follows:
(3.5)Ω=Ep,H(t,t′)=p2+m2-Ep+ΣR(t,t′;p),D(t,t′)=-iImΣR(t,t′;p),N(t,t′)=i2Σ(1)(t,t′;p),
where the quantities on the right-hand side correspond to the original field theory system and the quantities on the left-hand side correspond to the equivalent QBM model. The correspondence is valid no matter the renormalisation scheme chosen to fix 𝒵p and Rp.This representation provides thus a first rough interpretation for the retarded and Hadamard self-energies,ΣR(t,t′;p) and Σ(1)(t,t′;p), respectively. The retarded self-energy corresponds to the dissipation kernel, so that it determines the dissipative properties of the system, and it is independent of the state of the equivalent environment (though not independent of the state of the original environment). The Hadamard self-energy corresponds to the noise kernel, and thus it is basically related to fluctuations. Although the description in terms of the noise and dissipation kernels is often sufficient, the equivalent QBM system can be alternatively described in terms of the linear coupling constant to the effective environment, g, the distribution of frequencies of the environment ℐ(ω) (see Appendix B), and the occupation number of the modes of the effective environment n(p)=Tr[ρ̂e(eff)âp†âp]. The product g2ℐ(ω), which determines the coupling strength to the the ω-mode of the effective environment, can be obtained from the imaginary part of the self-energy:
(3.6a)ImΣR(ω;p)=iD(ω)=-g22ωℐ(ω).
This last equation implies that the equivalent coupling depends on the state of the real environment, since the retarded self-energy is state dependent in general. The occupation numbers of the effective environment n(p) can be reproduced from the Hadamard self-energy:
(3.6b)Σ(1)(ω;p)=-2iN(ω)=-2ig2|ω|ℐ(ω)[12+n(|ω|)].
The knowledge of the occupation numbers fully determines a Gaussian stationary state for the equivalent environment. These results follow from (B.14) and (B.16).Note that the analogy depends only on the Gaussian approximation, so that it can be extended to all orders in perturbation theory. Notice also that the results derived from the QBM interpretation are exact for all those properties which depend only on the two-point correlation functions; for the properties which depend on higher-order correlation functions, it is a correct approximation depending on the validity of the Gaussian approximation, that is, depending on the relative importance of the connected parts of the correlation functions with respect to the disconnected parts.


## 4. Example Applications

### 4.1. Interpretation of the Self-Energy

As a first example application, let us make use of the QBM correspondence to analyse the physical significance of the self-energy in general backgrounds. Our findings will coincide with the result by Weldon [67] (which nowadays is a textbook result [37]) as far as the imaginary part of the self-energy is concerned. However, while Weldon's original analysis was only valid to first-order in perturbation theory, our technique is be valid to all orders. Moreover, our analysis will not be tied to any field theory model. Additionally we will also obtain an interpretation for the other components of the self-energy. In this case the QBM analogy is exact since no four-point correlation functions are involved.To start we consider the probability that an excitation of energyω decays into the one-dimensional environment, in the equivalent QBM system. The probability Γ- that an excitation of the Browian particle with positive energy ω decays into the environment is given by (see, e.g., [37])
(4.1)Γ-(ω)=12ω∫dk2π2|k|2πδ(ω-|k|)|ℳ|2[1+n(|k|)],
where ℳ is the amplitude of the transition, and n(|k|) is the occupation number of the environment states with energy |k|. The factor 1+n(|k|), which is due to the Bose-Einstein statistics, enhances the decay probability to those states which are already occupied. Since the equivalent QBM system is linear, the squared decay amplitude is simply given to first-order in the linear coupling constant g by
(4.2)|ℳ|2=g2ℐ(ω)ω2,
where ℐ(ω) is the distribution of frequencies of the effective environment. The factor ω2 is a consequence of the derivative coupling in the QBM model. The decay probability is therefore
(4.3)Γ-(ω)=12g2ℐ(ω)[1+n(ω)].
Likewise, the probability that an excitation of positive energy ω is created spontaneously from the environment is given by
(4.4)Γ+(ω)=12ω∫dk2π2|k|2πδ(ω-|k|)|ℳ|2n(|k|)=12g2ℐ(ω)n(ω).
In the original system, Γ- can be interpreted as the probability that a (possibly off-shell) excitation with energy ω decays into the environment, and Γ+ can be interpreted as the probability that an environment spontaneously creates an excitation with energy ω.Notice that the notion of decay rate in quantum mechanics is meaningful only when the excitations are long-lived, or, to put it differently, when the productg2ℐ(ω) is very small. There is an inherent uncertainty in the concept of decay rate, which can be traced to the time-energy uncertainty principle. Therefore it is sufficient to present results to first-order in g: the inherent uncertainty to the notion of decay rate in quantum mechanics is of the same order as the error done by neglecting higher powers of g. In any case, this does not mean that we are doing any perturbative expansion in the original system: the decay rate can be computed to any desired order in the original perturbative coupling constant (e.g., λ in the case of the λϕ4 theory).We next analyse the self-energy components in the equivalent QBM mode. Given thatImΣR(ω)=(i/2)[Σ21(ω)-Σ12(ω)] (see Appendix A), we start by analyzing Σ21(ω). Applying CTP Feynman rules we get
(4.5)-iΣ21(t,t′)=-(ig)2∂t∂t′∫dp2πℐ(p)Trenv[ρ̂eφ̂Ip(t)φ̂I(-p)(t′)],
where φ̂Ip is p-mode of the environment field in the interaction picture. We have exploited the fact that first-order perturbation theory yields exact results for the self-energy in linear systems. Introducing two resolutions of the identity in the basis of eigenstates of the environment Hamiltonian is a simple exercise [68] to show that the above equation can be developed to
(4.6)Σ21(t,t′)=ig2∂t∂t′∫0∞dpℐ(p)p∑nρp,n[(n+1)e-ip(t-t′)+neip(t-t′)],
where ρp,n is the diagonal value of the reduced density matrix of the p-mode of the equivalent environment: Trq≠p〈np|ρ̂e|mp〉=ρp,nδnm. Introducing the Fourier transform we get
(4.7)Σ21(ω)=ig2ω2∫0∞dpℐ(p)p∑nρp,n[(n+1)δ(ω-p)+nδ(ω+p)],
and restricting to positive energies
(4.8)Σ21(ω)=-ig2ωℐ(ω)∑nρω,n(n+1)=ig2ωℐ(ω)[1+n(ω)],ω>0,
where we recall that n(ω)=Tr[ρ̂e(eff)âω†âω]=∑nρω,nn is the occupation number of the ω-mode of the environment. By comparing with (4.3), we thus see that i is proportional to the decay rate:
(4.9)Σ21(ω)=2iωΓ-(ω),ω>0.Repeating the calculation forΣ21(ω) we similarly find
(4.10)Σ12(ω)=ig2ωℐ(ω)∑nρω,nn=ig2ωℐ(ω)n(ω),ω>0,
having the corresponding interpretation in terms of the creation rate
(4.11)Σ12(ω)=2iωΓ+(ω),ω>0.
When the energies are negative Σ21(ω) and Σ12(ω) exchange roles.The imaginary part of the retarded self-energy is therefore given by(4.12)ImΣR(ω)=i2[Σ21(ω)-Σ12(ω)]=-12g2ωℐ(ω),
and can be interpreted as the net decay rate for an excitation of energy ω, that is, decay rate minus creation rate:
(4.13)ImΣR(ω)=-ω[Γ-(ω)-Γ+(ω)].
We therefore recover Weldon's result [67].We can additionally get an interpretation for the Hadamard self-energy. It is given by(4.14)Σ(1)(ω)=-Σ21(ω)-Σ12(ω)=-ig2|ω|ℐ(ω)[1+2n(|ω|)],
and it is proportional to the probability of decay plus the probability of creation,
(4.15)Σ(1)(ω)=-2i|ω|[Γ-(ω)+Γ+(ω)].


### 4.2. Master Equation for Relativistic Quasiparticles

As another example application of the correspondence, let us study the master equation for the second-quantized relativistic quasiparticles. A quasiparticle is an elementary excitation propagating in some backgroung, characterized by a momentump, an energy Ep, and a decay rate γp verifying γp≪Ep (so that the excitation is long lived). The basic properties characterizing the long-time evolution of the relativistic quasiparticles can be extracted from the analysis of the retarded propagator corresponding to the field mode in a model independent way [43] (differently to the short-time behaviour, which is model dependent). In particular, the real and imaginary parts of the self-energy, when evaluated on-shell, correspond to the physical energy and the decay rate of the quasiparticle, respectively,
(4.16)Ep2=m2+p2+ReΣR(Ep,p),γp=-1EpImΣR(Ep,p).
Using the QBM correspondence, these two equations can be equivalently written as:
(4.17)ReH(Ω)=0,γp=-iΩD(Ω).
The value of Ep (or equivalently Ω) is fixed by requiring that Ep represents the physical energy of the quantum mode as would be measured by a particle detector. We refer to [43] (which makes use of the QBM correspondence) for further details on the field theory analysis of the quasiparticle excitations.Thus, let us study the master equation (B.30) for q={ϕΣ,ϕΔ}. To this end it will prove useful to express the on-shell values of the noise and dissipation kernels in terms of the energy Ω=Ep and the decay rate γ=γp:
(4.18)D(Ω)=-H(Ω)=iγ,N(Ω)=γΩ(12+n),
where n=n(|Ω|) is the equilibrium occupation number of the mode in question. Notice that since the decay rate is necessarily small, the coupling constant to the effective environment is also small. Therefore the mode can be considered to be weakly coupled to the effective environment.In the weak coupling regime, the master equation coefficients are given by (B.31), and are usually divergent for short times [60]. These divergences are associated to the fact that the assumption of factorized initial conditions for the inital state is unphysical, as commented in Appendix B. Since we are interested in the long-time behaviour of the quasiparticles, we study the asymptotic value of these coefficients, which is free from the short-time divergences.It is easy to show that the frequency shift vanishes:(4.19)δΩ2=-2∫ti∞D(s,ti)cosΩ(s-ti)ds=12Re∫-∞∞H(s,ti)eiΩ(s-ti)ds=ReH(Ω)=0.
Taking into account that the dissipation kernel is antisymmetric, the evaluation of the dissipative factor is also straightforward:
(4.20)Γ=1Ω∫ti∞D(s,t)sinΩ(s-ti)ds=12Ωi∫-∞∞D(s,t)eiΩ(s-ti)ds=12ΩiD(ω)=γp2.
The first diffusion factor is also easily computed, this time recalling that the noise kernel is symmetric:
(4.21)Γh=∫ti∞N(s,ti)cosΩ(s-ti)ds=12∫-∞∞N(s,ti)eiΩ(s-ti)ds=N(Ω)2=Ωγ(12+n).
The second diffusion factor requires some extra work:
(4.22)Γf=1Ω∫ti∞N(s,ti)sinΩ(s-ti)ds=1ΩIm∫-∞∞N(s,ti)θ(s-ti)eiΩ(s-ti)ds=1ΩIm∫dω2πiN(ω)Ω-ω+iϵ=-1ΩP∫dω2πiD(ω)sign(ω)Ω-ω[1+2n(|ω|)].
The integrand in the last equality is only signifficantly different from zero when ω~Ω. Therefore as a first approximation we may write
(4.23)Γf≈-1ΩP∫dω2π-iD(ω)Ω-ω[1+2n(|Ω|)]=12ΩReH(Ω)(1+2n)=0.Thus the master equation can be written as(4.24)i∂∂tρs(q,q′,t)=[-12(∂2∂q2-∂2∂q′2)+12Ω2(q2-q′2)-i2γ(q-q′)(∂∂q-∂∂q′)-iΩγ(12+n)(q-q′)2]ρs(q,q′,t),
or,   equivalently, in terms of the Wigner function, as
(4.25)∂Ws∂t=-p∂Ws∂q+Ω2∂Ws∂p+γ∂pWs∂p+Ωγ4(1+2n)∂2Ws∂p2.
In this equation p is the canonical momentum associated to the variable q and has nothing to do with the physical quasiparticle momentum p. This is identical to the master equation found by Caldeira et al. [69], corresponding to weak couplings. Fleming et al.  [70] found logarithmically divergent results for the master equation coefficients in the ohmic dissipation case, even in the asymptotic regime and in the weak coupling limit. This divergences are associated to the unphysical ultraviolet behaviour of the ohmic dissipation model. Our approximate treatment overlooks the possible divergences associated to the high energy limit of the noise and dissipation kernel.The reduced density matrix and the Wigner function are to be interpreted in terms of the field modes, and not directly in terms of particles or quasiparticles themselves. Qualitatively, for Gaussian initial conditions the solution of the master equation can be described in terms of the cumulants of a Gaussian distribution in the following way [48, 70]. The expectation value of the field 〈q〉 follows the trajectory of a classical underdamped harmonic oscillator (q̈+γq̇+Ω2q=0), namely, it oscillates and slowly decays towards the origin at a rate γ/2. However, for single quasiparticle excitations the expectation value of the field is always vanishing. Therefore, it is more appealing to consider the dynamics of the second-order cumulants, and in particular the dynamics of the energy of the mode, (1/2)〈q̈+Ω2q〉. When perturbed by the introduction of a quasiparticle, the energy of the mode slowly decays at a rate γ towards its equilibrium value Ω(1/2+n), with n being the original occupation number of the mode [43].


## 4.1. Interpretation of the Self-Energy

As a first example application, let us make use of the QBM correspondence to analyse the physical significance of the self-energy in general backgrounds. Our findings will coincide with the result by Weldon [67] (which nowadays is a textbook result [37]) as far as the imaginary part of the self-energy is concerned. However, while Weldon's original analysis was only valid to first-order in perturbation theory, our technique is be valid to all orders. Moreover, our analysis will not be tied to any field theory model. Additionally we will also obtain an interpretation for the other components of the self-energy. In this case the QBM analogy is exact since no four-point correlation functions are involved.To start we consider the probability that an excitation of energyω decays into the one-dimensional environment, in the equivalent QBM system. The probability Γ- that an excitation of the Browian particle with positive energy ω decays into the environment is given by (see, e.g., [37])
(4.1)Γ-(ω)=12ω∫dk2π2|k|2πδ(ω-|k|)|ℳ|2[1+n(|k|)],
where ℳ is the amplitude of the transition, and n(|k|) is the occupation number of the environment states with energy |k|. The factor 1+n(|k|), which is due to the Bose-Einstein statistics, enhances the decay probability to those states which are already occupied. Since the equivalent QBM system is linear, the squared decay amplitude is simply given to first-order in the linear coupling constant g by
(4.2)|ℳ|2=g2ℐ(ω)ω2,
where ℐ(ω) is the distribution of frequencies of the effective environment. The factor ω2 is a consequence of the derivative coupling in the QBM model. The decay probability is therefore
(4.3)Γ-(ω)=12g2ℐ(ω)[1+n(ω)].
Likewise, the probability that an excitation of positive energy ω is created spontaneously from the environment is given by
(4.4)Γ+(ω)=12ω∫dk2π2|k|2πδ(ω-|k|)|ℳ|2n(|k|)=12g2ℐ(ω)n(ω).
In the original system, Γ- can be interpreted as the probability that a (possibly off-shell) excitation with energy ω decays into the environment, and Γ+ can be interpreted as the probability that an environment spontaneously creates an excitation with energy ω.Notice that the notion of decay rate in quantum mechanics is meaningful only when the excitations are long-lived, or, to put it differently, when the productg2ℐ(ω) is very small. There is an inherent uncertainty in the concept of decay rate, which can be traced to the time-energy uncertainty principle. Therefore it is sufficient to present results to first-order in g: the inherent uncertainty to the notion of decay rate in quantum mechanics is of the same order as the error done by neglecting higher powers of g. In any case, this does not mean that we are doing any perturbative expansion in the original system: the decay rate can be computed to any desired order in the original perturbative coupling constant (e.g., λ in the case of the λϕ4 theory).We next analyse the self-energy components in the equivalent QBM mode. Given thatImΣR(ω)=(i/2)[Σ21(ω)-Σ12(ω)] (see Appendix A), we start by analyzing Σ21(ω). Applying CTP Feynman rules we get
(4.5)-iΣ21(t,t′)=-(ig)2∂t∂t′∫dp2πℐ(p)Trenv[ρ̂eφ̂Ip(t)φ̂I(-p)(t′)],
where φ̂Ip is p-mode of the environment field in the interaction picture. We have exploited the fact that first-order perturbation theory yields exact results for the self-energy in linear systems. Introducing two resolutions of the identity in the basis of eigenstates of the environment Hamiltonian is a simple exercise [68] to show that the above equation can be developed to
(4.6)Σ21(t,t′)=ig2∂t∂t′∫0∞dpℐ(p)p∑nρp,n[(n+1)e-ip(t-t′)+neip(t-t′)],
where ρp,n is the diagonal value of the reduced density matrix of the p-mode of the equivalent environment: Trq≠p〈np|ρ̂e|mp〉=ρp,nδnm. Introducing the Fourier transform we get
(4.7)Σ21(ω)=ig2ω2∫0∞dpℐ(p)p∑nρp,n[(n+1)δ(ω-p)+nδ(ω+p)],
and restricting to positive energies
(4.8)Σ21(ω)=-ig2ωℐ(ω)∑nρω,n(n+1)=ig2ωℐ(ω)[1+n(ω)],ω>0,
where we recall that n(ω)=Tr[ρ̂e(eff)âω†âω]=∑nρω,nn is the occupation number of the ω-mode of the environment. By comparing with (4.3), we thus see that i is proportional to the decay rate:
(4.9)Σ21(ω)=2iωΓ-(ω),ω>0.Repeating the calculation forΣ21(ω) we similarly find
(4.10)Σ12(ω)=ig2ωℐ(ω)∑nρω,nn=ig2ωℐ(ω)n(ω),ω>0,
having the corresponding interpretation in terms of the creation rate
(4.11)Σ12(ω)=2iωΓ+(ω),ω>0.
When the energies are negative Σ21(ω) and Σ12(ω) exchange roles.The imaginary part of the retarded self-energy is therefore given by(4.12)ImΣR(ω)=i2[Σ21(ω)-Σ12(ω)]=-12g2ωℐ(ω),
and can be interpreted as the net decay rate for an excitation of energy ω, that is, decay rate minus creation rate:
(4.13)ImΣR(ω)=-ω[Γ-(ω)-Γ+(ω)].
We therefore recover Weldon's result [67].We can additionally get an interpretation for the Hadamard self-energy. It is given by(4.14)Σ(1)(ω)=-Σ21(ω)-Σ12(ω)=-ig2|ω|ℐ(ω)[1+2n(|ω|)],
and it is proportional to the probability of decay plus the probability of creation,
(4.15)Σ(1)(ω)=-2i|ω|[Γ-(ω)+Γ+(ω)].


## 4.2. Master Equation for Relativistic Quasiparticles

As another example application of the correspondence, let us study the master equation for the second-quantized relativistic quasiparticles. A quasiparticle is an elementary excitation propagating in some backgroung, characterized by a momentump, an energy Ep, and a decay rate γp verifying γp≪Ep (so that the excitation is long lived). The basic properties characterizing the long-time evolution of the relativistic quasiparticles can be extracted from the analysis of the retarded propagator corresponding to the field mode in a model independent way [43] (differently to the short-time behaviour, which is model dependent). In particular, the real and imaginary parts of the self-energy, when evaluated on-shell, correspond to the physical energy and the decay rate of the quasiparticle, respectively,
(4.16)Ep2=m2+p2+ReΣR(Ep,p),γp=-1EpImΣR(Ep,p).
Using the QBM correspondence, these two equations can be equivalently written as:
(4.17)ReH(Ω)=0,γp=-iΩD(Ω).
The value of Ep (or equivalently Ω) is fixed by requiring that Ep represents the physical energy of the quantum mode as would be measured by a particle detector. We refer to [43] (which makes use of the QBM correspondence) for further details on the field theory analysis of the quasiparticle excitations.Thus, let us study the master equation (B.30) for q={ϕΣ,ϕΔ}. To this end it will prove useful to express the on-shell values of the noise and dissipation kernels in terms of the energy Ω=Ep and the decay rate γ=γp:
(4.18)D(Ω)=-H(Ω)=iγ,N(Ω)=γΩ(12+n),
where n=n(|Ω|) is the equilibrium occupation number of the mode in question. Notice that since the decay rate is necessarily small, the coupling constant to the effective environment is also small. Therefore the mode can be considered to be weakly coupled to the effective environment.In the weak coupling regime, the master equation coefficients are given by (B.31), and are usually divergent for short times [60]. These divergences are associated to the fact that the assumption of factorized initial conditions for the inital state is unphysical, as commented in Appendix B. Since we are interested in the long-time behaviour of the quasiparticles, we study the asymptotic value of these coefficients, which is free from the short-time divergences.It is easy to show that the frequency shift vanishes:(4.19)δΩ2=-2∫ti∞D(s,ti)cosΩ(s-ti)ds=12Re∫-∞∞H(s,ti)eiΩ(s-ti)ds=ReH(Ω)=0.
Taking into account that the dissipation kernel is antisymmetric, the evaluation of the dissipative factor is also straightforward:
(4.20)Γ=1Ω∫ti∞D(s,t)sinΩ(s-ti)ds=12Ωi∫-∞∞D(s,t)eiΩ(s-ti)ds=12ΩiD(ω)=γp2.
The first diffusion factor is also easily computed, this time recalling that the noise kernel is symmetric:
(4.21)Γh=∫ti∞N(s,ti)cosΩ(s-ti)ds=12∫-∞∞N(s,ti)eiΩ(s-ti)ds=N(Ω)2=Ωγ(12+n).
The second diffusion factor requires some extra work:
(4.22)Γf=1Ω∫ti∞N(s,ti)sinΩ(s-ti)ds=1ΩIm∫-∞∞N(s,ti)θ(s-ti)eiΩ(s-ti)ds=1ΩIm∫dω2πiN(ω)Ω-ω+iϵ=-1ΩP∫dω2πiD(ω)sign(ω)Ω-ω[1+2n(|ω|)].
The integrand in the last equality is only signifficantly different from zero when ω~Ω. Therefore as a first approximation we may write
(4.23)Γf≈-1ΩP∫dω2π-iD(ω)Ω-ω[1+2n(|Ω|)]=12ΩReH(Ω)(1+2n)=0.Thus the master equation can be written as(4.24)i∂∂tρs(q,q′,t)=[-12(∂2∂q2-∂2∂q′2)+12Ω2(q2-q′2)-i2γ(q-q′)(∂∂q-∂∂q′)-iΩγ(12+n)(q-q′)2]ρs(q,q′,t),
or,   equivalently, in terms of the Wigner function, as
(4.25)∂Ws∂t=-p∂Ws∂q+Ω2∂Ws∂p+γ∂pWs∂p+Ωγ4(1+2n)∂2Ws∂p2.
In this equation p is the canonical momentum associated to the variable q and has nothing to do with the physical quasiparticle momentum p. This is identical to the master equation found by Caldeira et al. [69], corresponding to weak couplings. Fleming et al.  [70] found logarithmically divergent results for the master equation coefficients in the ohmic dissipation case, even in the asymptotic regime and in the weak coupling limit. This divergences are associated to the unphysical ultraviolet behaviour of the ohmic dissipation model. Our approximate treatment overlooks the possible divergences associated to the high energy limit of the noise and dissipation kernel.The reduced density matrix and the Wigner function are to be interpreted in terms of the field modes, and not directly in terms of particles or quasiparticles themselves. Qualitatively, for Gaussian initial conditions the solution of the master equation can be described in terms of the cumulants of a Gaussian distribution in the following way [48, 70]. The expectation value of the field 〈q〉 follows the trajectory of a classical underdamped harmonic oscillator (q̈+γq̇+Ω2q=0), namely, it oscillates and slowly decays towards the origin at a rate γ/2. However, for single quasiparticle excitations the expectation value of the field is always vanishing. Therefore, it is more appealing to consider the dynamics of the second-order cumulants, and in particular the dynamics of the energy of the mode, (1/2)〈q̈+Ω2q〉. When perturbed by the introduction of a quasiparticle, the energy of the mode slowly decays at a rate γ towards its equilibrium value Ω(1/2+n), with n being the original occupation number of the mode [43].


## 5. Summary and Discussion

In this paper we have explored the open quantum system viewpoint for a pair of field modes of opposite momentum, which are the relevant degrees of freedom for the analysis of particle-like excitations in field theory. The main results have been, on the one hand, showing that in any interacting field theory, assuming homogeneity, stationarity, and Gaussianity, this open quantum system can be equivalently represented by two identical quantum Brownian particles interacting linearly with an effective environment, and, on the other hand, exploring the details of the equivalence, which are expressed in (3.5), (3.6a), and (3.6b) and in Table 1.The Brownian motion equivalence is based on three simple well-known observations. First, the fact that the structure of the two-point correlation functions, which we re-examined in the appendices, is identical regardless of the nature of the interactions. A byproduct of this fact has been establishing a link between the self-energy and the noise and dissipation kernels, thus connecting the open quantum system and field theory notation and languages. Second, the recognition that a Gaussian truncation leads to a generating functional which coincides with that of a linear theory. Finally, the observation that when assuming homogeneity, stationarity, and isotropy the two-mode pair behaves as  two copies of a single degree of freedom.If the system is Gaussian, stationary, and homogeneous, the equivalence is exact. If the system is nonGaussian, the analogy is also exact as far as two-point correlation functions are concerned, but only approximate for higher-order correlation functions. If the system is nonstationary or nonhomogeneous, there will be corrections to the results of the order ofLp, where L is the characteristic inhomogeneity time or length scale, and p is the relevant energy or momentum scale. Therefore the analogy is perfectly valid in nonhomogeneous backgrounds as long as we consider modes whose characteristic wavelengths are much smaller than the inhomogeneity scale.It must be noted that for the equivalence to be useful one still needs the field theoretic computation of the two-point correlation functions (or, equivalently, the self-energy). The correspondence does not help in this calculation but must be regarded as a tool useful for interpreting and analyzing the dynamics of the two-mode system. In this sense, there are three basic characteristics which make the QBM equivalence appealing.First, the equivalence isuniversal, in the sense that it provides the most general description of the dynamics of the scalar two-mode pair within the Gaussian approximation. The description in terms of a linear open quantum system allows a unified description of many different quantum field theory systems: the details of the quantum field theory are unimportant once the equivalent noise and dissipation kernels are known. Different field theory models can be thus classified in the same equivalence class if they lead to the same QBM equivalent model.Second, the equivalence provides asimple  characterization of the interaction, given that an arbitrarily complicated coupling with any number of fields is reduced to a linear interaction with a one-dimensional field. Linearly interacting systems have been thoroughly studied in the literature (see, e.g., references given in the introduction and Appendix B) and exhibit many intersting properties, among which one can cite the fact that they are exactly solvable. Therefore, within the regime of validity of the Gaussian approximation, all the methods developed in the literature for linear quantum systems can be applied to analyse the dynamics of a two-mode pair.Finally, the equivalence isnonperturbative, since it does not rely on any perturbative development in the coupling constant of the original field theory. It can be therefore extended to all orders in perturbation theory and applied to strongly interacting systems.The QBM correspondence, which has also been implicitly used in [71, 72], proves a useful tool when analyzing the properties of particle-like excitations in general backgrounds from a field theory perspective [43, 44]. The utility of the QBM analogy is already highlighted by the first two immediate applications presented in this paper.As a first application, we reexplored the well-known result that the imaginary part of the retarded self-energy corresponds to the net decay rate of the particle excitations. This was a rather simple application which however improved the conventional textbook derivation by taking profit of the three properties we have remarked above: in the first place, the universality of the QBM analogy, given that the derivation presented in this paper was not tied to any specific field theory model; in the second place, its simplicity, since the calculation essentially reduced to doing trivial perturbative expansions in the equivalent linear system, and, finally, its nonperturbative character, since the derivation avoided any perturbative expansion in the original system.As a second application, we presented the relevant master equation for the dynamics of the modes corresponding to the quasiparticle momentum. This was a slighly more involved application, relying on the results of [43] (and also complementing them), which highligted the fact that the linear open quantum system machinery can readly be exported to quantum field theory. The presentation and analysis of the master equation, as done in the paper, must be understood as a first approximation to the problem, which might be sufficiently interesting by itself to deserve further work.



---

*Source: 278759-2009-08-02.xml*