MSSM-inspired multifield inflation

Despite the fact that experimentally with a high degree of statistical significance only a single Standard Model-like Higgs boson is discovered at the LHC, extended Higgs sectors with multiple scalar fields not excluded by combined fits of the data are more preferable theoretically for internally consistent realistic models of particle physics. We analyze the inflationary scenarios which could be induced by the two-Higgs-doublet potential of the Minimal Supersymmetric Standard Model (MSSM) where five scalar fields have non-minimal couplings to gravity. Observables following from such MSSM-inspired multifield inflation are calculated and a number of consistent inflationary scenarios are constructed. Cosmological evolution with different initial conditions for the multifield system leads to consequences fully compatible with observational data on the spectral index and the tensor-to-scalar ratio. It is demonstrated that the strong coupling approximation is precise enough to describe such inflationary scenarios.


Introduction
The inflationary models, which solve successfully the horizon, flatness and relic problems [1][2][3] and generate the primordial density perturbations finally initiating the formation of galaxies and large-scale structure [4][5][6], are the most reasonable models for the evolution of the early Universe.In simplest case inflation is controlled by a single scalar field (the inflaton) with an effective potential that plays a role of the cosmological constant during inflation.
A fundamental step towards the unification of physics at all energy scales could be the possibility to describe the inflation using particle physics models.In numerous models (for a review see [7]) the role of the inflaton has been performed by the Standard Model (SM) Higgs boson [8][9][10][11][12][13][14][15] or a boson in Grand Unified Theories (GUTs) [16,17] or a scalar boson in supersymmetric (SUSY) models [18][19][20] (see [21,22] as reviews).A number of advantages of simplified SUSY GUTs in comparison with nonsupersymmetric GUTs such as naturally longer period of exponential expansion and better stability of the effective Higgs potential with respect to radiative corrections due to cancelation of loop diagrams have been noted quite long ago [23].
Thus, the implementation of inflationary scenario within a well-defined model of particle physics consistent with collider phenomenology where the inflaton is unambiguously identified is a longstanding problem.The only candidate on the role of the inflaton in the SM is the Higgs boson.The Higgs-driven inflation [9][10][11][12][13][14][15] was originally proposed as a single-field model based on the SM in the unitary gauge.This minimal model uses the Higgs isodoublet Φ interaction with gravity of the form ξRΦ † Φ (R is the scalar curvature and ξ is a positive constant).The Higgs-driven inflation leads to the spectral index value n s = 0.967 and the tensor-to-scalar ratio r = 3 • 10 −3 which are in agreement with the Planck Collaboration data [24][25][26][27].However, the effects of Goldstone bosons should be included at an energy scale relevant to inflation in the model which is actually multifield.
For the Higgs-driven inflation it was found [28] that the multifield effects are negligibly small during inflation and do not influence the observable quantities, such as the spectral index of primordial perturbations and the ratio of squared amplitudes for the tensor and the scalar perturbations (tensor-to-scalar ratio).This fact and the known considerations about the need to extend the SM, which could be an effective limit of GUTs, supersymmetry, supergravity or other beyond the SM theories, lead to the belief that inflation that is compatible with recent observations [24][25][26][27] might have been generated by several fields.It has been shown [29,30] that there is a class of inflationary models with two scalar fields non-minimally coupled to gravity that provides good agreement with the Planck data.Hybrid inflation proposed in [31][32][33] which involves the potential of two scalar fields ensures inflationary expansion, explains the observed spectrum of density fluctuations not requiring the unnatural scalar field amplitudes at the Planck scale.At the same time, a greater degree of uncertainty arises in the theory.Identification of the two scalar fields in the framework of a gauge theory model is not simplified, the initial conditions are not unambiguously fixed [34] theoretically in the two-dimensional field space and their tuning is needed to ensure adequate phenomenological consequences.
It should be noted that some tension is observed between the experimental data and predictions of the minimal model.In to order to explain cosmic microwave background observables in the Higgs-driven inflationary scenario the parameter of non-minimal coupling should be very large (ξ ∼ 10 4 ).So large value of ξ is not satisfactory from general theoretical backgrounds because it leads to violation of perturbative unitarity at the scale M P l /ξ which is smaller than the expected inflationary range above the M P l / √ ξ (M P l denotes the Planck mass).In order to restore unitarity above the scale M P l /ξ, "new physics" (new particles interacting with the SM ones) should be introduced which modify the SM Higgs potential.The more serious problem of a large ξ value in the SM is the renormalization group evolution (RGE) of meaningful parameters which demonstrates unsatisfactory matching with the measured Higgs boson and top quark masses [35] as soon as the inflationary range of the order of M P l / √ ξ or above is concerned.In order to reduce the value of ξ, an extremely small value of the effective quartic coupling λ ef f (µ) near the Planck scale is needed.At the Higgs boson mass m h = 125 GeV such value of λ ef f (µ) can be achieved at the top quark mass which is more than 2σ below its observable central value [36][37][38].In this case, the value of ξ necessary for a satisfactory inflationary scenario decreases thus allowing to avoid the problem of perturbative unitarity violation below the inflationary scale.Note that the GUT motivated inflationary model [39] predicts the same order of the parameter for a non-minimal coupling.However, there are cosmological models (see, for example [40]) with the same function of the non-minimal coupling and even polynomial potential of the fourth order that could provide a suitable inflationary parameters at small values of ξ.
Apparent tensions arising in connection with parameter matching of the Higgs-driven inflation model increase the popularity of models with new physics at the TeV and multi-TeV scales.New particles consistent with restrictions on the new physics imposed by the LHC data provide extensive opportunities to improve significantly the Higgs-driven inflationary model.Analogue of this single-field model for the multifield scenarios is based on an observation that redefined fields in the Einstein frame practically coincide with primary fields in the Jordan frame at the low energy scale of the order of superpartners mass scale M SU SY , reproducing the MSSM potential, while at the scale higher than the GUT scale the potential in the redefined fields can be slowly changing respecting the slow-roll approximation of an inflationary scenarios.This observation is sufficiently general.
Recent analyses [41][42][43][44][45][46][47][48] of multifield models showed that unlike the single-field models they generically provide density (entropy) perturbations which can induce the curvature perturbation to evolve beyond the cosmological horizon in the process of inflation [49].Evolution of density perturbations in multifield models should be studied in order to analyze new features in the observables such as non-Gaussianities [50] which are absent in the single-field inflationary models.Deviations of observable power spectrum calculated in multifield models from predictions of the single-field models could take place in the power of one of the three criteria [48,51]: (i) noncanonical kinetic terms; (ii) violation of slow-roll approximation; (iii) nonstandard initial ground state (different from Bunch-Davies vacuum).The main feature of the multifield models which leads to nonstandard primordial spectrum is the ability of trajectories of slow-roll fields to rotate in the field space, that occurs due to the presence of bumps and ridges in the effective multifield potential.When the slow-roll field trajectories turn in the field space, nonstandard contributions to primordial spectra can be amplified enough to be detectable in the microwave background [52][53][54].
In this paper we analyze a multifield extension of the standard Higgs-driven inflation inspired by the MSSM.A few general observations let us make first.In the framework of a sypersymmetric model the natural class of cosmological models are those with local supersymmetry (supergravity models).For the case when interactions at an energy scale below M P l are described by an effective N = 1 supergravity, the general form for the effective potential of scalar fields in the Einstein frame was derived in [55].In the notation of [56] the Lagrangian can be written as where G k = ∂G/∂ϕ k , ĝ is the gauge coupling constant, T are generators of the groups, D µ are covariant with respect to gravity and gauge group, (z i , χ i ) is the chiral supermultiplet.
The Kähler potential G can be written in terms of the function φ which transforms as a real vector superfield and the superpotential g s in the following form: (1.2) Equation (1.1) is a consequence of a Lagrangian written in terms of chiral superfields Φ: here R is the superspace curvature and E is chiral density connected with local superspace basis (see [57]).One can show [58] that minimal coupling to gravity which takes place for φ = z * i z i − 3 in the Kähler potential can be modified to a non-minimal coupling R → R + p( Φ)R instead of the first term in eq.(1.3) by the replacement φ = z * i z i −3−3(p(z)+h.c.)/2 for a given polynomial form p. 1  In the case of the MSSM natural choice is p = ξ Φ1 Φ2 , where Φ1 and Φ2 are chiral Higgs-Higgsino superfields.This choice of p in combination with the general form for the superpotential g s = Λ + µ Φ1 Φ2 (Λ and µ are real constants) leads to problems of achieving a suitable inflationary scenario, see [58,59].For ξ parameter large enough when different regimes for the 'flat direction' tan β parameter are taken, either there is no slow roll or the potential takes negative values.Not referring here to the possibility of the MSSM extension with a gauge singlet (non-minimal MSSM or NMSSM) where unsuitable behavior can be cured, we introduce non-minimal couplings in the non-holomorphic form R (here H 1 and H 2 are SU (2) spinors and R is the Ricci scalar) that have no counterparts in supergravity.So only small electroweak quartic couplings g 1 and g 2 in the D-terms of eq.(1.1) which then appear in the tree-level scalar potential at the SUSY scale provide grounds to speak about 'MSSM-inspired' inflationary scenarios.This sort of model is not a direct extension of models associated with the MSSM which include scalar fields minimally coupled to gravity [62][63][64].The inflaton fields are identified as Higgs sector fields, thus, in this case one is talking about the multifield extension of the SM single-field Higgs inflation.Note that other realizations of the inflationary scenario in the MSSM are possible, when the inflaton is a combination of squark and slepton fields [65], while the process of inflation is controlled by flat directions of the MSSM potential which are lifted by non-renormalizeable superpotential terms and soft supersymmetry breaking terms.It is assumed that D-terms in eq.(1.1) vanish in the hidden sector.A number of other options of the MSSM-inspired inflation can be found in [18-23, 62-64, 66].
The model which is considered in the following sections includes two Higgs doublets coupled with gravity non-minimally.We focus on the two-Higgs doublet MSSM potential in the mass basis of scalar fields that has been analyzed starting from 1975 [67].This potential includes three massless Goldstone bosons and five massive Higgs bosons.Working in the physical gauge, in this paper we do not take Goldstone bosons into account and consider inflationary scenarios that include Higgs bosons only.We show that inflationary scenarios with suitable parameters n s and r are possible at the scale corresponding to the Hubble parameter H ∼ 10 −5 M P l .By this way a MSSM-inspired extension of the original Higgsdriven inflation is constructed.
The structure of the paper is as follows.In section 2 we define the MSSM two-Higgsdoublet potential in the basis of mass eigenstates for the five Higgs bosons at the superparticle mass scale.The mixing angles of the SU (2) field eigenstates are chosen in 1 Taking frequent in the literature point of view that the main qualitative features at the scale of the order of M P l are valid despite the loop effects of gravity, there is an opinion that in a simplest case for a single field any polynomial form p(x) can be adjusted by taking dgs/dx = (3 + p(x) exp(G(x))) /2 (where x is Re z, (zi, χi) is the chiral supermultiplet) with the following extension of a solution in the form of series expansion to complex z [60,61].
the form which is acceptable for the low-energy Higgs phenomenology.In section 3 the MSSM-inspired potential taken in the Jordan frame with the polynomial form of the nonminimal coupling function is transformed to the Einstein frame.Equations of motion in the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric are described in section 4. Numerical integration of the equations of motion with the initial conditions which are adjusted in a way suitable for reproduction of the observable values for the spectral index n s and the tensor-to-scalar ratio r is preformed in section 5.In section 6 we discuss briefly the strong coupling (SC) approximation for the MSSM-inspired potential under consideration.Results are summarized in section 7.

The MSSM-inspired Higgs potential
Two Higgs doublets of the MSSM can be parameterized using the SU (2) states where ω + 1,2 are complex scalar fields, η 1,2 and χ 1,2 are real fields, the vacuum expectation values v 1 and v 2 are usually redefined in (v, tan β) parametrization: Two doublets Φ 1 and Φ 2 can be used to form the SU (2) × U (1) invariant and renormalizable effective potential which breaks gauge symmetry.
The most general two-doublet effective potential can be written as [68]: Let us consider the action in the Jordan frame where g is the determinant of metric tensor gµν , and R is the scalar curvature.The factor in front of the kinetic term is not dependent on fields, so the case of Brans-Dicke gravity-like models are beyond our analysis.However, δ ab in front of the kinetic term is not narrowing the generality of consideration, see details in appendix A. Variation of action with respect to metric tensor gµν and isodoublets Φ a of the fields leads to the following equations ) where a = 1, 2, ∇ µ is a covariant derivative and the d'Alembert operator acting on the scalar fields is denoted by In the following we are using notations and normalization conventions for the potential V (Φ 1 , Φ 2 ) in the generic basis (with λ 6,7 terms) from [69,70], where the mass eigenstates for scalars were constructed.
Note that the potential in eq. ( 2.3) explicitly violates CP invariance if parameters µ 12 , λ 5 , λ 6 , or λ 7 are complex-valued.For simplicity, we are not considering such possibility in the following.At the tree-level λ 5 , λ 6 , and λ 7 are equal to zero in the MSSM two-doublet potential.Nonzero parameters λ 5,6,7 of the effective Higgs potential can be generated by radiative corrections coming from the sector of soft supersymmetry breaking terms, where scalars couple to quark superpartners.To simplify the analysis we will not consider this possibility remaining with the tree-level potential at the M SU SY scale.It is well-known that radiative corrections are large and in the context of this simplification (when the upper limit of the light CP-even state mass m h does not exceed the Z-boson mass m Z = 91.2GeV ) it is impossible to describe adequately the spectrum of Higgs boson masses.However, precision fitting of the collider data is not the primary purpose at this stage of consideration.
The mass basis of scalars is constructed in a standard way.The SU (2) eigenstates (ω ± a , η a and χ a , a = 1, 2) are expressed through mass eigenstates of the Higgs bosons h, H 0 , A and H ± and the Goldstone bosons G 0 , G ± by means of two orthogonal rotations where the rotation matrix Masses of the CP-even scalars h and H 0 are m h and m H 0 , the charged scalar mass is m H ± and the CP-odd scalar mass is m A .At the superpartners mass scale M SU SY the m A and tan β can be chosen as the input parameters which fix the dimension-two parameters µ 2 1 , µ 2 2 and µ 2 12 of the Higgs potential, while the dimensionless factors λ i (i =1,...,7) at the tree level are expressed, using the SU (2) and U (1) gauge couplings g 2 and g 1 , as follows [71,72] (2.9) The dimension-two parameters µ 2 1 , µ 2 2 and µ 2 12 are fixed using the minimization conditions: where m Z = v g 2 1 + g 2 2 /2.Then the potential given by eq. ( 2.3) can be rewritten in the mass basis of scalar bosons, which are massless Goldstone bosons G 0 , G + , G − and massive Higgs bosons h, H 0 , A, H + , H − : where ) Explicit forms of the interaction terms I 3 and I 4 are presented in the appendix B. The mixing angles α and β at the M SU SY scale are connected by the following equation The scalar resonance with mass 125 GeV which is experimentally observed at the LHC [73][74][75] has properties consistent with the SM.However, MSSM identifications are still possible with limited experimental statistics.Experimental data of the LHC Run I demonstrates the SM-like couplings of observed Higgs boson to fermions and vector bosons at the level of statistical significance only on the level slightly better than 2σ [76].In the following consideration, the CP-even state h of the MSSM, when it is overridden mass, which is determined by the radiation corrections from the squark sector, will be identified as the 125 GeV resonance.In the presence of other scalars H 0 , A, H + and H − , which are not experimentally observed, such identification is possible for the two specific features in the MSSM parameter space: (i) the decoupling regime [77] and/or (ii) the alignment limit [78,79].In the decoupling regime masses of scalars H 0 , A, and H ± are very large (they are at multi-TeV scale where also the lightest superpartners can be found), so their contributions to the observables at the top quark scale are strongly suppressed, while in the alignment limit H 0 , A, and H ± are not necessarily extremely heavy.The alignment limit will be used in the following consideration.In this limit β − α ≈ π/2 and the potential in eq. ( 2.3) can be simplified by a special choice of mixing angles α and β.After rotation of scalar isodoublets to so-called Higgs basis [71] and the choice of mixing angles β = π/2 and α = 0, the SU (2) components of isodoublets and the vacuum expectation values are (2.17) So, in the unitary gauge G 0 = G ± = 0, we get and the isodoublet convolutions are given by where and Ω 2 ± = 2 H + H − .The kinetic terms have canonical form (2.21) It follows from eq. (2.10) that µ 2 12 = 0 and the potential in eq.(2.11) becomes where .
The potential in eq.(2.22) qualitatively corresponds to the MSSM potential at the scale M SU SY .It is invariant under two-dimensional rotations in (H 0 , A) space and (H + , H − ) space, what is the consequence of the specific choice of the mixing angles α and β in the alignment limit.This property allows reducing the number of five physically significant fields h, H 0 , A, H + and H − to the three field combinations, h 2 v , Ω 2 0 and Ω 2 ± .Note that at h v = 0 the potential given by eq.(2.22) is invariant under rotations in the fourdimensional field space.Tree-level quartic couplings λ i , eq. (2.9), are expressed through the gauge couplings g 1,2 which are fixed by collider data, since the gauge boson masses at tree level m Z = v g 2 1 + g 2 2 /2, m W = v g 2 /2 and cross sections of W ± , Z production are precisely measured GeV and m W = 80.4 GeV , we obtain g 1 = 0.36 and g 2 = 0.65, which are used in numerical calculations of section 5.

The MSSM-inspired model with non-minimal interaction
Generic action which is dependent on N scalar fields φ I , I = 1, ..., N with the standard kinetic term and non-minimal coupling to gravity can be written as where tilde denominates the metric tensor and curvature in the Jordan frame.In our case V (φ I ) depends on five real scalar fields This action can be transformed to the following action in the Einstein frame [80] (see also [28,48]): where the reduced Planck mass M P l ≡ 1/ √ 8πG, f ,I = ∂f /∂φ I .Metric tensors in the Jordan and the Einstein frames are related by the equation In the single-field Higgs-driven inflation the function f has been chosen as a sum of the Hilbert-Einstein term and the induced gravity term.We choose the function f in an analogous form: where ξ 1 and ξ 2 are positive dimensionless constants.This form of function f follows from the requirement of renormalizability for quantum field theories in curved space-time [81][82][83][84], where non-minimal couplings appear as renormalization counterterms for scalar fields.We also assume that vacuum expectation values for scalar fields are negligibly small in comparison with M P l .Note that non-minimal interaction in the form of eq.(3.4) was considered [85] in the framework of the (nonsuperymmetric) two-Higgs-doublet model, when the boundary condition eq.(2.9) is not used and the Higgs potential includes seven quartic couplings.Arbitrariness of the choice of λ i is constrained imposing exact or approximate Z 2 symmetry (discrete symmetry whose breaking results in the appearance of the axion) on the generic two-Higgsdoublet potential which takes a specific functional form different from the 'MSSM-inspired', eq.(2.22).It is assumed that the Higgs doublets Φ 1 and Φ 2 in this simplified potential are , what happens if the fields ω 1,2 , η 1,2 and χ 1,2 in eq.(2.1) are taken to be zero, so the MSSM mass eigenstates h, H, A and H ± are not specified.In our case the function f depends on the five scalar fields:

Properties of the equations of motion in the FLRW metric
Let us consider a spatially flat FLRW universe with metric interval where a(t) is the scale factor.Varying the action in eq.(3.3) with respect to g µν and fields we get the following equations for the FLRW metric where the Hubble parameter H = ȧ/a, σ2 = G IJ φI φJ , and dots mean the time derivatives.
Field equations have the following form [48] φI + 3H φI + Γ where Γ I JK is the Christoffel symbol for the field-space manifold, calculated in terms of G IJ , W ,K = ∂W/∂φ K .Hereafter, primes denote derivatives with respect to the fields.Due to the relationship of inflationary evolution in the Jordan and the Einstein frames, eqs.(4.1) and (4.2) are equivalent to eqs.(2.5) and (2.6) after transformation of the latter to Einstein frame.
During inflation the Hubble parameter is positive and the scalar factor is a monotonically increasing function.To describe the evolution of scalar fields during inflation we use the number of e-foldings N e = ln(a/a e ), where a e is the value of the scalar factor at the end of inflation, as a new measure of time.The notation N * e = −N e will be also used for convenience.
Using d/dt = H d/dN e one can write eqs.(4.1) and (4.2) in the form where (σ ) 2 = H 2 ( σ) 2 .After substitution of eqs.(4.3) and (4.4) a system defined by eqs.(4.5) and (4.6) includes ten first order equations which are suitable for numeric integration.Integration was performed by means of built-in subroutines of several computer algebra systems with cross-checks of results.Note that so far in this section and in the previous section 3 we have not made any approximations.In order to calculate the observables, spectral index n s and tensor-to-scalar ratio r, slow-roll parameters are introduced analogously to the single-field inflation where σ I = φI / σ is the unit vector in the field space and D denotes a covariant derivative with respect to the field-space metric, D I φ J = ∂ I φ J + Γ J IK φ K .Then the spectral index n s and tensor-to-scalar ratio r at the time when a characteristic scale (50-65 e-foldings before the end of inflation) is of the order of the Hubble radius in the course of inflation, can be calculated using the single-field equations valid to lowest order in slow-roll parameters [48,49] n s = 1 − 6 + 2η σσ , r = 16 .(4.9)

Numerical solutions of the equations of motion
The isosurfaces for the potential W (h v , Ω 0 , Ω ± ) (one of the three variables is fixed) are shown in figure 1.At the fixed value of h v , Ω 0 or Ω ± of the order of 0.1 (in Planck units) the saddle configuration of the surface is observed, as shown in figure 1(a).A characteristic feature of W which demonstrates ridges and gullies is shown in figure 1(b).In the gullies evolution of the field system looks as an infinite expansion at the constant Hubble parameter.One can see that the slow-roll inflation is possible if the initial condition for h v or Ω 0 is chosen in the vicinity of zero, which is equivalent to four nonzero values of the fields (A, H 0 , H ± ) or three nonzero values of the fields (h v , H ± ).Initial conditions for a number of successful inflationary scenarios of this sort are presented in table 1.The evolution of fields superimposed on the Einstein-frame potential for the inflationary scenarios A and B, see table 1, is shown in figure 2, where the dashed fragments of field trajectories correspond to the inflationary stage when 0 ≤ N * e ≤ 65.If we assume that the number of e-foldings during inflation N * e = 65, then we get the initial conditions for inflationary trajectories presented in tables 2 and 3.For scenarios A 1 , A 2 , A 3 , B 3 , and B 4 (see table 4) the value of the Hubble parameter in the beginning of inflation is H init < 3.6 • 10 −5 M P l .Note that this value of the Hubble parameter is found to be in good agreement with the observational data [26,27].One can see that for the type A inflationary scenarios the field system rolls slowly down to the potential minimum (see also figures 3(a) and 3(b)), while for the B type inflationary scenarios, except B 4 , all nonzero fields demonstrate rapidly damped oscillations going to zero h v for the number of e-foldings before the end of inflation N * e 65, see figures 3(c) and 3(d).At the same time, significant nonzero value of h v in the initial field configuration is suitable for inflation in the case B 4 when ξ 1 = ξ 2 .Note that inflationary scenarios with initial conditions denoted by A and B in tables 2, 3 demonstrate remarkable stability of slow-roll parameters , η σσ and observables r and n s .In different cases with the Hubble parameter H ∼ 10 −5 M P l the values of n s and r coincide up to three digits.Such "attractor behavior" when over a wide range of initial conditions the system evolves along the same trajectory in the course of inflation is known for single-field models [86], but it is not an obvious observation, generally speaking, for multifield models.In this sense the phenomenological stability inherent to the single-field Higgs inflation is preserved for the multifield MSSM-inspired model under consideration.
The problem of perturbative unitarity violation at a large values of ξ parameters mentioned in the Introduction may persist in the MSSM although an order of magnitude smaller values of ξ appear in comparison with the SM Higgs inflation (except A 4 scenario, see Table 1).While in the SM for the Higgs inflation a simple unitarity bound can be derived E < M P l /ξ on the general basis of power-counting formalism for effective theory (for example, [87]), in the MSSM-inspired models with several fields such a simple criteria is not reliable and the situation with partial wave unitarity is much more difficult.Recent analysis [88] for the case of a general two-Higgs-doublet model without any discrete symmetry imposed on the scalar potential leads to non-trivial constraints on the masses and mixings which may depend on the scenario of new physics at a high energy scale.-14 -It follows configurations shown in tables 2 and 3 that quite different initial values of scalar fields and parameters ξ 1 and ξ 2 appear in all cases in combination with a very small value of h v .In all cases, but A 4 , the inflationary parameters practically coincide, see figure 4 and table 4. In this section, we show that such a pattern can be explained in the framework of the so-called "strong coupling approximation".It has been shown in a large number of analyses [89][90][91][92][93][94][95][96][97] that there are several classes of the single-field inflationary models such that within a given class all models predict the same values of observable parameters n s and r in the leading 1/N e approximation.These classes are known as cosmological attractors.
A similar analysis of two-field inflationary models has been made in [29,98].For stringmotivated supergravity theory in which both the field-space metric and the potential usually have poles at the same points, the inflationary dynamic and the corresponding attractor have been studied [99].The idea of a cosmological attractor is based on an observation that the kinetic term in Jordan frame practically does not affect the slow-roll parameters if the "strong coupling regime" is respected during inflation.In the case of multifield models the field system is in the SC regime if the following inequality is respected: In the approximation of eq. ( 6.1) the action given in eq. ( 3.3) can be written as and rewritten in the equivalent form where f 0 is a positive constant with the same dimension as f .The role of inflaton in the strong coupling approximation is performed by the "effective field" in terms of which the action S E includes the standard kinetic term of Θ and does not include kinetic terms of any other scalar fields which can be interpreted as model parameters.This circumstance allows one to calculate the inflationary parameters in the SC approximation using the single-field model.If we adjust Θ in such a way that Θ = 0 corresponds to Ω 0 = 0 and Ω ± = 0, then f 0 = M 2 P l /2.The single-field model consistent with the above-mentioned scenarios A and B (see tables 2 and 3) can be easily defined.In the scenario A we set φ 5 = h v = 0 during inflation, while in the scenario B (except B 4 ) one can observe that inflation starts when We do not consider the case B here.In all other cases we can neglect h v and write the potential in the form The function f is approximated by and thereby so the Einstein frame potential can be written as follows Using m2 2 ξ 1 M 2 P l ν 1 we get from eq. (6.6) The slow-roll parameters are With these analytic expressions for the slow-roll parameters in the SC approximation the inflationary parameters can be easily calculated.It is convenient to express the inflationary parameters as a functions of f sc Straightforward numerical cross-checks demonstrate that the ratio is less than 7×10 −5 in the scenario A and 2×10 −4 in the scenario B, so the SC approximation is meaningful.It is demonstrated in table 5 that in all cases the values of inflationary parameters r and n s calculated using eq.(6.8) are close to the parameter values that have been found numerically in section 5. Note in this connection that the primordial non-Gaussianities which do not arise in the single-field inflationary models should be very small in the case under consideration as soon as the reduction to a single-field scenario is precise enough.It should be mentioned that f in /M 2 P l close to 44 is not a sufficient condition for an inflationary scenario with suitable values of n s and r.For a large number of initial data with such values of f in /M

Summary
In this paper, we constructed a MSSM-inspired extension of the original Higgs-driven inflation [9][10][11][12][13][14][15] using the two-Higgs-doublet potential of the MSSM which is simplified in a way suitable for calculation of transparent symbolic and numerical results for the main observables, the spectral index n s and the tensor-to-scalar ratio r.The shape of the MSSM potential surface in the Einstein frame where ridges and bumps influence the trajectory in the fields space is different from the usual form in models of hybrid inflation.The model under consideration incorporates multiple non-minimally coupled scalar fields and non-canonical kinetic terms in the Einstein frame which are induced by the curvature of the field-space manifold.For these reasons, the evolution of fields is generically different from slow-roll, at least for some time interval during inflation.The analysis of the background inflation dynamics demonstrated that after setting up the initial conditions for the five-dimensional field configuration such simplified MSSMinspired model successfully describes the Higgs-driven inflation consistently with the observations of the Planck and BICEP2 collaborations.Two types of consistent inflationary scenarios are found with the initial conditions denoted as A and B, see table 1, which demonstrate the remarkable stability of the observables with respect to the shift of the initial field system configuration.The main difference between these two cases is the presence of rapid field oscillations in the initial phase of case B before the beginning of inflation, while oscillations are absent in case A. During the period of cosmological evolution which determines the observables, h v field is negligibly small so the value of ξ 2 parameter practically does not influence the result and in the MSSM-inspired model degenerate values of ξ 1 and ξ 2 are always meaningful.Inflation occurs for field values much smaller than the Planck scale, although no suitable expansion scenario was found for initial state when h v , Ω 0 and Ω ± are very small at the same time.In all cases trajectories of the system do not turn steeply in the field space, so specific features of the potential like bumps and ridges trajectory rapidly (after ∼ 0.05 e-foldings) rolls into the gully h 2 v = Ω 2 0 (see figure 1(b)).This direction is not absolutely flat (the case when critical points are degenerate and not isolated [100]), but so close to flat that cannot be analysed by numerical methods.Simple estimate with fin = 45 and ξ1 + ξ2 ∼ 2 • 10 3 gives an extremely long slow-roll with the number of e-foldings of the order of 10 12 .
are not expected to induce primordial non-Gaussianities with a magnitude large enough to be detectable in the cosmic microwave background.
Multifield model under consideration demonstrates rather strong attractor behavior and can be to the single-field model with the effective inflaton field defined by eq. ( 6.3).Such models share very close results for the spectral index and the tensor-toscalar ratio in combination with negligible non-Gaussianity, which are in good agreement with the latest experimental data.
In conclusion let us also note that an important point beyond our analysis is the stability of results with respect to radiative corrections.The flatness of the effective potential in the region of the field amplitudes of the order of M P l is an essential property for a suitable slow-roll.While the quantum gravity corrections are expected to be rather small of the order of V /M 4 P l ∼ g 2 p /ξ 2 , the corrections induced by the SM fields and the superpartner fields involved in the F and D soft supersymmetry breaking Lagrangian terms require careful analysis which is dependent on the MSSM parametric scenario under consideration.For example, in the "natural MSSM scenario" which is used for LHC analyses the superpartners of quarks show up at the multi-TeV scale, while gauginos decouple.At the one-loop resummed level the superpartner threshold corrections to the two-doublet MSSM potential are expressed by Coleman-Weinberg terms ∆V = 1/(64π 2 )Sp[(V (φ)) 2 (log(V (φ)/µ 2 ) − 3/2)], where second derivatives taken at the local minimum of Higgs potential are equal to masses of scalars.Nontrivial significant contributions are provided in the higher orders of perturbation theory by nonrenormalizable operators [101].Fermionic and bosonic loops give contributions of different signs which could partially compensate each other.Contributions of the SM vector bosons and fermions are smaller than the MSSM ones because of small gauge and Yukawa couplings, so main corrections from the third generation of quark superpartners interacting with Higgs isodoublets must not spoil a small slope of the potential.Important correction can be provided also by the renormalization group (RG) evolution of ξ non-minimal couplings from the top quark scale to the M P l scale.RG evolution gives at least a factor of two for the value of ξ in the framework of SM Higgs-driven inflation, but moderate changes of the order of ten percent in the inflationary region.Models which are described by the RG-improved effective action [102][103][104] should provide an improved precision for observables.Careful MSSM evaluations which are beyond our analysis are appropriate in order to ensure stability of results.

A General action for non-minimal Higgs interactions in the MSSM
In the general case one can write the action for non-minimal interaction of the MSSM Higgs doublets with gravity in the (here we redefine f ( Φ1 , Φ2 ) = One can find some (may be non-unitary) transformation Φa → Φ a = U ab Φb to diagonalize G IJ → δ IJ , so U ab is After such transformation the action can be written as Thus one can always start with the action in the form (A.2) or (2.4) without loss of generality.

B Higgs potential in the mass basis
The potential given in eq. ( 2.3) can be written in terms of the mass eigenstates, which are massless Goldstone fields G 0 , G + , G − and massive Higgs bosons h, H 0 , A, H + , H − ,3 in the following form ), tensor-to-scalar ratio r and spectral index n s for successful inflationary scenarios at N * e = 65, m A = 200 GeV .

Table 1 .
Initial conditions (in units of M P l ) for trajectories with successful inflationary scenarios, CP-odd Higgs boson mass m A = 200 GeV .The SU (2) and U (1) gauge couplings are g 1 = 0.36 and g 2 = 0.65.

Table 2 .
Initial conditions (fields in units of M P l ) at N * e = 65 in the scenarios of type A (see table1).

Table 3 .
Initial conditions (fields in units of M P l ) at N * e = 65 in the scenarios of type B (see table1).

Table 4 .
The Hubble parameter H (in units of M P l

Table 5 .
2P l , but beyond the abovementioned A and B type scenarios, acceptable inflationary evolution is not observed 2 .The inflationary parameters in the strong coupling approximation calculated at h v = 0, m Z = 0 and m A = 0.