{'fraction_non_alphanumeric': 0.03914373088685015, 'fraction_numerical': 0.009206502494769033, 'mean_word_length': 4.247635135135135, 'pattern_counts': {'":': 0, '<': 8, '<?xml version=': 0, '>': 11, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Rest frame vacua of massive Klein–Gordon fields on spatially flat FLRW spacetimes

We propose a method of projecting the quantum states from a state space of a given geometry into another state space generated by a different geometry, taking care of the correct normalization which is crucial in interpreting the quantum theory. Thanks to this method we can define on any spatially flat FLRW spacetime states in which genuine Minkowskian parameters are measured. We use these Minkowskian states for separating the frequencies in the rest frames of the massive scalar particles defining thus the scalar rest frame vacuum. We show that this vacuum is stable on the de Sitter expanding universe where the energy is conserved. In contrast, on a spatially flat FLRW spacetime with a Milne-type scale factor this vacuum is found to be dynamic, corresponding to a time-dependent rest energy interpreted as an effective mass. This dynamic vacuum gives rise to cosmological particle creation which is significant only in the early Milne-type universe considered here. Some interesting features of this new effect are pointed out in a brief analysis.

The main task here is just the criterion of separating the frequencies defining the particle and antiparticle modes and, implicitly, the current vacuum at a given time. The principal method used so far is to focus mainly on the asymptotic states whose behavior is similar to the usual Minkowskian particle and antiparticle mode functions. In this manner one may choose in and out states whose frequencies are separated as in the flat case defining thus the adiabatic vacua as, for example, the Bunch-Davies one [24], used in many applications.
Recently we proposed the rest frame vacuum (r.f.v.) of the massive Dirac field on (1 + 3)-dimensional spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) spacetimes [43]. We started with the observation that in the rest frame, where the particle momentum vanishes, the solutions of the Dirac equation on any FLRW spacetime have a Minkowskian behavior regardless of the time evolution of the background. Thus we can separate the frequencies as in special relativity obtaining a stable time-independent vacuum on any FLRW spacetime. This vacuum is different from the Bunch-Davies one, which has largely been used since a long time ago including in our previous papers [44,45]. Note that the r.f.v. can be defined only for massive particles since the massless ones do not have rest frames.
The next step might be the generalization of the r.f.v. to the Klein-Gordon and Proca fields seen as perturbations on the mentioned manifolds. Unfortunately, here we are faced with a serious difficulty since, in contrast with the Dirac field, the rest mode functions of these bosonic fields do not have the desired Minkowskian forms we need for defining the r.f.v. in a natural manner. Nevertheless, since our concept of particle and antiparticle comes from the Minkowskian quantum field theory, we are forced to impose the Minkowskian forms to the rest mode functions even though this is possible only at a given time, assuming that the time-dependent rest energy represents a dynamical effective mass.
Another challenge is to solve the ambiguity related to the normalization of these Minkowskian states, which can be done in two manners, either with respect to the scalar product of the curved manifold or by using the Minkowskian scalar product. This difficulty can be avoided since the state spaces are separable Hilbert spaces which are isometric among themselves such that we may look for an appropriate method of mapping the state spaces produced in different geometries into a unique one where we may calculate transition coefficient between these states. In what follows we would like to concentrate on these problems proposing a method of defining well-normalized Minkowskian rest mode functions on any spatially flat FLRW spacetime. These will help us to define at any moment the bosonic time-dependent r.f.v. associated to a time-dependent dynamical mass. When this is time-independent we say that the r.f.v. is stable. For concrete calculations we restrict ourselves to the massive Klein-Gordon field minimally coupled to the background gravity.
In our approach the r.f.v. of the scalar field are stable only on the FLRW manifolds where the energy is conserved as eigenvalue of a conserved energy operator given by a timelike Killing vector field, as in the case of the de Sitter expanding universe. In other FRLW manifolds we meet c.p.c. processes that may be studied deriving the Bogolyubov coefficients between states whose vacua are defined at two arbitrary moments, obtaining thus information about the time behavior of the c.p.c. in any FLRW geometry. For illustrating how our method works we give two examples, the stable r.f.v. on the de Sitter expanding universe and, for the first time, we present an example of time-dependent r.f.v. on a spatially flat FLRW with a Milne-type scale factor where only the r.f.v. of the Dirac field was studied [43]. On this last manifold we study the c.p.c. at finite times obtaining probabilities and rates which depend exclusively on the moments when the particle is prepared and then measured. Note that our results are different from other attempts of studying c.p.c. at finite times as in a series of papers reviewed recently [38,39] or in Ref. [40] where the vacuum depends, in addition, on kinetic parameters (as for example the momentum components).
We start in the next section presenting our basic assumptions concerning the scalar quantum modes prepared or measured by a global apparatus on a curved manifold and showing how the state space can be mapped into a Minkowskian one. The next section is devoted to the spatially flat FLRW spacetimes where we propose a concrete method of defining Minkowskian rest states, correctly normalized at a given arbitrary time, regardless the time evolution of the background geometry. By using such states we define the r.f.v. showing that these vacua are stable only on the FLRW manifolds where the energy is conserved. Section 4 is devoted to a general method of setting the r.f.v. deriving the Bogolyubov coefficients between two bases of mode functions whose frequency separation in rest frames was performed at two different moments. In the Section 5, the mentioned examples are studied, namely the stable r.f.v. on de Sitter expanding universe and the time-dependent one on the spatiallyflat Milne-type universe. In the last part of this section we present the c.p.c. produced by the vacuum instability discussing some physical consequences in a short analytical and graphical analysis. In the last section we present our concluding remarks.

Minkowskian scalar modes
Respecting ad litteram the principles of the quantum theory, we assume that the quantum states on any local Minkowskian manifold, (M, g), are prepared or measured by a global apparatus represented by the algebra of the quantum observables, i.e. the Hermitian operators defined globally as vector fields on the whole manifold or on a portion with an independent physical meaning, as in the case of the de Sitter expanding universe. The operators proportional to the Killing vector fields are conserved, commuting with the operator of the field equation. Our global apparatus prepares quantum modes whose mode functions are common eigenfunctions of a system of commuting conserved operators (s.c.c.o.) {E K G , A, B, ...}, which includes the operator of the field equation E K G . In addition, these mode functions are supposed to be normalized with respect to a specific relativistic scalar product on (M, g).
In general, the s.c.c.o. determining the quantum modes are not complete such that the mode functions remain with some integration constants which depend on the separation of the positive and negative frequencies defining the vacuum. Another possible manner of setting these constants is by defining the modes on (M, g) in which one measures the parameters corresponding to another geometry (M,ĝ), according to the method we present in what follows.
We start with the (1 + 3)-dimensional local Minkowskian manifold (M, g) where we consider a local chart {x} of coordinates x μ (labeled by natural indices α, ..., μ, ... = 0, 1, 2, 3) with x 0 = t and arbitrary space coordinates. The scalar field, Φ : M → C, of mass m, minimally coupled to the gravity of (M, g), satisfies the Klein-Gordon equation E K G Φ = m 2 Φ whose operator is defined as where the particle, a, a † , and antiparticle, a c , a c † , field operators must satisfy the canonical bosonic commutation relations [46]. The mode functions, f ∈ K , behave as tempered distributions or square integrable functions with respect to the indefinite Hermitian form written with the notation f From the physical point of view the mode functions of K ± have positive/negative frequencies, while those of K 0 do not have a physical meaning. For any f ∈ K + we have In fact, H is a Krein space [47] while K ± are the spaces of tempered distributions of the rigged Hilbert spaces (or Gel'fand triplets) [48,49] associated to the Hilbert spaces H ± equipped with the scalar products ± , M .
A complete system of orthonormal mode functions, { f α } α∈I ⊂ K + , forms a (generalized) basis of positive frequencies in K + related to the negative frequencies one, In this manner one defines a frequency separation associated to a specific vacuum state of the Fock space. It is well known that two different bases define different vacuum states when these are related among themselves through a non-trivial Bogolyubov transformation that mixes the positive and negative frequency modes. Otherwise the vacuum state remains stable.
Furthermore, we consider another manifold (M,ĝ) whose local chart {x} is defined on the same domain of the flat model as the chart {x} of (M, g). This means that there exists a coordinate transformationx = χ(x) allowing us to relate the set K discussed above to the setK of the scalar mode functions on (M,ĝ) equipped with the Hermitian form , M , defined as in Eq. (3). We observe that the physical parts of the setŝ K and K are separable Hilbert spaces between which we can define the isometry μ : H + →Ĥ + which satisfies Then for any normalized mode functions f α ∈ H + andf β ∈ H + which obey we can construct the amplitude which, in general, depends on time. This gives the quantity | α|β t | 2 , which can be interpreted as the probability of measuring at the time t the parameters β in the state α prepared on (M, g) or, reversely, as the probability of measuring the parameters α in the state β prepared on (M,ĝ). For this rea- The isometry μ is complicated since this involves the coordinate transformationx = χ(x) but which can be eliminated by choosing the same coordinates for the both manifolds under consideration by taking χ = id →x = x. Note that this is possible, since we assumed that the local charts of (M, g) and (M,ĝ) are included in the same domain of the flat model. With this choice the isometry takes the simple form which can be used in applications. An important particular case is when (M,ĝ) is just the Minkowski spacetime, which is the flat model of (M, g). Then we can set at any time χ = id and, in addition, we get the opportunity of defining in (M, g) states in which one measures exclusively Minkowskian parameters at a given time t 0 . Thus for any normalized mode functionf ∈K on the Minkowski spacetime we may define the corresponding Mikowskian state on (M, g) whose normalized mode function f ∈ K is defined such that the functions andf have a contact of order k at the time t 0 , satisfying the system of k + 1 algebraic equations able to give all the integration constants of f in terms of the Minkowskian parameters of the functionf we chose. Obviously, the number k+1 of equations we may use depends on the number of the undetermined integration constants, left by the incomplete s.c.c.o., or other parameters we need to find as, for example, the rest energy. With this method we can apply the definitions of Minkowskian particles or antiparticles to any manifold (M, g) but only at a given time since, in general, these states are evolving in time.

Rest frame vacua
Let us consider now the family of (1 + 3)-dimensional spatially flat FLRW spacetimes for which we use the same coordinates of the FLRW chart, {t, x}, i.e. the proper (or cosmic) time t ∈ D t and the Cartesian space coordinates We denote by M the spacetime whose line element depends on the scale factor a(t) which is assumed to be a smooth function on D t giving for the conformal time of the conformal chart {t c , x}. The line elements of these charts are where we denoted a(t c ) = a[t (t c )]. The Minkowski spacetime, denoted from now on simply asM, is the particular case when a(t) = 1 and t c = t.
In the chart {t, x} the massive scalar field Φ : M → C of mass m satisfies the Klein-Gordon equation which allows for a system of plane wave solutions, i.e. eigenfunctions of the momentum operators P i = −i∂ i corresponding to the eigenvalues ( p 1 , p 2 , p 3 ), representing the components of the conserved momentum p. These mode functions can be written as in terms of the time modulation functions F p : D t → C which depend on p = |p| satisfying the equation where This equation does not determine completely the form of the functions F p , there remaining integration constants which have to be fixed by supplemental assumptions defining the vacuum.
The fundamental solutions (14) form an orthonormal basis with respect to the scalar product (3) that now reads allowing us to impose the normalization condition requiring the time modulation functions to satisfy Then the Klein-Gordon field can be expanded as in terms of the particle a, a † and antiparticle a c , a c † field operators which satisfy the canonical commutation relations In the particular case of the Minkowski spacetimeM the mode functions of positive frequencies of a scalar field of massm, are eigenfunctions of the energy operator i∂ t depending on the conserved energy E = p 2 +m 2 and satisfying the orthonormalization condition with respect to the scalar product On the other hand, we have shown that in any FLRW spacetime there exists an energy operator that in the FLRW chart, {t, x}, has the form [50,51] In general, this operator does not commute with the momentum p but in the rest frames (where p = 0) this coincides with the Minkowski one, H 0 = i∂ t , suggesting us to determine the integration constants of the solutions (14) by separating the frequencies just in rest frames by using the Minkowskian rest states on M defined in the previous section. Thus we may set the r.f.v. of the Klein-Gordon field on the FLRW manifolds under consideration. Without introducing new notations we suppose that now the mode functions (14) are the Minkowskian states in which one measures in the rest frame, at the time t 0 , the parameters of the mode functions (23) for p → 0 but with another rest energy,m = m, which we call the dynamical mass. Therefore, we may consider the system (10) with k = 2 giving the following equations: which suffice for separating the frequencies in the rest frame and finding the dynamical massm(t 0 ). Thus the first two equations give the normalized integration constants corresponding to the r.f.v., while the third one helps us to find the associated dynamical mass in the rest frame. All these quantities may depend on the time t 0 when we impose the Minkowskian form of the mode functions in the rest frame. This means that, in general, the r.f.v. is dynamic, being associated with a time-dependent dynamical mass. Nevertheless, this vacuum becomes stable on the FLRW manifolds where the energy operator is conserved, i.e. the Minkowski and de Sitter spacetimes, since then the energy operator in the rest frame commutes with P completing thus the s.c.c.o. but only in the rest frame.

Applications
For solving concrete examples we may start with a time modulation function of the general form where φ p is a particular solution satisfying The normalized solutions of positive frequency, f p ∈ K + , have time modulation functions which satisfy The problem is to derive these constants and the dynamical mass from the system (26).
In the rest frame (where p = 0) we denote simply φ = φ p | p=0 such that the system (26) can be written as The first two equations give the normalized integration constants corresponding to the r.f.v., while the third one gives us the associated dynamical mass in the rest frame, sinceφ = −Ω 2 φ as in Eq. (15). Thus we find that a particle prepared in r.f.v. at the time t 0 has the mode function whose time modulation function depends on the integration constants (33) and (34) which comply with the normalization condition The set { f p,t 0 |p ∈ R 3 } forms a basis in K + , while the set { f * p,t 0 |p ∈ R 3 } is the corresponding basis of K − in the r.f.v. prepared at t = t 0 .
As mentioned, the r.f.v. becomes stable only in the de Sitter and Minkowski spacetimes where the energy operators are conserved, satisfying [H 0 , Ω] = i∂ t Ω = 0. On other FLRW manifolds, the r.f.v. is dynamic, being associated with a timedependent dynamical massm(t) = Ω(t) ∈ R. The time domain D t = D + t ∪ D − t is split into the tardyonic part D + t = {t|Ω(t) 2 > 0} and the tachyonic one, D − t = {t|Ω(t) 2 < 0}. All the tachyonic states with Ω(t) = i|Ω(t)| are eliminated, having null norms. Thus in r.f.v. the scalar field survives only on D + t . When the r.f.v. is dynamic this gives rise to c.p.c., which can be pointed out studying how the states prepared in the r.f.v. at the time t 0 ∈ D + t are measured in the same vacuum at the time t > t 0 . For this purpose we have to calculate the Bogolyubov transformation between the bases { f p (t 0 )} and { f p (t)} which gives the following coefficients: since the system of particular solutions {φ p , φ * p } is supposed to be orthonormal. Then the density of the new particles or antiparticles created between t 0 and t is proportional to In addition, we observe that the rate of c.p.c. can also be estimated as Thus we obtain densities and rates independent on the kinetic parameters as the momentum components. This result is different from those obtained according to the method of Hamiltonian diagonalization by using Bogolyubov transformations [3,4,[9][10][11]14]. This is because in our case the Hamiltonian is diagonal only in the rest frame, while in the mentioned approach the Hamiltonian is put in diagonal form in any frame. It is remarkable that the integration constants of the corresponding mode functions have the same form as in Eqs. (33) and (34) but with φ p and p instead of φ and . This means that the Hamiltonian diagonalization may be performed only by the Minkowskian states as defined above. Thus it is obvious that the difference between these two approaches comes from the fact that the Hamiltonian diagonalization considers exclusively Minkowskian states in all the frames, for any momenta, while in the r.f.v. the only Minkowskian states are those of the rest frames, where p = 0.

Examples
Now we focus on two typical examples, namely the stable r.f.v. of the de Sitter expanding universe in contrast with the dynamical r.f.v. of a spatially flat FLRW manifold with a Milne-type scale factor where the behavior of the scalar field was never studied.

de Sitter expanding universe
Let us consider first the stable r.f.v. on the expanding portion of the de Sitter spacetime, M, having the scale factor a(t) = e 2ωt (where ω is the Hubble de Sitter constant in our notation) defined for t ∈ (−∞, ∞), giving the conformal time t c and the function a(t c ) as In the conformal chart the Klein-Gordon equation is analytically solvable giving the mode functions of the momentum basis of the form (14) having the time modulation functions where By using Eq. (A.4) we find that the set {φ p , φ * p } is orthonormal such that the normalization condition (19) is fulfilled only if we take We assume first that m > 3 2 ω solving the system (26) in the conformal chart {t c , x} where the de Sitter time modulation function has the form (44) with ν = iκ, while the Minkowski one (23) readŝ Moreover, since in this case the limit to p → 0 is sensitive, we solve first this system for p = 0 and then we evaluate this limit. From the first two equations we obtain the integration constants while from the last one giving the expected dynamical mass, related to the well-known rest energy [51]. Then, for p → 0 we obtain the constants which have the absolute values where the general phase factor was introduced for ensuring the correct limit for p → 0 as given by Eq. (A.6). These functions are correctly normalized since the integration constants (52) and (53) satisfy the condition (46). Note that these results can be rewritten in terms of the cosmic time t according to Eq. (43). Furthermore, we consider the case of m < 3 2 ω applying the same method for fixing the r.f.v.. We solve the first two equations of the system (26) for p = 0 and an arbitrary time t c obtaining From the third equation we find the expected condition giving the tachyonic dynamical massm = ±iνω. Moreover, we find that in the rest frame we have which means that if we set the r.f.v. then the particles with m < 3 2 ω cannot survive on the de Sitter expanding portion. The above results can be now gathered in the synthetic form of the mode functions of positive frequency in the conformal chart, that hold for any real or imaginary value of ν, given by Eq. (45). In the tachyonic case, when ν takes real values, the squared norm of f p vanishes, since then we have Similar results can be obtained for other couplings which modify only the dynamical mass and implicitly the limit between the tachyonic and tardyonic domains.
Thus we have shown that the scalar r.f.v. on the de Sitter expanding universe is stable corresponding to a timeindependent dynamical mass (51), which does make sense only when m > 3 2 ω. In other words, the frequency separation in the rest frames can be done only for the scalar fields which satisfy this condition. Otherwise we have either to eliminate the scalar fields with m < 3 2 ω or to resort to another vacuum as the adiabatic Bunch-Davies one [46] which can be set for particles of any mass by taking c 1 = 1 and c 2 = 0.
The principal consequence of the vacuum stability is that on the de Sitter expanding universe the c.p.c. is inhibited when we set the r.f.v.. Then the integration constants are independent on time such that the density (41) and the rate (42) of the created particles vanish. Note that this is a specific property of the r.f.v. since there are other situations when the c.p.c. may be active as, for example, in the case of the states which diagonalize the Hamiltonian in any frame [25,26,34].

Milne-type spatially flat FLRW spacetime
Let us consider now an example of manifold M where we do not have adiabatic vacua remaining only with an unstable r.f.v. corresponding to a time-dependent dynamical mass. This is the (1+3)-dimensional spatially flat FLRW manifold with the scale factor a(t) = ωt determining the conformal time as The constant ω, introduced because of dimensional considerations, is a useful free parameter which in the case of the genuine Milne universe (of negative space curvature) must be fixed to ω = 1 for eliminating the gravitational sources [46]. This spacetime M is produced by isotropic gravitational sources, i.e. the density ρ and pressure p, evolving in time as and vanishing for t → ∞. These sources govern the expansion of M that can be better observed in the chart {t,x}, of 'physical' space coordinatesx i = ωt x i , where the line element lays out an expanding horizon at |x| = t and tends to the Minkowski spacetime when t → ∞ and the gravitational sources vanish.
In the FLRW chart {t, x} of this spacetime the Klein-Gordon equation is analytically solvable, the fundamental solutions having the time modulation functions which can take either real or pure imaginary values. We must specify that here we cannot speak of adiabatic vacua as long as the functions (63) are singular in t = 0. Therefore, we must focus only on the r.f.v. for which the time-dependent integration constants, result from Eqs. (33) and (34). The corresponding dynamical mass readŝ The functions (65) and (66) are singular in t = 0 and t = t m ≡ √ 3 2m whenm(t) vanishes (as in Fig. 1). From Eq. (67) we see that a particle of mass m shows a tachyonic behavior in the domain D − t = (0, t m ) and a tardyonic one only if t ∈ D + t = (t m , ∞). As in the general case, we can verify that showing that on the tachyonic domain the wave function is of null norm, having thus no physical meaning. This means that the scalar particles can be prepared only in the tardyonic domain t > t m wherem(t) increases with t such that, for t → ∞, when M becomes just the Minkowski spacetime, it tends to m. Moreover, in this limit we recover the usual Minkowski scalar modes since the functions K behave as in Eq. (A.5) such that All these results can be encapsulated in the definitive form of the mode functions of positive frequency, prepared at the time t 0 > t m and defined for t > t 0 , which read where ν depends on p as in Eq. (64). The instability of r.f.v. on this expanding manifold gives rise to a c.p.c. that can be analyzed thanks to our previous results that hold for any t > t 0 . We can study how the particles created at t 0 can be measured at any moment t > t 0 calculating the density (41) and the rate (42), which depend now on the integration constants (65) and (66). We plot these functions in Figs. 2 and 3 as functions of ct instead of t for avoiding too small numbers. Furthermore, we observe that the dynamic r.f.v. tends to stability when the time is increasing, since then as we deduce from Eqs. (69). Moreover, we see that the dynamical effect is visible only for the very old particles, prepared at t 0 < 5 t m , since the function c 2 (t 0 ) decreases rapidly to zero when t 0 increases andm(t 0 ) → m. Thus for the younger particles, prepared at t 0 > 5 − 10 t m , the dynamical effect is suppressed leaving us with an apparently stable r.f.v. of the Bunch-Davies-type (with c 1 = 1 and c 2 = 0), in which the mode functions can be approximated as independent on the moment t 0 when the particle was prepared.
Finally we must specify that, in general, the dynamic effect discussed above is very fast, during an extremely short period of time, even at the quantum scale, since by definition t m = √ 3 2m (or √ 3 2 mc 2 in SI units) is very small. For example, if we take m to be just the electron mass m e , then t m e ∼ 1.1 10 −21 s such that for the particles born at cosmic times t 0 > 10 −20 s the r.f.v. is apparently stable. Only the particles prepared at t 0 < 10 −20 s present this effect, which may be of interest only at quantum scale in the cosmology of the very early spatially flat Milne-type universe.

Concluding remarks
We proposed here a method of projecting the quantum states from a state space of a given geometry into another state space generated by a different geometry, keeping the correct normalization which is crucial in interpreting the quantum quantities (probabilities, expectation values, transition amplitudes, etc.). This method helped us to define, on any spatially flat FLRW spacetime, the Minkowskian states we need for setting the r.f.v. of the massive scalar field which, in contrast to the Dirac one, does not show a Minkowskian behavior in the rest frames of the FLRW manifolds. In this manner, we obtained a stable r.f.v. on the de Sitter expanding universe and, for the first time, we found a dynamical vacuum, corresponding to a time-dependent dynamical mass on a Milne-type spacetime. In this last case, the dynamic r.f.v. gives rise to a very fast c.p.c., which could be of interest but only in a very early Milne-type universe. It is remarkable that in r.f.v. all the possible tachyonic behaviors (e.g. for m < 3 2 ω in the de Sitter case and t < t m in the Milne-type universe) are eliminated in a natural manner, the corresponding mode functions having null norms. These results may improve the study of the c.p.c. on the FLRW manifolds combining the r.f.v. with the other vacua proposed so far.
On the other hand, we must stress that the r.f.v. cannot be defined for the massless fields, which do not have rest frames. In the case of the Maxwell and massless Dirac fields this is not an impediment, since the neutrino and Maxwell equations are conformally covariant such that in the conformal charts of the FLRW spacetimes one may take over the frequency separation from the flat case. The only problem which remains partially unsolved is the vacuum of the massless scalar field whose equation is no longer covariant under conformal transformations. This sensitive case is being revisited time and again in the hope of finding a convenient interpretation [52].
Another approach is the quantum theory of interacting fields on curved manifolds in which the amplitudes of the quantum transitions can be calculated by using perturbations in terms of free fields [53][54][55][56][57][58][59][60] as in our recent de Sitter QED [45,61,62]. Even though in this framework only adiabatic vacua were considered so far, we have now the opportunity of using many types of vacua for improving the calculation of the transition amplitudes. Thus, for example, in a collision process we may take the incident beam in the adiabatic vacuum and the target in the r.f.v.. Moreover, for the internal lines of the Feynman diagrams the r.f.v. is the favorite candidate since this can be defined naturally for the massive fields on any spatially flat FLRW spacetime. Thus, by using many well-defined vacua, we could combine the methods of c.p.c. with those of the perturbative quantum field theory for analyzing various quantum effects in evolving universes.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The paper is theoretical having no data to be deposited.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .