nonclassical light generation has attracted much attention over the last 40 years , partly because it not only provides a powerful probe for addressing fundamental issues of quantum theory such as complementarity , hidden variables , and other aspects central to the foundations of quantum mechanics @xcite ; but also holds promise for many potential applications to quantum information processing @xcite , quantum computation and communication @xcite , quantum cryptography @xcite , quantum imaging @xcite , quantum lithography @xcite , and quantum metrology @xcite . in particular , entangled photon pairs
have already been established as a standard research tool in the field of quantum optics .
traditionally , paired photons are produced from spontaneous parametric down conversion ( spdc ) @xcite , in which a strong pump laser drives the atomic oscillators in a noncentrosymmetric crystal into a nonlinear regime , and then two down - converted beams are radiated by these oscillators .
the two down - converted photons from such a nonlinear process usually have very broad bandwidth ( typically in the thz range ) and very short coherence time ( typically around few picoseconds ) so that their waveforms are not resolvable by existing single - photon detectors ( which have resolution around the nanosecond range ) .
in addition , the short coherence length of the spdc biphotons , roughly about @xmath0 @xmath1 m , also limits their application in long - distance quantum communication .
the optical properties of these wide - band biphotons have been well studied both experimentally and theoretically for more than two decades @xcite .
narrow - bandwidth biphotons are ideal for a number of recently proposed protocols for long - distance quantum communication based on coherent interaction between single photons and atomic ensembles @xcite , which require efficiently absorbing biphotons and storing the entanglement .
other applications include the long - distance quantum state teleportation which requires long temporal coherence time of paired photons @xcite . using cavity - enhanced spdc in nonlinear crystals , paired photons with a bandwidth of about 10 mhz and coherence time of about @xmath2 ns
have been reported @xcite .
however , the spdc process is usually tuned to occur in the region of far - off - resonant atomic transitions , which limits the conversion efficiency and requires a large pump power . moreover , in the far - off - resonant pumping the nonlinear parametric coupling coefficient can be treated as a constant and the information about the material internal transition structure is unresolvable from the two - photon correlation measurement .
the two most recently developed advanced technologies make the near- and on - resonance nonlinear processes possible .
the first is laser cooling and trapping of neutral atoms @xcite . because cold atoms with a temperature of about 100 @xmath3 have negligible doppler broadening and very long lifetime , their atomic hyperfine structures
can be resolved without need of a doppler - free setup .
the second is the electromagnetically induced transparency ( eit ) @xcite which not only can eliminate absorption on resonance but also enhance nonlinear interactions at low light level @xcite . following the theoretical proposal in @xcite and early experiments by lukin and kimble groups @xcite , generation of nonclassical correlated photon pairs by exploring the  writing - reading "
procedure has attracted much research interest .
however , the created paired photons are not maximally time - frequency entangled due to two separate operation processes .
using cold atomic - gas media and eit - assisted four - wave mixing ( fwm ) , harris group at stanford university was the first to generate time - frequency entangled narrow - band biphotons with cw driving fields in a double-@xmath4 system @xcite . by controlling the eit at a sufficiently high optical depth ( od ) around 50
, biphotons with subnatural linewidth have been produced in a recent experiment , where the coherence time is up to about 1 @xmath1s in a two - dimensional magneto - optical trap ( mot ) @xcite .
vuletic and his colleagues have also produced biphotons with a temporal coherence length of 100 ns by trapping cold atoms inside an optical cavity @xcite .
, by driving a two - level system into the nonlinear region , have also created paired photons with correlation time of same order of magnitude @xcite .    to describe the physics behind narrow - band biphoton generation from the eit - based multilevel atomic system ,
two different approaches have been presented in the literature .
one approach is to use the heisenberg - langevin method by solving the coupled field operator equations in the heisenberg picture , and the results agree well with the experimental data in the two - photon correlation measurement @xcite .
the advantage of this method is that it provides a sophisticated calculation of the effects from the langevin noises , the medium absorption and optical gain on the biphoton wavepacket .
working in the interaction picture , the other approach developed by wen _
@xcite focuses on the description of the two - photon state vector after the nonlinear medium using perturbation theory .
this state vector theory not only properly describes the experimental results but also offers a clear physical picture on biphoton generation mechanism .
it is the intent of this paper to review and generalize this latter approach by taking into account both linear loss and gain , and provide more insights on the roles of nonlinear interaction and linear response to the biphoton amplitude . using this state vector picture
, we show that the two - photon wave function is a convolution of the nonlinear and linear optical responses .
according to this convolution , the two - photon temporal correlation is considered in two regimes , damped rabi oscillation and group delay .
the paper is organized as the following . in sec .
[ sec : biphotonstate ] , using perturbation theory we lay out a general formulism of the two - photon state and biphoton wavepacket ( or amplitude ) generated from the fwm parametric process . in sec .
[ sec : opticalresponse ] , taking a four - level double-@xmath4 atomic system as an example , we show that the linear and nonlinear optical responses to the generated fields play an important role in determining the two - photon amplitude . by looking at the nonlinear susceptibility , we illustrate the mechanism of biphoton generation in such a system . in sec .
[ sec : singleatomregime ] we describe the damped rabi oscillation regime where the biphoton amplitude is mainly characterized by the structure of the third - order nonlinear susceptibility .
we show that the two - photon temporal correlation can exhibit damped and over - damped rabi oscillaitons . in sec .
[ sec : groupdelayregime ] , we describe the group - delay regime where the biphoton wavepacket is mainly determined by the phase matching due to the linear optical response . in sec .
[ sec : interference ] , we examine the two - photon interference using biphotons generated from such a cold atomic - gas medium .
we show that by manipulating the linear optical response , it is possible to switch the two - photon anti - bunching - like effect to bunching - like effect . finally , we draw the conclusions and outlook in sec .  [
sec : conclusion ] .
( color online ) biphoton generation via a four - level double-@xmath4 atomic system . ( a ) the level structure , where in the presence of a cw pump ( @xmath5 ) and coupling ( @xmath6 ) beams , paired stokes ( @xmath7 ) and anti - stokes ( @xmath8 ) photons are spontaneously created from the four - wave mixing processes in the low - gain regime .
( b ) the forward generation configuration .
( c ) the backward generation geometry .
, width=302 ]    a schematic of biphoton generation via a four - level double-@xmath4 atomic system is shown in fig .
[ fig : fourlevelconfiguration ] , where in the presence of a cw pump ( @xmath5 ) and coupling ( @xmath6 ) lasers , phase - matched , paired stokes ( @xmath7 ) and anti - stokes ( @xmath8 ) photons are spontaneously produced from the fwm process in the low - gain limit . as shown in fig .
[ fig : fourlevelconfiguration](a ) , the strong coupling laser forms a standard three - level @xmath4 eit scheme with the generated anti - stokes field
. therefore , the role of the coupling laser here is that it not only assists the fwm nonlinear process but also creates a transparency window for the anti - stokes photons by the slow - light effect . with the use of eit , it is possible to generate the anti - stokes photons near atomic resonance @xmath9 , which greatly enhances the efficiency of fwm process .
the phase - matching condition allows forward and backward two generation configurations , as shown in fig .
[ fig : fourlevelconfiguration](b ) and ( c ) , respectively . in the backward generation configuration paired stokes and anti - stokes photons can propagate collinearly with the pump and coupling beams @xcite , and can also propagate in a right - angle geometry @xcite .
these flexible generation setups depend on whether the relationship @xmath10 is well satisfied , where @xmath11 are wave vectors of the input pump and coupling fields .
one advantage for the backward generation geometry is that it allows us to easily separate the generated weak fields from the strong inputs .    in the following discussions , we assume that the cold atomic - gas medium consists of identical four - level atoms prepared in their ground state @xmath12 [ see fig .
[ fig : fourlevelconfiguration](a ) ] .
the idealized atoms or molecules are confined within a long , thin cylindrical volume with a length @xmath13 and atomic density is @xmath14 . for simplicity ,
in this paper we will not take into account the doppler broadening and polarization effects .
in the two - photon limit , the quantum langevin noise introduces unpaired photons which are not of interest here and so is ignored .
in addition , we concentrate on the two - photon temporal correlation .
we start with the geometries shown in fig .
[ fig : fourlevelconfiguration](b ) and ( c ) , where stokes - anti - stokes photon pairs propagate along the @xmath15 axis .
we denote the electric field as @xmath16=\frac{1}{2}[e^{(+)}+c.c.]$ ] , where @xmath17 and @xmath18 are the positive- and negative - frequency parts .
the input pump and coupling beams are assumed to be strong classical fields , so that @xmath19},\nonumber \\
e^{(+)}_c(z , t)&=&e_ce^{i[k_c z-\omega_c t]},\label{eq : pcfields}\end{aligned}\ ] ] where @xmath20 are pump and coupling field wave numbers .
the @xmath21 sign in eq .
( [ eq : pcfields ] ) stands for the forward and backward propagation geometries for the pump laser .
the single - transverse - mode stokes and anti - stokes fields are taken as quantized , @xmath22 } , \nonumber \\
\hat{e}^{(+)}_{as}(z , t)&=&\frac{1}{\sqrt{2\pi}}\int d\omega \sqrt{\frac{2 \hbar \omega}{c \varepsilon_0 a}}\hat{a}_{as}(\omega)e^{i[k_{as}(\omega)z-\omega t ] } , \label{eq : fieldoperators}\end{aligned}\ ] ] where @xmath23 are wave numbers of stokes and anti - stokes photons , and @xmath24 is the single - mode cross section area .. the annihilation operators @xmath25 and @xmath26 satisfy the commutation relation @xmath27=[\hat{a}_{as}(\omega),\hat{a}^\dagger_{as}(\omega')]= \delta(\omega-\omega ' ) .
\label{eq : fcommutation}\end{aligned}\ ] ] in the interaction picture , the effective interaction hamiltonian for the four - wave mixing parametric process takes the form @xcite @xmath28 where @xmath29 means the hermitian conjugate , and @xmath30 is the third - order nonlinear susceptibility to the stokes ( or anti - stokes ) field defined by the nonlinear polarizability @xmath31 @xcite .
with use of eqs .
( [ eq : pcfields ] ) and ( [ eq : fieldoperators ] ) , after the @xmath15 integration we rewrite the hamiltonian ( [ eq : hamiltonian0 ] ) as @xmath32 where @xmath33 is the phase mismatching for the forward ( @xmath34 ) and backward ( @xmath35 ) configurations .
when the pump and coupling fields are incident with an angle respect to the @xmath15 axis , the phase mismatching can be generalized as @xmath36 , where @xmath37 is the unit vector along the @xmath38 axis .
the integral over the length @xmath13 gives a sinc function @xmath39 , which determines the two - photon natural spectrum width .
@xmath40 is the nonlinear parametric coupling coefficient , where @xmath41 and @xmath42 are central frequencies of the anti - stokes and stokes fields , respectively .
perturbation theory gives the photon state at the output surface(s ) approximately as a linear superposition of @xmath43 , where @xmath44 is the vacuum state .
the two - photon ( biphoton ) state @xmath45 can be expressed as @xmath46 ) in eq .
( [ eq : biphotonstate0 ] ) , the time integral yields a @xmath47 function , @xmath48 , which expresses the energy conservation for the process .
if @xmath13 is infinite , then the sinc function @xmath49 in ( [ eq : hamiltonian ] ) becomes a @xmath47 function , @xmath50 . in this case , the conditions @xmath51 both hold and the phase matching is said to be perfect .
the phase matching condition arises from the fact that the fwm process is a coherent process .
now the two - photon state ( [ eq : biphotonstate0 ] ) becomes @xmath52 as seen from eq .
( [ eq : biphotonstate ] ) , the two - photon state is entangled in frequency and wave number , but is not entangled in polarization . in frequency space ,
the entanglement is the result of the frequency phase - matching condition , which implies that the detection of a photon at frequency @xmath8 requires the detection of the other photon at frequency @xmath53 .
the frequency correlation has interesting consequences for the temporal behavior of the pair @xcite , as we shall see .
the state is also entangled with respect to the wave number since the sinc function @xmath39 can not be factorized as a function of @xmath54 times a function of @xmath55 . in the general noncollinear case , the wave - number entanglement has implication for the spatial correlation of photon pairs @xcite .
the polarization - entanglement generation has been discussed in @xcite .    to study the optical properties of paired stokes and anti - stokes photons generated from such a four - level system , let us look at the simple experiments of two - photon joint - detection measurement as shown in fig .
[ fig : fourlevelconfiguration](b ) and ( c ) , where anti - stokes photons go to detector d@xmath56 and stokes photons go to d@xmath57 .
the two - photon temporal correlation is the quantity of primary interest in this paper .
the field annihilation operators in time domain on the surfaces are expressed as @xmath58 } , \nonumber \\
\hat{a}_{as}(t)&=&\frac{1}{\sqrt{2\pi}}\int d\omega \hat{a}_{as}(\omega)e^{i[k_{as}(\omega)l/2-\omega t]}.   \label{eq : operatorintime}\end{aligned}\ ] ] it can be shown that they obey the following commutation relation @xmath59=[\hat{a}_{as}(t),\hat{a}^\dagger_{as}(t')]=\delta(t - t').\label{eq : tcommutation}\end{aligned}\ ] ] the joint - detection probability is defined by glauber s theory , @xmath60 in general , @xmath61 is referred to as the _ two - photon amplitude _ , or _
biphoton wavepacket_. with the help of eqs .
( [ eq : biphotonstate ] ) and ( [ eq : operatorintime ] ) , the two - photon wave function ( or amplitude ) now is @xmath62 where @xmath63 with the relative time delay @xmath64 and @xmath65 . in eq .
( [ eq : biphotonwavefunction1 ] ) , @xmath66 is defined as the _ longitudinal detuning function _ ,
@xmath67 from eq .
( [ eq : biphotonwavefunction1 ] ) we see that the biphoton wave function is determined by both the nonlinear coupling coefficient ( @xmath68 ) and the longitudinal detuning function ( @xmath69 ) . in the time domain , @xmath70 is a convolution ( @xmath71 ) of @xmath72 and @xmath73 @xmath74=l\int d\tau'\tilde{\kappa}(\tau')\tilde{\phi}(\tau-\tau ' ) .
\label{eq : convolution}\end{aligned}\ ] ] here , @xmath75 is fourier transform of the nonlinear coupling coefficient @xmath76 , @xmath77 and similarly , @xmath78 is fourier transform of the longitudinal detuning function @xmath66 , @xmath79 the photon pair emission rate @xmath80 can be evaluated , by using eqs .
( [ eq : biphotonwavefunction1 ] ) , as @xmath81    now let us consider a photon counting experiment .
let the detectors be at equal distances from the source so that the relative time delay caused by the free - space propagation of electric fields are equal .
suppose the stokes photon is detected by d@xmath57 at time @xmath82 , and after a controllable delay ( @xmath83 ) detector d@xmath56 is opened to detect the anti - stokes photon .
assuming perfect detection efficiency , the averaged two - photon coincidence counting rate is determined by @xmath84 @xmath85 is the coincidence window function within the time bin width @xmath86 : @xmath87 for @xmath88 ; otherwise , @xmath89 .
when the two - photon coherence time is much larger than the detector time bin width @xmath86 , eq .
( [ eq : rcc ] ) becomes @xmath90 , which is a measure of the two - photon correlation function .
this is the essence of narrow - band paired photons measurement discussed in this paper
.    equation ( [ eq : biphotonwavefunction ] ) is the starting point for the discussions in the following sections . by extending wave numbers as complex variables
, the effects of loss and gain inside the material can be properly included in the above theory , as we shall see .
it is found that to be consistent with the heisenberg - langevin theory @xcite in the low - gain limit , the argument in @xmath69 should be replaced by @xmath91 , where @xmath92 is the conjugate of @xmath93 .    before proceeding the discussion let us make a short summary .
in this section , we have laid out the general formulism of the two - photon state and wave function generated in a fwm medium pumped by two cw classical fields using perturbation theory .
the two - photon state ( [ eq : biphotonstate ] ) is a coherent superposition and it is entangled in frequency and wave numbers . the biphoton wave function , as shown in eq .
( [ eq : convolution ] ) , is a convolution of @xmath72 and @xmath73 .
alternatively , eq .  ( [ eq : biphotonwavefunction1 ] ) states that the two - photon spectrum is determined by both the nonlinear susceptibility and the phase matching .
this implies two regimes for the temporal correlation measurement , as discussed in sec .
[ sec : singleatomregime ] and sec .
[ sec : groupdelayregime ] .
the methodology explored here is equivalent to the heisenberg - langevin theory @xcite at the two - photon level .
moreover , this method offers a clear physical picture of the paired photons , which makes it easier to understand the two - photon interference experiment @xcite . using
the above theory new phenomena have been analyzed @xcite which are not easy to be interpreted from the heisenberg - langevin theory . for example , the two - photon correlation at a high optical depth @xcite , the simulation by numerically solving the heisenberg - langevin coupled - operator equations shows a good agreement in the leading edge spike with the experiment data , but the nature of optical precursors is much easier to be understood in the present biphoton state picture @xcite .
as shown in eq .
( [ eq : biphotonwavefunction ] ) , the two - photon wavefunction is determined by both the nonlinear susceptibility and the longitudinal detuning function .
therefore , in this section we look at the linear and nonlinear optical responses to the generated fields . on the one hand ,
the importance of the nonlinear susceptibility is that it not only characterizes the strength of the nonlinear parametric process but also is the mechanism of biphoton generation . as a consequence , it determines the feature of the two - photon temporal correlation in one regime . on the other hand
, the longitudinal detuning function ( related to the linear susceptibilities ) sets the natural spectral width of biphotons , and governs the pattern of the temporal correlation in another regime .
the optical response of a four - level double-@xmath4 eit system has been discussed in @xcite .
et al_. have further systematically studied the optical responses in a three - level double-@xmath4 scheme @xcite and a two - level system @xcite .
recently , the responses of a four - level inverted - y system have been analyzed in @xcite .    for simplicity
, here let us still choose paired stokes - anti - stokes generation shown in fig .
[ fig : fourlevelconfiguration ] as an example to illustrate the physics behind fwm .
assuming that the intensity of the off - resonant pump laser is chosen so that the atomic population remains primarily in the ground level @xmath12 , the third - order nonlinear susceptibility for the generated anti - stokes field is @xcite @xmath94}\nonumber \\
= \frac{-n\mu_{13}\mu_{32}\mu_{24}\mu_{41}/(\varepsilon_0\hbar^3)}{4(\delta_{p}+i\gamma_{14})(\omega-\omega_e/2+i \gamma_{e})(\omega+\omega_e/2+i\gamma_{e})},\label{eq : chi3}\end{aligned}\ ] ] in eq .
( [ eq : chi3 ] ) , @xmath95 are the electric dipole matrix elements , @xmath96 is the coupling rabi frequency , and @xmath97 are dephasing rates , respectively .
@xmath98 is the pump detuning from the atomic transition @xmath99 .
@xmath100 is the detuning of the anti - stokes photons from the transition @xmath9 , and we will take @xmath101 as the anti - stokes central frequency .
@xmath102 is the effective coupling rabi frequency .
@xmath103 is the effective dephasing rate .
the third - order nonlinear susceptibility has two resonances separated by @xmath104 and each is associated with a linewidth of @xmath105 .
the two - photon coherence time is inversely proportional to @xmath105 .
more importantly , the two resonances here indicate two types of fwm processes .
alternatively , two types of paired stokes and anti - stokes photons can be generated from these two fwm processes in the spontaneous emission region .
one fwm process happens when the central frequency of the anti - stokes photons is @xmath106 and the central frequency of the correlated stokes photons is @xmath107 .
the other fwm occurs as the anti - stokes photons have the central frequency at @xmath108 while the paired stokes photons have the frequency at @xmath109 .
this is the generation mechanism of biphotons produced from the fwm in atomic - gas media .
consequently , the interference between these two types of biphotons will appear in the two - photon temporal correlation , as shall be discussed in sec .
[ sec : singleatomregime ] .
the results obtained here agree well with the dressed - state picture @xcite .
the linear susceptibilities at the stokes and anti - stokes frequencies are , respectively , @xmath110 where @xmath111 is the pump rabi frequency .
the complex wave numbers of stokes and anti - stokes photons are obtained from the relations @xmath112 and @xmath113 , where the imaginary parts stand for the raman gain and eit loss , respectively . the linear susceptibility ( [ eq : chilinear ] ) to the anti - stokes field
takes the form of standard eit .
as expected , the coupling beam @xmath114 not only assists the propagation of the anti - stokes photons by creating an eit window , but also manipulates the group velocity using the slow - light effect .
another role played by eit here is that it allows the nonlinear optics occurring near atomic resonance with small absorption and hence enhances the efficiency of the nonlinear interaction .
it is instructive to examine the linear responses under some simplifying approximation . taking @xmath115 , eqs .
( [ eq : chis ] ) and ( [ eq : chilinear ] ) give @xmath116 and @xmath117 so that the wave - number mismatching is @xmath118 . here
@xmath119 is the central wave number of the anti - stokes field , @xmath120 is its group velocity , and @xmath121 , the imaginary part of the anti - stokes wave number , characterizes the eit loss , respectively . now eq .  ( [ eq : detuningfunction ] ) can be approximated as @xmath122 which sets a coherence time about @xmath123 . in eq .
( [ eq : phasematchingapp ] ) , @xmath124 where @xmath125 is the on - resonance absorption cross section in the transition @xmath9 .
as discussed in this section , there are three important characteristic frequencies that determine the shape of the biphoton wave function .
the first is the coupling effective rabi frequency @xmath104 , which determines the two - resonance spectrum of the nonlinear susceptibility .
the second is the linewidth @xmath105 of the two resonances in the nonlinear susceptibility .
the third is the full - width - at - half - maximum ( fwhm ) phase - matched bandwidth determined by the sinc function , @xmath126 , where @xmath127 is the anti - stokes group delay time . in time domain , they correspond to the rabi time @xmath128 , nonlinear coherence time @xmath129 and group delay time @xmath127 .
the group delay time can be estimated from @xmath130 with the optical depth defined as @xmath131 .
the competition between @xmath132 , @xmath133 and @xmath127 will determine which effect plays a dominant role in governing the feature of the two - photon correlation .
therefore , in the following we will discuss the two - photon joint - detection measurement in two regimes , damped rabi oscillation and group delay .
let us first look at the damped rabi oscillation regime , where optical properties of the two - photon amplitude ( [ eq : biphotonstate ] ) are mainly determined by the nonlinear coupling coefficient .
this regime requires that the effective coupling rabi frequency @xmath104 and linewidths @xmath134 be smaller than the phase - matching bandwidth @xmath135 or @xmath136 and @xmath137 , so that we can treat the longitudinal detuning function as @xmath138 and @xmath139 .
the biphoton spectral generation rate is proportional to @xmath140 .
hence both biphoton spectrum intensity and emission rate are proportional to @xmath141 .
if @xmath142 , the two - photon temporal correlation clearly exhibits a damped rabi oscillation with oscillation period @xmath143 ; if @xmath144 , the second and following oscillations are suppressed because of the short dephasing time , as experimentally demonstrated in @xcite .
let us consider the case of @xmath145 first , which implies a real effective coupling rabi frequency @xmath104 .
the two - photon wave function now is determined by fourier transform of the nonlinear coupling coefficient ( [ eq : nonlinearresponse ] ) , which is @xmath146\theta(\tau ) .
\label{eq : f}\end{aligned}\ ] ] and @xmath147 .
the heaviside step function here is defined as @xmath148 for @xmath149 , and @xmath150 for @xmath151 . combining with eq .
( [ eq : biphotonwavefunction ] ) , the biphoton wavepacket becomes @xmath152\theta(\tau),\label{eq : total}\end{aligned}\ ] ] where @xmath153 and @xmath101 have been used .
the physics of eq .
( [ eq : total ] ) is understood as follows .
because the two - photon state is entangled , it can not be factorized into a function of @xmath154 times a function of @xmath155 .
@xmath156 , depends only on the relative time delay @xmath83 , which implies that the pair is randomly generated at any time .
if the wave packets of the pump and the coupling beams are not plane waves , the term @xmath157 would become a wave packet with a coherence length determined by the pump and coupling . equation ( [ eq : total ] ) shows also destructive interference between the two biphoton channels . the first term in the bracket on the rhs of eq .
( [ eq : total ] ) represents the two - photon amplitude of paired anti - stokes at @xmath158 and stokes at @xmath159 ; while the second term is the two - photon amplitude of paired anti - stokes at @xmath160 and stokes at @xmath161 .    to further see the interference ,
let us look at the two - photon glauber correlation function , @xmath162\theta(\tau),\label{eq : rabig2}\end{aligned}\ ] ] the physics of eq .
( [ eq : rabig2 ] ) shows the beating ( or interference ) between two types of paired photons generated from the two fwm processes .
the two - photon correlation ( [ eq : rabig2 ] ) exhibits damped rabi oscillations of period @xmath143 .
the damping rate is determined by the resonant linewidth @xmath105 in the doublet . at @xmath163
, @xmath164 vanishes and as @xmath165 , @xmath164 also approaches to zero .
this indicates photon anti - bunching - like effect [ @xmath166 . in the counter - propagation geometry , if @xmath167 , @xmath168 becomes a symmetric distribution about @xmath163 , because spontaneously emitted paired photons can trigger both detectors d@xmath56 and d@xmath57 as experimentally verified in @xcite .
figure  [ fig : rabioscillation](a ) shows the first experimental demonstration of this type of biphoton source by bali _
_ in a cold @xmath169rb gas @xcite .
the solid lines in fig .
[ fig : rabioscillation](a1 ) and ( a2 ) are theoretical curves obtained from the heisenberg - langevin theory . in fig .
[ fig : rabioscillation](b1 ) and ( b2 ) , the curves are theoretical simulations using eq .
( [ eq : rabig2 ] ) with the same experimental parameters in ( a1 ) and ( a2 ) .
comparing with the heisenberg - langevin theory , the state vector picture also gives an agreement with the experimental data . from eq .
( [ eq : rabig2 ] ) , fig .
[ fig : rabioscillation](a1 ) and ( b1 ) correspond to the case of @xmath170 where the oscillations can be clearly observed before damping out ; while fig .
[ fig : rabioscillation](a2 ) and ( b2 ) correspond to the case of @xmath171 where only the first oscillation is observable and other disappear due to the fast dephasing rate @xmath172 .
although in @xcite the lack of other oscillations in fig .
[ fig : rabioscillation](a2 ) were interpreted due to the group delay , we believe that both cases shown in fig .
[ fig : rabioscillation](a1)and ( a2 ) are in the damped rabi oscillation regime and the correlation time is well fitted with @xmath133 and @xmath132 .
this can be further verified by reducing the optical depth by a factor of 1/2 , which is corresponding to a reduction of the group delay by a factor of 1/2 .
however the correlation time widthes in fig .
[ fig : rabioscillation](a2 ) hold the same .    if @xmath173 , the effective coupling rabi frequency becomes purely imaginary , @xmath174 .
the biphoton wave function becomes @xmath175 at the weak coupling limit , @xmath176 and eq .  ( [ eq : overdamped1 ] ) takes the form @xmath177 this can be understood as follows .
since the coupling beam is very weak , the atomic - level structure is not changed due to the input power .
two types of biphotons are generated with the same central frequency but are associated with different linewidths , @xmath178 and @xmath179 , respectively .
the interference between these two types of biphotons gives a manifested exponential decay .
in such a limit , the eit disappears and the medium absorption decreases the biphoton generation efficiency .
this case might be observable in the dilute atomic - gas medium with very low optical depth ( @xmath180 ) .
suggested by balic _ _ et al.__@xcite and demonstrated by du _ _ et al.__@xcite , the group delay regime is defined as @xmath181 and the ( eit ) slow light effect can be used to dynamically control the biphoton temporal correlation time .
the group delay condition is equivalent to @xmath182 , i.e. , the biphoton bandwidth is determined by the phase - matching spectrum .
therefore , in this section we treat the third - order nonlinear susceptibility as a constant over the phase - matching spectrum . as a consequence , the mechanism of double resonance of biphoton generation
is washed out by the longer group - delay time .
the two - photon correlation in this regime looks more like that with conventional spdc photons .
the coincidences tends to be a square - like pattern .
the biphoton wave function can be approximated as @xmath183 where @xmath184 is the on - resonance nonlinear coupling constant .
we find that in this regime , the eit profile modifies the biphoton wave function through transparency window and slow light .
when the eit bandwidth ( @xmath185 ) is larger than the phase - matching bandwidth , the anti - stokes loss can be ignored and the two - photon wave packet approaches a rectangular shape . however , when the eit loss is significant , the biphoton wave packet follows an exponential decay .
when @xmath186 , we can ignore the loss in eq .
( [ eq : phasematchingapp ] ) and rewrite it as @xmath187 its fourier transform gives the biphoton wave function @xmath188 the rectangular function @xmath189 , ranging from @xmath163 to @xmath190 , shows that the anti - stokes photon is always delayed with respect to its paired stokes photon because of the slow light effect .
it is clear that the two - photon correlation time is determined by the group delay @xmath191 .
the photon pair generation rate can be calculated from eq .
( [ eq : pairrate ] ) as @xmath192 thus even though the spectral generation rate , @xmath193 , scales as @xmath141 , the total rate of paired counts scales linearly as @xmath194 because the bandwidth reduces linearly with optical depth @xcite . as explained by rubin _
et al_. @xcite , the rectangular waveform may be understood as follows : the paired photons are always produced from the same , but unknown , space point in the nonlinear medium .
if emitted from the front surface , the anti - stokes photon has no delay , and both photons arrive at detectors simultaneously .
however , if emitted from the back surface , the anti - stokes photon is delayed relative to the stokes photon by @xmath127 . for conventional spdc photons ,
the rectangular - shaped biphoton wave packet has been well known and its sub - ps correlation time has been directly measured in @xcite .
interestingly , we find that to observe the rectangular shape , the condition @xmath186 sets a lower bound for the optical depth . using @xmath185 and @xmath195
, we obtain @xmath196 .
it can be shown that such an optical depth is necessary to have the eit delay - bandwidth product greater than two @xcite .
( color online ) two - photon correlation function in the group delay regime .
( a ) rectangular - shape correlation for an ideal spontaneous parametric down converter with a relative group delay @xmath190 and pair emission rate @xmath80 . ( b )
experimental data ( red ) and theoretical curve ( solid blue line ) obtained from the heisenberg - langevin theory in ref .
( c ) propagation delay of an anti - stokes pulse ( red ) with respect to a reference pulse ( black ) .
data are taken from ref .
( d ) theoretical simulation using the perturbation theory fitted with same experimental parameters in ( b).,width=302 ]    figure  [ fig : groupdelay](a ) shows an ideal rectangular - shape correlation with a group delay of @xmath191 . in fig .
[ fig : groupdelay](b ) the experimental demonstration of a near - rectangular - shape correlation is obtained at the optical depth around 53 @xcite , where experimental data are denoted by  + " .
the solid line is theoretical fitting using the heisenberg - langevin theory . as seen from fig .
[ fig : groupdelay](b ) , the correlation time of about 300 ns is achieved and is consistent with the anti - stokes pulse delay measurement shown in the inset ( c ) , where the smaller curve represents the delayed signal and the higher one is the reference pulse .
the experimental parameters used here are @xmath197 mhz , @xmath198 , @xmath199 , @xmath200 , @xmath201 , and @xmath202 , respectively .
the eit transparency width is estimated around @xmath203 mhz and the phase - matching spectral width about @xmath204 mhz .
the theoretical simulation using the perturbation theory [ eq .
( [ eq : biphotonwavefunction1 ] ) ] , plotted in fig .
[ fig : groupdelay](c ) , also fits very well with the reported experimental data .
it is found that the exponential - decay behavior at the tail is due to the finite eit loss , which alters the correlation function away from the ideal rectangular shape . in the experiment @xcite ,
_ have also experimentally verified that the total coincidence counts are linearly dependent on the optical depth .
the sharp peak shown in the leading edge of the two - photon coincidence counts in fig .
[ fig : groupdelay](b ) is the first observed sommerfeld - brillouin precursor @xcite at the biphoton level @xcite .
mathematically , the integration in the two - photon amplitude ( [ eq : biphotonwavefunction1 ] ) is similar with the classical pulse propagation . in the stationary - phase approximation , starting from eq .
( [ eq : biphotonwavefunction1 ] ) one can show that the sharp peak is resulted from the interference of biphotons which are generated out of the eit opacity window .
the detailed analysis of the precursor generation at the two - photon level has been presented in @xcite by using the theory laid out in sec .
[ sec : biphotonstate ] .
although by solving the coupled - operator equations the heisenberg - langevin theory predicts the correct pattern , the physics behind the precursor is not easy to be retrieved from this theory .
suppose now the eit loss is not negligible , i.e , @xmath205 .
in such a case , the longitudinal detuning function ( [ eq : phasematchingapp ] ) can be approximated as @xmath206 the two - photon wave function becomes @xmath207 the correlation function @xmath208 has an exponential decay time of @xmath209 .
this means that the anti - stokes photon detected at @xmath83 is generated at the depth of @xmath210 such that the amplitude decreases by @xmath211 .
in sec .  [ sec : singleatomregime ] and sec .  [
sec : groupdelayregime ] , we have discussed the two - photon temporal correlation in two regimes . in general , the two - photon wave function is a convolution of both the nonlinear and linear responses as shown in eq .
( [ eq : convolution ] ) . in both regimes ,
the correlation functions have non - zero values only at @xmath212 . the observed anti - bunching - like effect in the damped rabi oscillation regime
is due to the destructive interference between two types of biphoton generation . in principle , if the detection system is frequency sensitive enough to resolve these two types of biphotons , the damped rabi oscillations will disappear and the square - wave - like pattern manifested by exponential decay will be observed .
( color online ) biphoton generation in a three - level system .
( a ) right - angle experimental configuration .
( b ) three - level double-@xmath4 energy level diagram.,width=302 ]    in this section , we are going to briefly show that it is possible to switch the two - photon coincidence counting rate from photon anti - bunching - like effect to bunching - like effect , by manipulating the linear optical responses to both stokes and anti - stokes fields @xcite .
we consider a case that stokes and anti - stokes photons are degenerate at both spectrum and polarization so that both stokes and anti - stokes experience slow light propagation . as shown in fig .
[ fig : feyman](a ) and ( b ) , we consider paired stokes and anti - stokes photons produced from a three - level double-@xmath4 scheme in a right angle configuration . in the counter - propagation geometry ,
frequency - phase matching allows both stokes and anti - stokes fields to reach both detectors .
the two generated fields both may experience slow - light effect inside the medium due to the linear optical responses .
the biphoton wave function now is expressed as @xmath213\ast\pi(\tau ;- l / v_g , l / v_g).\label{eq : convolution2}\end{aligned}\ ] ] even though there is no temporal overlap between @xmath75 and @xmath214 , their convolution with @xmath215 mixes them together in the time domain . by manipulating the group delay , the two - photon interference in eq .
( [ eq : convolution2 ] ) can be switched from the destructive to the constructive .
alternatively , the coincidence counting rate is changed from photon anti - bunching - like effect to bunching - like effect .
this can be achieved by making the indistinguishability of their propagation pathways which gives rise to the biphoton interference .
the interference comes from two indistinguishable paths for paired stokes and anti - stokes photons .
in one possible propagation pathway , the pair is created at @xmath216 with stokes photon going to the left and anti - stokes photon to the right . in the other propagation pathway ,
the pair is created at @xmath217 but with stokes photon going to the right and anti - stokes photon to the left .
the interference between the two paths then modifies the two - photon correlation by further altering the group delay .
the details of this two - photon interference are presented in @xcite .
recently , wen _
_ has proposed a scheme to switch the two - photon interference from photon anti - bunching - like effect to bunching - like effect by manipulating the linear optical responses to the paired photons in a strongly driven two - level system @xcite .
in conclusion , using perturbation theory , we have discussed the nature of the two - photon state and its waveform for narrow - band biphoton generation near atomic resonance via fwm parametric process using eit and slow light . in general , the biphoton waveform is determined by a convolution of the nonlinear response ( determined by the third - order nonlinear susceptibility ) and the linear response ( determined by the longitudinal detuning function ) .
therefore , the two - photon wave function can be analytically solved in two regimes , the damped rabi oscillation and the group delay . in the damped rabi oscillation regime , where the biphoton wave function is mainly determined by the third - order nonlinear susceptibility , the two - photon temporal correlation exhibits a damped rabi oscillation because of the interference between two types of fwms . in the group delay regime , where the two - photon amplitude is mainly determined by the phase - matching condition
, the correlation approaches a rectangular shape manifested by an exponential decay .
the theoretical simulations agree very well with recent experimental results @xcite and the heisenberg - langevin theory @xcite .
although the discusses are based on a four - level double-@xmath4 eit system , the analysis can be easily extended to three - level and two - level systems .    before concluding , we note that a number of other interesting avenues are presently being explored that are closely related with this new type of biphoton source . by manipulating the linear optical responses to the generated fields
, bell inequality can be tested by two - photon beatings using biphotons generated from a two - level system @xcite .
it is interesting to note that the correspondence principle fails in describing paired - photon generation in a two - level system . as mentioned in @xcite , it is not generally valid to simply replace the coupled - field maxwell s equations by coupled field - operator equations without changing the coefficients .
but the underlying physics is not clear at this moment
.    it will be interesting to explore quantum networks with these narrow - band paired photons . for example , the feasibility of an eit - based quantum memory can be easily demonstrated using current biphoton source by looking at the two - photon interference , as suggested in @xcite . it may be interesting to look at the squeezing properties after an eit system if the input squeezed light is initially generated from the atomic - gas medium @xcite .
as proposed in @xcite experimental demonstration of polarization - entangled narrow - band biphoton generation would be useful for quantum computation and communication .
another problem is that it is not clear whether the maximally time - frequency entangled photon pairs can be generated in the room - temperature atomic - gas medium .
since all the experimental and theoretical work now are focused on cold atomic ensembles , it would be interesting to know if the doppler broadening becomes dominant in the process , what kind of feature would be observed in the two - photon correlation measurement ? by optimizing the system ,
can we still observe the damped rabi oscillations ?    apart from temporal correlation , a possible application of these narrow - band biphotons may involve the transverse correlation @xcite .
it may be worth studying whether this new type of source is useful for quantum imaging and lithography .
s. du acknowledges the start - up financial support from the department of physics at the hong kong university of science and technology . j .-
wen and m. h. rubin were supported in part by the u.s .
army research office under muri grant w911nf-05 - 1 - 0197 .
a. n. boto , p. kok , d. s. abrams , s. l. braunstein , c. p. williams , and j. p. dowling ,  quantum interferometric optical lithography : exploiting entanglement to beat the diffraction limit , " phys .
lett . * 85 * , 2733 - 2736 ( 2000 ) .
a. migdall , r. datla , a. v. segienko , j. s. orszak , and y .- h .
shih ,  measuring absolute infrared spectral radiance with correlated visible photons : technique verification and measurement uncertainty , " appl
. opt . * 37 * , 3455 - 3463 ( 1998 ) .
a. kuzmich , w. p. bowen , a. d. boozer , a. boca , c. w. chou , l .- m .
duan , and h. j. kimble ,  generation of nonclassical photon pairs for scalable quantum communication with atomic ensembles , " nature ( london ) * 423 * , 731 - 734 ( 2003 ) .
wen and m. h. rubin ,  transverse effects in paired - photon generation via an electromagnetically induced transparency medium .
i. perturbation theory , " phys .
a * 74 * , 023808 ( 2006 ) ;  transverse effects in paired - photon generation via an electromagnetically induced transparency medium .
ii . beyond perturbation theory , " phys . rev . *
74 * , 023809 ( 2006 ) .
wen , m. h. rubin , and s. du ,  two - photon beatings using biphotons generated from a two - level system , " _ slow light and fast light conference _
( utah , 2007 ) ; j .-
wen , s. du , m. h. rubin , and e. oh ,  two - photon beatings using biphotons generated from a two - level system , " opt .
express ( to be published ) .


observations of the cmb , supernovae , lensing and baryon acoustic oscillations have led to the cosmological standard model which requires an accelerated expansion of the late universe , driven by dark energy but despite many years of research its origin has not yet been identified .
there are two major explanations for the origin and properties of dark energy .
the first solution consists of introducing a cosmological constant @xmath0 with a constant energy density causing an effective repulsive force between cosmological objects at large distances . from the particle physics perspective the cosmological constant could correspond to the vacuum energy density .
the theoretical expectations for the vacuum energy density caused by fluctuating quantum fields , however , exceeds the observational bounds on @xmath0 by up to 120 orders of magnitude .
this discrepancy remains for more almost a century one of the most challenging puzzles in physics , @xcite .
alternatively , the acceleration of the universe can be explained by introducing new dynamical degrees of freedom , either by invoking new fluids @xmath1 with negative pressure or by changing the geometrical part of einstein s equations . in particular , weakening gravity on cosmological scales could not only tackle the cosmological constant problem , but would also come hand in hand with new degrees of freedom which might be responsible for a late - time speed - up of the hubble expansion .
such scenarios could arise in massive gravity or in higher - dimensional frameworks .    in the higher dimensional picture ,
the dvali - gabadadze - porrati ( dgp ) model is one of the important large scale modified theories of gravity @xcite . in this braneworld model
our universe is confined to a three - brane embedded in a five - dimensional bulk . on small scales ,
four - dimensional gravity is recovered due to an intrinsic einstein hilbert term sourced by the brane curvature , whereas on larger , cosmological scales gravity is systematically weaker as the graviton acquires a soft mass @xmath2 which limits its effective range . being a fundamentally higher dimensional theory ,
the effective four - dimensional graviton on the brane carries five degrees of freedom , namely the usual helicity-2 modes , two helicity-1 modes and one helicity-0 mode .
whilst the helicity-1 mods typically decouple , the helicity-0 one can mediate an extra fifth force . in the limit @xmath3 , one recovers general relativity ( gr ) through the vainshtein mechanism : the basic idea is to decouple the additional modes from the gravitational dynamics via nonlinear interactions of the helicity-0 mode of the graviton , @xcite . as a result , at the vicinity of matter , the non - linear interaction for the helicity-0 mode become large and hence suppresses its coupling to matter .
this decoupling of the nonlinear helicity-0 mode is manifest in the limit where @xmath4 and @xmath5 while the strong coupling scale @xmath6 is kept fixed .
this limit enables a linear treatment of the usual helicity-2 mode of gravity while the helicity-0 mode @xmath7 is described non - linearly , which is the so - called decoupling limit .
one of the successes of the dgp model is the existence of a self - accelerating solution , where the acceleration of the universe is sourced by the graviton own degrees of freedom ( more precisely its helicity-0 mode ) .
unfortunately that branch of solution seems to be plagued by ghost - like instabilities @xcite , in the dgp model , but this issue could be avoided in more sophisticated setups , for instance including gauss - bonnet terms in the bulk @xcite .
more recently , it has been shown that the decoupling limit of dgp could be extended to more general galilean invariant interactions @xcite .
this galileon model relies strongly on the symmetry of the helicity-0 mode @xmath7 : invariance under internal galilean and shift transformations , which in induced gravity braneworld models can be regarded as residuals of the 5-dimensional poincar invariance .
these symmetries and the postulate of ghost - absence restrict the construction of the effective @xmath7 lagrangian .
there exist only five derivative interactions which fulfill these conditions . from the five dimensional point of view
these galilean invariant interactions are consequences of lovelock invariants in the bulk of generalized braneworld models , @xcite . since their inception there has been a flurry of investigations related to self - accelerating de sitter solutions without ghosts @xcite , galileon cosmology and its observations @xcite , inflation @xcite , lensing @xcite , superluminalities arising in spherically symmetric solutions around compact sources @xcite , k - mouflage @xcite , kinetic gravity braiding @xcite , etc ... .
furthermore , there has been some effort in generalizing the galileon to a non - flat background .
the first attempt was then to covariantize directly the decoupling limit and to study its resulting cosmology @xcite .
in particular it was shown in @xcite that the naive covariantization would yield ghost - like terms at the level of equation of motion but a given unique nonminimal coupling between @xmath7 and the curvature can remove these terms resulting in second order of equations of motion @xcite , which are also consistent with a higher - dimensional construction @xcite . in this paper , we will pursue the same strategy when constructing our proxy theory .
the outcome of this kind of covariantization method was explored further in @xcite .
while this covariantization is ghost - free , the galileon symmetry is broken explicitly in curved backgrounds .
only recently there has been a successful generalization to the ( anti- ) de sitter background and ultimately to maximally symmetric backgrounds where it has been discovered that the de sitter galileon interactions acquire additional , potential - like terms being functions of @xmath7 and mixtures of @xmath7 and gradients of @xmath7 , all fulfilling a generalized galileon symmetry @xcite .
there exists a parallel to theories centered on a massive graviton : galileon - type interaction terms naturally arise in gravitational theories using a massive spin-2 particle as an exchange particle , which has , in addition , been constructed to be ghost - free be it in three dimensions , @xcite or for an generalized fierz - pauli action in four dimensions @xcite .
not only is the existence of a graviton mass a fundamental question from a theoretical perspective , it could also have important consequences both in cosmology and in solar system physics , @xcite .
although solar system observations have confirmed gr to high accuracy and placed bounds on the graviton mass to be smaller than a few @xmath8ev , even such a small mass would become relevant at the hubble scale which corresponds to the graviton compton wavelength .
in particular , it has been successfully shown that this massive gravity theory exhibits a stable self - accelerated solution in the decoupling limit since the scalar mode can generate a constant negative pressure density . in the decoupling limit
, the expansion history of the universe in this self - accelerating branch was found to be indistinguishable from @xmath0cdm , @xcite .
while the self - accelerating solutions in the above models yield viable expansion histories including late - time acceleration , they do not address the cosmological constant problem , the giant mismatch between the theoretically computed high energy density of the vacuum and the low observed value .
a possible answer comes from the idea of degravitation , which asserts that the energy density could be as large as the theoretically expected value , but would not bear a large effect on the geometry .
technically , gravity is less strong on large scales ( ir - limit ) and could act as a high - pass filter suppressing the gravitational effect of a potentially large vacuum energy . since such modifications of gravity in the ir naturally arise in models of massive gravity , they logically provide a possible mechanism to degravitate the vacuum energy density , @xcite , which was observed in bi - galileon models @xcite as well as in the decoupling limit of massive gravity , @xcite .
analogously , the dgp braneworld model can be extended to higher dimensions to tackle the cosmological constant problem as well , @xcite .    in this paper , we focus on the covariantization of a ghost - less extension of fierz - pauli massive gravity recently proposed in @xcite and
show that this proxy model allows for a stable self - accelerating solution .
hereby , we have performed the covariantization in the jordan as well as in the einstein frame .
we discuss the differences between the two approaches and the consequences of choosing the wrong starting frame . in the well defined proxy theory
we are able to tackle the puzzle about the self - acceleration of the universe but not the one about the cosmological constant problem .
furthermore , we study the perturbations around the self - accelerating background and provide the expression for the modified evolution equation for these density perturbations .
as expected , the @xmath7 field enhances the gravitational clustering resulting in a rapid growth of structures .
this result is quite different from the one we had obtained in @xcite , where the self - accelerating solution was indistinguishable from a @xmath0cdm , but we emphasize that the theory used in this paper , is a proxy model distinct from massive gravity . at early times ,
the enhancement of clustering is restrained since the galileon self - interactions are the dominant ones suppressing their energy density relative to that of matter or radiation .
once the matter density has dropped sufficiently , the galileon become an important contribution to the dark sector of the universe driving cosmic expansion .
+ the paper is organized as follows : we start in section [ sec : mg1 ] with a summary of the non - linear theory of massive gravity @xcite which we use for constructing our proxy theory , where we discuss the successful implementation of ghost - free massive gravity in the decoupling limit .
we then move to our proxy theory in section [ sec : proxy ] where we first give some motivations on why we use a proxy theory rather than the exact non - linear theory .
we then work out the consequences of this covariantization for cosmology and specifically for late - time acceleration in section [ sec : sa ] where we study also the stability conditions and the reason why a degravitating solution can not be constructed in our covariantized theory .
we then move onto more general cosmology in section [ sec : cosmology ] and follow the helicity-0 mode contribution to the universe throughout its evolution , before quickly presenting the consequences for structure formation in section [ sec : structure ] .
we use the same notation as @xcite , with @xmath9 , @xmath10 and @xmath11 where the covariant derivative is taken w.r.t .
@xmath12 and square brackets @xmath13 $ ] represent the trace of a tensor @xmath14=\pi_{\mu\nu}\pi^{\mu\nu}$ ] and @xmath15 ^ 2=\pi^\mu_\mu\pi^\nu_\nu$ ] .
the first theory of massive gravity was proposed by fierz and pauli in 1939 @xcite , but was shown to be unstable by boulware - deser ( bd ) @xcite , due to the non - propagation of the hamiltonian constraint at the non - linear level .
this instability can also be seen in the language @xcite where in the decoupling limit , the helicity-0 mode typically has equations of motion with more than two derivatives , and hence does not possess a well defined cauchy surface .
however it was shown in @xcite that the graviton potential could be built in such way as to remove any higher derivative term in the equations of motion , and obtain a galileon - type of action for the additional helicity-0 mode . to review this ,
let us start with a graviton of mass @xmath2 described by l=(r - u(g , h))[eq : full ] with the potential @xmath16 , where the tensor @xmath17 is constructed in terms of the metric @xmath18 and the four fields @xmath19 by @xmath20 .
we can then split the fields into helicity-1 and -0 contributions , but the helicity-1 mode decouples in the decoupling limit and can be consistently set to zero .
we therefore focus on the helicity-0 mode @xmath7 and write @xmath21 such that @xmath22 . defining the quantity @xmath23 the most generic potential that bears no ghosts in the decoupling limit is u(g , h)=-4@xmath24[eq : fullu ] where @xmath25 are two free parameters and _ 2&=&[]^2-[^2 ] + _ 3&=&[]^3 - 3 [ ] [ ^2]+2[^3 ] + _ 4&=&[]^4 - 6[^2][]^2 + 8[^3 ] [ ] + 3[^2]^2 - 6[^4 ] .
notice that @xmath26 can be expressed in terms of @xmath27 and the tadpole @xmath28 $ ] , @xcite .
it has then been shown that ghost - like pathologies in this theory of massive gravity theory disappear to all orders in the decoupling limit and at least up to quartic order beyond the decoupling limit , as well as completely non - linearly in some specific cases . while this theory is completely covariant , studying the cosmology as well as other non - trivial curved geometries can be extremely complicated . instead
, we will here focus on its version in the decoupling limit and covariantize the theory directly from this limit .
we emphasize that the resulting proxy theory will be distinct from the theory of massive gravity presented above , but presents nevertheless some interesting features for cosmology .
in the decoupling limit , taking the scales @xmath29 and @xmath3 , while keeping the strong coupling scale @xmath30 fixed , the above lagrangian reduces to the more compact expression l=-12h^e^_h_+h^x^(1)_+h^x^(2)_+h^x^(3)_+h^t_[eq : massivegravity ] where @xmath31 stands for the helicity-2 mode canonically normalized , @xmath32 is the lichnerowicz operator , the coefficients @xmath33 are related to the free parameters @xmath25 and @xmath34 denote the interactions with the helicity-0 mode @xcite x^(1 ) _ & = & g_- _ [ eq : x1 ] + x^(2 ) _ & = & ^2_-_-12([^2]-[]^2)g _
+ x^(3 ) _ & = & 6 ^ 3_-6[]^2_+3([]^2-[^2])_-g_([]^3 - 3[^2][]+2[^3]).[eq : x3 ] it is worth to mention that these interaction terms are all transverse and at most second order in time derivatives to ensure the absence of ghost .
being ghost - free , these interactions are closely related to the galileon interactions and fulfill the same internal symmetry . in the next section we will covariantize these interaction terms @xmath35 and discuss their physical properties . + in the case of external sources
there is a coupling between the metric @xmath31 and the stress energy tensor via @xmath36 but there is not such a direct coupling with the field @xmath7 .
nevertheless , if one diagonalizes the first interaction term @xmath37 by a change of variables of the form @xmath38 , then the coupling between @xmath7 and external sources will become transparent , @xmath39 .
instead of studying the cosmology in the exact non - linear covariantized theory which can be extremely hard , we use the alternative approach of covariantizing the lagrangian in the decoupling limit , and use the resulting theory as a proxy .
we claim that the covariantized version of the above lagrangian is simply given by [ eq : totalaction ] s=(^2 r+^(,g)+l^matter(,g ) ) , where @xmath40 is the lagrangian for the matter fields @xmath41 living on the geometry , the lagrangian for @xmath7 is @xcite [ eq : covjor ] ^=@xmath42 .
and the tensor @xmath43 stands for the dual riemann tensor l^=2r^+2(r^g^+r^g^-r^g^-r^g^ ) + r(g^g^-g^g^ ) .
this form of tensor structure has been first discussed by horndeski @xcite in the context of the most general scalar - tensor theory , and more recently in @xcite .
however we point out here that these interactions come as a direct outcome of massive gravity .
we can show explicitly the following correspondences , h^x^(1 ) _ & & -r + h^x^(2 ) _ & & -__g^ + h^x^(3 ) _ & & - _ _ _ l^. [ eq : corres ]    the einstein equation is given by [ eq.:einstein ] g_=t^_+t^matter _ with t^_=t^(1)_-t^(2)_-t^(3 ) _ and the structure of the einstein and riemann dual tensor ensure that @xmath7 enters at most with two derivatives in the stress - energy tensor , t^(1)_&=&x^(1)_+g _ [ eq : t1 ] + t^(2)_&=&x^(2)_+12l_^^+12g_()^2 [ eq : t2 ] + t^(3)_&=&x^(3)_+32l_^()^2 [ eq : t3 ] where @xmath44 are defined in ( [ eq : x1 ] - [ eq : x3 ] ) .
furthermore , the fact that @xmath45 , @xmath46 , @xmath47 and @xmath48 have at most one time - derivative guarantees the propagation of constraints .    since we are not in the einstein frame , these stress - energy tensor are only transverse on - shell , and satisfy the relation , @xmath49 where @xmath50 is the equation of motion with respect to @xmath7 . since both the einstein tensor and the riemann dual tensor are transverse , this equation of motion is also at most second order in derivative , _ & = & + & = & -r - g^-l^(__+r^___)=0,[eq : pij ] where we have used the fact that d_d_d_l^&=&-r^__l^=14 ^l_gb with @xmath51 . in the rest of this paper
, we study the resulting cosmology in this proxy theory , starting with the existence of self - accelerating solutions .
in what follows , we focus on the cosmology of the covariantized theory ( [ eq : totalaction ] , [ eq : covjor ] ) , and focus for that on a frw background with scale factor @xmath52 and hubble parameter @xmath53 . the resulting effective energy density and pressure for the field @xmath7 are then ^&=&(6h+6h^2-h^2 ^ 2-h^3 ^ 3)[eq : effdens ] + p^&=&3,[eq : effpress ] and the equation of motion for @xmath7 ( [ eq : pij ] ) in the frw space - time is equivalent to @xmath54+@xmath55=r .
similarly as in @xcite , this expression can be rewritten more compactly + 3h - r=0 if we define @xmath56 as = h^2@xmath57 .
now , we would like to study the self - acceleration solution with @xmath58const and @xmath59 .
for the @xmath7 field we make the ansatz @xmath60 .
furthermore , we assume that we are in a regime where @xmath61 so that we can neglect terms proportional to @xmath7 and consider only the terms including @xmath62 or @xmath63 .
thus , the friedmann and field equations can be recast in h^2=(6q-9a_2 q^2 - 30a_3 q^3)[eq : modfrid ] + h^2(18a_2 q+54a_3 q^2 - 12)=0 .
assuming @xmath64 , the field equation then imposes , q= while the friedmann equation ( [ eq : modfrid ] ) sets the hubble constant of the self - accelerated solution .
similar to what has been found in @xcite our proxy theory admits a self - accelerated solution , with the hubble parameter set by the graviton mass . for the stability condition of this self - accelerating solution
the first constraint we have is to demand @xmath65 .
the other constrain comes from the stability condition for perturbations on the background which we discuss in the following subsection .      in the last subsection
we have shown explicitly that our proxy theory exhibits a self - accelerating solution with @xmath66 .
now , we would like to study whether the perturbations on this background are stable and what the constraints are .
for this purpose , consider perturbations on the background solution of the following form = _ 0(t)+(t , x , y , z)[eq : perturbations ] the second order action for the perturbations is l=-__g^-__^(0)_l^-__0_d_d_l^[eq : perlag ] , which can be written in the form l = k_tt(^2-()^2 ) where k_tt=-(a_2 + ) and c_s^2=13 ( 2 + ) .
the condition for the stability is then given by @xmath67 , @xmath68 and @xmath69 , which are fulfilled if [ eq : stability ] a_2>0 0>a_3>-18 a_2 ^ 2 . to compare this result with the condition obtained in the decoupling limit , @xcite
, we first mention that @xmath70 and @xmath71 . in terms of the parameters used there ,
we need to compare our conditions @xmath72 , @xmath73 to the conditions @xmath74 and @xmath75 .
we see that our theory is less constraining but is still within the parameter space derived in @xcite .
it is not surprising that the stability condition in the decoupling limit and in our covariantized theory do not coincide totally as we have explicitly broken the symmetry when getting the proxy and in particular , our solution does now spontaneously break lorentz invariance , which was not the case in the decoupling limit , @xcite .    it is also worth pointing out that the self - accelerating solution by itself does not propagate any superluminal mode , since @xmath76 .
we emphasize as well that the constant @xmath77 solution is a dynamical attractor .
for this we just consider time dependent perturbations @xmath78 which is a special case of fulfilling the same stability conditions .
the equation of motion for perturbations simplifies to _ t(a^3 ) = 0 the solution for @xmath79 is given by ( t)~a^-3 .
thus , these perturbations @xmath80 redshift away exponentially compared to the @xmath81 self - accelerating solution .
therefore , the self - accelerating solution is an attractor .
more interestingly , one can wonder whether degravitation can be exhibited in these class of solutions .
if one take @xmath82 and @xmath83 , it is straightforward to see that we obtain @xmath84 , so the field has absolutely no effect and can not help the background to degravitate .
the situation is however different when setting @xmath85 , which was possible in the decoupling limit . in this proxy theory
, such a behaviour will involve explicit space dependances in the equations of motion which should be dealt with specifically .
interestingly , the interactions considered here are precisely of the same form as that studied recently in @xcite . there as well , in the absence of spatial curvature @xmath86 , the contribution from the scalar field vanishes if @xmath83 . comparing with @xcite
, we can hence wonder whether the addition of spatial curvature @xmath87 in our proxy theory could help achieving degravitation , but relying strongly on spatial curvature brings concerns over instabilities which are beyond the scope of this study .
in the following we would like to discuss in more detail the interplay of all the constituents of the universe .
we assume that matter , radiation and the scalar field @xmath7 contribute to the total energy density of the universe .
+ ^mat ) consider the scalar field @xmath7 as a perfect fluid with the effective energy density and pressure given by ( [ eq : effdens ] , [ eq : effpress ] ) .
thus , the equation of state parameter of this new field would be _ & = & + & & + . at this point
one should mention that the energy density for the @xmath7-field is not conserved but rather given by @xmath49 ( where @xmath50 is the equation of motion for @xmath7 ) , which is not surprising since @xmath7 is non - minimally coupled to gravity in the jordan frame .
therefore , we can have @xmath88 .
+ in the following we will first assume that at early times in the evolution history of the universe we can neglect the extra density coming from the helicity-0 @xmath89 .
we will then check this assumption by plugging the solution for @xmath53 back in the equation of motion for @xmath7 .
if we assume that at early times the radiation density dominates , we simply have h^2=_0^rada^-4 a~t^1/2=1/3[eq : raddom ] during the radiation era , the dominant terms in the equation of motion for @xmath7 are then @xmath90 which can be solved assuming the previous expression for @xmath53 ( [ eq : raddom ] ) _
rad~t^1.75 ^_rad~t^-1/4    , @xmath91 and @xmath89 during the epochs of radiation , matter and @xmath0-domination normalised to today @xmath89 . during the radiation domination
the energy density for @xmath7 goes as @xmath92 and during matter dominations as @xmath93 and is constant for later times @xmath94.,scaledwidth=80.0% ]     + at later times when the matter dominated epoch starts we have h^2=_0^mata^-3 a~t^2/3=0
now the dominant terms in the equation of motion for @xmath7 are @xmath95 .
we get for @xmath7 this time ^mat~c_2t+ ^_mat = c_2 t^-1 + .
summarizing , during radiation domination the effective energy density for the @xmath7-field goes like @xmath96 while during matter domination as @xmath97 and approaches a constant at late time . as shown in the figure @xmath98 can be neglected at early times where @xmath99 and @xmath100 .
we end the cosmological analysis by looking at the evolution of matter density perturbations .
the density perturbations follow the evolution _
m+2h_m= where @xmath41 is the newtonian potential .
the effects of @xmath7 are all encoded in its contribution to the poisson equation .
consider perturbations of @xmath7 around its cosmological background solution @xmath101 .
@xmath102 gives a contribution to the newtonian potential of the form @xmath103 . in the newtonian approximation
we have @xmath104 .
the equation of motion for the scalar field in first order in @xmath102 is -g^ d_d_- l_(4_0^ d^d^+2r^__0^)=r which is equivalent to ( neglecting @xmath105 ) = r and last but not least we need the trace of einstein equation , ( [ eq.:einstein]-[eq : t3 ] ) . perturbing the trace to first order
, we get -r&= & + & & + to reach that point , we have neglected the perturbations of the curvature of the form @xmath106 as they are negligible compared to @xmath107 since we work in the regime where @xmath108 .
we have also ignored terms of the form @xmath109 , which could a priori be relevant , but their inclusion would require solving the full einstein equation , which is beyond the scope of this study .
as a first approximation , such terms are hence ignored .
the perturbations for the source is just given by @xmath110 for non - relativistic sources , thus we have = where @xmath111 stands for & & q1-(2 h _
0+_0+(2 h + 3 h^2 ) ) + & & -@xmath112 .
finally , the modified evolution equation for density perturbations is _ m+2h_m=(1 + 13q ) . knowing the background configuration it is then relatively straightforward to derive the effect on structure formation .
we recover the usual result that when the field is screened , @xmath113 , the extra force coming from the helicity-0 is negligible and the formation of structure is similar as in @xmath0cdm .
after having studied the cosmology and the structure formation of our proxy theory in the jordan frame , the natural question is whether we would expect similar results in a different frame . instead of covariantizing our lagrangian in the jordan frame ,
it is on an equal footing to go to the einstein frame where the ricci scalar is not multiplied by the scalar field @xmath7 and covariantize the theory at that stage .
since it is unclear at first sight which frame is the physical frame we will consider both frames and discuss their differences .
our starting lagrangian was l=-12h^e^_h_+h^x^(1)_+h^x^(2)_+h^x^(3)_+h^t _ now , when we do the following change of variables h_=|h_+_[eq : ef1 ] we can diagonalize the first mixed term @xmath37 such that the lagrangian takes the following form l=-12|h^e^_|h_+32+|h^x^(2)_-32()^2 + + |h^x^(3 ) _
-2()^2([^2]-^2 ) + ( |h_+_)t^ covariantizing this action is straightforward .
we use again the correspondences in and it has been shown explicitly that the covariant equivalence to @xmath114 and @xmath115-\box\pi^2)$ ] are given by @xmath114 and @xmath116-\frac14(\partial\pi)^2r ) $ ] respectively , which do not yield any ghostlike instabilities ( @xcite , @xcite ) . thus , the covariantized action in the einstein frame is simply given by & = & ^2 r+32-__g^-32()^2 + & & - _ _ _ l^ + 2()^2(^2-[^2]-14()^2r ) [ eq : einstein ] + & & + l_m ( , ( 1+)g_)similarly as before , the properties of @xmath117 and @xmath43 ensure that their equations of motion lead at most to second order derivative terms . to find a self - accelerating solution we set again a pure de sitter metric , with @xmath118 . the friedmann and the field equations
then take the form & & ^2h^2=3a_2 ^ 3q^2 + 10a_3 ^ 3q^3 - 12q^2 + 3a_2q^3 + 15a_3q^4 + & & ^3(-1 + 3q(a_2 + 4a_3q))+2(a_2 + 3a_3q)h^2=0[eq : frein ] when comparing the above friedmann and the field equations with the one we had in the jordan frame , we see significant differences coming from the extra terms which were not there in the jordan frame .
these terms yield friedmann and field equations proportional to @xmath119 and @xmath120 which are more difficult to solve .
+ for fairness , we should compare both actions in the same frame .
we do so by performing a conformal transformation on the action : g_=^2 g _ ^2=(1- ) . for simplicity
we consider the case for which @xmath121 , so under this conformal transformation the covariantized action becomes _ j= ^2r-32 ^ -4()^2-(^^g_+32()^2 + 54()^4 ) in the limit where @xmath122 we have then finally the following expression _
j= ^2 r+32-__g^-32()^2 , which coincides with the theory obtained from the einstein frame , with @xmath121 .
so within the regime of validity of our results , our conclusions are independent of the choice of frame .
however beyond this regime the theory originally constructed from the jordan frame could violate the null energy condition from the term proportional to @xmath123 .
common to generic massive gravity models is the presence of at least one additional helicity-0 degree of freedom which originates from a spin - two graviton field .
in contrast to an arbitrary scalar field , its dynamics is not driven by a potential but rather by specific derivative terms fixed by symmetries .
such a degree of freedom is experimentally testable . on solar system and galactic scales gravity
is very well compatible with gr and correctional terms from a modified gravity model are excluded . on these scales , however , the effects of massive gravity might be cloaked by the vainshtein mechanism , where the helicity-0 mode interactions become appreciable to freeze out the field fluctuations .
a number of screening mechanisms have been devised such as the chameleon and symmetron mechanisms .
contrary to those the success of the vainshtein mechanism does not rely on some specific potential but instead on derivative interactions , which cause the helicity-0 mode to become decoupled from matter and light on short distances , which nevertheless could have observational signatures on larger scales in cosmic structure formation .
a further benefit of the vainshtein mechanism would be a natural solution to the problem that the matter density and the cosmological constant are of similar magnitudes today .
the cosmological application of the vainshtein mechanism works such that at early times , the galileon interactions are dominated by self - interaction , which suppresses their energy density relative to that of matter or radiation . if the matter density has decreased sufficiently by cosmic expansion , the galileon constitutes an important contribution to the energy density of the universe and drives cosmic expansion .    although the precise cosmology of the theory of massive gravity proposed in @xcite has not been derived , most models of massive gravity give rise to alterations of the hubble function similar to dgp - gravity .
experimental tests of the expansion dynamics of the universe include the distance - redshift relation of supernovae , and measurements of the angular diameter distance as a function of redshift , as in the case of the cosmic microwave background and the baryon acoustic oscillations .
apart from these geometrical tests , the time - dependence of evolving cosmic structures can be investigated , and the influence of the gravitational theory on the geodesics of relativistic ( photons ) and nonrelativistic ( dark matter ) test particles .
the first category includes gravitational lensing and the sachs - wolfe effects , which have been shown to differ from their gr - expectation in some modified gravity theories , and a similar result can be expected in massive gravity .
the second category includes the homogeneous growth of the cosmic structure , and the formation of galaxies and clusters of galaxies by gravitational collapse .
again , massive gravity would influence the time sequence of gravitational clustering and the evolution of peculiar velocities , as well as the number density of collapsed objects .
in particular , massive gravity would enhance gravitational clustering since they tend to lower the collapse threshold for density fluctuations in the large - scale structure , leading to a higher comoving number density of galaxies and clusters of galaxies .
naturally , these changes are degenerate with a different choice of cosmological parameters and with introducing non - gaussian initial conditions , which would be very interesting to quantify .
recent discrepancies of @xmath0cdm with observational data on large scales include the number of very massive clusters , the strong lensing cross section , anomalous multipole moments of the cmb , the axis of evil , and large peculiar velocities .
it is beyond the scope of this paper to address these issues using massive gravity , but we propose how to proceed from constructing a massive gravity theory to providing observationally testable quantities . from our point of view
it is advisable to focus on probes of large scales , due to difficulties related to nonlinear structure formation and the influence of baryons on small scales .
natural questions concern the homogeneous dynamics of the universe , the formation of structures and the shape of geodesics of relativistic and nonrelativistic particles .
basically , using our proxy theory one should be able to make predictions concerning these four issues and the combination of the four should give insight into the nature of the gravitational sector .
in our proxy theory , from the modified field equation the hubble - function can be derived easily , which allows the definition of distance measures , needed in the interpretation of supernova data .
cosmic structure growth tests the newtonian limit for slowly - moving particles and describes the clustering of galaxies and the growth of structures which are investigated by e.g. gravitational light deflection .
lensing , in turn , makes use of the geodesic equation for relativistic particles , and measures the correlation function of the matter density , weighted with the lensing efficiency function , which in turn is derived from distance measures .
it would be quite interesting to study these observational consequences and constraints of our proxy theory in a future work .
we would like to thank clare burrage , paul de fromont , gregory gabadadze , parvin moyassari and andrew tolley for useful discussion .
furthermore , special thanks for robert caldwell ! this work is supported by the snf .
g.  r.  dvali , g.  gabadadze and m.  porrati , phys .
b * 485 * , 208 ( 2000 ) [ arxiv : hep - th/0005016 ] ; g.  r.  dvali and g.  gabadadze , phys .
d * 63 * , 065007 ( 2001 ) [ arxiv : hep - th/0008054 ] .
a.  i.  vainshtein , phys .
b * 39 * ( 1972 ) 393 .    c.  deffayet , g.  r.  dvali and g.  gabadadze , phys .
d * 65 * , 044023 ( 2002 ) [ arxiv : astro - ph/0105068 ] .
k.  koyama , phys .
d * 72 * , 123511 ( 2005 ) [ arxiv : hep - th/0503191 ] . c.  charmousis , r.  gregory , n.  kaloper , a.  padilla , jhep * 0610 * , 066 ( 2006 ) .
[ hep - th/0604086 ] .
a.  nicolis , r.  rattazzi and e.  trincherini , phys .
d * 79 * , 064036 ( 2009 ) [ arxiv:0811.2197 [ hep - th ] ] . c.  de rham and a.  j.  tolley , jcap * 1005 * , 015 ( 2010 ) [ arxiv:1003.5917 [ hep - th ] ] .
f.  p.  silva and k.  koyama , phys .
d * 80 * , 121301 ( 2009 ) [ arxiv:0909.4538 [ astro-ph.co ] ] .
n.  chow and j.  khoury , phys .
d * 80 * , 024037 ( 2009 ) [ arxiv:0905.1325 [ hep - th ] ] .
j.  khoury , m.  wyman , phys .
rev .   * d80 * , 064023 ( 2009 ) .
[ arxiv:0903.1292 [ astro-ph.co ] ] ; t.  kobayashi , h.  tashiro and d.  suzuki , phys .
d * 81 * , 063513 ( 2010 ) [ arxiv:0912.4641 [ astro-ph.co ] ] ; t.  kobayashi , arxiv:1003.3281 [ astro-ph.co ] ; m.  wyman , j.  khoury , phys .  rev .   * d82 * , 044032 ( 2010 ) . [ arxiv:1004.2046 [ astro-ph.co ] ] ; r.  gannouji and m.  sami , phys .
d * 82 * , 024011 ( 2010 ) [ arxiv:1004.2808 [ gr - qc ] ] ; a.  de felice , s.  mukohyama and s.  tsujikawa , phys .
d * 82 * , 023524 ( 2010 ) [ arxiv:1006.0281 [ astro-ph.co ] ] .
a.  ali , r.  gannouji and m.  sami , phys .
d * 82 * , 103015 ( 2010 ) [ arxiv:1008.1588 [ astro-ph.co ] ] ; a.  de felice and s.  tsujikawa , arxiv:1008.4236 [ hep - th ] ; a.  de felice , r.  kase and s.  tsujikawa , phys .
d * 83 * , 043515 ( 2011 ) [ arxiv:1011.6132 [ astro-ph.co ] ] ; k.  hirano , z.  komiya , h.  shirai , [ arxiv:1103.6133 [ astro-ph.co ] ] ; k.  hirano , [ arxiv:1012.5451 [ astro-ph.co ] ] ; p.  brax , c.  burrage , a.  -c .
davis , [ arxiv:1106.1573 [ hep - ph ] ] ; x.  gao , [ arxiv:1106.0292 [ astro-ph.co ] ] .
p.  creminelli , a.  nicolis and e.  trincherini , jcap * 1011 * , 021 ( 2010 ) [ arxiv:1007.0027 [ hep - th ] ] ; l.  p.
levasseur , r.  brandenberger , a.  -c .
davis , [ arxiv:1105.5649 [ astro-ph.co ] ] ; t.  kobayashi , m.  yamaguchi and j.  yokoyama , phys .  rev
.  lett .   * 105 * , 231302 ( 2010 ) [ arxiv:1008.0603 [ hep - th ] ] ; t.  kobayashi , m.  yamaguchi and j.  yokoyama , arxiv:1103.1740 [ hep - th ] .    c.  burrage , c.  de rham , d.  seery and a.  j.  tolley , jcap * 1101 * , 014 ( 2011 ) [ arxiv:1009.2497 [ hep - th ] ] ; p.  creminelli , g.  damico , m.  musso , j.  norena and e.  trincherini , jcap * 1102 * , 006 ( 2011 ) [ arxiv:1011.3004 [ hep - th ] ] . s.  mizuno , k.  koyama , phys .  rev .
* d82 * , 103518 ( 2010 ) .
[ arxiv:1009.0677 [ hep - th ] ] ; s.  renaux - petel , [ arxiv:1105.6366 [ astro-ph.co ] ] .
e.  babichev , c.  deffayet and r.  ziour , int .
j.  mod .
d * 18 * , 2147 ( 2009 ) [ arxiv:0905.2943 [ hep - th ] ] ; c.  deffayet , o.  pujolas , i.  sawicki and a.  vikman , jcap * 1010 * , 026 ( 2010 ) [ arxiv:1008.0048 [ hep - th ] ] ; r.  kimura and k.  yamamoto , arxiv:1011.2006 [ astro-ph.co ] .    c.  deffayet , g.  esposito - farese and a.  vikman , phys .
d * 79 * , 084003 ( 2009 ) [ arxiv:0901.1314 [ hep - th ] ] ; c.  deffayet , s.  deser and g.  esposito - farese , phys .
d * 80 * , 064015 ( 2009 ) [ arxiv:0906.1967 [ gr - qc ] ] ; c.  deffayet , x.  gao , d.  a.  steer , g.  zahariade , [ arxiv:1103.3260 [ hep - th ] ] .
a.  de felice and s.  tsujikawa , phys .
lett .   * 105 * , 111301 ( 2010 ) [ arxiv:1007.2700 [ astro-ph.co ] ] ; s.  nesseris , a.  de felice and s.  tsujikawa , phys .
d * 82 * , 124054 ( 2010 ) [ arxiv:1010.0407 [ astro-ph.co ] ] .
g.  goon , k.  hinterbichler , m.  trodden , [ arxiv:1103.5745 [ hep - th ] ] ; g.  goon , k.  hinterbichler , m.  trodden , [ arxiv:1103.6029 [ hep - th ] ] ; c.  burrage , c.  de rham , l.  heisenberg , jcap * 1105 * , 025 ( 2011 ) .
[ arxiv:1104.0155 [ hep - th ] ] .
g.  gabadadze , phys .
b * 681 * , 89 ( 2009 ) [ arxiv:0908.1112 [ hep - th ] ] ; c.  de rham , phys .  lett .
b * 688 * , 137 ( 2010 ) [ arxiv:0910.5474 [ hep - th ] ] . c.  de rham and g.  gabadadze , phys .
b * 693 * ( 2010 ) 334 [ arxiv:1006.4367 [ hep - th ] ] .
s.  f.  hassan , r.  a.  rosen , [ arxiv:1104.1373 [ hep - th ] ] .
l.  berezhiani , m.  mirbabayi , [ arxiv:1104.5279 [ hep - th ] ] .
g.  dvali , g.  gabadadze and m.  shifman , _ published in * minneapolis 2002 , continuous advances in qcd * 566 - 581 _ ; arxiv : hep - th/0208096 ; g.  dvali , g.  gabadadze and m.  shifman , phys .
044020 ( 2003 ) [ arxiv : hep - th/0202174 ] . n.  arkani - hamed , s.  dimopoulos , g.  dvali and g.  gabadadze , arxiv : hep - th/0209227 . g.  dvali , s.  hofmann and j.  khoury , phys .
d * 76 * , 084006 ( 2007 ) [ arxiv : hep - th/0703027 ] .
g.  gabadadze and m.  shifman , phys .
d * 69 * , 124032 ( 2004 ) [ arxiv : hep - th/0312289 ] . c.  de rham , g.  dvali , s.  hofmann , j.  khoury , o.  pujolas , m.  redi and a.  j.  tolley , phys .
lett .   * 100 * , 251603 ( 2008 ) [ arxiv:0711.2072 [ hep - th ] ] ; c.  de rham , s.  hofmann , j.  khoury and a.  j.  tolley , jcap * 0802 * , 011 ( 2008 ) [ arxiv:0712.2821 [ hep - th ] ] ; c.  de rham , j.  khoury and a.  j.  tolley , phys .
lett .   * 103 * , 161601 ( 2009 ) [ arxiv:0907.0473 [ hep - th ] ] ; c.  de rham , j.  khoury and a.  j.  tolley , arxiv:1002.1075 [ hep - th ] .