partly motivated by our increased capability in recent years to perform accurate observations on the sky , we have turned our attention to investigating the impact of quantum effects on the evolution of the universe . in pursuit of this goal , quantum field theory in cosmological spacetimes continues to be an increasingly important field of research . among other things
, it investigates how quantum fluctuations affect the background spacetime in perturbative quantum gravity , a process also known as quantum backreaction .
a line of research deals for example with the backreaction of quantum fields whose spectrum is nearly flat .
examples of these fields are the minimally coupled massless scalar and the graviton .
consequently , these fields are expected to yield a substantial backreaction @xcite . also , fermions in curved spacetimes and in particular maximally symmetric spacetimes were investigated ( see e.g. : @xcite ) .    another approach to this question
is concerned with calculating the possible effect of the trace anomaly on the background spacetime .
some authors argue that this effect , particularly in relation with the cosmological constant problem , could be significant @xcite , however also see @xcite .
an essential element to study quantum effects in curved spaces is the propagator . due to its high degree of symmetry
, de sitter spacetime proves to be ideally suited for calculating various quantum effects .
moreover , results are immediately applicable to inflation , cosmologists favourite paradigm for a brief exponentially fast expansion in the early universe .
chernikov and tagirov calculated the scalar propagator in de sitter spacetime @xcite , also see @xcite .
the vector propagator was constructed by @xcite and the graviton propagator received contributions from @xcite . candelas and raine calculated the fermion propagator @xcite .
however , de sitter spacetime suffers from several drawbacks .
firstly , pure de sitter spacetime is never realised in nature .
since de sitter spacetime corresponds to a globally constant hubble parameter @xmath4 , this is in reality never attained . secondly , de sitter spacetime is non - dynamical .
it is therefore inconsistent to study backreaction effects while at the same time assuming a constant background .
we should allow the background to change to encompass all possible backreaction effects consistently .
therefore , it is necessary to consider propagators in more general spacetimes .
the scalar and graviton propagator in quasi de sitter spacetime were calculated by @xcite and subsequently generalised to friedmann - lematre - robertson - walker or flrw spacetimes with constant deceleration by @xcite . in this paper
we calculate the fermion propagator in flrw spacetimes with constant deceleration .
we should touch upon an important issue . in the massive case
we need two additional constraints to be satisfied in order to solve for the propagator :    [ assumptions ] @xmath5    recall that all interesting epochs in the evolution of our universe satisfy the first constraint . in the matter era
we have @xmath6 , in the radiation dominated epoch we find @xmath7 .
it also serves as an approximation to both inflation and the current dark energy dominated epoch when @xmath8 . note that @xmath9 coincides with the slow - roll parameter of inflation and it is straightforwardly related to the somewhat more familiar deceleration parameter @xmath10 .
hence ( [ epsilon ] ) is equivalent to requiring a constant deceleration parameter .    the second assumption ( [ mhassumption ] ) is required to analytically derive the propagator in the massive case . in yukawa theory , the mass of the fermion is generated by a scalar field for which we can achieve @xmath11 in several cases to which we will turn shortly .
we then calculate the one loop effective action induced by fermions using our propagator .
the one loop backreaction arises from integrating out a free , quadratic fermion field and , using dimensional regularisation , this generates a correction to the ( classical ) friedmann equations .
the hubble parameter occurring in this effective action is now a dynamical quantity .
consequently , the backreaction of these fermions on the background spacetime can thus be analysed dynamically .
also , when the fermions are coupled to a scalar field , this opens up the possibility to study the impact of the fermions on the evolution of scalar fields .
this paper is outlined as follows . in section
[ fermions in flrw spacetimes ] we review the basic theory required to study fermions in curved spaces and we establish our notation . sections [ flrw fermion propagator : the massless case ] and [ flrw fermion propagator : the massive case ] are devoted to deriving the fermionic propagator in flrw spacetimes in the massless and massive case , respectively .
finally , in section [ one loop effective action ] we apply our propagator to calculate the one loop effective action .
fermions are in general @xmath12 dimensional curved spacetimes described by the action @xcite : @xmath13 - m \bar{\psi}\psi\right\ } , \ ] ] where the dirac matrices @xmath14 satisfy the following anti - commutation relations : @xmath15 variation with respect to @xmath16 yields the equation of motion the fermion field @xmath17 satisfies : @xmath18 we will make use of the vierbein formalism which can be thought of as a transformation of the metric tensor to a locally flat minkowski metric : @xmath19 where @xmath20 is the minkowski metric .
we use greek letters to run over spacetime indices whereas we use latin letters either for the tangent space indices or for spinor indices
. it will be clear from the context which interpretation we are using .
the minkowski metric and the ( flat space ) dirac matrices @xmath21 are spacetime independent .
spacetime and tangent space indices are raised or lowered by making use of the full metric or minkowski metric , respectively .
the covariant derivative acting on a dirac spinor is defined as : @xmath22 where the spin connection is given by : @xmath23 \,,\ ] ] such that @xmath24 .
let us specialise to flat flrw spacetimes in which the metric is given by @xmath25 . here , @xmath26 is the scale factor of the universe in conformal time defined by @xmath27 . in flrw spacetimes the vierbeins are a function of conformal time only :    @xmath28    we thus find : @xmath29 this useful identity relates covariant derivatives to partial derivatives in the tangent space .
since we work in the chirality representation , the flat space dirac matrices are given by : @xmath30 where @xmath31 is the @xmath32 identity matrix and : @xmath33 and where the @xmath34 denote the pauli matrices .    ) assuming @xmath35 .
near the minimum of the potential ( [ phiproptoh ] ) , assumption ( [ mhassumption ] ) is satisfied .
[ fig : potential ] ]    _ _    in yukawa theory the mass of the fermion appearing in equation ( [ fermionaction ] ) is time dependent and generated by a scalar field . let us thus briefly consider the complex system of a fermion whose mass is generated by a homogeneous scalar field @xmath36 , both of which are coupled to gravity .
the action is : @xmath37 \right\}\label{actiontotal } .
\phantom{1}\end{aligned}\ ] ] here @xmath38 the yukawa coupling constant .
by comparing equations ( [ actiontotal ] ) and ( [ fermionaction ] ) we can thus identify the mass of the fermion as : @xmath39 furthermore , the ( classical ) evolution of the homogeneous background field @xmath36 is governed by : @xmath40 where @xmath41 is the potential : @xmath42 here , @xmath43 denotes some non - minimal coupling to gravity and @xmath44 is the ricci scalar . in flrw spacetimes in @xmath45 we have @xmath46 , where the hubble parameter is given by @xmath47 .
dots denote derivatives with respect to cosmic time and dashes correspond to conformal time derivatives . assuming @xmath48 , the solution of ( [ eomphi ] ) is : @xmath49}\,.\ ] ] clearly , @xmath11 and assumption ( [ mhassumption ] ) is satisfied .
moreover , observe we need @xmath50 .
the potential is depicted in figure [ fig : potential ] . given a potential of the form ( [ potential ] ) , we identify some interpretations leading to this scenario . in the early universe the higgs field satisfies @xmath51 and hence ( [ phiproptoh ] ) .
a quintessence field has @xmath52 and ( [ phiproptoh ] ) is also satisfied in the late universe .
finally , we can infer that a massless minimally coupled scalar field is also allowed in the regime where @xmath53 or @xmath54 .
we work in the approximation where @xmath55 , which is not equivalent to the standard slow - roll approximation of inflation , where one introduces the slow - roll parameters @xmath9 and @xmath56 .
these parameters are allowed to vary slowly in time such that @xmath57 and @xmath58 are higher order in the slow - roll parameters and can be neglected when one is interested in the leading order behaviour only .
when @xmath9 is allowed to vary in time , equation ( [ phiproptoh ] ) generalises to @xmath59 , where : @xmath60 } + \mathcal{o}\left(\dot{\alpha}\right ) \,.\ ] ] the leading order slow - roll approximation to this equation corresponds to : @xmath61 } \,.\ ] ] note the equation above depends on @xmath9 , whereas it does not depend on @xmath56 .      depending on which pole prescription one uses in the spirit of the schwinger - keldysh formalism ( see e.g. @xcite )
, one can obtain different propagators such as the anti - time ordered , time ordered or wightman propagators .
we are primarily interested in the time ordered or feynman fermion propagator , defined by : @xmath62 the symbols @xmath63 and @xmath64 are the spinor indices and @xmath65 denotes the state of the system .
also note the minus sign in front of the second @xmath66-function .
the fermion propagator at tree level satisfies : @xmath67 is_{\mathrm{f}}^{ab } ( x,\tilde{x } ) = i \delta^{\scriptscriptstyle{d}}\!(x-\tilde{x } ) \mathbb{i}^{ab}\,,\ ] ] and a likewise equation when @xmath68 and @xmath69 are interchanged .
we already discussed that flrw spacetimes can be characterised by the scale factor @xmath63 as a function of either cosmic or conformal time .
there are however a few more important properties of this spacetime that are worth mentioning here .
an important relation that we will use throughout this manuscript is : @xmath70 this equation is equivalent to assumption ( [ epsilon ] ) , see @xcite .
let us define some relevant geometrical functions :    [ y ] @xmath71    which all vanish near the lightcone . here
, @xmath72 in ( [ y ] ) refers to the feynman or time ordered pole prescription of the feynman pole prescription is unrelated to the `` slow roll '' parameter @xmath9 in flrw spacetimes . ] . in de sitter spacetime ,
i.e. : when @xmath73 , the de sitter invariant function @xmath74 is related to the geodesic length @xmath75 as @xmath76 .
massless fermions in flrw spacetimes are not difficult to deal with . by making use of equation ( [ covariantderivative2 ] ) the propagator for massless ( conformal ) fermions should satisfy : @xmath77 let us recall the following identity : @xmath78 we can
thus immediately infer the solution of equation ( [ fermionpropmassless ] ) : @xmath79 here , @xmath80 .
this completes the calculation for the propagator of massless fermions in flrw spacetimes in @xmath12 dimensions .
this propagator is valid in any flrw spacetime and assumption ( [ epsilon ] ) can be relaxed . because massless fermions are conformal in any dimension , their propagator is much easier to calculate than the massless scalar propagator @xcite .
a massive fermion is not conformal and hence its propagator in general contains some complicated mass dependence ) is much more complicated than for scalar fields . ] .
we firstly solve for the fermionic mode functions .
we generalise the approach outlined in @xcite to incorporate flrw spacetimes with constant @xmath9 . using these mode functions
, we then return to position space to construct the feynman propagator .
let us firstly define the rescaled fermionic spinor : @xmath81 keeping an eye on equation ( [ covariantderivative2 ] ) one can easily check that this factor is chosen such to conveniently transform the covariant derivative into a partial derivative .
equation of motion ( [ fermioneom ] ) thus reads : @xmath82 the canonical momentum associated with @xmath83 follows from promoting the variation of ( [ fermionaction ] ) with respect to @xmath84 to an operator . according to the usual moves in quantum field theory
, we impose anti - commutation relations between @xmath83 and its associated canonical momentum : @xmath85 with the other anti - commutators vanishing . by making use of the rescaling in equation ( [ fermionrescaled ] ) ,
appreciate that the anti - commutation relations above simplify :    [ anticommutation2 ] @xmath86      therefore we expand the rescaled spinors @xmath87 and @xmath88 in creation and annihilation operators as follows :    [ creationannihilationexpansion ] @xmath89    the expansion above merits a few remarks .
firstly , @xmath90 and @xmath91 are the fermion and anti - fermion annihilation operators of helicity @xmath92 , respectively , in the usual sense : @xmath93 .
the helicity @xmath92 , i.e. : the spin in the direction of motion , can be either @xmath94 or @xmath95 in units of @xmath96 .
let us for simplicity return to the @xmath45 setting we are familiar with .
we will generalise the following considerations shortly to arbitrary dimensions .    in the equation above @xmath97
is a 4-spinor of helicity @xmath92 .
we decompose the 4-spinor @xmath98 into a direct product of chirality and helicity 2-spinors : @xmath99 here , @xmath100 and @xmath101 are left- and right - handed 1-spinors of helicity @xmath92 , respectively .
furthermore , @xmath102 is the helicity 2-eigenspinor : @xmath103 where @xmath104 is a shorthand notation for the pauli matrices .
furthermore note that @xmath105 .    for future convenience ,
let us explicitly derive the two vectors @xmath102 .
we write : @xmath106 and keep in mind that : @xmath107 it is a trivial exercise to solve for the two eigenvectors of @xmath108 , which correctly normalised to unity read :    [ helicity1 ] @xmath109    moreover , the helicity eigenstates are mutually orthogonal : @xmath110 finally , we have tacitly ignored the expansion of the anti - particle contribution @xmath111 in chirality and helicity spinors in equation ( [ creationannihilationexpansion ] ) .
we will return to this subtlety shortly .
our chirality and helicity decomposition of the spinor degrees of freedom in equation ( [ fermiondecompose ] ) also works in @xmath12 spacetime dimensions .
once a particular representation of the gamma matrices has been found , the projection operator : @xmath112 splits the @xmath113-spinor degrees of freedom in two equal contributions of definite chirality , i.e. : a left- and right - handed @xmath114-spinor ( see e.g. @xcite ) . here : @xmath115 where @xmath116 is fixed by requiring that @xmath117 is a proper projection operator : @xmath118 this yields : @xmath119 \,.\ ] ] in @xmath45 , we recognise the familiar @xmath120 .
an helicity operator that commutes with the projection operator can be defined in fourier space as : @xmath121 where @xmath122 is the levi - civita tensor in @xmath123 dimensions .
the exponential phase ensures that @xmath108 is hermitian .
finally , the gamma function has been inserted to account for the number of different permutations of gamma matrices .
hence , @xmath124 and @xmath125 are .
note the definition of the helicity operator is consistent with ( [ helicity ] ) .
most importantly , the chirality and helicity decomposition as proposed in equation ( [ fermiondecompose ] ) carries through for any dimension .    when an explicit form of @xmath102 is required , we return to @xmath45 and use equation ( [ helicity1 ] ) .
we subsequently generalise the result of such calculations by analytically extending it to arbitrary @xmath12 .      in textbooks on quantum field theory , such as @xcite ,
one expands the fermion in momentum space in spin eigenstates , whereas we expand in helicity eigenstates .
this results in a different normalisation requirement on the spinors in momentum space and consequently the standard textbook results can not be straightforwardly copied .
the reason for expanding in helicity eigenstates is that in curved spacetimes helicity is more convenient to work with rather than spin .
spin is a conserved quantity , i.e. : an appropriate quantum number , in static backgrounds , whereas helicity is conserved in time dependent but spatially homogeneous backgrounds @xcite .
one can easily show that the helicity operator constructed above in ( [ helicityd ] ) commutes with the kinetic operator in wigner space for spatially homogeneous correlators .
we normalise the spinors using two conditions : a consistent canonical quantisation and orthogonality of particle and anti - particle spinors .
we thus impose the usual anti - commutation relations between creation and annihilation operators :    [ anticommutation3 ] @xmath126    and all other anti - commutators vanish identically .
now ( [ anticommutation3 ] ) is consistent with ( [ anticommutation2 ] ) if and only if : @xmath127 note that this is a condition on @xmath128 matrix elements in spinor space , i.e. : for each @xmath63 and @xmath64 the above equality has to be satisfied .    secondly , we require that the particle and anti - particle spinors are mutually orthogonal : @xmath129 this restricts 4 matrix elements in helicity space , i.e. : for each @xmath130 the above equality needs to be obeyed .    peskin and schroeder @xcite
impose normalisation conditions on both particle and anti - particle spinors in spin space .
the corresponding conditions in helicity space are only partially satisfied .
both particle and anti - particle spinors of different helicity are trivially orthogonal by construction ( [ kvectorsorthogonal ] ) .
however , we do not relate the amplitude of equal helicity states for particles and anti - particles to the mass of the spinor because helicity eigenstates for a particle at rest can not be constructed .
we will now solve the dirac equation .
let us insert equation ( [ creationannihilationexpansion ] ) into the dirac equation ( [ fermioneom2 ] ) .
making use of our chirality and helicity decomposition ( [ fermiondecompose ] ) transforms the dirac equation to a first order coupled system of differential equations :    [ diraceq1 ] @xmath131    the arguments of the functions in the equation above have been omitted for notational convenience .
the left- and right - handed spinors however are now a function of the magnitude of the fourier mode @xmath132 only , because the differential equation is invariant under @xmath133 .
standard model fermions are chiral , in the sense that @xmath134 . for simplicity
, we take one mass @xmath135 only .
the following derivation can easily be generalised when standard model fermions are considered by making use of the @xmath117 projectors in equation ( [ p+- ] ) .
we now return to the expansion in chirality and helicity spinors of the anti - particle contribution @xmath111 , which we did not need so far .
since the dirac equation is a linear differential equation , the solution for @xmath111 can not contain any new degrees of freedom once we completely solved for @xmath98 .
in other words : the form of @xmath111 is dictated once we have completely specified @xmath98 .
we expand @xmath111 slightly differently : @xmath136 note we flipped the position of the left- and right - handed spinors compared to ( [ fermiondecompose ] ) .
consequently , the resulting equations of motion for both @xmath137 and @xmath138 are identical to ( [ diraceq1 ] ) .
alternatively , we could have expanded differently and flipped helicity , i.e. : send @xmath139 in equation ( [ fermiondecompose ] ) . because the mass enforces mixing between different chirality states , labelling these states differently has no physical relevance .
we return to equation of motion ( [ diraceq1 ] ) and define : @xmath140 this transforms equation ( [ diraceq1 ] ) to :    [ diraceq2 ] @xmath141    we exploit assumption ( [ epsilon ] ) , which is equivalent to ( [ epsilonconstant ] ) . inserting this expression into the equation above yields :    [ diraceq3 ] @xmath142    recall that the mass of a fermion in yukawa theory is time dependent in flrw spacetimes
. there does not exist a rescaling of the functions @xmath143 that removes the mass dependence in the equation of motion above .
even when one assumes that the mass @xmath144 is time independent , one can not find such a rescaling .
it is particularly simple to choose @xmath145 constant as following from equation ( [ phiproptoh ] ) and for convenience we define : @xmath146 we can proceed and decouple the two linear differential equations to find : @xmath147 the solutions are given by : @xmath148 where @xmath149 are two normalisation constants that still need to be determined and @xmath150 is the hankel function of the first kind of order : @xmath151 this implies :    @xmath152    the hankel function of the second kind also solves equation ( [ diraceq4 ] ) . in the infinite asymptotic past ,
equivalent to the deep uv , the fermions become effectively massless for @xmath153 and the distinction between @xmath154 and @xmath155 vanishes as a consequence .
therefore we require that the vacuum mode functions in this regime equal the standard conformal vacuum solutions :    [ uvlimit ] @xmath156    this excludes the hankel function of the second kind for the moment from contributing . away from the vacuum
, we should of course also incorporate the second solution to ( [ diraceq4 ] ) to allow for mode mixing .
note that the uv limits above yield the familiar uv behaviour for the left- and right - handed spinors ( e.g. : @xcite ) :    [ uvlimit2 ] @xmath157    when solutions ( [ diraceq5sols ] ) are substituted in ( [ diraceq3a ] ) one finds : @xmath158 we are free to fix the remaining coefficient @xmath159 by the following condition : @xmath160 we postpone the motivation for this requirement to section [ normalisation of the mode functions ] . for the moment ,
let us simply normalise accordingly for convenience .
we thus find : @xmath161 where we have chosen the unobservable phase factor in accordance with the uv limit ( [ uvlimit ] ) .
we use some well - known identities involving hankel functions which we list for convenience in appendix [ properties of hankel functions ] .
concluding , we find that the vacuum solutions to equation ( [ diraceq3 ] ) neglecting derivatives of @xmath9 and treating the quotient @xmath145 as time independent are given by :    [ diraceq6sols ] @xmath162    as we send @xmath163 and approach de sitter space , our solution is in accordance with @xcite , as can be seen from ( [ order ] ) . in this solution we have tacitly assumed that @xmath164 such that @xmath165 , equivalent to an accelerated expansion .
if however @xmath166 or @xmath167 , which corresponds to a decelerating universe , one should simply use relations ( [ hankelprops8][hankelprops9 ] ) in appendix [ properties of hankel functions ] .
the case @xmath168 requires special attention , see e.g. @xcite for the analogous scalar propagator case .
let us now turn our attention to finding the solution for the second mode function corresponding to the excited state .
the procedure is completely analogous : we again let @xmath169 and match to the appropriate uv expansion involving the second solution @xmath170 .
we normalise the solution according to ( [ normalisation2 ] ) and choose the phase factor in accordance with the corresponding uv limit
. we can thus construct the total solution by linear superposition :    [ diraceqtotalsol1 ] @xmath171    where :    [ diraceqtotalsol2 ] @xmath172    note that we have omitted all @xmath132 dependence in ( [ diraceqtotalsol1 ] ) for notational convenience .
note finally that the two fundamental solutions ( [ diraceqtotalsol2 ] ) in @xmath173 are not complex conjugates of each other , unlike for example for the scalar field case .
however , the second solution @xmath174 can be obtained from complex conjugation of the fundamental solutions comprising @xmath175 , as can clearly be seen in equation ( [ diraceqtotalsol1 ] ) above .    for the moment
, we leave the normalisation constants @xmath176 and @xmath177 undetermined .
appreciate we are nevertheless free to normalise the fundamental solutions to unity .
of course we will determine the normalisation constants shortly by applying our analysis performed in subsection [ spinorial normalisation conditions ] . but
first , we need to solve for the anti - particle mode functions .
let us stress again that the other spinor @xmath178 can not contain any new degrees of freedom because dirac s equation is linear .
we only have to determine the connection with @xmath97 . when we expand according to ( [ fermiondecompose2 ] ) the coupled linear differential equations @xmath137 and @xmath138 obey are identical to ( [ diraceq1 ] ) .
therefore , we can analogously define : @xmath179 we are rewarded for our extensive discussion of chirality and helicity .
the solutions of @xmath180 are identical to ( [ diraceqtotalsol1 ] ) , where we should only denote the coefficients in this solution differently :    [ diraceqtotalsol5 ] @xmath181    the functions appearing in this equation are given by ( [ diraceqtotalsol2 ] ) . we need to determine how the coefficients @xmath176 and @xmath177 of @xmath182 relate to @xmath183 and @xmath184 of @xmath180 .
moreover , we need to find the correct normalisation condition for @xmath176 and @xmath177 .      inspired by our discussion in [ spinorial normalisation conditions ] , we impose two normalisation conditions on the spinors : a consistent canonical quantisation ( [ consistentquantisation ] ) and orthogonality between particle and anti - particle spinors ( [ orthogonality ] ) .    after quite some work the , naively , 10 constraint equations in @xmath45 in ( [ consistentquantisation ] ) yield the following 3 conditions on our coefficients :    [ consistentquantisation2 ] @xmath185 the orthogonality condition ( [ orthogonality ] ) yields one more condition : @xmath186    the solution is :    [ consistentquantisation3 ] @xmath187 the last line relates the phases of our normalisation constants , e.g. : @xmath188 denotes the phase of @xmath183 . moreover , from ( [ consistentquantisation2 ] ) we can derive : @xmath189 ) results in a change of the right - hand side of equations ( [ consistentquantisation2a ] ) , ( [ consistentquantisation3d ] ) and ( [ consistentquantisation3e ] ) . normalising the fundamental solutions to unity as in ( [ normalisation2 ] )
is particularly convenient as it allows for a particle interpretation of @xmath190 and @xmath191 .
the normalisation constants provide us with three physical degrees of freedom .
of course @xmath192 is physical which determines @xmath193 through ( [ consistentquantisation3d ] ) .
hence , @xmath194 is physical just as in the scalar field case . since @xmath195 is a phase that can be removed without physical consequences , the phase relation ( [ consistentquantisation3c ] ) determines @xmath188 in terms of the third physical phase @xmath196 .
note finally that one of the two phase relations ( [ consistentquantisation2c ] ) and ( [ consistentquantisation2d ] ) is redundant in this part of the analysis .
a final remark is in order . to derive ( [ consistentquantisation2 ] ) we needed the explicit form of @xmath102 in equation ( [ helicity1 ] ) , which is only valid in @xmath45 .
however , we analytically continue this result to arbitrary dimension and one naturally finds ( [ consistentquantisation2 ] ) .      using the solutions for @xmath197 and @xmath198 given in ( [ diraceqtotalsol1 ] ) and ( [ diraceqtotalsol5 ] ) respectively , we recall relations ( [ uplusminus ] ) and ( [ nuplusminus ] ) to find :    [ summary1 ] @xmath199 \phantom{1111 } \label{summary1a } \\
\chi_{\mathrm{\scriptscriptstyle{r}},h}(k,\eta ) & = & \frac{1}{\sqrt{2}}\left [ c\ { \mathcal{h } + h \tilde{\mathcal{h}}^{\ast } \ } + d\ { \tilde{\mathcal{h } } - h \mathcal{h}^{\ast } \ } \right ] \label{summary1b } , \phantom{1111}\end{aligned}\ ] ] and moreover : @xmath200\phantom{1111 } \label{summary1c } \\
\nu_{\mathrm{\scriptscriptstyle{r}},h}(k,\eta ) & = & \frac{1}{\sqrt{2}}\left [ f\ { \mathcal{h } + h \tilde{\mathcal{h}}^{\ast } \ } + g\ { \tilde{\mathcal{h } } - h \mathcal{h}^{\ast } \}\right ] , \phantom{1111}\label{summary1d}\end{aligned}\ ] ]    where the fundamental solution are given in equation ( [ diraceqtotalsol2 ] ) .
the normalisation conditions ( [ consistentquantisation3 ] ) determine @xmath183 and @xmath184 in terms of @xmath176 and @xmath177 .
moreover , we have @xmath201
.      having discussed the fermionic mode functions in some detail , we can turn our attention to solving the feynman propagator for fermions .
we insert the rescaled spinors , expanded in terms of creation and annihilation operators , into the formal definition of the feynman propagator ( [ feynmanpropagator ] ) , which yields : @xmath202 where @xmath203 and @xmath80 .
the reader can easily see that we have to make use of the solutions ( [ summary1 ] ) above to calculate the 16 matrix elements of the propagator .
let us outline this rather cumbersome calculation in seven steps .
we recall the form of e.g. the massless fermionic propagator ( [ fermionpropmassless2 ] ) and realise we have to extract an operator acting on @xmath68 , rather than on @xmath69 .
the reader can easily see that the 16 matrix elements of the propagator contain products of two hankel functions .
some of these products involve hankel functions at the same order , but others encompass hankel functions where the order differs . since we can only easily integrate products of hankel functions of equal order we need to transform precisely the latter products by making use of identities ( [ hankelprops ] ) , e.g. : @xmath204 where @xmath205 or @xmath206 applies when @xmath144 equals 1 or 2 , respectively . clearly , the `` extracted operator '' is a function of @xmath56 .
we rewrite the hankel functions in terms of macdonald functions . by making use of @xcite
we find :    [ hankelmacdonald ] @xmath207    where @xmath208 is a macdonald function .
the emerging matrix structure of equation ( [ feynmanpropagator2 ] ) is most easily uncovered by evaluating its @xmath32-matrix constituents successively and writing the result in terms of pauli matrices .
recalling equations ( [ kvector1 ] ) and ( [ kvector2 ] ) we derive : @xmath209 where the @xmath206 or @xmath205 sign applies depending on whether one deals with the @xmath210 or @xmath211 contribution in ( [ feynmanpropagator2 ] ) , respectively .
the last step is possible because only the exponent depends on @xmath212 .
note finally that the index @xmath213 in equation ( [ convertingderivatives ] ) is supposed to run over ( flat ) spatial indices exclusively .
the @xmath214 independent contribution to the propagator can thus be pulled out of the fourier integral .
these matrices simplify considerably when we consider the following identities : @xmath215 clearly , we recognise the complex conjugate of the dirac operator and the @xmath216 structure we are familiar with from the de sitter propagator .
let us present an intermediate result that already captures much of the final structure of our propagator . for brevity ,
let us only consider the contribution to the propagator arising from the first @xmath66-function in the propagator , i.e. : when @xmath217 .
this part of the feynman propagator reads : @xmath218 \nonumber , \end{aligned}\ ] ] where :    @xmath219    and where : @xmath220    the contribution to the propagator when @xmath221 is given by the same expression , where we only need to replace the normalisation constants @xmath176 and @xmath177 by @xmath183 and @xmath184 , respectively .
the complicated structure of the propagator in spinor space has dramatically simplified by extracting the appropriate operators in position space , which could already be expected from the propagator in de sitter spacetime .
we need to perform the fourier integrals over the functions @xmath222 and @xmath223 next .
we use : @xmath224 which is valid for any function @xmath225 that depends solely on @xmath226 . here , @xmath227 is a bessel function of the first kind .
the reader can easily verify that all the fourier integrals we need to perform are of the following form : @xmath228 where the arguments of the macdonald functions are purely imaginary :    [ fourierintegral2 ] @xmath229    we perform these integrals in appendix [ fourier transforming hankel functions ] by making use of @xcite .
the result of this integral can be expressed in terms of the gauss hypergeometric function @xmath230 .
the @xmath66-functions in the propagator in combination with the arguments of the hypergeometric functions allow us to rewrite these arguments in terms of the @xmath231-functions ( [ y ] ) we introduced in section [ properties of flrw spacetimes ] .    a well - known fact in minkowski or de
sitter spacetime extends to our analysis in flrw spacetimes : the small real contribution needed to make the fourier integrals convergent determines the @xmath232 pole prescription of the various contributions to the propagator .
the final result for the constant @xmath9 and constant @xmath233 flrw fermion propagator reads : @xmath234 \,,\ ] ] where : @xmath235 where the @xmath231-functions are given in equation ( [ y ] ) and @xmath233 in ( [ zeta ] ) .
moreover , to write the propagator in terms of the @xmath231-functions we introduced the following notation : @xmath236    one can think of @xmath237 as an antipodal coordinate ( see e.g. : @xcite ) . note that we extract an operator in ( [ fermionpropagator ] ) that does not depend on @xmath69 .
however , there is nothing special about the @xmath68 leg of the propagator . keeping remark 1 above in mind
, we know that we could have equally well written the propagator ( [ fermionpropagator ] ) in terms of an operator acting on @xmath69 by replacing : @xmath238 this propagator would then have satisfied the dirac equation ( [ fermionprop ] ) in terms of @xmath69 .
as a first check of this result , note that if we send @xmath239 and @xmath240 , the massive fermionic flrw propagator ( [ fermionpropagator ] ) correctly yields the massless fermionic flrw propagator earlier derived in equation ( [ fermionpropmassless2 ] ) .
again , in order to derive this result we needed the explicit form of @xmath102 given by ( [ helicity1 ] ) . the operator in ( [ fermionpropagator ] ) has thus been constructed in @xmath45 .
however , we can analytically continue again to arbitrary dimensions .
finally , let us study the ir behaviour of our propagator .
note that the coincidence propagator ( in fourier space ) is ir finite due to the pauli exclusion principle . in position space ,
note that the ir limit of all the y - functions in equation ( [ y ] ) is : @xmath241 where the plus or minus sign applies for infinite spacelike or timelike separation , respectively .
if we recall the following well - known relation :    @xmath242    we can easily extract the leading order ir behaviour of the propagator in equation ( [ fermionpropagator2 ] ) : @xmath243    note that for the spacelike ir , all y - functions behave identically and the result above can be simplified further . for the timelike ir ,
the @xmath232-prescription in equation ( [ y ] ) dictates how to take the powers of the y - functions .
of course , one has to insert equation ( [ fermionpropagator2ir ] ) above into ( [ fermionpropagator ] ) to obtain the full propagator in the ir limit .
the fermionic propagator ( [ fermionpropagator ] ) is analogous to the scalar and graviton propagator in universes with constant deceleration parameter @xmath0 derived in @xcite . unlike in the case of scalars and gravitons , where quite generically the bunch - davies vacuum is ir divergent , the fermionic propagator is ir finite due to the pauli exclusion principle which forbids an accumulation of fermions in the deep ir .
the ir behaviour for the scalar and graviton are important , as it can comprise secular effects ( that grow as a power in time ) in for example @xmath244 ( see : @xcite ) and @xmath245 ( see : @xcite ) . mathematically speaking ,
pauli blocking is enforced by the imaginary @xmath246 in equation ( [ diraceq4 ] ) , which ultimately carries through in the complex indices of the hypergeometric functions in the propagator ( [ fermionpropagator2 ] ) .
having derived the fermionic flrw propagator , let us verify whether it correctly reduces to the fermionic de sitter propagator , which is a well - known result in the literature .
the fermion propagator in de sitter spacetime reads @xcite : @xmath247 \nonumber,\end{aligned}\ ] ] where the superscript @xmath248 denotes de sitter spacetime and where : @xmath249 here , @xmath231 is shorthand for @xmath250 .
if we send @xmath163 and set @xmath251 and @xmath252 , equation ( [ fermionpropagator ] ) above correctly reproduces the de sitter result .
moreover , using : @xmath253 note that the massless limit of ( [ fermionpropmassdesitter ] ) indeed corresponds to ( [ fermionpropmassless2 ] ) in de sitter space .
we conclude that the fermionic propagator in flrw spacetimes correctly reduces to the known cases in existing literature .
as a simple application for the propagator derived in section [ the feynman propagator ] , we calculate the one loop contribution to the effective action .
this is important to study the impact of fermions on the evolution of the background spacetime and of scalar fields coupled to these fermions . we renormalise using the minimal subtraction dimensional renormalisation technique .
the one loop effective action formally reads : @xmath254 \\ & = &   \int^{m}\ !
\mathrm{d}\bar{m } \ , \mathrm{tr } \left [ \sqrt{-g } \ , i s_{\mathrm{f}}^{ab}(x,\tilde{x})\right ] \nonumber\,.\end{aligned}\ ] ] here ,
the trace is both over spatial and spinorial indices .
we have been able to evaluate the logarithm at the expense of losing all mass independent contributions to the effective action .    the effective action for a massless fermion is completely fixed by the trace anomaly which this field is known to exhibit ( see for example @xcite ) .
the trace anomaly can be generated from a finite non - local effective action ( see e.g. : @xcite ) .
it can alternatively be generated from an infinite but local effective action @xcite .
we will henceforth neglect all contributions from the trace anomaly to the effective action in this paper because it has already been extensively discussed in the literature .
we treat the covariant derivative according to ( [ covariantderivative2 ] ) and we can easily expand the hypergeometric functions in equation ( [ fermionpropagator2 ] ) at coincidence . in dimensional regularisation , the @xmath12 dependent powers in this expansion do not contribute .
these are generated for the @xmath255 and @xmath256 contributions to ( [ fermionpropagator2 ] ) by the first term in : @xmath257 moreover , when we trace over spacetime indices note there is no contribution at all from the derivative term hitting one of the hypergeometric functions : @xmath258 here , @xmath259 is the integer valued coefficient of the taylor expansion of a hypergeometric function .
a similar identity holds for @xmath256 .
likewise , for @xmath260 or @xmath261 , we have e.g. : @xmath262 note that both @xmath260 and @xmath261 depend on the antipodal coordinate @xmath237 in equation ( [ fermionpropagator ] ) . at coincidence
, we also use : @xmath263 the result is :    @xmath264    let us compare this result with the known calculations in the literature in de sitter spacetime ( @xmath265 ) and in the vacuum ( @xmath266 ) .
note the two omissions in the candelas and raine effective action @xcite .
when tracing over spinor indices , we have : @xmath267 and this trace does not , as is apparent from @xcite , equal @xmath268 because a spinor in @xmath12 spacetime dimensions has @xmath269 degrees of freedom ( see e.g. : @xcite ) .
indeed , this has already been noted by @xcite .
moreover , the effective action in @xcite misses a factor of 2 stemming from the separate particle and anti - particle contributions , which has been corrected for in e.g. @xcite .
the reader can easily verify that ( [ effectiveactionfullresult ] ) , taking these errors into account agrees with @xcite in de sitter spacetime and in the vacuum .      in the spirit of dimensional regularisation ,
we now expand around @xmath45 .
we can use the familiar result : @xmath270 + \mathcal{o}(\tilde{\epsilon}^{2})\,,\ ] ] where @xmath271 and where @xmath272 is the digamma function defined by @xmath273 . anticipating the form of the counterterms we will add shortly ,
we introduce a scale @xmath274 by : @xmath275 \nonumber , \end{aligned}\ ] ] we we neglected @xmath276 contributions .
finally , we recall the familiar result : @xmath277 where @xmath278 is the euler  mascheroni constant .
we thus arrive at : @xmath279 \right ) \label{effectiveactionexpansiond4 } \\ & & \ ! + \
! \int^{m}\!\!\mathrm{d}\tilde{m}\ !
\left(\tilde{m}^{3}\!+\!\tilde{m}h^{2}(1\!-\!\epsilon)^{2 } \right)\!\left [ \psi(2\!-\!i\tilde{\zeta})\!+\!\psi(2\!+\!i\tilde{\zeta})\ ! \right]\ ! \bigg]\ ! \nonumber \\ & & \
! + 2 h^{3}(1-\epsilon)^{3 } \mathrm{im}\left(c d^{\ast}\right ) \ !
\mathrm{d}\tilde{m } \ , \gamma(2+i \tilde{\zeta } ) \gamma(2-i\tilde{\zeta})\bigg\ } \nonumber \!.\end{aligned}\ ] ] here , @xmath280 denotes the imaginary part of @xmath281 .
the one loop backreaction contains a divergence when @xmath45 .
however , appreciate that the contribution to the effective action multiplying the mixed coefficients is finite in @xmath45 .
anticipating the renormalisation procedure in the next section , we need the tree level , i.e. : classical , friedmann equations of motion responsible for driving the expansion of the universe .
the system of interest is thus given by : @xmath282 \right\ } \ , .
\nonumber\end{aligned}\ ] ] here , @xmath283 represents the rescaled newton constant and @xmath248 denotes the cosmological constant .
the reader can easily recognise the usual einstein - hilbert action on the first line .
we introduce a new scalar field @xmath284 , responsible for the dynamics of the universe . by an appropriate choice of the potential @xmath285
, a scalar field can mimic any mixture of fluids relevant for the evolution of our universe ( see e.g. @xcite ) .
the last two lines in equation ( [ actiontotalgravity ] ) above are identical to ( [ actiontotal ] ) and contain the fermion field and the scalar field @xmath286 generating the mass of the fermion through a yukawa coupling , as in ( [ fermionmass ] ) .
we write :    [ backgroundfields1 ] @xmath287    where we assume the background fields @xmath288 and @xmath289 to be homogeneous .
moreover , we assume that , such that @xmath290 drives the dynamics of the universe as stated above .
the classical friedmann equations of motion are :    [ treelevelfriedmann1 ] @xmath291 and moreover , the classical scalar field equation of motion reads : @xmath292    we can thus derive the following identities :    [ treelevelfriedmann2 ] @xmath293    here , we have used assumption ( [ epsilon ] ) .      in order to derive equation ( [ effectiveactionexpansiond4 ] ) ,
we have assumed a constant deceleration @xmath9 , as in equation ( [ epsilon ] ) .
we now promote this constant to a dynamical quantity : @xmath294 we will motivate this step shortly . the effective action ( [ effectiveactionexpansiond4 ] ) now contains divergences in both @xmath295 and @xmath296 that we ought to cancel by an appropriate counterterm action .
it is not possible to identify local covariant counterterms of curvature invariants only that remove the singularities in both of these quantities at the level of the effective action .
we can use the tree level equations of motion derived above in ( [ treelevelfriedmann2 ] ) to renormalise the effective action , which is first order in @xmath96 . by making use of an infinitesimal field redefinition
, one can show the following : substituting the zeroth order equation of motion into the first order contribution to the action yields , up to one loop , the same equation of motion as directly varying the action .
we thus add the following counterterm action : @xmath297 the form of the counterterm action is unique , which can be shown by making use of dimensional analysis and by requiring that only two time derivatives can act on @xmath288 or @xmath298 for stability .
these conditions limit us to the counterterms above and possibly @xmath299 .
however , one can easily verify that the latter term does not have the correct form required to cancel the divergences in ( [ effectiveactionexpansiond4 ] ) .
we recall relation ( [ fermionmass ] ) .
we need to expand the various terms in the counterterm action above around @xmath45 , by making use of @xmath300 and equation ( [ treelevelfriedmann2 ] ) :    [ counterterms1 ] @xmath301    where we neglected all @xmath302 contributions in the equations above .
we can now easily solve for the coefficients in the counterterm action :    [ counterterms2 ] @xmath303    where the divergent coefficients in @xmath45 are fixed to cancel the divergences occurring in ( [ effectiveactionexpansiond4 ] ) , and where e.g. @xmath304 is a finite but arbitrary coefficient of the counterterm action .
the renormalised one loop effective action thus reads :    @xmath305 \right ) - \frac{1}{12}m^{2}h^{2}(7 - 5\epsilon )   \nonumber \\ & & \qquad\qquad\qquad\qquad\qquad\quad   + \int^{m}\!\mathrm{d}\tilde{m } \left(\tilde{m}h^{2}(1-\epsilon)^{2 } + \tilde{m}^{3 } \right ) \left [ \psi(2-i\tilde{\zeta})+\psi(2+i\tilde{\zeta } ) \right ] \bigg ]   \label{effectiveactionrenormalised } \\ & & \qquad\qquad\quad + 2 h^{3}(1-\epsilon)^{3 } \mathrm{im } \left ( c d^{\ast } \right ) \int^{m}\ !
\mathrm{d}\tilde{m } \
, \gamma(2+i \tilde{\zeta } ) \gamma(2-i\tilde{\zeta})\bigg\ } \nonumber \\ & & + \int \mathrm{d}^{4}\!x \sqrt{-g } \big(c_{1}^{\mathrm{f } } \phi^{2 } r + c_{2}^{\mathrm{f } } \phi^{2}\frac{\partial^{2}v } { \partial\varphi^{2}}(\varphi_{0 } ) - c_{3}^{\mathrm{f } } \kappa \phi^{2 } g^{\mu\nu}\partial_{\mu}\varphi \partial_{\nu}\varphi + c_{4}^{\mathrm{f } } \phi^{4 } \big ) \nonumber\,.\end{aligned}\ ] ]
the effective action above comprises the one loop backreaction of fermions .
this term gives rise to additional physical corrections to the ( classical ) friedmann equations , thus obtaining the so - called quantum corrected friedmann equations .
we will study its effect on the evolution of our universe and its impact on the evolution of the scalar field to which the fermions are coupled in a future publication . in this equation
we have neglected all contributions to the one loop effective action arising from the trace anomaly .
if we expand equation ( [ effectiveactionrenormalised ] ) to first order in @xmath9 , we can find out how large the first order correction to the de sitter result is .
we have :    @xmath306\right\ } - \frac{1}{2}m^{4 } \label{effectiveactionrenormalisedfirstorder } \\ & & \quad + \int^{m}\!\mathrm{d}\tilde{m } i \left(\tilde{m}^{2}h + \frac{\tilde{m}^{4}}{h } \right ) \big [ \psi'(2+i\tilde{m}/h)-\psi'(2-i\tilde{m}/h ) \big]-2\tilde{m}h^{2}\big [ \psi(2+i\tilde{m}/h)+\psi(2-i\tilde{m}/h ) \big ] \bigg ]    \nonumber \\ & & \quad + 2 h^{3 } \mathrm{im } \left ( c d^{\ast } \right ) \int^{m}\ ! \mathrm{d}\tilde{m } \
, \gamma(2+i \tilde{m}/h ) \gamma(2-i\tilde{m}/h ) \left[i\frac{\tilde{m}}{h}\big ( \psi(2+i\tilde{m}/h)-\psi(2-i\tilde{m}/h ) \big)-3 \right ] \bigg\ } \nonumber + \mathcal{o}\left(\epsilon^{2}\right),\end{aligned}\ ] ]    where the prime denotes a derivative with respect to the argument , and the superscript @xmath248 represent the well - known de sitter result , which trivially follows from ( [ effectiveactionrenormalised ] ) by letting @xmath265 .
let us at this point however already consider the @xmath307 limit , to extract the physical behaviour of the effective potential resulting from ( [ effectiveactionrenormalised ] ) .
note firstly that : @xmath308 where @xmath309 denotes the one loop effective potential for @xmath286 . in order to evaluate the integrals over @xmath310 occurring in equation ( [ effectiveactionrenormalised ] )
, we need the following asymptotic expansions ( @xmath311 ) : @xmath312 and also : @xmath313 & = & \ln [ \zeta^{2 } ] + \frac{13}{6}\frac{1}{\zeta^{2 } } \\ & & - \frac{119}{60}\frac{1}{\zeta^{4 } } + \mathcal{o}\left(\frac{1}{\zeta^{6}}\right ) \nonumber \,.\end{aligned}\ ] ] keeping the leading order terms of the effective potential ( [ effectivepotential1 ] ) in this limit only , we find from ( [ effectiveactionrenormalised ] ) and ( [ limitbigmass1][limitbigmass2 ] ) : @xmath314 \\ & &   \!+   \frac{\mathrm{im}(c d^{\ast})}{\pi^{2 } }   g_{\mathrm{\scriptscriptstyle{y}}}^{3 } \phi^{3 } h(1\!-\!\epsilon ) \exp\left[{-\pi \frac{g_{\mathrm{\scriptscriptstyle{y } } }    equation ( [ fermionmass ] ) has been used . without studying the dynamics resulting from this asymptotic effective potential , let us briefly consider its stability properties .
clearly , this potential is unstable in the vacuum ( @xmath266 ) because the effective potential is unbounded from below in this case .
this wild behaviour has already been recognised in de sitter spacetime @xcite .
an interesting phenomenon can be observed precisely when @xmath315 .
this is realised for example in a thermal state with a temperature @xmath316 much larger than the uv cutoff of the theory . in this case , the asymptotic effective potential is exponentially suppressed and its ( small ) contributions can either be positive or negative .
let us make a few general remarks .
consider quantum field theories in de sitter spacetime .
the hubble parameter @xmath4 in de sitter spacetime is a constant .
an effective action for fermions or scalars in de sitter spacetime contains divergences for example proportional to a power of @xmath4 ( see e.g. : @xcite ) .
the same statement holds for divergences in the stress energy tensor in for example the trace anomaly literature ( see e.g. : @xcite ) . at this point
, one can interpret the results differently :    1 .
the hubble parameter can quite conservatively be treated as a _
constant_. the argument for this case is simple : throughout the calculation one simply assumed a constant hubble parameter .
the hubble parameter is promoted to a _
dynamical _ quantity : @xmath317 a reason for supporting this case would be the following : when @xmath8 one would find the de sitter result as a leading order contribution in each time interval where the hubble parameter varies only adiabatically slowly , a spacetime also known as locally de sitter spacetime .
depending on which of the two interpretations one adheres to , one would renormalise differently . if one advocates the first point of view , all terms of e.g. an effective action merely contribute to the renormalised cosmological constant . if , however , one prefers the second interpretation , one adds ( dynamical ) curvature invariants in order to cancel these divergences in a covariant manner .
we observe the following . in the effective action ( [ effectiveactionrenormalised ] ) above , the hubble parameter is a dynamical quantity . hence , covariant counterterms added to cancel uv divergences can and should be treated dynamically .
this allows us to study the backreaction of fermions in a dynamical manner .
albeit strictly speaking not allowed to relax the assumption @xmath318 in de sitter spacetime , we conclude that promoting the hubble parameter to a dynamical quantity @xmath319 has been a useful approach .
when in our calculation we keep @xmath9 constant as assumed , result ( [ effectiveactionfullresult ] ) favours the second interpretation above .
let us now return to our calculation to see what the above considerations can tell us about our effective action ( [ effectiveactionfullresult ] ) .
we face the following two options :    1 .
we treat @xmath9 as a _
constant_. 2 .
we promote @xmath9 to a _
dynamical _ quantity , as already stated in ( [ epsilonpromotion ] ) . we could call this a spacetime with an almost constant deceleration .    despite the fact that throughout the calculation we have assumed that @xmath9 is a constant
, it will come as no surprise that we argue in favour of the latter option for two reasons : if one would allow @xmath9 to change , one would probably find our main result for the effective action before renormalisation ( [ effectiveactionfullresult ] ) in each time interval where @xmath9 changes only very slowly .
secondly , the latter option is the generalisation of the confirmed de sitter case above
.    an obvious disadvantage of our choice is that we can only renormalise on - shell .
consider the divergences in ( [ effectiveactionexpansiond4 ] ) .
as already mentioned , it is not possible to identify local covariant counterterms of curvature invariants only that remove all singularities in @xmath295 and @xmath296 at the level of the effective action .
therefore , we relied on an on - shell renormalisation technique and inserted the zeroth order ( tree level ) friedmann equations of motion into the effective action at first order in @xmath96 .    finally , let us try to comfort the reader who would rather keep @xmath9 constant , but who would simultaneously like to understand how to improve on his final result for the renormalised effective action . in this case , the only appropriate counterterms to renormalise ( [ effectiveactionexpansiond4 ] ) available to us are the following : @xmath320 we can then renormalise the theory as usual would have to be considered separately , for @xmath321 in a radiation dominated universe .
] , with :    [ countertermscoefficientsalt ] @xmath322    we could now study mode mixing . in their paper , tsamis and woodard study mode mixing for a massless minimally coupled scalar field in constant @xmath9 flrw spacetimes .
the basic strategy is as follows : in each  possibly infinitesimally small  time interval @xmath9 can be well approximated by a constant .
the mode functions in two adjacent time intervals can be related by imposing two conditions : ( i ) continuity of the mode functions and ( ii ) continuity of their first derivatives .
hence , the mode functions at some later time @xmath56 can be obtained from the initial condition at @xmath323 by means of a transfer matrix . ] .
one could simply allow @xmath9 and @xmath145 to vary in time and choose ( if necessary infinitesimally small ) time intervals in which @xmath9 and @xmath145 vary only slowly .
subsequently , one can match the solutions in two neighbouring intervals and study mode mixing in the spirit of @xcite . in this case
, this procedure would imply that we would have to renormalise for each time interval separately , for the coefficients in equation ( [ countertermscoefficientsalt ] ) depend on @xmath9 .
this is not surprising : a sudden jump , even infinitesimally small , in for example @xmath145 is not a physical phenomenon and would call for non - local counterterms ( see e.g. : @xcite ) .
in this paper , we have constructed the fermionic propagator in realistic flrw spacetimes with constant deceleration parameter @xmath0 . moreover
, we assumed that @xmath145 is also constant , which can be realised in yukawa theory by means of a scalar field for which @xmath11 .
we have derived the propagator both in the ( trivial ) massless case ( [ fermionpropmassless2 ] ) and in the massive case ( [ fermionpropagator ] ) .
two pillars support our derivation .
we split the fermionic degrees of freedom into a direct product of chirality and helicity eigenspinors .
we normalise the spinors using a consistent canonical quantisation and we require orthogonality of particle and anti - particle spinors .
moreover , we verify that upon sending @xmath9 to zero , we recover the known vacuum - to - vacuum de sitter results in the literature .
we calculate the one loop effective action induced by fermions using our propagator .
the one loop backreaction arises from integrating out a free , quadratic fermion field and , using dimensional regularisation , this generates a correction to the ( classical ) friedmann equations .
this effective action allows for the first time for a dynamical interpretation of @xmath4 such that the study of dynamical backreaction either on the background spacetime or on scalar fields coupled to these fermions is within reach .
* acknowledgements *    jfk thanks gerben stavenga , shun - pei miao and timothy budd for useful suggestions .
the authors thank gerard t hooft , gleb arutyunov and tomas janssen and acknowledge financial support from fom grant 07pr2522 and by utrecht university .
in section [ particle mode functions ] we used the following well - known identities involving hankel functions :    [ hankelprops ] @xmath324 the wronskian of two hankel functions reads : @xmath325 & = & -\frac{4 i}{\pi z }   \label{hankelprops5}\,.\end{aligned}\ ] ] moreover , we made use of the following recurrence relation : @xmath326
consider the following integral equation : @xmath327 here , @xmath328 is a bessel function of the first kind and @xmath208 is a macdonald function .
note that for all integrals in ( [ fourierintegral ] ) both @xmath63 and @xmath64 are purely imaginary .
the requirements from @xcite are such that :          we choose @xmath331 to make the integral convergent .
we now use equation ( 6.578.10 ) from @xcite to find : @xmath332 where : @xmath333 here , @xmath334 represents the associated legendre function which we can rewrite in terms of a gauss hypergeometric function @xmath335 by means of equations ( 8.702 ) and ( 9.131.1 ) of @xcite : @xmath336 upon taking @xmath337 we find : @xmath338 \,.\ ] ] finally note that we can further simplify the argument of the hypergeometric function in equation ( [ i6 ] ) to : @xmath339 let us finally stress that we have introduced @xmath331 to make the integral convergent .
it is however of significant physical importance for it is the same @xmath340 that dictates the @xmath232 pole prescription in the propagator ( see equations ( [ fermionpropagator][fermionpropagator2 ] ) and ( [ y ] ) ) .
this should not come as a surprise : we observe the same behaviour in e.g. minkowski and de sitter spacetimes .
99 l.  h.  ford , quantum instability of de sitter space - time , phys .  rev .
d * 31 * ( 1985 ) 710 .
i.  antoniadis and e.  mottola , graviton fluctuations in de sitter space , j.  math .
phys .   * 32 * ( 1991 ) 1037 . n.  c.  tsamis and r.  p.  woodard , the structure of perturbative quantum gravity on a de sitter background , commun .
* 162 * ( 1994 ) 217 .
n.  c.  tsamis and r.  p.  woodard , quantum gravity slows inflation , nucl .
b * 474 * ( 1996 ) 235 [ arxiv : hep - ph/9602315 ] .
n.  c.  tsamis and r.  p.  woodard , one loop graviton self - energy in a locally de sitter background , phys .
d * 54 * ( 1996 ) 2621 [ arxiv : hep - ph/9602317 ] .
l.  r.  w.  abramo , r.  h.  brandenberger and v.  f.  mukhanov , the energy - momentum tensor for cosmological perturbations , phys .
d * 56 * ( 1997 ) 3248 [ arxiv : gr - qc/9704037 ] .
f.  finelli , g.  marozzi , g.  p.  vacca and g.  venturi , adiabatic regularization of the graviton stress - energy tensor in de sitter space - time , phys .
d * 71 * ( 2005 ) 023522 [ arxiv : gr - qc/0407101 ] .
t.  prokopec , a solution to the cosmological constant problem , arxiv : gr - qc/0603088 . t.  prokopec , n.  c.  tsamis and r.  p.  woodard , stochastic inflationary scalar electrodynamics , annals phys .
* 323 * , 1324 ( 2008 ) [ arxiv:0707.0847 [ gr - qc ] ] .
a.  bilandzic and t.  prokopec , quantum radiative corrections to slow - roll inflation , phys .
d * 76 * ( 2007 ) 103507 [ arxiv:0704.1905 [ astro - ph ] ] .
t.  janssen and t.  prokopec , a graviton propagator for inflation , class .
* 25 * ( 2008 ) 055007 [ arxiv:0707.3919 [ gr - qc ] ] .
t.  janssen , s.  p.  miao and t.  prokopec , graviton one - loop effective action and inflationary dynamics , arxiv:0807.0439 [ gr - qc ] .
t.  janssen and t.  prokopec , implications of the graviton one - loop effective action on the dynamics of the universe , arxiv:0807.0447 [ gr - qc ] .
t.  m.  janssen , s.  p.  miao , t.  prokopec and r.  p.  woodard , infrared propagator corrections for constant deceleration , class .
* 25 * ( 2008 ) 245013 [ arxiv:0808.2449 [ gr - qc ] ] .
p.  candelas and d.  j.  raine , general relativistic quantum field theory - an exactly soluble model , phys .
d * 12 * ( 1975 ) 965 .
b.  allen and c.  a.  lutken , spinor two point functions in maximally symmetric spaces , commun .
* 106 * ( 1986 ) 201 .
r.  camporesi , the spinor heat kernel in maximally symmetric spaces , commun .
* 148 * ( 1992 ) 283 .
miao and r.  p.  woodard , the fermion self - energy during inflation , class .
* 23 * ( 2006 ) 1721 [ arxiv : gr - qc/0511140 ] .
s.  p.  miao and r.  p.  woodard , leading log solution for inflationary yukawa , phys .
d * 74 * ( 2006 ) 044019 [ arxiv : gr - qc/0602110 ] .
e.  t.  tomboulis , dynamically adjusted cosmological constant and conformal anomalies , nucl .
b * 329 * ( 1990 ) 410 .
i.  antoniadis and e.  mottola , 4-d quantum gravity in the conformal sector , phys .
d * 45 * ( 1992 ) 2013 .
i.  antoniadis , dynamics of the conformal factor in 4-d gravity , [ arxiv : hep - th/9211055 ] .
i.  antoniadis , p.  o.  mazur and e.  mottola , fractal geometry of quantum spacetime at large scales , phys .
b * 444 * ( 1998 ) 284 [ arxiv : hep - th/9808070 ] .
h.  salehi and y.  bisabr , conformal anomaly and large scale gravitational coupling , int .
j.  theor .
* 39 * ( 2000 ) 1241 [ arxiv : hep - th/0001095 ] .
i.  antoniadis , p.  o.  mazur and e.  mottola , cosmological dark energy : prospects for a dynamical theory , new j.  phys .
* 9 * ( 2007 ) 11 [ arxiv : gr - qc/0612068 ] . j.  f.  koksma and t.  prokopec , effect of the trace anomaly on the cosmological constant , phys .
d * 78 * ( 2008 ) 023508 [ arxiv:0803.4000 [ gr - qc ] ] .
n.  a.  chernikov and e.  a.  tagirov , quantum theory of scalar fields in de sitter space - time , annales poincare phys .
a * 9 * ( 1968 ) 109 .
b.  allen , vacuum states in de sitter space , phys .
d * 32 * ( 1985 ) 3136 .
b.  allen and a.  folacci , the massless minimally coupled scalar field in de sitter space , phys
.  rev .
d * 35 * ( 1987 ) 3771 .
b.  allen and t.  jacobson , vector two point functions in maximally symmetric spaces , commun .
* 103 * ( 1986 ) 669 .
n.  c.  tsamis and r.  p.  woodard , a maximally symmetric vector propagator , j.  math .
* 48 * ( 2007 ) 052306 [ arxiv : gr - qc/0608069 ] .
b.  allen , the graviton propagator in de sitter space , phys .
d * 34 * ( 1986 ) 3670 .
b.  allen and m.  turyn , an evaluation of the graviton propagator in de sitter space , nucl .
b * 292 * ( 1987 ) 813 .
e.  g.  floratos , j.  iliopoulos and t.  n.  tomaras , tree level scattering amplitudes in de sitter space diverge , phys .
b * 197 * ( 1987 ) 373 .
a.  higuchi and s.  s.  kouris , the covariant graviton propagator in de sitter spacetime , class .
grav .   * 18 * ( 2001 ) 4317 [ arxiv : gr - qc/0107036 ] .
n.  d.  birrell and p.  c.  w.  davies , quantum fields in curved space , cambridge monographs on mathematical physics , cambridge university press ( 1982 ) .
t.  prokopec and e.  puchwein , photon mass generation during inflation : de sitter invariant case , jcap * 0404 * ( 2004 ) 007 [ arxiv : astro - ph/0312274 ] .
s.  weinberg , quantum contributions to cosmological correlations , phys .
d * 72 * ( 2005 ) 043514 [ arxiv : hep - th/0506236 ] .
j.  f.  koksma , t.  prokopec and g.  i.  rigopoulos , the scalar field kernel in cosmological spaces , class .
* 25 * ( 2008 ) 125009 [ arxiv:0712.3685 [ gr - qc ] ] .
i.  i.  cotaescu , polarized dirac fermions in de sitter spacetime , phys .  rev .
d * 65 * ( 2002 ) 084008 [ arxiv : hep - th/0109199 ] . b.  garbrecht and t.  prokopec , fermion mass generation in de sitter space , phys .
d * 73 * ( 2006 ) 064036 [ arxiv : gr - qc/0602011 ] .
b.  de wit and j.  smith , field theory in particle physics .
1 , north - holland personal library ( 1986 ) .
m.  e.  peskin and d.  v.  schroeder , an introduction to quantum field theory , addison - wesley ( 1995 ) . c.  itzykson and j.  b.  zuber , quantum field theory , mcgraw - hill ( 1980 ) .
t.  prokopec , m.  g.  schmidt and s.  weinstock , transport equations for chiral fermions to order h - bar and electroweak baryogenesis , annals phys .
* 314 * ( 2004 ) 208 [ arxiv : hep - ph/0312110 ] .
i.  s.  gradshteyn and i.  m.  ryzhik , table of integrals , series , and products , academic press ( 1965 ) .
a.  p.  prudnikov , yu .
a.  brychkov and o.  i.  marichev , integrals and series , volume 2 , special functions , gordon and breach science publishers ( 1986 ) .
e.  mottola , particle creation in de sitter space , phys .
d * 31 * ( 1985 ) 754 .
k.  goldstein and d.  a.  lowe , a note on alpha - vacua and interacting field theory in de sitter space , nucl .
b * 669 * , 325 ( 2003 ) [ arxiv : hep - th/0302050 ] . m.  b.  einhorn and f.  larsen , squeezed states in the de sitter vacuum , phys .
d * 68 * ( 2003 ) 064002 [ arxiv : hep - th/0305056 ] . t.  m.  janssen , s.  p.  miao , t.  prokopec and r.  p.  woodard , the hubble effective potential , arxiv:0904.1151 [ gr - qc ] . t.  m.  janssen , s.  p.  miao , t.  prokopec and r.  p.  woodard , paper in preparation .
t.  inagaki , s.  mukaigawa and t.  muta , a soluble model of four fermion interactions in de sitter space , phys .
d * 52 * ( 1995 ) 4267 [ arxiv : hep - th/9505058 ] .
t.  inagaki , t.  muta and s.  d.  odintsov , dynamical symmetry breaking in curved spacetime : four - fermion interactions ,  prog .
* 127 * ( 1997 ) 93 [ arxiv : hep - th/9711084 ] .
t.  padmanabhan , accelerated expansion of the universe driven by tachyonic matter , phys .
d * 66 * ( 2002 ) 021301 [ arxiv : hep - th/0204150 ] .
n.  c.  tsamis and r.  p.  woodard , plane waves in a general robertson - walker background , class .
grav .   * 20 * ( 2003 ) 5205 [ arxiv : astro - ph/0206010 ] .
a.  vilenkin and l.  h.  ford , gravitational effects upon cosmological phase transitions , phys .
d * 26 * ( 1982 ) 1231 .