inflation @xcite is the leading paradigm for the very early universe cosmology .
inflation has been proposed to explain the horizon , flatness and monopole problems in the standard hot big bang cosmology , and almost all the predictions of the simplest inflation model have now been tested .
the observational tests of inflation includes    * coherent and nearly scale invariant power spectrum of density perturbations .
the power spectrum of the simplest slow roll inflation is @xcite @xmath0 * a small tilt of the scalar power spectrum .
@xmath1 where @xmath2 is the slow roll parameter defined from expansion .
now @xmath3 is ruled out under the assumptions of simplest inflation models . *
nearly gaussian density fluctuations .
the non - gaussianities of the density fluctuations are tightly constrained at @xmath4 for the local shape and equilateral shape non - gaussianities respectively .
those numbers indicate that , non - gaussian components of the primordial fluctuations , even if exist , have to be at least 3@xmath54 orders of magnitudes smaller than the gaussian component . * gravitational waves .
the recent bicep2 experiment reports an over @xmath6 detection of gravitational waves @xcite , with tensor to scalar ratio @xmath7 this corresponds to a gravitational wave fluctuation amplitude @xmath8    despite of the great success of inflation , there are still a few outstanding challenges for theorists and experimentalists .    on the theoretical side , large field inflation is now favored .
however , large field inflation is hard to construct from the effective field theory , and stringy uv completion points of view .
the uv completion of inflation has long suffered from an @xmath9-problem @xcite , in which the mass of the inflaton is theoretically too large to allow enough e - folds of inflation .
however , with the current data , a more serious @xmath10-problem emerges  the observed energy scale of inflation is too high for an effective field theory or stringy model building to be under control . for single field
slow roll inflation , at every e - fold , the inflaton rolls a distance of order @xmath11 . in perturbative string theory , this field motion per e - fold is comparable with , or greater than the string scale @xmath12 . as a result , one can no longer safely globally expand the inflaton field and ignore non - renormalizable terms .
more discussions and a local reconstruction of the inflationary potential can be found in @xcite    on the observational side , there is yet another ( and maybe the last in the foreseeable future , unless nature is so kind as to imprint other relics on the cmb sky or in the large scale structure ) test for inflation which is possible in light of bicep2 , but not yet achieved
the tilt of the tensor power spectrum .
the simplest inflation models predict a consistency relation between @xmath13 and @xmath14 as @xmath15 currently the data has not been good enough to test @xmath13 precisely .
however , there are a lot of ongoing and upcoming experiments in the near future @xcite , measuring @xmath14 at different scales , which provides a possibility for a precise measurement of @xmath13 .    in this paper
, we shall explore the possibility of blue @xmath13 . in section
[ sec : hints - blue - tensor ] , the bound on @xmath13 is derived the bicep2 @xcite and polarbear @xcite data .
the string gas cosmology , null energy condition ( nec ) violating inflation , and general initial condition are addressed in sections [ sec : string - gas - cosmology ] , [ sec : infl - viol - nec ] and [ sec : infl - gener - init ] respectively . in section [ sec : other - possibilities ] , a few other possibilities are discussed , including external sources for tensor modes , modified gravity and matter bounce .
we conclude in section [ sec : concl - disc ] .
in this section , we show the constraints on @xmath14 and @xmath13 from bicep2 , _ planck _ and polarbear data .
the public codes cosmomc @xcite and class @xcite are used in the calculation . here
we vary @xmath14 and @xmath13 in the boltzmann codes , and other cosmological parameters @xmath16 are set at the best fit value from planck results @xcite where not explicitly mentioned .
the relation between @xmath14 and @xmath13 is illustrated in fig .
[ fig : rnt_sample ] . at @xmath17 ,
the flat tensor spectrum gives a good fit to the data , which is consistent with the result of the bicep2 @xcite . with the fixed ratio @xmath14 at pivot scale @xmath18 , by increasing or decreasing the tilt @xmath13 by @xmath19 , it starts to be ruled out by the current data . if we consider smaller tensor - to - scalar ratio , indicated from the right panel of fig .
[ fig : rnt_sample ] , a blue tensor spectrum is required by the bicep2 data .
one can observe from the right panel of fig .
[ fig : rnt_sample ] that very blue @xmath13 fits all data points within @xmath20 .
this is different from the @xmath21 case , where some data points are too high and may be considered as outliers .
a more precise correlation of @xmath13 and @xmath14 is illustrated in fig .
[ fig : rnt_contour ] , where the pivot scale at @xmath22 is the same as the analysis of bicep2 paper  @xcite .     power spectrum with different @xmath23 and @xmath13 .
left panel : @xmath23 is fixed and variation of @xmath13 is illustrated .
right panel : some sets of @xmath23 and @xmath13 are chosen to get good fit against data . in both panels @xmath14 is calculated at pivot scale @xmath24 mpc@xmath25 . ]
power spectrum with different @xmath23 and @xmath13 . left panel : @xmath23 is fixed and variation of @xmath13 is illustrated .
right panel : some sets of @xmath23 and @xmath13 are chosen to get good fit against data . in both panels @xmath14 is calculated at pivot scale @xmath24 mpc@xmath25 . ]     vs @xmath13 to fit against the data .
left panel : fit the data of the bicep2 and polarbear .
right panel : fit the first five data bins of bicep2 and treat the others as upper bound .
note that the strong correlation between @xmath26 and @xmath23 is an artifact from the choice of pivot scale .
we have chosen @xmath27 to match with bicep2 conventions . however , the bicep2 experiments measures @xmath28 .
thus given a similar tensor spectrum at @xmath28 , the large tilt modifies the @xmath23 value significantly . ]
vs @xmath13 to fit against the data .
left panel : fit the data of the bicep2 and polarbear .
right panel : fit the first five data bins of bicep2 and treat the others as upper bound .
note that the strong correlation between @xmath26 and @xmath23 is an artifact from the choice of pivot scale .
we have chosen @xmath27 to match with bicep2 conventions . however , the bicep2 experiments measures @xmath28 .
thus given a similar tensor spectrum at @xmath28 , the large tilt modifies the @xmath23 value significantly . ]    however , as one can observe from fig .
[ fig : rnt_contour ] , @xmath14 and @xmath13 are strongly correlated when the pivot scale is chosen at @xmath22 .
this makes detailed analysis inconvenient .
the strong correlation comes from the fact that the bicep2 experiment is actually not observing the b - mode signals at their pivot scale @xmath22 , but instead the observation is made at about @xmath29 .
this does not matter for the bicep2 data analysis itself because @xmath21 is chosen in their data analysis .
however , in our case , when @xmath13 is allowed as a free parameter , we had better pay more attention to the pivot scale . in the remainder of this section , we shall choose the pivot scale at @xmath29 and calculate @xmath14 at this pivot scale .
in table  [ tab : chi2 ] , the best fit parameters and their likelihood are summarized .
one can observe that from the bicep2 data only , including @xmath13 improves @xmath30 by 3.161 .
considering that one additional parameter is induced , it is useful compute the akaike information criterion ( aic ) or the bayesian information criterion ( bic ) .
the bicep2 e - mode polarization , the wmap polarization and the high @xmath31 part of the planck temperature data has a large amount of data points which are not sensitive to the change of @xmath13 .
thus we choose aic to compare between models . ] to see how much improvements @xmath13 would bring .
aic can be calculated as @xmath32 where @xmath33 is the number of model parameters . thus including @xmath13 results in @xmath34 .
the number not very significant , though showing hint that blue @xmath13 would fit data slightly better . to explore the goodness - of - the - fit
, we can also look at the relative @xmath30 and p - value of the bicep2 bb correlation .
including bicep2 bb correlation and polarbear , the relative @xmath35 , p - value is @xmath36 in the case of @xmath13 as a free parameter and @xmath37 , p - value is @xmath38 for @xmath21 .
both of @xmath39 are close to @xmath40 meaning both of them are good fit , but the one with the parameter @xmath13 is slightly better .
the improvement mainly comes from the small scale excess of the b - mode power spectrum .
keck ( preliminary ) cross correlation @xcite .
however , in the bicep2 @xmath41 keck ( preliminary ) cross correlation , the large scale data points are lower than the theoretical prediction , again showing some evidence of blue @xmath13 . ]
however , the significance of blue @xmath13 increase dramatically once planck data is considered . after allowing to vary @xmath13 , @xmath30 is improved by @xmath42 for this single parameter .
namely , the fitting of @xmath43 part of planck temperature data is improved by @xmath44 .
the fitting of @xmath45 part of planck temperature data is improved by @xmath46 .
the fitting of bicep2 data is improved by 4.945 .
on the other hand , blue @xmath13 actually fit slightly worse for wmap polarization , by @xmath47 . as a result
, one can calculate that @xmath48 , showing strong evidence of blue @xmath13 .
it is straightforward to understand why inclusion of planck data strongly prefers blue @xmath13 .
note that among the likelihoods listed in table  [ tab : chi2 ] , the most significant improvement comes from the @xmath43 part of planck temperature data .
this is because , there is already an anomalous temperature power deficit in the small @xmath31 part of the planck data , by an amount of @xmath49 , at the statistical significance of @xmath50 .
on the other hand , the existence of tensor mode , when assuming @xmath21 , adds to the theoretically predicted temperature power by another @xmath49 , which fits data even worse .    by inclusion of @xmath13 , the above tension between planck and bicep is reconciled .
as one can compare to the last column of table  [ tab : chi2 ] , after adding blue @xmath13 , the likelihood of planck low @xmath31 , planck high @xmath31 and wmap polarization are all comparable with the best fit values before bicep2 .
thus the tension between planck and bicep2 can be reconciled by blue @xmath13 .    to further get constraints of @xmath14 and @xmath13 , we run mcmc chains to get samples according to the likelihood code of bicep2 , planck and wmap polarization .
the @xmath51 and @xmath20 contours for bicep2 only ( 9 bins ) is plotted in fig .
[ fig : b9plot ] , and bicep2+planck+wp is plotted in fig .
[ fig : b9pwplot ] .
again one observes that blue @xmath13 is preferred .
the mean value and standard derivation of @xmath14 and @xmath13 from bicep2 ( 9 bins ) is @xmath52 the constraints from bicep2 ( 9 bins ) + planck ( 2013 ) + wmap polarization is @xmath53    on the other hand , it is worth noticing that if one only take the first 5 bins from bicep2 data , there is no evidence of blue @xmath13 .
this is plotted in fig .
[ fig : b5plot ] for comparison purpose . because as we have discussed , the two major sources of blue @xmath13 are the observed small scale excess of the b - mode power spectrum of bicep2 , and the observed large scale suppression of the temperature spectrum of planck .
none of those are present in the first 5 bins of bicep2 data .
the constraints from those 5 bins are @xmath54    .[tab : chi2 ] the best fit values for @xmath14 and @xmath13 , likelihood and aic .
the 5 columns corresponds to : [ column 1 ] bicep2 only ( 9 bins ) with varying @xmath14 only , [ column 2 ] bicep2 only ( 9 bins ) with varying @xmath14 and @xmath13 , [ column 3 ] bicep2 ( 9 bins ) + planck ( 2013 ) + wmap polarization with varying @xmath14 only , [ column 4 ] bicep2 ( 9 bins ) + planck ( 2013 ) + wmap polarization with varying @xmath14 and @xmath13 , and [ column 5 ] planck ( 2013 ) + wmap polarization with @xmath55 ( for comparison purpose ) , respectively . note that in the first two columns ( bicep2 data only ) the other cosmological parameters @xmath16 are fixed because bicep2 alone does not provide reasonable constraints on those parameters . in the minimization of the last two columns ( combining with planck and wp ) those additional cosmological parameters and nuisance parameters are also varying . [
cols="^,^,^,^,^,^",options="header " , ]      the @xmath56 contour ( left panel ) , and the likelihood for @xmath14 ( middle panel ) and @xmath13 ( right panel ) from bicep2 ( 9 bins ) . ]     the @xmath56 contour ( left panel ) , and the likelihood for @xmath14 ( middle panel ) and @xmath13 ( right panel ) from bicep2 ( 9 bins ) + planck ( 2013 ) + wmap polarization . ]     the @xmath56 contour ( left panel ) , and the likelihood for @xmath14 ( middle panel ) and @xmath13 ( right panel ) from bicep2 ( 5 bins ) . ]    note that when @xmath57 , the primordial tensor contribution to the b - mode polarization becomes more important than lensing signal even at small scales ( @xmath58 ) . in this case , the polarbear detection of  lensing b - mode " becomes a constraint on blue @xmath13 .
the @xmath51 and @xmath20 exclusion curves from polarbear is plotted in fig .
[ fig : combined ] . on the other hand ,
the cross - correlation - detection of b - mode from sptpol does not put constraint on @xmath13 , because the primordial b - mode does not cross correlate with the density perturbation in the simplest inflation models .
another rather theoretical constraint for blue @xmath13 is that , very blue tensor tilt can not last long . assuming the running of @xmath13 is not significant , then for @xmath59
, it takes about 23 e - folds to bring the tensor mode to be non - perturbative ( @xmath60 ) .
those non - perturbative tensor modes forms primordial black holes ( pbh ) , the induced scalar perturbation also becomes of order one , which forms primordial black holes .
] , which is constrained from current observations . for @xmath61
, it takes about 12 e - folds . thus assuming @xmath62 we get @xmath63 for 60 e - folds of inflation and
@xmath64 for 50 e - folds of inflation .
those constraints are plotted in fig .
[ fig : combined ] . nevertheless ,
if @xmath13 is large , @xmath65 may be large as well . if @xmath65 is negative @xcite , it relaxes the above constraint .    finally , although the current bound on @xmath13 is far from testing the consistency relation of inflation , it is already informative to disfavor scenarios with very blue or red spectrum .
for example , a sharp pulse of gravitational waves with rapid decaying tail towards both ends ( @xmath66 or @xmath67 ) are disfavored by the current data .
combined constraints on @xmath14 and @xmath13 .
the blue contours are the @xmath51 and @xmath20 constraints from bicep2 ( 9 bins ) .
the red contours are the constraints from bicep2 + planck ( 2013 ) + wmap polarization .
the blue shaded region are the @xmath51 and @xmath20 exclusion curves from polarbear .
the gray shaded region is the exclusion curves from pbh constraints , assuming the blue @xmath13 is also present in smaller scales than cmb observations.,scaledwidth=60.0% ]
before a survey of inflationary possibilities for the blue tensor spectra , let us mention alternatives to inflation .
this is because the first prediction of a slightly blue @xmath13 comes from the string gas cosmology @xcite ( see @xcite for a recent discussion in light of bicep2 data ) .. this is consistent with our @xmath13 bound if all the bicep2 data points are fitted .
however , the amplitude of the tensor mode is exponentially suppressed on cosmological scales , which can not be @xmath68 at @xmath69 .
there is a similar issue for the bouncing galileons @xcite .
also , there is recently another model with blue tensor tilt derived from the hagedorn phase @xcite .
] on the contrary , for the other models that we shall show below , though they have the possibility to tune the tensor spectrum to be blue , the blueness is not a firm prediction .    in string gas cosmology , the universe was in a string hagedorn phase before expansion starts .
it is conjectured that due to t - duality , and the huge specific heat in the string hagedorn phase , the universe should stay at a nearly constant temperature ( hagedorn temperature @xmath70 ) for a long time , until the string winding modes decays and allow the expansion of the universe ( which is also a possible explanation of three large spatial dimensions ) .
string gas cosmology .
the universe starts from a string hagedorn phase and a phase transition brings the universe into radiation dominated.,scaledwidth=80.0% ]    the density and tensor fluctuations in string gas cosmology are thermal fluctuations .
those fluctuations were in causal contact in the string hagedorn phase , and get frozen on super - hubble scales because of the rapid shrinking of the hubble horizon ( fig .
[ fig : sgc ] ) .
the scalar and tensor fluctuations can be calculated as @xcite @xmath71 where @xmath72 is the newtonian potential .
the scalar spectral index is @xmath73 from the string gas picture , one expects that @xmath74 is an decreasing function of @xmath33 and thus string gas cosmology has red scalar spectrum .
the tensor power spectrum is @xmath75~.\end{aligned}\ ] ] the tensor to scalar ratio is @xmath76~.\end{aligned}\ ] ] on the other hand , the tensor to scalar ratio measured by bicep2 is between @xmath77 and the comoving curvature perturbation @xmath78 .
this results in a conversion factor . on super - hubble scales , @xmath79 where @xmath80 is
the comoving hubble parameter and prime denotes derivative with respect to the comoving time . assuming the hagedorn phase is fast enough to the radiation dominated era , the factor can be calculated as @xmath81 thus the tensor to scalar ratio using @xmath78 is @xmath82~.\end{aligned}\ ] ] inserting the cobe normalization @xmath83 and the bicep2 central value @xmath68 , the string scale is derived as @xmath84 .
the tensor spectral index is @xmath85 thus string gas cosmology predicts blue tensor spectra .
however , note that is suppressed by slow roll parameters .
thus it would not explain the hint of order one @xmath13 shown in data .
it is also worth to mention that string gas cosmology produces highly gaussian density perturbations ( unless the string scale is near tev scale , which is not preferred from current data ) .
the @xmath86 estimator for string gas cosmology can be calculated as @xcite @xmath87 where @xmath88 corresponds to the scale of our present observable universe . with the detection of @xmath14 ,
the string scale is known in string gas cosmology . as a result , @xmath89 .
this is orders - of - magnitude smaller than the observational bound .
as we shall see , the inflationary candidates with blue @xmath13 , as far as we consider , produce considerable amount of non - gaussianities . thus in case
a blue tensor tilt is detected , non - gaussianity should be the next test ( and maybe already much better constrained than now by the time of @xmath13 measurement ) to distinguish between string gas and inflationary models .
in standard einstein - hilbert gravity , the amplitude of the tensor modes is determined by the energy scale of inflation . and thus for the case of blue @xmath13 , super inflation is needed , with null energy condition ( nec ) violation @xcite .
it has been debated for a long time that if the cosmological ( inflation , alternatives or dark energy ) background violates nec , the scalar sector of the inflationary perturbations should become a ghost . however , this situation has been changed since galileons are introduced to inflation .
here we shall explore the parameter space of the simplest g - inflation scenario , with slight generalization from the explicit model with an exponential potential in @xcite into a general slow roll functional dependence .
we start from the action of g - inflation @xmath90 where @xmath91 . here
we restrict our attention to the models with @xmath92    a special class of de sitter solution is worked out in @xcite : @xmath93 we shall consider small derivation of the de sitter solution . for this purpose ,
define small parameters @xmath94 when expanding around the de sitter solution , the @xmath10 , @xmath95 , and @xmath96 parameters should indeed be small .
the consistency for the smallness of the @xmath97 parameter shall be checked later .    in terms of those slow roll parameters , and near the de sitter solution , the order of magnitude of the following quantities can be estimated : @xmath98    the second order action of the scalar perturbations can be written as @xmath99~,\end{aligned}\ ] ] where with the above slow roll approximation , the @xmath100 and @xmath101 functions can be written as ( for unapproximated definition , see @xcite ) @xmath102 thus the scalar perturbations are stable in general for @xmath103 satisfying the above slow roll conditions .
the power spectrum and spectral index of the scalar sector can be calculated as @xmath104 where @xmath105 note that @xmath106 and @xmath107 near the de sitter solution .
thus @xmath108 .
the gravity sector is not modified .
thus the tensor mode has the conventional spectrum @xmath109 the tensor spectral index is @xmath110 and the tensor to scalar ratio in the slow - roll approximation is written as @xmath111 here the tensor tilt can be blue , but the blueness is suppressed by a slow roll parameter . also note that for the tensor mode to be blue , we need @xmath112 .
this is indeed possible .
however , it is worth noting that from , @xmath113 in other words , near the de sitter solution , the scalar and tensor modes tilt towards the same direction  both red spectra or both blue spectra .
note that a blue scalar spectrum is not favored by observations .
however , in more general cosmological models , for example , with a different number of neutrinos , a blue scalar spectrum is not yet ruled out .    on the other hand
, we are not sure at this point if the same direction of tilt for scalar and tensor spectra is a general feature of g - inflation ( and generalized galileons @xcite ) , or there exists unexplored models with different tilt of two sectors .    also , one has to note the sound speed of the scalar sector of perturbations is @xmath114 this would induce an equilateral non - gaussianity @xmath115 .
this is consistent with the current _ planck _ bound , and on the other hand accessible in the future .
it is worth mentioning that the small @xmath116 and large @xmath117 are not a coincidence in this particular toy model , but should be a rather model - independent statement ( unless fine tuned ) . to violate nec at the background level while keeping the perturbations stable , the relevant galileons lagrangian at the energy scale of inflation should be highly non - linear .
non - linear self - coupling introduces equilateral non - gaussianity .
experiments in the near future would reduce the bound for @xmath117 ( or detection ) and examine the possibility of galileons and nec violation .
finally , one can check that when @xmath118 , @xmath119 .
thus the smallness of @xmath97 is consistent with data .
it has been an open question if the inflationary ( scalar and tensor ) perturbations originates from a bunch - davies ( bd ) vacuum , or special care needs to be taken for the non - bd initial conditions .
one motivation for choosing non - bd initial conditions , for example , is the transplanckian problem of inflation @xcite ( see @xcite for a recent review ) .
inflation requires uv completion . during inflation ,
the perturbations originates from the uv completion scale , for example , planck scale , before they expand and cross the hubble horizon .
note that the planck scale could have been replaced by some lower scales , for example , the string scale . with those energy scales for new physics
, we can no longer make sure that the inflationary perturbations were in their lowest energy states before they got stretched by the cosmic expansion .
this situation is illustrated in fig .
[ fig : trans ] .
the transplanckian problem .
the inflationary perturbations are initialized on sub - planckian scales ( or above other uv completion scales).,scaledwidth=80.0% ]    in the context of tensor perturbations @xcite ( for a related recent study , see also @xcite ) , the action of tensor mode is @xmath120 decomposing the tensor modes by polarization , we have @xmath121~,\end{aligned}\ ] ] as a result , the action contains two copies of modes @xmath122~.\end{aligned}\ ] ] for each mode , one can impose a different initial condition from bd .
the @xmath123 ( collectively denote @xmath124 and @xmath125 ) field can be quantized as @xmath126 where @xmath127 the quantization condition requires @xmath128 the lowest energy state has @xmath129 , which corresponds to the bd state .
the tensor power spectrum for each polarization mode is @xmath130    @xmath131 ~.   \end{aligned}\ ] ]    note that the @xmath33-dependence of @xmath132 is determined by high energy physics and is unknown at low energy scales .
thus this @xmath33-dependence could in principle tilt the tensor spectrum to the blue end and overwrite the red tilt from decreasing @xmath133 .
this requires @xmath134\right\ } > 0~.\end{aligned}\ ] ] and the tensor tilt is @xmath135}{d \ln k}~. \end{aligned}\ ] ]    we would like to make a few comments before closing up this section :    * the non - bd state of tensor modes also opens up possibilities for small field inflation to be consistent with the observed large tensor to scalar ratio @xmath68 .
however , note that a reason is needed for the scalar sector and the tensor sector to achieve different non - bd coefficients . *
non - bd initial condition generically introduce non - gaussianities of the folded shape @xcite .
the size of the non - gaussianity is proportional to @xmath136 $ ] and peaked at the folded limit @xmath137 . in the squeezed limit @xmath138
, the folded non - gaussianity also blows up more quickly than local shape ( however , in the regime under theoretical control , the squeezed limit remains small @xcite ) .
thus the tensor modes generated by non - bd coefficients are non - gaussian . despite of the non - gaussianity in the tensor sector ,
it is also reasonable to expect a folded non - gaussianity for the scalar sector . because if the non - bd issue exists , it is likely not tensor - only unless a reason is provided
there is a huge landscape of inflation models , and also quite a few alternatives to inflation .
some of them may also be able to produce blue @xmath13 .
we list a few possibilities here .
besides gravitational waves from quantum fluctuations during inflation , there is a possibility that the dominant gravitational waves are generated by particles or strings @xcite , or particle states @xcite produced in the period inflation .
the kinetic energy of the inflaton @xmath139 dumps to the particles , and it is large enough to produce visible tensor spectrum .
the emission rate is related to the square of the coupling constant , i.e , for the particles with energy @xmath140 , @xmath141 .
one can derive a bound on @xmath10 from those class of particle production mechanism of sourcing tensor mode .
one can estimate the energy density of the classically sourced gravitational waves energy density @xmath142 near horizon crossing : @xmath143 where @xmath144 is the kinetic energy density of the inflaton . on the other hand , near horizon
crossing , the physical wave length of the tensor modes is of order hubble scale .
thus @xmath145 thus the tensor power spectrum @xmath146 from the observed tensor to scalar ratio @xmath147 we get @xmath148 from cobe normalization , we get @xmath149 below those scales , no classical source could generate primordial tensor modes .    in some simple models proposed in @xcite ,
the produced particles have a time - dependent mass .
@xmath150 due to its coupling to the rolling of inflaton . compared with the fact that inflation has the decreasing hubble constant , which predict a red tensor spectrum
, the increasing mass will afford an explanation of the blue tensor spectrum .
it is also important to note that non - gaussianities are also produced from the particle production process @xcite .
the constraints from non - gaussianity of particle production remains to be tested .
based on the slow roll expansion @xcite , it is recently pointed out that violation of slow roll can generate blue tensor spectrum @xcite . in @xcite , the tensor spectral index
is calculated up to second order in slow roll as @xmath151 thus if the slow roll expansion is extrapolated to @xmath152 , blue tensor spectrum is obtained .
for example , the possible blue tensor spectrum in the ultra - slow - roll inflation @xcite is discussed in @xcite .
it would be interesting to further examine this possibility without slow roll restrictions , and see if the introduced large running is consistent with data .
the tensor modes come from the gravity sector .
thus there is no surprise that if gravity is modified , the tensor spectrum could change .    among other possibilities ,
massive gravity is one example to achieve this goal
. massive gravity on a time dependent background is a challenge @xcite . nevertheless , there exist viable models @xcite and application to inflation @xcite .
we may expect tensor modes with blue spectra in massive gravity : if @xmath153 for the graviton during inflation , the tensor modes which exit the horizon earlier have more time to roll back to the origin of their mass potential , and thus are suppressed more . as a result
, a blue tensor tilt can be generated .
however , the massive gravity in the late universe does not help for the tensor modes .
this is because the mass of the graviton is too small .
also , the return - to - horizon effect actually generates a redder spectrum for the gravitational waves @xcite .      as another interesting class of possibilities ,
some inflation models come with space - like condensates of fields .
for example , space - like gradient of scalars in solid inflation @xcite and so(3 ) massive gravity @xcite , or space - like vector fields in chromo - natural inflation @xcite .
typically , so(3 ) spatial rotational symmetry is imposed to keep isotropy . in this class of models
, there can be additional tensor components from decomposition of the so(3 ) group , other than the graviton .
those components can also contribute to the b - mode power spectrum .
for the additional components , there is no nec type restriction against blue spectra .
for example , solid inflation predicts a slightly blue tensor spectrum @xcite .
note that the blue tilt in solid inflation is on the one hand suppressed by @xmath10 but on the other hand boosted by sound speed . for @xmath68
one finds @xmath154 .      with the assumption of isotropy
, scalar modes can not source the tensor modes at the linear level .
however , beyond the linear level , tensor modes can be sourced by scalar modes @xcite .
those second order effects are typically unobservably small because they are suppressed by the inflationary power spectrum .
however , the second order contribution to tensor mode can be boosted by a small sound speed of an isocurvature scalar sector @xcite .
the generated gravitational waves are highly non - gaussian .
the tilt of tensor spectrum coming from those contributions are not constrained by the null energy condition , thus a blue tilt is possible .      as
another alternative to inflation scenario , in matter bounce , the universe was in a contracting phase before bouncing back and heating up ( see @xcite for a review ) . in matter
bounce , the tensor modes are nearly scale invariant but the amplitudes have been too high .
recently , a two - field bounce model is introduced and the tensor modes can be tuned as a parameter @xcite .
it is thus interesting to see if this model ( among other possibilities ) fit the current observations and has a potential for blue tensor spectrum .
to conclude , we fit the tensor - to - scalar ratio and the tensor spectral tilt with data .
the current data is not good enough to test the inflationary consistency relation but nevertheless blue tensor spectra are favored when all 9 bins of bicep2 data , in combination with the polarbear data , are used .    from theoretical aspects ,
string gas cosmology predicts blue tensor spectra .
however , the tilt is small , at the same order - of - magnitude of scalar tilt .
thus the future experiments targeting to test the inflationary consistency relation can also test this prediction of string gas cosmology . on the other hand ,
string gas cosmology predicts highly gaussian density and tensor perturbations .
this is unlike the other inflationary mechanisms , where a blue tilt also implies non - gaussianities .
the simplest model of g - inflation , on the other hand , tilts scalar and tensor power spectra in the same way .
thus with a red scalar spectral tilt , the tensor spectra are also red . in the parameter regime where both scale and tensor spectra are blue , the blueness of tensor perturbation is suppressed by slow roll parameter @xmath10 .
the equilateral non - gaussianity of g - inflation , at @xmath68 , is about @xmath155 and close to the current observational bound .
the non - gaussianities should be a model independent feature for super - inflation type models which generate blue tensor spectrum , because the perturbations have to behave differently from the background to avoid ghosts .
generalized initial conditions of inflation is left largely unconstrained .
however , the generalized initial conditions are also sources of non - gaussianities .
thus non - gaussianities in the tensor sector would be a test of those class of models .
inflationary particle production is another possible source of tensor modes .
we derived a bound of inflationary hubble scale and energy density for this mechanism to work .
non - gaussianities are also present in the case of particle productions .
we hope to investigate the particle production mechanism and its relation between the blue tensor spectra in a future work .
the current data shows hint of blue @xmath156 .
the models with slow roll suppressed @xmath13 would not be enough to explain such a hint .
on the other hand , from modified initial conditions , external sources , inflation beyond slow roll , modified gravity and some models of matter bounce , blue and large @xmath13 may be produced .
it would be interesting to study those models in more details to see how the fitting of data is improved in those models .
* note added : *    two related works @xcite appeared on arxiv on the same day as ours .
@xcite ( see also @xcite ) overlaps with our section [ sec : hints - blue - tensor ] , and @xcite overlaps with section [ sec : infl - gener - init ] .    the data analysis of tensor tilt is also investigated by @xcite , @xcite and @xcite , within a few days before / after our paper . among those papers , @xcite and @xcite reports a nearly zero central value , with small errorbars @xmath157 ( bicep2 only ) and @xmath158 ( bicep2 + planck + wp , with running of scalar spectral index ) . while the following works prefer blue tilt with considerably larger errorbars : @xcite ( bicep2 only ) , our result ( bicep2 , bicep2 + planck + wp , polarbear ) , @xcite ( bicep2 + planck + wp and bicep2 x keck + planck + wp ) and @xcite ( bicep2 + planck + wmap + bao ) .
yw is supported by a starting grant of the european research council ( erc stg grant 279617 ) , and the stephen hawking advanced fellowship .
this work was undertaken on the cosmos shared memory system at damtp , university of cambridge operated on behalf of the stfc dirac hpc facility .
this equipment is funded by bis national e - infrastructure capital grant st / j005673/1 and stfc grants st / h008586/1 , st / k00333x/1 .
p.  a.  r.  ade _ et al .
_ [ bicep2 collaboration ] , arxiv:1403.3985 [ astro-ph.co ] .
m.  j.  mortonson and u.  seljak , arxiv:1405.5857 [ astro-ph.co ] .
r.  flauger , j.  c.  hill and d.  n.  spergel , jcap * 1408 * , 039 ( 2014 ) [ arxiv:1405.7351 [ astro-ph.co ] ] .
e.  j.  copeland , a.  r.  liddle , d.  h.  lyth , e.  d.  stewart and d.  wands , phys .
d * 49 * , 6410 ( 1994 ) [ astro - ph/9401011 ] .
y.  -z .
ma and y.  wang , arxiv:1403.4585 [ astro-ph.co ] .
p.  a.  p.  ade et al .
, planck 2013 results xxii , arxiv : 1303.5082          a.  lewis and s.  bridle , phys .
d * 66 * , 103511 ( 2002 ) [ astro - ph/0205436 ] .
d.  blas , j.  lesgourgues and t.  tram , jcap * 1107 * ( 2011 ) 034 [ arxiv:1104.2933 [ astro-ph.co ] ] .
g.  hinshaw _ et al . _
[ wmap collaboration ] , astrophys .
j.  suppl .
* 208 * ( 2013 ) 19 [ arxiv:1212.5226 [ astro-ph.co ] ] .
a.  ashoorioon , k.  dimopoulos , m.  m.  sheikh - jabbari and g.  shiu , arxiv:1403.6099 [ hep - th ] .
p.  a.  r.  ade _ et al .
_ [ planck collaboration ] , arxiv:1303.5075 [ astro-ph.co ] .
d.  hanson _ et al . _
[ sptpol collaboration ] , phys .
lett .   * 111 * , 141301 ( 2013 ) [ arxiv:1307.5830 [ astro-ph.co ] ] .
y.  -f .
cai and y.  wang , arxiv:1404.6672 [ astro-ph.co ] .
r.  h.  brandenberger and c.  vafa , nucl .
b * 316 * , 391 ( 1989 ) .
a.  nayeri , r.  h.  brandenberger and c.  vafa , phys .
lett .   *
97 * , 021302 ( 2006 ) [ hep - th/0511140 ] .
r.  h.  brandenberger , a.  nayeri , s.  p.  patil and c.  vafa , phys .
lett .   * 98 * , 231302 ( 2007 ) [ hep - th/0604126 ] .
r.  h.  brandenberger , a.  nayeri , s.  p.  patil and c.  vafa , int .
j.  mod .
a * 22 * , 3621 ( 2007 ) [ hep - th/0608121 ] .
r.  h.  brandenberger , a.  nayeri and s.  p.  patil , arxiv:1403.4927 [ astro-ph.co ] .
j.  khoury , b.  a.  ovrut , p.  j.  steinhardt and n.  turok , phys .
d * 64 * , 123522 ( 2001 ) [ hep - th/0103239 ] .
t.  biswas , t.  koivisto and a.  mazumdar , jhep * 1408 * , 116 ( 2014 ) [ arxiv:1403.7163 [ hep - th ] ] .
l.  a.  boyle , p.  j.  steinhardt and n.  turok , phys .
d * 69 * , 127302 ( 2004 ) [ hep - th/0307170 ] .
t.  qiu , j.  evslin , y.  -f .
cai , m.  li and x.  zhang , jcap * 1110 * , 036 ( 2011 ) [ arxiv:1108.0593 [ hep - th ] ] . b.  chen , y.  wang , w.  xue and r.  brandenberger , arxiv:0712.2477 [ hep - th ] . y.  -s .
piao and y.  -z .
zhang , phys .
d * 70 * , 063513 ( 2004 ) [ astro - ph/0401231 ] .
t.  kobayashi , m.  yamaguchi and j.  i .
yokoyama , phys .
lett .   * 105 * , 231302 ( 2010 ) [ arxiv:1008.0603 [ 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 and g.  zahariade , phys .
d * 84 * , 064039 ( 2011 ) [ arxiv:1103.3260 [ hep - th ] ] .
t.  kobayashi , m.  yamaguchi and j.  i .
yokoyama , prog .
phys .   * 126 * , 511 ( 2011 ) [ arxiv:1105.5723 [ hep - th ] ] .
r.  h.  brandenberger and j.  martin , mod .
a * 16 * , 999 ( 2001 ) [ astro - ph/0005432 ] .
j.  martin and r.  h.  brandenberger , phys .  rev .
d * 63 * , 123501 ( 2001 ) [ hep - th/0005209 ] .
r.  h.  brandenberger and j.  martin , class .
grav .   * 30 * , 113001 ( 2013 ) [ arxiv:1211.6753 [ astro-ph.co ] ] .
t.  biswas , r.  brandenberger , t.  koivisto and a.  mazumdar , phys .
d * 88 * , no . 2 , 023517 ( 2013 ) [ arxiv:1302.6463 [ astro-ph.co ] ] .
a.  ashoorioon , k.  dimopoulos , m.  m.  sheikh - jabbari and g.  shiu , jcap * 1402 * , 025 ( 2014 ) [ arxiv:1306.4914 [ hep - th ] ] .
h.  collins , r.  holman and t.  vardanyan , arxiv:1403.4592 [ hep - th ] .
x.  chen , m.  -x .
huang , s.  kachru and g.  shiu , jcap * 0701 * , 002 ( 2007 ) [ hep - th/0605045 ] .
w.  xue and b.  chen , phys .
d * 79 * ( 2009 ) 043518 [ arxiv:0806.4109 [ hep - th ] ] .
r.  flauger , d.  green and r.  a.  porto , jcap * 1308 * , 032 ( 2013 ) [ arxiv:1303.1430 [ hep - th ] ] .
a.  aravind , d.  lorshbough and s.  paban , jhep * 1307 * , 076 ( 2013 ) [ arxiv:1303.1440 [ hep - th ] ] .
l.  senatore , e.  silverstein and m.  zaldarriaga , arxiv:1109.0542 [ hep - th ] . j.  l.  cook and l.  sorbo , phys .
d * 85 * , 023534 ( 2012 ) [ erratum - ibid .
d * 86 * , 069901 ( 2012 ) ] [ arxiv:1109.0022 [ astro-ph.co ] ] .
d.  carney , w.  fischler , e.  d.  kovetz , d.  lorshbough and s.  paban , jhep * 1211 * , 042 ( 2012 ) [ arxiv:1209.3848 [ hep - th ] ] .
n.  barnaby , j.  moxon , r.  namba , m.  peloso , g.  shiu and p.  zhou , phys .
d * 86 * , 103508 ( 2012 ) [ arxiv:1206.6117 [ astro-ph.co ] ] .
j.  -o .
gong , class .
grav .   * 21 * , 5555 ( 2004 ) [ gr - qc/0408039 ] .
j.  -o .
gong , arxiv:1403.5163 [ astro-ph.co ] . n.  c.  tsamis and r.  p.  woodard , phys .
d * 69 * , 084005 ( 2004 ) [ astro - ph/0307463 ] .
x.  chen , h.  firouzjahi , m.  h.  namjoo and m.  sasaki , europhys .
lett .   * 102 * , 59001 ( 2013 ) [ arxiv:1301.5699 [ hep - th ] ] .
g.  damico , c.  de rham , s.  dubovsky , g.  gabadadze , d.  pirtskhalava and a.  j.  tolley , phys .
d * 84 * , 124046 ( 2011 ) [ arxiv:1108.5231 [ hep - th ] ] .
a.  e.  gumrukcuoglu , c.  lin and s.  mukohyama , jcap * 1111 * , 030 ( 2011 ) [ arxiv:1109.3845 [ hep - th ] ] .
a.  de felice and s.  mukohyama , phys .
b * 728 * , 622 ( 2014 ) [ arxiv:1306.5502 [ hep - th ] ] .
a.  de felice , a.  e.  gumrukcuoglu and s.  mukohyama , phys .
d * 88 * , 124006 ( 2013 ) [ arxiv:1309.3162 [ hep - th ] ] .
s.  dubovsky , r.  flauger , a.  starobinsky and i.  tkachev , phys .
d * 81 * , 023523 ( 2010 ) [ arxiv:0907.1658 [ astro-ph.co ] ] . c.  lin , arxiv:1307.2574 .
a.  e.  gumrukcuoglu , s.  kuroyanagi , c.  lin , s.  mukohyama and n.  tanahashi , class .
* 29 * , 235026 ( 2012 ) [ arxiv:1208.5975 [ hep - th ] ] .
s.  endlich , a.  nicolis and j.  wang , jcap * 1310 * , 011 ( 2013 ) [ arxiv:1210.0569 [ hep - th ] ] .
m.  akhshik , r.  emami , h.  firouzjahi and y.  wang , arxiv:1405.4179 [ astro-ph.co ] .
p.  adshead and m.  wyman , phys .
lett .   * 108 * , 261302 ( 2012 ) [ arxiv:1202.2366 [ hep - th ] ] .
x.  chen , r.  emami , h.  firouzjahi and y.  wang , jcap * 1408 * , 027 ( 2014 ) [ arxiv:1404.4083 [ astro-ph.co ] ] .
s.  matarrese , s.  mollerach and m.  bruni , phys .
d * 58 * , 043504 ( 1998 ) [ astro - ph/9707278 ] .
m.  biagetti , m.  fasiello and a.  riotto , phys .  rev .
d * 88 * , no .
10 , 103518 ( 2013 ) [ arxiv:1305.7241 [ astro-ph.co ] ] .
r.  h.  brandenberger , int .
j.  mod .
01 * , 67 ( 2011 ) [ arxiv:0902.4731 [ hep - th ] ] .
y.  -f .
cai , e.  mcdonough , f.  duplessis and r.  h.  brandenberger , jcap * 1310 * , 024 ( 2013 ) [ arxiv:1305.5259 [ hep - th ] ] .
m.  gerbino , a.  marchini , l.  pagano , l.  salvati , e.  di valentino and a.  melchiorri , arxiv:1403.5732 [ astro-ph.co ] .
antony lewis , http://cosmocoffee.info/viewtopic.php?t=2302 .