Research Article  Open Access
Quantum Brownian Representation for the Quantum Field Modes
Abstract
When analyzing the particlelike 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 twomode system can be equivalently represented in terms of a pair of quantum Brownian oscillators under a Gaussian approximation. In other words, the twomode 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 particlelike excitation. As immediate applications of the equivalence, we reanalyze the interpretation of the selfenergy in terms of decay rates in a general background state and present the master equation for the field mode corresponding to the particle momentum.
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 socalled EulerHeisenberg 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 particlelike excitation over a given field background is another situation in which there is a natural systemenvironment 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énLehmann 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 particlelike 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 particlelike 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 systemenvironment division directly, since the modedecomposed field is not real, but is a complex quantity obeying the contraint and acting on both the Hilbert space sector with momentum and the sector with momentum . 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 and , and the other modes , with , form the environment. A similar systemenvironment division has been used in [45–47] in an inflationary context.
The main goal of this paper is to show that this twomode 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), or inin, method, which was originally proposed by Schwinger [62] and Keldysh [63]. The most characteristic feature of the CTP method, in contrast to the conventional inout method, is the doubling of the number of degrees of freedom.
The paper is organised as follows. In Section 2, the relevant systemenvironment 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 twopoint propagators. In Section 4, we illustrate the utility of the QBM representation with twoexample applications. First, we rederive the interpretation of the imaginary part of the selfenergy. 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 selfcontained. In Appendix A, we briefly present the CTP approach to field theory and apply it to the analysis of the structure of the twopoint propagators, and in Appendix B, we introduce the theory of linear open quantum systems, focusing also in the structure of the twopoint 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 , denote quantum mechanical operators with a hat, and use a volumedependent 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 systemenvironment separation. For concreteness, we consider a selfinteracting 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 where is the volume of the space, a formally infinite quantity. The factor in the definition of is chosen so that the propagators verify
As stated in the introduction, given a particular momentum , the system is composed by the two modes and , and the environment is composed by the other modes of the field, , with . 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 , where in turn . 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 , and the state of the environment is . Generally speaking, the state for the entire system is not a factorised product state (i.e., ).
The action can be decomposed as , where is the renormalized system action, is the appropriate counterterm action, (with being the bare mass of the field), and and 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 and for an arbitrary frequency of the twomode system (which needs not be necessarily of the form ).
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 and vice 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 onshell renormalisation scheme; in nonvacuum situations, it is investigated in [43], and will be briefly discussed in Section 4. Notice that the form of and 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 Lorentzbreaking contributions. Anyway for the purposes of most of this paper the values of and are not relevant and will be left unspecified (some comments will be made in Section 4.2 though).
The system variables and are complex quantities verifying . We can construct real degrees of freedom by introducing the following change of variables: In terms of the real variables, the system action can be reexpressed as 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 This expression is somewhat cumbersome and needs some clarification. Latin indices 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 and and make reference to the two field modes and . An Einstein summation convention is used both for Latin and Greek indices. The propagator is the 2point propagator connecting CTP indices and , whose external legs correspond to particles with momenta and . When the state is translationinvariant and isotropic, momentum conservation imposes: In turn, is the connected part of the fourpoint correlation function having external legs with momenta , , and . For translationinvariant 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 point correlation functions, with being arbitrarily large. However, in many situations one is interested in properties which only depend on the twopoint correlation functions. In other cases one is doing a perturbative expansion of the generating functional, and connected higherorder correlation functions are usually also of higherorder in the expansion parameter. Finally there are situations in which one only has access to the twopoint correlation functions and expects (or simply hopes) that connected higherorder correlation functions are subdominant. In any of this situations one can be tempted to approximate the generating functional by the following Gaussian expression: This Gaussian approximation can be controlled in the framework or the large expansion (where 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 or equivalently, in terms of the real variables and and their corresponding classical sources, or Equation (3.4) shows that for translationinvariant states within the Gaussian approximation the reduced twomode 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 twomode in terms of the twopoint 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 twopoint 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 twomode pair treated under the Gaussian approximation. In other words, within the Gaussian approximation, every twomode 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 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.

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 kernel —or its closely related counterpart ; see (B.13b)—and the noise kernel . 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 twopoint correlation functions are fully characterised by the frequency and the retarded and Hadamard selfenergies, and respectively. By comparing again (B.25) and (B.26) with (A.20) we realise that the precise analogy goes as follows: where the quantities on the righthand side correspond to the original field theory system and the quantities on the lefthand side correspond to the equivalent QBM model. The correspondence is valid no matter the renormalisation scheme chosen to fix and .
This representation provides thus a first rough interpretation for the retarded and Hadamard selfenergies, and respectively. The retarded selfenergy 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 selfenergy 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, , the distribution of frequencies of the environment (see Appendix B), and the occupation number of the modes of the effective environment . The product , which determines the coupling strength to the the mode of the effective environment, can be obtained from the imaginary part of the selfenergy: This last equation implies that the equivalent coupling depends on the state of the real environment, since the retarded selfenergy is state dependent in general. The occupation numbers of the effective environment can be reproduced from the Hadamard selfenergy: 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 twopoint correlation functions; for the properties which depend on higherorder 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 SelfEnergy
As a first example application, let us make use of the QBM correspondence to analyse the physical significance of the selfenergy 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 selfenergy is concerned. However, while Weldon's original analysis was only valid to firstorder 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 selfenergy. In this case the QBM analogy is exact since no fourpoint correlation functions are involved.
To start we consider the probability that an excitation of energy decays into the onedimensional 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]) where is the amplitude of the transition, and is the occupation number of the environment states with energy . The factor , which is due to the BoseEinstein 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 firstorder in the linear coupling constant by where is the distribution of frequencies of the effective environment. The factor is a consequence of the derivative coupling in the QBM model. The decay probability is therefore Likewise, the probability that an excitation of positive energy is created spontaneously from the environment is given by In the original system, can be interpreted as the probability that a (possibly offshell) 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 longlived, or, to put it differently, when the product is very small. There is an inherent uncertainty in the concept of decay rate, which can be traced to the timeenergy uncertainty principle. Therefore it is sufficient to present results to firstorder in : 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 . 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 theory).
We next analyse the selfenergy components in the equivalent QBM mode. Given that (see Appendix A), we start by analyzing . Applying CTP Feynman rules we get where is mode of the environment field in the interaction picture. We have exploited the fact that firstorder perturbation theory yields exact results for the selfenergy 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 where is the diagonal value of the reduced density matrix of the mode of the equivalent environment: . Introducing the Fourier transform we get and restricting to positive energies where we recall that is the occupation number of the mode of the environment. By comparing with (4.3), we thus see that is proportional to the decay rate:
Repeating the calculation for we similarly find having the corresponding interpretation in terms of the creation rate When the energies are negative and exchange roles.
The imaginary part of the retarded selfenergy is therefore given by and can be interpreted as the net decay rate for an excitation of energy , that is, decay rate minus creation rate: We therefore recover Weldon's result [67].
We can additionally get an interpretation for the Hadamard selfenergy. It is given by and it is proportional to the probability of decay plus the probability of creation,
4.2. Master Equation for Relativistic Quasiparticles
As another example application of the correspondence, let us study the master equation for the secondquantized relativistic quasiparticles. A quasiparticle is an elementary excitation propagating in some backgroung, characterized by a momentum , an energy , and a decay rate verifying (so that the excitation is long lived). The basic properties characterizing the longtime 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 shorttime behaviour, which is model dependent). In particular, the real and imaginary parts of the selfenergy, when evaluated onshell, correspond to the physical energy and the decay rate of the quasiparticle, respectively, Using the QBM correspondence, these two equations can be equivalently written as: The value of (or equivalently ) is fixed by requiring that 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 . To this end it will prove useful to express the onshell values of the noise and dissipation kernels in terms of the energy and the decay rate : where 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 longtime behaviour of the quasiparticles, we study the asymptotic value of these coefficients, which is free from the shorttime divergences.
It is easy to show that the frequency shift vanishes: Taking into account that the dissipation kernel is antisymmetric, the evaluation of the dissipative factor is also straightforward: The first diffusion factor is also easily computed, this time recalling that the noise kernel is symmetric: The second diffusion factor requires some extra work: The integrand in the last equality is only signifficantly different from zero when . Therefore as a first approximation we may write
Thus the master equation can be written as or, equivalently, in terms of the Wigner function, as In this equation is the canonical momentum associated to the variable and has nothing to do with the physical quasiparticle momentum . 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 follows the trajectory of a classical underdamped harmonic oscillator (), namely, it oscillates and slowly decays towards the origin at a rate . 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 secondorder cumulants, and in particular the dynamics of the energy of the mode, . When perturbed by the introduction of a quasiparticle, the energy of the mode slowly decays at a rate towards its equilibrium value , with 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 particlelike 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 wellknown observations. First, the fact that the structure of the twopoint correlation functions, which we reexamined in the appendices, is identical regardless of the nature of the interactions. A byproduct of this fact has been establishing a link between the selfenergy 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 twomode 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 twopoint correlation functions are concerned, but only approximate for higherorder correlation functions. If the system is nonstationary or nonhomogeneous, there will be corrections to the results of the order of , where is the characteristic inhomogeneity time or length scale, and 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 twopoint correlation functions (or, equivalently, the selfenergy). The correspondence does not help in this calculation but must be regarded as a tool useful for interpreting and analyzing the dynamics of the twomode system. In this sense, there are three basic characteristics which make the QBM equivalence appealing.
First, the equivalence is universal, in the sense that it provides the most general description of the dynamics of the scalar twomode 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 a simple characterization of the interaction, given that an arbitrarily complicated coupling with any number of fields is reduced to a linear interaction with a onedimensional 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 twomode pair.
Finally, the equivalence is nonperturbative, 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 particlelike 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 wellknown result that the imaginary part of the retarded selfenergy 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.
Appendices
A. The Closed Time Path Method and the TwoPoint Propagators in Field Theory
In this appendix we give a brief introduction to the closed time path (CTP) method (also called inin method, in contrast to the conventional inout method), stressing those aspects relevant for this paper, and apply it to the analysis of the twopoint propagators. We address the reader to [73–77] for further details on the CTP method, and to [68, 78–80] for further details on the structure of the twopoint functions. For the purposes of this appendix we will consider a free or an interacting scalar field , although results also apply for a single quantum mechanical degree of freedom.
The pathordered generating functional is defined as where is the field operator in the Heisenberg picture, is a certain path in the complex plane, means time ordering along this path, and is a classical external source. By functional differentiation of the generating functional with respect to , pathordered correlation functions can be obtained Introducing a complete basis of eigenstates of the field operator in the Heisenberg picture, , as a representation of the identity, the generating functional can be expressed as The functional measures and go over all field configurations of the fields at fixed initial time . If the path begins and ends at the same point , then the transition element of the evolution operator can be computed via a path integral: where is the Lagrangian density of the scalar field.
Let us consider the time path shown in Figure 1. If we define and for , then the generating functional can be reexpressed as In the following it will prove useful to use a condensed notation where neither the boundary conditions of the path integral nor the integrals over the initial and final states are explicit. With this simplified notation the above equation becomes An operator representation is also possible:
By functionally differentiating the generating functional, the Feynman and Dyson propagators and the Whightman functions can be obtained These four propagators can be conveniently organised in a matrix, the socalled direct matrix, Lowercase roman indices may acquire the values 1 and 2 are raised and lowered with the “CTP metric’’ . Higherorder correlation functions can be obtained in a similar way.
We may also consider the PauliJordan or commutator propagator, and the Hadamard or anticonmutator function For linear systems (systems whose Heisemberg equations of motion are linear) the PauliJordan propagator is independent of the state and carries information about the system dynamics. Finally, one can also consider the retarded and advanced propagators The kernels and are called, respectively, dissipation and noise kernels. The frequency shift is a formally divergent quantity given by
The dissipation and noise kernels can be computed following (B.10) and (B.13a)–(B.13e). The value of the dissipation kernel in the frequency space is [68] The dissipation kernel is closely related to the kernel (see (B.13a) and (B.13b), which is also stateindependent and given by In contrast, the noise kernel (B.13d) is state dependent. For a general Gaussian stationary environments, characterised by the occupation numbers , the noise kernel in Fourier space is given by The dissipation and noise kernels are related through For the particular case of an environment in thermal equilibrium at a temperature the occupation numbers are given by and the above equation becomes the fluctuationdissipation theorem:
By considering an arbitrary distribution of frequencies and an arbitrary Gaussian state for the environment the dissipation and noise kernels may adopt almost any value. In the rest of the appendix we will try to express all results in terms of the dissipation and noise kernels. To this end, it will prove useful to reexpress (B.13b) in Fourier space: The kernel can be decomposed in its real and imaginary parts as We have used the property . Notice that . The real and imaginary parts of the kernel in frequency space correspond, respectively, to the even and odd parts in time space. Notice also the KramersKronig relation between the real and imaginary parts of the kernel : The frequency shift can be always absorbed in . From now on, and in the main body of the paper, we will assume that such absorption has been carried out.
For a Gaussian environment and asymptotic initial boundary conditions the generating functional can be expressed as [82] where , and is the retarded propagator of the kernel that is, the kernel which verifies
In the general case there would be a prefactor taking into account the initial conditions for the system. However, since the system has a dissipative dynamics and the initial conditions are given in the remote past initial conditions for the system turn out to be completely irrelevant.
Explicit expressions are most easily obtained in the Fourier space, in which the retarded propagator adopts the form It can be checked that , besides being the retarded propagator of the kernel , in the sense of (B.24), is also the retarded propagator of the quantum mechanical system, in the sense of (A.12a).
Differentiating the CTP generating functional, according to (A.8), we obtain the different correlation functions in terms of the noise and dissipation kernels: The generating functional can be alternatively expressed in terms of these correlation functions as
If desired, the dynamics of the Brownian oscillator can be analyzed in terms of a Langevin equation where is a stochastic Gaussian field defined by the correlation functions with meaning stochastic average. The stochastic correlation functions derived from the Langevin equation correspond to a subset of the quantum correlation functions [3, 82].
It is also possible to study the master equation for the reduced density matrix of the system , which is given by [60, 82, 86]
where , , and are a frequency shift, a dissipative factor, and two dispersive factors, respectively, which in the weak coupling limit are given by [60] The corresponding expressions valid for arbitrary couplings can be found in [60, 82]. Recall that the time is the time at which the density matrix is assumed to factorize.
Alternatively, it is also possible to introduce the reduced Wigner function in terms of which the master equation adopts a FokkerPlanck form The Wigner function has some similarities with a classical distribution function, although it cannot be interpreted as a true probability density, because, among other reasons, it can adopt negative values [3].
Acknowledgments
The author is very grateful to Albert Roura for several helpful discussions and suggestions, and to Renaud Parentani and Enric Verdaguer for their comments on the original manuscript. This work is partially supported by the Research Projects MEC FPA200766665C0202 and DURSI 2005SGR00082.
References
 E. B. Davies, Quantum Theory of Open Systems, Academic Press, London, UK, 1976. View at: Zentralblatt MATH  MathSciNet
 H.P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press, New York, NY, USA, 2002. View at: Zentralblatt MATH  MathSciNet
 C. W. Gardiner and P. Zoller, Quantum Noise, Springer Series in Synergetics, Springer, Berlin, Germany, 2nd edition, 2000. View at: Zentralblatt MATH  MathSciNet
 U. Weiss, Quantum Dissipative Systems, World Scientific, River Edge, NJ, USA, 1999.
 C. Itzykson and J. B. Zuber, Quantum Field Theory, International Series in Pure and Applied Physics, McGrawHill, New York, NY, USA, 1980. View at: MathSciNet
 P. M. V. B. Barone and A. O. Caldeira, “Quantum mechanics of radiation damping,” Physical Review A, vol. 43, no. 1, pp. 57–63, 1991. View at: Publisher Site  Google Scholar
 C. Anastopoulos, “nparticle sector of field theory as a quantum open system,” Physical Review D, vol. 56, no. 10, pp. 6702–6705, 1997. View at: Google Scholar
 C. Anastopoulos and A. Zoupas, “Nonequilibrium quantum electrodynamics,” Physical Review D, vol. 58, no. 10, Article ID 105006, 12 pages, 1998. View at: Google Scholar  MathSciNet
 J. A. Wheeler and R. P. Feynman, “Classical electrodynamics in terms of direct interparticle action,” Reviews of Modern Physics, vol. 21, no. 3, pp. 425–433, 1949. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. Calzetta and B. L. Hu, “Noise and fluctuations in semiclassical gravity,” Physical Review D, vol. 49, no. 12, pp. 6636–6655, 1994. View at: Publisher Site  Google Scholar
 B. L. Hu and A. Matacz, “Back reaction in semiclassical gravity: the EinsteinLangevin equation,” Physical Review D, vol. 51, no. 4, pp. 1577–1586, 1995. View at: Publisher Site  Google Scholar
 B. L. Hu and S. Sinha, “Fluctuationdissipation relation for semiclassical cosmology,” Physical Review D, vol. 51, no. 4, pp. 1587–1606, 1995. View at: Google Scholar  MathSciNet
 A. Campos and E. Verdaguer, “Stochastic semiclassical equations for weakly inhomogeneous cosmologies,” Physical Review D, vol. 53, no. 4, pp. 1927–1937, 1996. View at: Google Scholar  MathSciNet
 R. Martín and E. Verdaguer, “Stochastic semiclassical fluctuations in Minkowski spacetime,” Physical Review D, vol. 61, no. 12, Article ID 124024, 26 pages, 2000. View at: Google Scholar  MathSciNet
 E. A. Calzetta, B. L. Hu, and F. D. Mazzitelli, “Coarsegrained effective action and renormalization group theory in semiclassical gravity and cosmology,” Physics Reports A, vol. 352, no. 46, pp. 459–520, 2001. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 B. L. Hu and E. Verdaguer, “Stochastic gravity: a primer with applications,” Classical and Quantum Gravity, vol. 20, no. 6, pp. R1–R42, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 B. L. Hu and E. Verdaguer, “Stochastic gravity: theory and applications,” Living Reviews in Relativity, vol. 7, p. 3, 2004. View at: Google Scholar
 E. Calzetta, A. Campos, and E. Verdaguer, “Stochastic semiclassical cosmological models,” Physical Review D, vol. 56, no. 4, pp. 2163–2172, 1997. View at: Google Scholar
 F. C. Lombardo and F. D. Mazzitelli, “EinsteinLangevin equations from running coupling constants,” Physical Review D, vol. 55, no. 6, pp. 3889–3892, 1997. View at: Google Scholar
 R. Martín and E. Verdaguer, “Stochastic semiclassical gravity,” Physical Review D, vol. 60, no. 8, Article ID 084008, 24 pages, 1999. View at: Google Scholar  MathSciNet
 R. Martín and E. Verdaguer, “On the semiclassical EinsteinLangevin equation,” Physics Letters B, vol. 465, no. 1–4, pp. 113–118, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. Calzetta, A. Roura, and E. Verdaguer, “Vacuum decay in quantum field theory,” Physical Review D, vol. 64, no. 10, Article ID 105008, 2001. View at: Publisher Site  Google Scholar
 E. Calzetta, A. Roura, and E. Verdaguer, “Dissipation, noise, and vacuum decay in quantum field theory,” Physical Review Letters, vol. 88, no. 1, Article ID 010403, 4 pages, 2002. View at: Google Scholar
 F. Lombardo and F. D. Mazzitelli, “Coarse graining and decoherence in quantum field theory,” Physical Review D, vol. 53, no. 4, pp. 2001–2011, 1996. View at: Google Scholar
 C. Greiner and B. Müller, “Classical fields near thermal equilibrium,” Physical Review D, vol. 55, no. 2, pp. 1026–1046, 1997. View at: Google Scholar
 F. C. Lombardo, F. D. Mazzitelli, and R. J. Rivers, “Decoherence in field theory: general couplings and slow quenches,” Nuclear Physics B, vol. 672, no. 3, pp. 462–486, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Matacz, “Inflation and the finetuning problem,” Physical Review D, vol. 56, no. 4, pp. R1836–R1840, 1997. View at: Google Scholar
 T. Tanaka and M.A. Sakagami, “Classical nature of the inflaton field with selfinteraction,” Progress of Theoretical Physics, vol. 100, no. 3, pp. 547–580, 1998. View at: Publisher Site  Google Scholar
 F. C. Lombardo and D. L. Nacir, “Decoherence during inflation: the generation of classical inhomogeneities,” Physical Review D, vol. 72, no. 6, Article ID 063506, 19 pages, 2005. View at: Publisher Site  Google Scholar
 J. Zanella and E. Calzetta, “A nonequilibrium renormalization group approach to turbulent reheating,” Journal of Physics A, vol. 40, no. 25, pp. 6927–6933, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory, AddisonWesley, Reading, Mass, USA, 1998.
 L. S. Brown, Quantum Field Theory, Cambridge University Press, Cambridge, UK, 1992. View at: MathSciNet
 S. Weinberg, The Quantum Theory of Fields, Cambridge University Press, Cambridge, UK, 1995.
 J. I. Latorre, P. Pascual, and R. Tarrach, “Speed of light in nontrivial vacua,” Nuclear Physics B, vol. 437, no. 1, pp. 60–82, 1995. View at: Publisher Site  Google Scholar
 D. Arteaga, R. Parentani, and E. Verdaguer, “Propagation in a thermal graviton background,” Physical Review D, vol. 70, no. 4, Article ID 044019, 15 pages, 2004. View at: Google Scholar  MathSciNet
 D. Arteaga, R. Parentani, and E. Verdaguer, “Retarded green functions and modified dispersion relations,” International Journal of Theoretical Physics, vol. 44, no. 10, pp. 1665–1689, 2005. View at: Publisher Site  Google Scholar
 M. Le Bellac, Thermal Field Theory, Cambridge Monographs on Mathematical Physics, Cambridge University Press, Cambridge, UK, 1996. View at: MathSciNet
 J. Bros, H. Epstein, and U. Moschella, “The lifetime of a massive particle in a de Sitter universe,” Journal of Cosmology and Astroparticle Physics, vol. 2008, no. 2, 8 pages, 2008. View at: Publisher Site  Google Scholar
 I. T. Drummond and S. J. Hathrell, “QED vacuum polarization in a background gravitational field and its effect on the velocity of photons,” Physical Review D, vol. 22, no. 2, pp. 343–355, 1980. View at: Google Scholar  MathSciNet
 G. M. Shore, “Quantum gravitational optics,” Contemporary Physics, vol. 44, no. 6, pp. 503–521, 2003. View at: Publisher Site  Google Scholar
 D. Arteaga, “Particle propagation in cosmological backgrounds,” Journal of Physics A, vol. 40, no. 25, pp. 6901–6906, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Arteaga, R. Parentani, and E. Verdaguer, “Particle propagation in cosmological backgrounds,” International Journal of Theoretical Physics, vol. 46, no. 9, pp. 2227–2241, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. Arteaga, “Quasiparticle excitations in relativistic quantum field theory,” Annals of Physics, vol. 324, no. 4, pp. 920–954, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Arteaga, “A field theory characterization of interacting adiabatic particles in cosmology,” Classical and Quantum Gravity, vol. 25, no. 15, Article ID 154014, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Campo and R. Parentani, “Inflationary spectra and violations of Bell inequalities,” Physical Review D, vol. 74, no. 2, Article ID 025001, 2006. View at: Publisher Site  Google Scholar
 D. Campo and R. Parentani, “Inflationary spectra and partially decohered distributions,” Physical Review D, vol. 72, no. 4, Article ID 045015, 25 pages, 2005. View at: Publisher Site  Google Scholar
 D. Campo and R. Parentani, “Spacetime correlations in inflationary spectra: a wavepacket analysis,” Physical Review D, vol. 70, no. 10, Article ID 105020, 2004. View at: Publisher Site  Google Scholar
 W. G. Unruh and W. H. Zurek, “Reduction of a wave packet in quantum Brownian motion,” Physical Review D, vol. 40, no. 4, pp. 1071–1094, 1989. View at: Google Scholar  MathSciNet
 L. D. Romero and J. P. Paz, “Decoherence and initial correlations in quantum Brownian motion,” Physical Review A, vol. 55, no. 6, pp. 4070–4083, 1997. View at: Publisher Site  Google Scholar
 A. O. Caldeira and A. J. Leggett, “Quantum tunnelling in a dissipative system,” Annals of Physics, vol. 149, no. 2, pp. 374–456, 1983. View at: Publisher Site  Google Scholar
 D. Arteaga, E. Calzetta, A. Roura, and E. Verdaguer, “Activationlike processes at zero temperature,” International Journal of Theoretical Physics, vol. 42, no. 6, pp. 1257–1280, 2003. View at: Publisher Site  Google Scholar
 E. Calzetta and E. Verdaguer, “Realtime approach to tunnelling in open quantum systems: decoherence and anomalous diffusion,” Journal of Physics A, vol. 39, no. 30, pp. 9503–9532, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. O. Caldeira and A. J. Leggett, “Influence of dissipation on quantum tunneling in macroscopic systems,” Physical Review Letters, vol. 46, no. 4, pp. 211–214, 1981. View at: Publisher Site  Google Scholar
 S. Massar, R. Parentani, and R. Brout, “On the problem of the uniformly accelerated oscillator,” Classical and Quantum Gravity, vol. 10, no. 2, pp. 385–395, 1993. View at: Publisher Site  Google Scholar  MathSciNet
 B. L. Hu and A. Matacz, “Quantum Brownian motion in a bath of parametric oscillators: a model for systemfield interactions,” Physical Review D, vol. 49, no. 12, pp. 6612–6635, 1994. View at: Google Scholar  MathSciNet
 D. F. Walls and G. J. Milburn, Quantum Optics, Springer, Berlin, Germany, 1994.
 A. O. Caldeira and A. J. Leggett, “Path integral approach to quantum Brownian motion,” Physica A, vol. 121, no. 3, pp. 587–616, 1983. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. P. Feynman and F. L. Vernon Jr., “The theory of a general quantum system interacting with a linear dissipative system,” Annals of Physics, vol. 24, pp. 118–173, 1963. View at: Publisher Site  Google Scholar
 R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, McGrawHill, New York, NY, USA, 1965.
 B. L. Hu, J. P. Paz, and Y. Zhang, “Quantum Brownian motion in a general environment: exact master equation with nonlocal dissipation and colored noise,” Physical Review D, vol. 45, no. 8, pp. 2843–2861, 1992. View at: Google Scholar  MathSciNet
 B. L. Hu, J. P. Paz, and Y. Zhang, “Quantum Brownian motion in a general environment. II. Nonlinear coupling and perturbative approach,” Physical Review D, vol. 47, no. 4, pp. 1576–1594, 1993. View at: Google Scholar  MathSciNet
 J. Schwinger, “Brownian motion of a quantum oscillator,” Journal of Mathematical Physics, vol. 2, no. 3, pp. 407–432, 1961. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. V. Keldysh, “Diagram technique for nonequilibrium processes,” Zhurnal Eksperimental noi i Teoreticheskoi Fiziki, vol. 47, pp. 1515–1527, 1965, translated as Soviet Physics JETP, vol. 20, pp. 1018–1026, 1965. View at: Google Scholar  MathSciNet
 E. Tomboulis, “$1/N$ expansion and renormalization in quantum gravity,” Physics Letters B, vol. 70, no. 3, pp. 361–364, 1977. View at: Publisher Site  Google Scholar
 F. Cooper and E. Mottola, “Initialvalue problems in quantum field theory in the large$N$ approximation,” Physical Review D, vol. 36, no. 10, pp. 3114–3127, 1987. View at: Publisher Site  Google Scholar
 F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz, and P. R. Anderson, “Nonequilibrium quantum fields in the large$N$ expansion,” Physical Review D, vol. 50, no. 4, pp. 2848–2869, 1994. View at: Publisher Site  Google Scholar
 H. A. Weldon, “Simple rules for discontinuities in finitetemperature field theory,” Physical Review D, vol. 28, no. 8, pp. 2007–2015, 1983. View at: Publisher Site  Google Scholar
 D. Arteaga, Particle propagation in nontrivial backgrounds: a quantum field theory approach, Ph.D. thesis, Universitat de Barcelona, Barcelona, Spain, 2007.
 A. O. Caldeira, H. A. Cerdeira, and R. Ramaswamy, “Limits of weak damping of a quantum harmonic oscillator,” Physical Review A, vol. 40, no. 6, pp. 3438–3440, 1989. View at: Publisher Site  Google Scholar
 C. Fleming, B. Hu, and A. Roura, “Solutions to master equations of quantum Brownian motion in a general environment with external force,” Quantum Physics. In press. View at: Google Scholar
 R. Parentani, “Constructing QFT's wherein Lorentz Invariance is broken by dissipative effects in the UV,” http://rxiv.org/abs/0705.2766. View at: Google Scholar
 R. Parentani, “Confronting the transPlanckian question of inflationary cosmology with dissipative effects,” Classical and Quantum Gravity, vol. 25, no. 15, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. Chou, Z. Su, B. Hao, and L. Yu, “Equilibrium and nonequilibrium formalisms made unified,” Physics Reports, vol. 118, no. 12, pp. 1–131, 1985. View at: Google Scholar  MathSciNet
 E. Calzetta and B. L. Hu, “Closedtimepath functional formalism in curved spacetime: application to cosmological backreaction problems,” Physical Review D, vol. 35, no. 2, pp. 495–509, 1987. View at: Google Scholar  MathSciNet
 A. Campos and B. L. Hu, “Nonequilibrium dynamics of a thermal plasma in a gravitational field,” Physical Review D, vol. 58, no. 12, Article ID 125021, 1998. View at: Google Scholar
 A. Campos and E. Verdaguer, “Semiclassical equations for weakly inhomogeneous cosmologies,” Physical Review D, vol. 49, no. 4, pp. 1861–1880, 1994. View at: Publisher Site  Google Scholar
 S. Weinberg, “Quantum contributions to cosmological correlations,” Physical Review D, vol. 72, no. 4, Article ID 043514, 19 pages, 2005. View at: Google Scholar  MathSciNet
 A. Das, Finite Temperature Field Theory, World Scientific, Singapore, 1997.
 I. D. Lawrie, “Peturbative description of dissipation in nonequilibrium field theory,” Physical Review D, vol. 40, no. 10, pp. 3330–3342, 1989. View at: Google Scholar  MathSciNet
 I. D. Lawrie and D. B. McKernan, “Nonequilibrium perturbation theory for complex scalar fields,” Physical Review D, vol. 55, no. 4, pp. 2290–2297, 1997. View at: Google Scholar
 N. P. Landsman and Ch. G. van Weert, “Real and imaginarytime field theory at finite temperature and density,” Physics Reports, vol. 145, no. 34, pp. 141–249, 1987. View at: Google Scholar  MathSciNet
 E. Calzetta, A. Roura, and E. Verdaguer, “Stochastic description for open quantum systems,” Physica A, vol. 319, no. 1–4, pp. 188–212, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Grabert, P. Schramm, and G.L. Ingold, “Quantum Brownian motion: the functional integral approach,” Physics Reports, vol. 168, no. 3, pp. 115–207, 1988. View at: Google Scholar  MathSciNet
 J. R. Anglin, J. P. Paz, and W. H. Zurek, “Deconstructing decoherence,” Physical Review A, vol. 55, no. 6, pp. 4041–4053, 1997. View at: Publisher Site  Google Scholar
 A. Roura, , Ph.D. thesis, Universitat de Barcelona, Barcelona, Spain, 2001.
 J. J. Halliwell and T. Yu, “Alternative derivation of the HuPazZhang master equation of quantum Brownian motion,” Physical Review D, vol. 53, no. 4, pp. 2012–2019, 1996. View at: Google Scholar  MathSciNet
Copyright
Copyright © 2009 Daniel Arteaga. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.