The Early Universe with HighScale Supersymmetry
Abstract
A small tensortoscalar ratio may lead to distinctive phenomenology of highscale supersymmetry. Assuming the same origin of SUSY breaking between the inflation and visible sector, we show model independent features. The simplest hybrid inflation, together with a new linear term for the inflaton field which is induced by large gravitino mass, is shown to be consistent with all experimental data for of order . For superpartner masses far above the weak scale we find that the reheating temperature after inflation might be beneath the value required by thermal leptogenesis if the inflaton decays to its products perturbatively, but above it if nonperturbatively instead. Remarkably, the gravitino overproduction can be evaded in such highscale supersymmetry because of the kinematically blocking effect.
I Introduction
After the discovery Higgs () of standard model (SM) Higgs boson at the Large Hadron collider (LHC), lowscale supersymmetry (SUSY) susy () which is favored by the naturalness argument naturalness () has been extensively explored. These studies show the difficulties in the both theoretic explanation of 125 GeV Higgs mass and experimental fits to the LHC data. The second run of LHC will shed light on the prospect of such natural SUSY models. Given above consideration some efforts have been devoted to the study of highscale SUSY.
Even though highscale SUSY cannot be detected at the TeV LHC, they can be still studied via their effects on the evolution of early universe. Measurement on the tensortoscalar ratio via experiments such as WAMP, Plank and BICEP that are devoted to measure Cosmic Microwave Background (CMB) temperature anisotropy and polarization during inflation, may probe highscale SUSY with mass spectrum far above the weak scale. The measured value of reported by the Plank Collaboration is of order at 95 % CL 1502.02114 (); 1303.5076 (); 1303.5082 (), from which the energy scale of inflation can be directly inferred. Since the energy scale of inflation is proportional to , it is only mildly sensitive to . So the study of highscale SUSY remains well motivated as long as is not extremely small.
In this paper, we consider inflationary models with far below the Plank bound value . The motivation is mainly based on two facts. At first, the stability of the SM electroweak vacuum requires 1505.04825 (), where refers to the value at which the Higgs potential is maximal. For the central value of top quark pole mass GeV, which implies that should be smaller than Plank bound value GeV corresponding to 1504.08093 (). In this sense small is more favored to guarantee the electroweak vacuum stable against quantum fluctuation during inflationary epoch. Secondly, for , it can still generate SUSY mass spectrum large enough to escape the LHC constraints.
For simplicity we adopt the assumption that the inflation and visible (namely the minimal supersymmetric standard model (MSSM )) sector share the same origin of SUSY breaking. This assumption is rational, as it can be realized in model building. Moreover, it allows us to discuss reheating in the early universe after inflation, once the SUSY mass spectrum and the inflaton decay are identified explicitly.
The paper is organized as follows. In section 2 we reanalyze the model independent consequences from above assumption within the range . In section 3, we consider hybrid inflation as an example in the course of highscale SUSY breaking ^{1}^{1}1For earlier attempts to address this issue, see, .e.g, 1009.5340 (); 1211.0994 ().. We will show that a new linear term for inflaton field with a large coefficient proportional to affects the inflation significantly, and the simplest hybrid inflation is consistent with of order .
In the second part of this paper we discuss the reheating in the early universe after inflation in section 4. In particular, reheating temperature after inflation is estimated for superpartner mass spectrum above TeV. We find that might be beneath the value GeV required by thermal leptogenesis if inflaton decays to its products perturbatively, but above it if nonperturbatively instead. The gravitino overproduction in conventional highscale SUSY can be easily evaded because of kinematically blocking effect. Finally we conclude in section 5.
Ii Implications of the value of to Inflation
In this section we revise the model independent implications (together with Plank and 9year WAMP data) to singlefield inflation for . These results provide useful information on the model building of inflation, as a reliable inflation model should at least explain observable quantities as what follows.
, First of all the scale of energy density during inflation is directly related to as,
(1) 
where is the amplitude of the power spectrum of the curvature perturbation and in the context of singlefield inflation. Recall that , where is the value of when wavenumber Mpc crossed outside the horizon. Substituting the reported value by Plank Collaboration 1303.5076 () into Eq.(1) gives rises to,
(2) 
Eq.(2) is valid independent of the explicit form of inflation potential. If the MSSM and inflation share the same origin of SUSY breaking, as we have assumed in this paper, the highscale SUSY breaking scale in particle physics will be order . For example, is of order GeV for , and slightly reduced to be of order GeV for .
, In the context of slow roll inflation the spectral index (for scalar) and (for tensor) are given by,
(3) 
respectively. Here and , with subscript denoting derivative of over . The combination of Plank and 9year WAMP data measures the value of in high precision 1303.5082 (),
(4) 
For one finds that . This tight bound is crucial to constrain inflation model.
, The gravitino mass can be determined. The constant superpotential , which is required to cancel out positive term in the potential so as to explain the smallness of cc, gives rise to
(5) 
, Finally the number of efold that undergoes during inflation is given by,
(6) 
where . Subscript “in" and “end" corresponds to initial and end value of during inflation, respectively. For realistic inflation models, is bounded as . If doesn’t change significantly during inflation, Eq.(6) can be expressed as . This is known as Lyth bound Lythbound (), which shows the need of small field inflation for .
Iii The Simplest Hybrid Inflation
The section is devoted to the study of inflation building in the course of highscale SUSY. We take the simplest hybrid inflation as an explicit illustration. We will show that a new linear term due to the assumption adopted in this paper significantly affects the choice on initial condition. Also this assumption introduces new constraints on parameters in the model, which make the simplest hybrid inflation only possible with of order .
iii.1 Scalar Potential
The scalar potential in hybrid inflation is constructed from superpotential ,
(7) 
and Kahler potential ,
(8) 
Here denotes the inflaton superfield, with its lowest component inflaton field . is a singlet of standard model gauge groups . and denote waterfall superfields which are in the bifundamental representation of ^{2}^{2}2 Alternatively, can be extended to include a local symmetry so as to explain leptogenesis.. The Kahler potential in Eq.(8) takes into account the noncanonical term, with a real coefficient. The noncanonical term provides the inflaton mass term. GeV is the reduced Plank mass, while is assumed to be far below .
Substituting Eq.(7) and Eq.(8) into the SUGRA potential
(9) 
one obtains the scalar potential of hybrid inflation. Here is the Kahler metric and . It is well known that the history of inflation can be naturally divided into two periods. In the first one inflation usually starts from an initial value of order towards to which is a critical value separating the two periods. The vacuum in the first period corresponds to , from which the energy density reads from Eq.(9) 1007.5152 (),
(10) 
where we have defined . Here and inflaton mass is given by for negative ,
(11) 
The term in Eq.(10) represents the contribution due to mass splitting in waterfall fields 9406319 (), with the cutoff scale. The linear term with coefficient proportional to arises from a constant superpotential added to , which is needed to cancel out positive contribution to energy density due to SUSY breaking, and explain the smallness of cc.
When approaches to , becomes massless as shown from its mass squared . After the time when is below , starts to roll towards to its global minimum value from , which is known as the second period of inflation. The evaluation of field (including angular component ) and during each period is determined by their equations of motion,
(12)  
where the Hubble constant is subject to the Friedmann constraint
(13) 
The time for each period is controlled by the magnitude of or relative to Hubble constant . As pointed out in hybrid (), the second period is very short in compared with the first one for wide ranges of parameter choices. Substituting Eq.(13) into the last equation in Eq.(III.1), we obtain the constraint for such property,
(14) 
In the next subsection, we will discuss in more details the initial conditions on inflaton field and the field value of when inflation ends.
iii.2 Initial Conditions
The inflation usually begins at some field value near Plank scale. The choice on is subtle when the inflaton potential has either a few local minimums at s, or local maximum at s. If one adopts bigger than , inflaton is probably trapped at these local minimums of inflaton potential along the trajectory, which leads to inflation with insufficient efold number . In order to avoid this, one should choose . On the other hand, one wants that inflation proceeds with exactly decreasing . This is only allowed if is less than . In other words, we should impose the initial condition
(15) 
In Fig.1 we show how extremes in depend on and by evaluating . The sign of changes when for and for . This implies that for realistic inflation. Otherwise, , which is too small to provide enough efold number . This observation has been noted in 1007.5152 () for of order electroweak scale, and further verified for larger value of GeV. With initial value , the initial value can be chosen in wide range, as shown in Fig.1. The evaluation of from is the same as original hybrid model because of absence of linear term in the first equation in Eq.(III.1). In this sense, inflation mainly ends at field value .
The efold number produced during inflation and can be both estimated in terms of slow roll parameters and in the model, which are given by respectively,
(16) 
With and , and mainly depend on parameter .
Fig.2 shows the bound on for two choices of (red curve) and 0.5 (black curve) respectively. Note that is constrained from slow roll condition . Given the range shown in Eq.(14) for , should be below unity, which implies that large field inflation is excluded under our assumption. Fig.2 shows that for .
In Fig.3 we show how changes for two typical choices of subtracted from Fig. 2. It clearly indicates that for observed value of , is of order . The simplest hybrid inflation can provide large efolds number and small . Nevertheless, large and large can not be induced at the same time. In the next section, we focus on reheating after (or during) inflation.
Iv Reheating in HighScale SUSY
When inflation ends the conversion of energy to the MSSM matters from the inflaton begins immediately. The efficiency of energy transfer depends on how inflaton is coupled to the MSSM matters, the magnitude of their couplings, and the SUSY mass spectrum. In general, the ways of energy transfer include the perturbative and nonperturbative decay of inflaton ^{3}^{3}3For reviews, see, e.g., 0507632 (); 1001.0993 ().. The later way is known as preheating 9405187 (); 9704452 (); 9711360 (). The conditions between these two ways of energy transfer are rather different. In the later case, parameter resonance requires a quartic interaction term with large magnitude of . This only happens if one allows renormalizable superpotential term of mass dimension 4 0603244 (),
(17) 
where is the inflaton superfield, and are Higgs doublet superfields. denotes the parity of inflaton, which is useful to keep the dark matter stable. In contrast, in the SM the inflaton couples to SM chiral fermions and gauge bosons in terms of nonrenormalizable interactions of mass dimension 5. Quartic term above doesn’t exist in the SM, and therefore the way of energy transfer in the SM is perturbative decay. In what follows, we consider these two ways separately.
iv.1 Perturbative Decay
As briefly mentioned above, perturbative decay happens either when there is no renormalizable interaction in Eq.(17) or the quartic coupling constant is tiny. Instead, the inflaton only decays to SM matters via fivedimensional operators such as
(18) 
Here s refer to strengths of SM gauge fields, s refer to SM fermions and represents the mass scale appearing in the fivedimensional operators. The plasma will be MSSMlike if the reheating temperature is larger than the typical scale of superpartner mass, . Otherwise, the plasma is actually SMlike.
Now we calculate the reheating temperature. We organize the decay width of inflaton to SM particles as,
(19) 
We simply take but leave as a free parameter. The thermal equilibrium of relativistic plasma is dominated by , the rate for SM inelastic scatterings of processes perturbative () . is related to as ,
(20) 
where is the SM fine structure constant. Therefore the reheating temperature for the perturbative decay isn’t equal to the conventional one, defined as , denotes the number of relativistic number. Instead, in terms of Eq.(20) we have
(21) 
where denote the quantum numbers of SM gauge groups. Note that Eq.(21) is valid for , which implies that the reheating temperature can serve as the upper bound on .
In Fig.4 we show reheating temperature as function of inflaton mass. Note that captures the magnitude of coupling between inflaton and “mediate " field, which also couples to the SM matter and gauge fields. Here a few comments are in order. (1), For , is below the lower bound GeV required by thermal leptogenesis in the whole range of . (2), Since is the upper bound on superpartner mass spectrum , one finds that is upper bounded as TeV TeV for the case of perturbative decay. (3), As is far below the gravitino mass of Eq.(5), there is no overproduction problem of gravitino in highscale SUSY. Superheavy gravitino mass of order GeV kinetically blocks its production in the thermal bath.
iv.2 Nonperturbative Decay
If it admits renormalizable superpotential Eq.(17), there exists quartic interaction between inflaton and its decay products , Nonperturbative decay can happen in wide range of parameter space for potential of type^{4}^{4}4Assuming inflaton and MSSM matters share the same origin of SUSY breaking, inflaton mass is dynamical induced by SUSY breaking. In this sense, the mass term is a soft SUSYbreaking term other than arises from SUSY treelevel mass superpotential . Consequently, there is no cubic interaction in compared with earlier discussions in 0603244 (); 0512227 ().,
(22) 
where the inflaton mass term is included and is the mass for . We would like to mention that include soft SUSY breaking contribution of order and dynamical mass induced by VEV of flat direction negative () through the quartic interaction 0512227 (),
(23) 
where is the quartic coupling constant. The magnitude of is determined by the selfinteraction potential for flat direction .
Now we consider the potential for flat direction. includes soft breaking mass, Hubble parameter induced term and high dimensional operators,
(24) 
where is real coefficient. Since isn’t far beneath the Hubble constant at the beginning of inflation, VEV depends on the sign of , which can be either positive or negative negative (); positive (). In particular, for the case of either positive or negative but with . It implies that SM gauge symmetry is unbroken during the whole history of early universe. On the other hand, for negative but with . It implies that SM gauge symmetry is broken in the eary universe, with gauge boson mass of order , then restored after the epoch of reheating.
The potentials we define in Eq.(22) to Eq.(23) are rather general, which can be applied to both cases in Eq.(17). To discuss the condition for parameter resonance, one starts with the modified KleinGordon equation for Fourie modes . Whether WKB approximation is viable for the study can be analyzed in term of a quantity defined as 0507632 (),
(25) 
where dot refers to derivative over time and is the frequency, with the expansion factor and the momentum. If , the WKB approximation is valid, the produced particle number of doesn’t grow in this case. If instead, the WKB approximation isn’t valid, which leads to significant production of . In long wavelengths limit, this constraint is given by ^{5}^{5}5There is a coefficient of order one in front of for either or . Here we simply take it equal to unity.,
(26) 
where we have used Eq.(22) and Eq.(23). Moreover, in order to keep that the parameter resonance isn’t spoiled by expansion, an additional constraint must be imposed,
(27) 
where refers to the amplitude of inflaton oscillations. For more details, we refer to reader to 0507632 () and references therein.
Fig.5 shows reheating temperature for the case in which . In this case, SM gauge bosons are massive and the rate for thermal equilibrium 0512227 () depends on its magnitude relative to . For , depends on both and , whereas it mainly depends on for . In this figure, we take GeV and . The bounds on are due to a few considerations. The first one is that condition Eq.(27) from parameter resonance requires . The second one is that overproduction problem gravitino () of gravitino in highscale SUSY can be kinematically blocked if such that nonperturbatively induced mass during parameter resonance is beneath ^{6}^{6}6 Ref. 1404.1914 () provides an example how gravitino problem in highscale SUSY is evaded in the context of miniSplit SUSY. In comparison with 1404.1914 (), gravitino mass is far heavier in this paper, and kinematically blocking is the solution to the overproduction of gravitino.. The bound on arises from condition Eq.(26) which shows .
The panel in Fig.5 shows that in the range GeV GeV reheating temperature GeV in the allowed range of for . With modifying ,
(28) 
which is always above the value required by thermal leptogenesis. The panel in Fig.5 shows that in the range GeV GeV reheating temperature GeV GeV if , and changes similarly to Eq.(28) for modifying . This implies that is also always above the the value required by thermal leptogenesis. Note that is excluded by condition Eq.(26) from parameter resonance for .
Fig.6 shows reheating temperature for the case in which . In this case, SM gauge symmetries are unbroken in the epoch of reheating. Thermalization cannot occur before the inflaton decay has completed. Due to the scatterings are already efficient when . The reheat temperature in this case is given by the standard expression: . However, calculated via renormalizable couplings is rather different from Eq.(19) calculated via nonrenormazible couplings. This difference between Fig.4 and Fig.6 is obvious. Typically, we have GeV in nonperturbative decay into MSSM and GeV in perturbative decay into SM.
V Conclusions
In the light of both LHC data and Plank bound on , Highscale SUSY is more favored in compared with lowscale SUSY. In this paper, we discussed the implications of highscale SUSY to the early universe. In particular, we assumed that the inflation and visible sector share the same origin of SUSY breaking, and derived model independent consequences based on this assumption. We find that the reheating temperature for superpartner mass spectrum above TeV might be beneath the value required by thermal leptogenesis if inflaton decays to its products perturbatively but above it if nonperturbatively instead. We also observed that problem of gravitino overproduction can be evaded through kinematically blocking in a wide range of parameter space in the later way.
As an illustration for the model building of inflation in the course of highscale SUSY, in section 3 we revise the simplest hybrid inflation that includes a new linear term for inflaton with coefficient proportional to . It is shown that this term significantly affects the choices on initial condition of inflaton fields. We found that with the assumption the simplest hybrid inflation is consistent with present experimental data for of order .
Under our assumption only the dark matter is a light SUSY state with mass near the weak scale Zheng (), which is the target of LUX and Xenon experiments, etc. Hopefully, it can be addressed in the near further.
Acknowledgements.
We would like to thank J.h, Huang for discussion. This work is supported in part by the Fundamental Research Funds for the Central Universities under Grant No.CQDXWL2013015.References

(1)
G. Aad et al, [],
“Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC,”
Phys. Lett. B 716, 1 (2012),
arXiv:1207.7214 [hepex];
S. Chatrchyan et al, [], “Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC,” Phys. Lett. B 716, 30 (2012), arXiv:1207.7235 [hepex].  (2) See, e.g., S. P. Martin, “A Supersymmetry primer,” Adv. Ser. Direct. High Energy Phys. 21, 1 (2010). [hepph/9709356].
 (3) See, e.g., J. L. Feng, “Naturalness and the Status of Supersymmetry,” Ann. Rev. Nucl. Part. Sci. 63, 351 (2013), arXiv:1302.6587 [hepph], and references therein.
 (4) P. A. R. Ade et al. [Planck Collaboration], arXiv:1502.02114 [astroph.CO].
 (5) P. A. R. Ade et al. [ Collaboration], “Planck 2013 results. XVI. Cosmological parameters,” arXiv:1303.5076 [astroph.CO].
 (6) P. A. R. Ade et al. [ Collaboration], “Planck 2013 results. XXII. Constraints on inflation,” arXiv:1303.5082 [astroph.CO].
 (7) J. R. Espinosa, G. F. Giudice, E. Morgante, A. Riotto, L. Senatore, A. Strumia and N. Tetradis, “The cosmological Higgstory of the vacuum instability,” arXiv:1505.04825 [hepph].
 (8) See, e.g., S. Zheng, “Can Higgs Inflation be Saved with Highscale Supersymmetry ?,” arXiv:1504.08093 [hepph].
 (9) Q. Shafi and J. R. Wickman, Phys. Lett. B 696, 438 (2011), arXiv:1009.5340 [hepph].
 (10) T. Higaki, K. S. Jeong and F. Takahashi, JHEP 1212, 111 (2012), arXiv:1211.0994 [hepph].
 (11) D. H. Lyth, “What would we learn by detecting a gravitational wave signal in the cosmic microwave background anisotropy?,” Phys. Rev. Lett. 78, 1861 (1997), [hepph/9606387].
 (12) K. Nakayama, F. Takahashi and T. T. Yanagida, “Constraint on the gravitino mass in hybrid inflation,” JCAP 1012, 010 (2010), arXiv:1007.5152 [hepph].
 (13) G. R. Dvali, Q. Shafi and R. K. Schaefer, “Large scale structure and supersymmetric inflation without fine tuning,” Phys. Rev. Lett. 73, 1886 (1994). [hepph/9406319]
 (14) A. D. Linde, “Hybrid inflation,” Phys. Rev. D 49, 748 (1994), [astroph/9307002].
 (15) B. A. Bassett, S. Tsujikawa and D. Wands, “Inflation dynamics and reheating,” Rev. Mod. Phys. 78, 537 (2006), [astroph/0507632].
 (16) A. Mazumdar and J. Rocher, “Particle physics models of inflation and curvaton scenarios,” Phys. Rept. 497, 85 (2011), arXiv:1001.0993 [hepph].
 (17) L. Kofman, A. D. Linde and A. A. Starobinsky, Phys. Rev. Lett. 73, 3195 (1994), [hepth/9405187].
 (18) L. Kofman, A. D. Linde and A. A. Starobinsky, “Towards the theory of reheating after inflation,” Phys. Rev. D 56, 3258 (1997), [hepph/9704452].
 (19) J. GarciaBellido and A. D. Linde, “Preheating in hybrid inflation,” Phys. Rev. D 57, 6075 (1998), [hepph/9711360].
 (20) R. Allahverdi and A. Mazumdar, “Reheating in supersymmetric high scale inflation,” Phys. Rev. D 76, 103526 (2007), [hepph/0603244].

(21)
S. Davidson and S. Sarkar,
“Thermalization after inflation,”
JHEP 0011, 012 (2000),
[hepph/0009078];
R. Allahverdi and M. Drees, “Thermalization after inflation and production of massive stable particles,” Phys. Rev. D 66, 063513 (2002), [hepph/0205246].  (22) M. Dine, L. Randall and S. D. Thomas, “Supersymmetry breaking in the early universe,” Phys. Rev. Lett. 75, 398 (1995), [hepph/9503303]; “Baryogenesis from flat directions of the supersymmetric standard model,” Nucl. Phys. B 458, 291 (1996), [hepph/9507453].
 (23) M. K. Gaillard, H. Murayama and K. A. Olive, “Preserving flat directions during inflation,” Phys. Lett. B 355, 71 (1995), [hepph/9504307].
 (24) R. Allahverdi and A. Mazumdar, “Supersymmetric thermalization and quasithermal universe: Consequences for gravitinos and leptogenesis,” JCAP 0610, 008 (2006), [hepph/0512227].
 (25) M. Kawasaki, F. Takahashi and T. T. Yanagida, “Gravitino overproduction in inflaton decay,” Phys. Lett. B 638, 8 (2006), [hepph/0603265].
 (26) J. Fan, B. Jain and O. Ozsoy, “Heavy Gravitino and Split SUSY in the Light of BICEP2,” JHEP 07, 073 (2014), arXiv:1404.1914 [hepph].
 (27) S. Zheng, “On Dark Matter Selected HighScale Supersymmetry,” JHEP 1503, 062 (2015), arXiv:1409.2939 [hepph].