text:
# Affleck-Dine leptogenesis via right-handed sneutrino fields in a supersymmetric hybrid inflation model

## I introduction

The hybrid inflation model  provides an attractive possibility for solving a range of cosmological problems . It has an unique dynamical feature due to interplay between two scalar fields, inflaton $`\sigma `$ and an auxiliary field, so called orthogonal scalar $`\varphi `$ which is actually responsible for providing the potential energy throughout the inflationary phase. Quite remarkably, a proper form of the inflaton potential can be naturally realized in the context of supersymmetric theories . The simplest supersymmetric model for hybrid inflation is based on the following superpotential including the inflaton superfield $`S`$ and the auxiliary one $`\mathrm{\Phi }`$ :

$$W=\mathrm{\Lambda }^2S+\lambda S\mathrm{\Phi }^2,$$
(1)
where $`\mathrm{\Lambda }`$ and $`\lambda `$ are the model parameters. The latter are constrained by COBE normalization for the density perturbations as $`\mathrm{\Lambda }6.510^{16}ϵ^{1/4}`$ GeV , or equivalently

$$\mathrm{\Lambda }1.310^{15}|\eta |\lambda ^{1/2}\mathrm{GeV},$$
(2)
where $`ϵ,\eta 1`$ are the slow roll parameters.

The scalar components of these superfields have a potential

$$V(\sigma ,\varphi )=\mathrm{\Lambda }^4\lambda \mathrm{\Lambda }^2\varphi ^2+\frac{\lambda ^2}{4}\varphi ^4+\lambda ^2\sigma ^2\varphi ^2.$$
(3)
which has a global (supersymmetric) minima at $`\sigma =0`$ and $`\varphi ^2=2\mathrm{\Lambda }^2/\lambda `$. However, for $`\varphi =0`$ and $`\sigma >\sigma _c=\mathrm{\Lambda }/\sqrt{\lambda }`$, the potential is flat along the $`\sigma `$ axis with non-zero vacuum energy $`V(\sigma ,0)=\mathrm{\Lambda }^4`$. This flat direction is lifted by radiative corrections resulting from the supersymmetry breaking and perhaps also by some supergravity corrections  that give an appropriate slope to the inflaton potential on which $`\sigma `$ can slowly roll down and produce inflation. During inflation the orthogonal field is held in a false vacuum $`\varphi =0`$, while the inflaton must evolve from an initial value close to the Planck scale $`M_\mathrm{P}=1.2\times 10^{19}`$ GeV. The inflationary phase ends when the inflaton field approaches the critical value $`\sigma _c`$, after which both fields $`\sigma `$ and $`\varphi `$ begin to oscillate near their global minima and reheat the universe.

However, as any inflationary paradigm, hybrid inflation models also suffer from the initial condition problems. A generic problem concerns the difficulty to inflate an arbitrary space-time patch. In order to trigger inflation, the inflaton field should be extremely homogeneous in an initial patch of the Universe at scales larger than the Hubble radius of this epoch . This becomes a challenging task, since we know from the observations that the present inhomogeneity is primordial in nature and necessarily ties its origin with inflation. This obstacle can be evaded for a chaotic initial condition provided the fields take values close to the Planck scale.

In the case of hybrid inflation this generic problem is aggravated by the fact that strong fine tuning is required also for the initial values for the fields . This consists in the following. At the Planck epoch, one can expect that all scalar fields, not only the inflaton, take their initial values close to the $`M_\mathrm{P}`$. However, if the orthogonal field has a value $`\varphi M_\mathrm{P}`$, it does not provide the right condition for the onset of inflation. The reason is simply the following. While the field $`\sigma M_\mathrm{P}`$ provides large effective mass term to $`\varphi `$ through the last term in (3), the value $`\varphi M_\mathrm{P}`$ in turn would induce large mass term to $`\sigma `$. In this case both the auxiliary field and inflaton would merely oscillate at their local minima without inflating the Universe. In order to trigger inflation, the initial field configuration should be settled down along the $`\sigma `$-valley, with $`\sigma `$ having the value of order $`M_P`$ while $`\varphi `$ must be close to zero with extreme accuracy, $`\varphi <10^5M_\mathrm{P}`$ . This is precisely the point which is not very natural; why the orthogonal field should start so close to the false minimum of its potential? In other terms, in order to inflate the initial patch of space after the Planck epoch, it should not only be strongly homogeneous at distances much larger than the corresponding Hubble radius, but it also should have energy density $`\mathrm{\Lambda }^4`$, about 15 orders of magnitude smaller than $`M_\mathrm{P}`$.

Certainly, such stringent initial conditions can be accepted on purely anthropic grounds. However, it is always appealing to obtain a more natural solution of the problem. It has been shown recently in Ref. , that the supersymmetric hybrid model can allow a solution to this fine-tuning problem quite elegantly. The solution is simple as it assumes that all the fields take their initial value close to the Planck scale. However, the orthogonal field $`\varphi `$ must decay into the quanta of some extra fields which do not interact with the inflaton, and thereby do not contribute to the curvature of the inflaton potential. This can be easily achieved by adding to the superpotential (1) a term $`\kappa \mathrm{\Phi }\mathrm{\Psi }^2`$, where $`\mathrm{\Psi }`$ is the extra superfield. The field $`\varphi `$ oscillating near its false minimum, which one can consider as a matter dominated phase before the onset of inflation, decays into $`\mathrm{\Psi }`$ and settles down to zero quickly enough such that the value of the inflaton remains $`\sigma M_\mathrm{P}`$ which allows inflation to begin.

The graceful exit of inflation take place when the inflaton crosses the critical point on the valley which unfolds the strong positive curvature of the false vacuum of an auxiliary field to a smooth potential with a negative curvature which allows this field to roll down to its true minimum. When this happens both the inflaton and the auxiliary field begin oscillations around their respective global minima. One of the interesting aftermath of any inflationary dynamics is that the fields which take part in inflating the Universe, produce entropy and a thermal bath during their oscillations, the process known as the reheating of the Universe.

Shortly after, or possibly during reheating another crucial event, must occur: baryogenesis. As far as any preexisting baryon asymmetry has been exponentially diluted during inflation, one has to find a valid mechanism for generating the observed baryon asymmetry.

By now there are plenty of models for baryogenesis. One of the attractive possibilities would be to produce the baryon asymmetry just during the reheating, in an out of equilibrium decay of the inflaton particles themselves. Notice, that the departure from the thermal equilibrium is one of the three basic criterion for baryogenesis other than $`CP`$ and $`B`$ (or rather $`BL`$) violation . It is not easy, however, to naturally realize this possibility in the context of supersymmetric hybrid inflation.

Another interesting idea is related to leptogenesis , where first the non-zero lepton number $`L`$ is generated in some out-of-equilibrium processes, usually assumed as a right-handed neutrino , and it is partially reprocessed into baryon number $`B`$ via $`B+L`$ violating sphaleron effects which however preserve $`BL`$ .

Yet another approach, known as Affleck-Dine mechanism , takes advantage of the flat directions in supersymmetry which carry $`BL`$ global charge. The non-zero baryon number density can be produced in decay of these modes at later times in the Universe evolution.

Baryogenesis mechanism which we propose in this paper is a blend of these three ideas discussed above, and it can be naturally realized in the context of the supersymmetric hybrid inflation model .

Namely, a new step which we make here is to identify the extra superfield $`\mathrm{\Psi }`$, needed for alleviating the problem of initial conditions, as a right-handed (RH) neutrino. First of all, this proposal can be motivated by the observation that the inflationary energy scale $`\mathrm{\Lambda }10^{15}`$ GeV is also of interest as the RH neutrino mass scale in the context of the see saw mechanism . Indeed, in the global minimum the field $`\varphi `$ receives a vacuum expectation value (VEV) and thus the coupling $`\mathrm{\Phi }\mathrm{\Psi }^2`$ induces the mass $`M_\mathrm{\Psi }\mathrm{\Lambda }`$.

Second, in this case is the scalar component of this superfield, the RH sneutrino field $`\stackrel{~}{\mathrm{\Phi }}`$, carries the lepton number. The associated exact global symmetry, which is actually $`U(1)_{BL}`$, is violated by the VEV of $`\varphi `$. However, during the inflation the system is trapped in a false vacuum with $`\varphi =0`$, where the $`U(1)_{BL}`$ symmetry is restored and $`\stackrel{~}{\mathrm{\Psi }}`$ behaves as massless field. Along with $`\sigma `$ and $`\varphi `$, also $`\stackrel{~}{\mathrm{\Psi }}`$ must have initial value $`M_\mathrm{P}`$ and before inflation it oscillates around zero. However, once the field $`\varphi `$ decays, $`\stackrel{~}{\mathrm{\Psi }}`$ become essentially massless mode. Therefore, during the de Sitter phase its evolution is slow and finally at the end of inflation it has non-zero value of the order of the Hubble parameter $`H\mathrm{\Lambda }^2/M_\mathrm{P}`$. In postinflationary epoch this field becomes massive and starts to oscillate near the origin. As we show below, their evolution in this epoch epoch could generate, due to dynamical breaking of the associated global $`U(1)`$ symmetry, an adequate amount of $`BL`$ in the Universe, which can be converted into the baryon asymmetry via sphaleron transitions.

The paper is organized as follows. We begin with presenting our model. In next section we discuss the evolution patterns of the scalar fields before, during, and, after inflation and the reheating of the Universe. Then we calculate the lepton asymmetry generated by the RH sneutrino field and discuss see-saw phenomenology for neutrino masses. Finally, we conclude with a brief discussion of our findings.

## II The model

From the point of view of the particle physics, our model is nothing but the MSSM including the standard fermion superfields: leptons $`l=(\nu ,e)`$, $`e^c`$ and quarks (which we do not write explicitly), two Higgs doublets $`\phi _{1,2}`$, and in addition an extra superfield $`\mathrm{\Psi }`$ to be identified with the right-handed neutrino, which is a gauge singlet of the standard model. For simplicity, let us consider only one fermion generation.

The charged lepton gets mass from the Yukawa term $`hle^c\phi _1`$, while the similar coupling $`gl\mathrm{\Psi }\phi _2`$ induces the Dirac mass for the neutrino component. For implementing the hybrid inflation scenario, we consider the superpotential of the following form :

$$W=\mathrm{\Lambda }^2S+\lambda S\mathrm{\Phi }^2+\kappa \mathrm{\Phi }\mathrm{\Psi }^2,$$
(4)
which is the simplest modification of the original one (1). It is essential the term $`\kappa \mathrm{\Phi }\mathrm{\Psi }^2`$ which communicates the standard particle sector to the fields producing inflation. The above superpotential as well as the Yukawa terms respect the global $`R`$-symmetry with the fermion superfields $`\mathrm{\Psi },l,e^c,\mathrm{}`$ carrying the $`R`$ charge $`1/2`$ and the Higgs ones $`S`$ and $`(\mathrm{\Phi },\phi _{1,2})`$ carrying charges $`1`$ and $`0`$ respectively. Notice that this $`R`$ symmetry forbids the R-violating terms $`l\phi _2`$, $`lle^c`$ etc. - in other words, the theory has an automatic matter parity under which the all fermion superfields change sign while the Higgs ones remain invariant. Interestingly, $`R`$ symmetry forbids also the supersymmetric mass term $`\mu \phi _1\phi _2`$ which seems to be a good starting point for solving the hierarchy problem. Needless to say, we assume that the Higgsino mass can emerge in some ways in consequence of the supersymmetry breaking. Moreover, in global supersymmetry limit, the only mass scale in the theory is $`\mathrm{\Lambda }`$ in (4).

Let us take the scalar field components in a form:

$$S=\frac{\sigma }{\sqrt{2}},\mathrm{\Phi }=\frac{\varphi }{\sqrt{2}},\stackrel{~}{\mathrm{\Psi }}=\frac{\psi _1+i\psi _2}{\sqrt{2}},$$
(5)
where $`R`$-symmetry transformation has been used to make $`S`$ field real, and the imaginary part of $`\mathrm{\Phi }`$ has been put to zero for simplicity. (Sometimes it is convenient to present the last field in terms of its modulus and phase: $`\stackrel{~}{\mathrm{\Psi }}=(\rho /\sqrt{2})\mathrm{exp}(i\delta )`$.) Their potential reads:

$`V`$ $`=`$ $`V(\sigma ,\varphi )+\kappa ^2\varphi ^2(\psi _1^2+\psi _2^2)+{\displaystyle \frac{\kappa ^2}{4}}(\psi _1^2+\psi _2^2)^2`$ (7)
$`+\kappa \lambda \sigma \varphi (\psi _1^2\psi _2^2),`$
where $`V(\sigma ,\varphi )`$ is given by (3). It has a global minimum with a non-zero VEV $`\mathrm{\Phi }`$:

$$\sigma =0,\psi _{1,2}=0,\varphi =\varphi _0=\sqrt{2}\lambda ^{1/2}\mathrm{\Lambda },$$
(8)
At this minimum the superfields $`S`$ and $`\mathrm{\Phi }^{}=\mathrm{\Phi }\mathrm{\Phi }`$ have equal masses

$$M_\sigma =M_\varphi =2\lambda ^{1/2}\mathrm{\Lambda }\eta \times 2.510^{15}\mathrm{GeV},$$
(9)
while the last coupling in (4) induces the mass of $`\mathrm{\Psi }`$:

$$M_\mathrm{\Psi }=\frac{2\kappa \mathrm{\Lambda }}{\lambda ^{1/2}}\frac{\kappa \eta }{\lambda }\times 2.510^{15}\mathrm{GeV},$$
(10)
where (2) has been used for numerical estimations. Interestingly, this is just the right range for the RH neutrino mass. Then, by means of the seesaw mechanism, the ordinary neutrino gets small Majorana mass

$$m_\nu =\frac{g^2\phi _2^2}{M_\mathrm{\Psi }}\frac{g^2\lambda }{\kappa \eta }\times 1.210^2\mathrm{eV},$$
(11)
where for numerical estimation we have taken $`\phi _2170`$ GeV (i.e. not very small $`\mathrm{tan}\beta `$).

The superpotential (4) provides the same pattern for the inflation as the original one (1). For the field values $`\sigma >\sigma _c=\lambda ^{1/2}\mathrm{\Lambda }`$, the potential (7) has a flat direction along the $`\sigma `$ axis, in which the orthogonal field is trapped in false vacuum with the energy $`V(\sigma ,0)=\mathrm{\Lambda }^4`$. This flat direction is lifted by radiative corrections resulting from the supersymmetry breaking by the non-zero vacuum energy, and perhaps also by some other supergravity corrections . The precise form of these corrections is not important, and we can simply assume that they result in an effective $`\sigma `$ dependent potential, e.g. in the form of mass term, $`\overline{m}^2\sigma ^2`$, $`m\mathrm{\Lambda }`$, that gives an appropriate slope to the inflaton potential on which $`\sigma `$ can slowly roll down and produce inflation. For achieving this, the parameters

$`ϵ={\displaystyle \frac{M_\mathrm{P}^2}{16\pi }}\left({\displaystyle \frac{V^{}}{V}}\right)^2{\displaystyle \frac{\lambda ^2M_\mathrm{P}^2\sigma ^2}{\pi \mathrm{\Lambda }^4}}\left({\displaystyle \frac{\overline{m}}{M_\sigma }}\right)^4,`$ (12)
$`\eta ={\displaystyle \frac{M_\mathrm{P}^2}{8\pi }}\left({\displaystyle \frac{V^{\prime \prime }}{V}}\right){\displaystyle \frac{\lambda M_\mathrm{P}^2}{2\pi \mathrm{\Lambda }^2}}\left({\displaystyle \frac{\overline{m}}{M_\sigma }}\right)^2,`$ (13)
have to be much smaller than one. For producing inflation, the orthogonal field should be settled in a false vacuum $`\varphi =0`$ while the inflaton evolves from an initial value $`\sigma M_\mathrm{P}`$. The inflationary phase ends up when the inflaton approaches the critical value $`\sigma _c`$, after which both fields $`\sigma `$ and $`\varphi `$ begin to oscillate near their global minima and reheat the universe.

Let us now turn to the role of the superfield $`\mathrm{\Psi }`$. In the global minimum (8), its fermion component, the RH neutrino, gets large Majorana mass and thus violates the lepton number conservation. However, the Lagrangian of the scalar component $`\stackrel{~}{\mathrm{\Psi }}`$, the RH sneutrino, maintains the global $`U(1)`$ symmetry. Indeed, since $`\sigma =0`$, the scalar potential (7) becomes a function of the field $`\stackrel{~}{\mathrm{\Psi }}`$ modulus $`\rho =(\psi _1^1+\psi _2^2)^{1/2}`$ and does not depend on its phase $`\delta `$. Needless to say, that the Yukawa coupling to lepton and Higgsino, $`\stackrel{~}{\mathrm{\Psi }}l\stackrel{~}{\phi }_2+\mathrm{h}.\mathrm{c}.`$, are conserving the lepton number.

On the other hand, in the false minimum $`\varphi =0`$, the superfield $`\mathrm{\Psi }`$ is massless and completely decoupled form the inflaton $`\sigma `$. Now the scalar potential of $`\stackrel{~}{\mathrm{\Psi }}`$ contains only the last term $`\kappa ^2\rho ^4`$ in (7). Thus, the global $`U(1)`$ symmetry associated with the lepton number is restored.

Therefore, the $`U(1)`$ symmetry in the potential of RH sneutrino fields $`\stackrel{~}{\mathrm{\Psi }}`$ can be violated only by the last term in (7), when that both $`\varphi `$ and $`\sigma `$ have non-zero values. This can occur only during the epoch when the background fields $`\varphi ,\sigma `$ oscillate near their global minima. As we show below, at this phase the non-zero lepton number can be generated in the classical motion of the RH sneutrino fields, which will be transferred to the standard leptons via their decay $`\stackrel{~}{\mathrm{\Psi }}l\stackrel{~}{\phi }_2`$.

## III Dynamics of the fields

In this section we shall describe the evolution of the fields in our model. Starting our consideration from an earliest time when the Planck era ends and classical general relativity becomes applicable, we assume that at this moment all classical fields have initial values $`M_\mathrm{P}`$ and hence the energy density of the Universe is $`M_\mathrm{P}^4`$, quite a generic situation. We study the scalar field dynamics before the inflation onset, during exponential expansion, their postinfaltionary oscillations and, finally, we discuss the reheating of the Universe. This will provide us all necessary tools to calculate the lepton asymmetry.

### A Pre-inflation

The detailed analysis of the scalar field dynamics before inflation in the model with the superpotential (4) has been performed in Ref. . Here we briefly recall its main features.

All scalar fields $`\sigma `$, $`\varphi `$ and $`\psi _{1,2}`$ have initial values of close to the Planck scale, and with the Hubble expansion they evolve down due to interaction terms in (7). Therefore, the cosmological energy density, initially $`M_\mathrm{P}^4`$, decreases with the Hubble expansion and its evolution settles quickly in a pattern when the energy of the system is dominated by the regular oscillations of $`\varphi `$ around zero. From this moment the Universe expands as in a matter dominated era.

Along the flat direction the effective mass for $`\varphi `$ field is quite large, $`M_{\varphi ,\mathrm{eff}}^2=2\lambda ^2(\sigma ^2\sigma _c^2)`$, and it vanishes only when $`\sigma =\sigma _c`$. On the other hand, once $`\varphi `$ has a large amplitude, it induces the large curvature for the inflaton field $`\sigma `$ and the latter starts to quickly roll down to the origin. Therefore, in order to produce inflation, the amplitude of $`\varphi `$ must settle to zero before $`\sigma `$ will manage to substantially decrease its initial value $`M_\mathrm{P}`$. The friction provided by the Universe expansion itself cannot really help, since the Hubble parameter in this epoch is actually smaller than $`M_{\varphi ,\mathrm{eff}}`$ and thus the field $`\varphi `$ suffers many oscillations during the Hubble time.

However, the problem can be with help of the coupling $`\kappa \mathrm{\Phi }\mathrm{\Psi }^2`$, which permits the oscillating field $`\varphi `$ to decay into $`\mathrm{\Psi }`$ particles. The decay rate $`\mathrm{\Gamma }(\varphi \mathrm{\Psi }\mathrm{\Psi })=(\kappa ^2/8\pi )M_{\varphi ,\mathrm{eff}}`$ can be large as far as the coupling constant is reasonably large, $`\kappa >0.1`$ or so, and $`\sigma M_\mathrm{P}`$. In addition, as soon as $`\varphi `$ becomes zero, $`\sigma `$ stops to feel the classical fields $`\psi _{1,2}`$ as well as their quanta produced in the decay, large curvature of the inflaton potential disappears and its value freezes at values $`\sigma M_\mathrm{P}`$. By the time $`t1/\mathrm{\Gamma }`$ the cosmological energy density becomes dominated by the relativistic $`\mathrm{\Psi }`$ particles, and it fastly decreases with the universe expansion until it becomes dominated by the false vacuum energy. Since this moment Universe starts to expand exponentially while the inflaton proceeds its slow roll due to small curvature term.

The decay of $`\varphi 2\mathrm{\Psi }`$ can be understood as a process of preinflationary reheating after which the system from a generic initial state, with all fields having magnitudes $`M_\mathrm{P}`$, settles in false vacuum, Therefore, it provides a simple mechanism to obtain the correct initial condition for the inflation to take place. Once inflation commences it dilutes the produced $`\mathrm{\Psi }`$ quanta.

Let us turn now to the evolution of the classical sneutrino field $`\stackrel{~}{\mathrm{\Psi }}`$ (5). As far as we did not introduce any mass term for the RH neutrino, in the superpotential (4), the fields $`\psi _{1,2}`$ are intrinsically massless. Nevertheless, before the decay of $`\varphi `$ oscillations, they have field dependent effective mass terms $`m_{1,2}^2=\kappa ^2\varphi ^2\pm \kappa \lambda \sigma \varphi `$ and therefore during the matter dominated phase prior to inflationary phase they simply roll down from the initial values order $`M_\mathrm{P}`$ . However, when $`\varphi `$ settles to zero and inflation begins, these fields become massless and their oscillations are damped. By this moment these fields still have reasonably big values. However, they have further non-trivial evolution during and after inflation, which we shall discuss in subsequent sections.

### B Inflation

As was told in the above, the inflation proceeds when $`\varphi `$ is set to zero while $`\sigma `$ still has a large value $`M_\mathrm{P}`$. On the other side, $`\sigma `$ field has a small effective mass term $`\overline{m}^2\sigma ^2`$ which allows it to roll down the potential. The main contribution to the energy density comes from the false vacuum, $`V(\sigma ,0)\mathrm{\Lambda }^4`$, and the Universe expands exponentially up to the time $`t`$ when the inflaton field reaches the critical value $`\sigma _c`$, where on the Universe exits gracefully from almost de-Sitter expansion during which the cosmological scale has grown by a factor $`\mathrm{exp}(Ht)`$. The Hubble parameter during the inflation is given by

$$H\sqrt{\frac{8\pi }{3}}\frac{\mathrm{\Lambda }^2}{M_\mathrm{P}}\frac{\eta ^2}{\lambda }\times 410^{11}\mathrm{GeV}.$$
(14)
As soon as $`\varphi `$ settles to zero and inflation begins, $`\psi _{1,2}`$ become massless and the curvature of their potential is induced only by quartic self-couplings in (7). These fields are completely decoupled from the inflaton and thus they evolve almost independently during the inflationary era, with the following equations of motion:

$$\ddot{\psi }_{1,2}+3H\dot{\psi }_{1,2}=\kappa ^2(\psi _1^2+\psi _2^2)\psi _{1,2},$$
(15)
or, in terms of the modulus $`\rho `$, $`\ddot{\rho }+3H\dot{\rho }=\kappa ^2\rho ^3`$. As for the phase $`\delta `$, it essentially becomes a flat mode as far as $`\varphi =0`$ and the lepton number conservation is restored.

The classical field $`\rho `$ keeps rolling down until it approaches values $`H`$. After this its dynamics is almost frozen and the rest of evolution is determined by the equation $`3H\dot{\rho }\kappa ^2\rho ^3`$. From here one immediately obtains that by the time $`t`$ when the slow roll ends up, the field value will be $`\rho \kappa ^1(3H/t)^{1/2}`$. This can be rewritten as

$$\rho =\frac{\sqrt{C}}{\kappa }H,$$
(16)
where $`C=3/N_\mathrm{e}`$, with $`N_\mathrm{e}=Ht`$ being the total number of e-foldings. Notice also, that in hybrid inflation models, unlike chaotic models, the number of e-foldings can not be arbitrarily large and for our purpose we assume that maximum number of e-folding can at most be $`N_\mathrm{e}100`$.

So far in the evolution of $`\mathrm{\Psi }`$ fields we have not taken into consideration any kind of supergravity correction to their masses. In particular, for a generic Kähler potential the fields $`\stackrel{~}{\mathrm{\Psi }}`$ could get the mass term $`CH^2(\psi _1^2+\psi _2^2)`$, $`C`$ being order 1 coefficient. The sign of $`C`$ is model dependent and it can not be determined correctly. If the correction is positive, then during inflation $`\rho `$ evolves as $`e^{3Ht}`$ and it will essentially vanish at the inflation exit. However, if the mass correction turns out to be negative, then $`\stackrel{~}{\mathrm{\Psi }}`$ fields will have a false minimum with a non-zero VEV which breaks the lepton number: $`\rho ^2=CH^2/\kappa ^2`$. Therefore, during inflation these fields will fastly evolve down until being trapped in the false minimum, and remain stuck to it until the effective mass of $`\mathrm{\Psi }`$ field overtakes the expansion parameter $`H`$. This happens very soon after the end of inflation, because $`M_\mathrm{\Psi }^2\mathrm{\Lambda }^2H^2`$. Therefore, by the end of the inflation era $`\rho `$ has a non-vanishing magnitude which can be given still by (16), but this time $`C`$ being some unknown order 1 coefficient.

### C Post-inflation

The evolution of $`\sigma `$ and $`\varphi `$ after the inflation exit has been studied in Refs. . After the phase transition the fields oscillate with more or less similar amplitude near their global minima (8), the maximum amplitude attained by $`\sigma `$ field is $`\sigma _c=\varphi _0/\sqrt{2}`$, while $`\varphi `$ takes at most $`\varphi _0`$. The initial conditions for the oscillations are fixed by the inflationary dynamics and are given by :

$`\sigma _\mathrm{i}`$ $`=`$ $`\sigma _c\pm {\displaystyle \frac{H_\mathrm{i}}{2\pi }},\varphi _\mathrm{i}={\displaystyle \frac{H_\mathrm{i}}{2\pi }},`$ (17)
$`\dot{\sigma }_\mathrm{i}`$ $`=`$ $`{\displaystyle \frac{1}{3H_\mathrm{i}}}{\displaystyle \frac{V}{\sigma }},\dot{\varphi }_\mathrm{i}={\displaystyle \frac{1}{3H_\mathrm{i}}}{\displaystyle \frac{V}{\varphi }},`$ (18)
where the initial velocities are calculated at $`\varphi _\mathrm{i}`$ and $`\sigma _\mathrm{i}`$, and $`H_\mathrm{i}`$ is the Hubble parameter at the end of inflation.

The oscillations are essentially an-harmonic in nature. However, near the global minimum, there exists a particular solution which satisfies a straight line trajectory in $`\varphi \sigma `$ plane ;

$$\varphi =\sqrt{2}(\sigma _c\sigma ).$$
(19)
Near the bottom of the potential the oscillations are harmonic and their frequency is governed by the mass of the fields at their minima, see (9). It is possible to give an approximate analytical solution

$$\frac{\varphi (t)}{\varphi _0}1+\frac{A(t)}{3}\mathrm{cos}(m_\varphi t),$$
(20)
where $`A(t)1/t`$ is slowly time-varying amplitude of the oscillations, which depends crucially upon the ratio $`H_\mathrm{i}/M_\varphi `$. The smaller the ratio is, the larger is the number of oscillations of the fields before they feel the Hubble expansion. For smaller inflationary scales, such as $`\mathrm{\Lambda }10^3M_\mathrm{P}`$, the ratio $`H_\mathrm{i}/M_\varphi `$ is small, $`H_\mathrm{i}/M_\varphi 1.3\mathrm{\Lambda }/\lambda ^{1/2}M_\mathrm{P}1`$. This is an interesting characteristics of the supersymmetric hybrid inflation model, which tells us that there are many oscillations of the background fields with an almost constant amplitude. The effect of expansion is felt after many oscillations, and this has been verified numerically in Refs. .

Let us consider now post-inflationary dynamics of the classical fields $`\stackrel{~}{\mathrm{\Psi }}`$. Although they have been essentially the massless fields during the inflation, after inflation they become massive, with mass $`M_\mathrm{\Psi }\mathrm{\Lambda }`$. In other words, $`\psi _{1,2}`$ fields get quadratic terms in the potential and they start to oscillate near the origin with the frequencies governed by mass $`M_\mathrm{\Psi }`$ (10). Depending on the situation, the initial amplitudes for these oscillations are the field which they have at the exit of inflation, see (16). For keeping more generality, let us parameterize the magnitudes of these amplitudes as

$$\psi _{1(2)}^\mathrm{i}=\rho _\mathrm{i}\mathrm{cos}\delta _\mathrm{i}(\mathrm{sin}\delta _\mathrm{i}),\rho _\mathrm{i}^2=\frac{C}{\kappa ^2}H$$
(21)
where we assume that the phase $`\delta `$ is order one, and $`C=3/N_\mathrm{e}`$ or $`C1`$, depending on the situation whether at inflation stage these fields had the supergravity induced order $`H^2`$ mass terms or not. Therefore, the classical fields $`\psi _{1,2}`$ have the initial amplitudes $`H`$, while $`\sigma `$ and $`\varphi `$ fields have much larger amplitudes $`\mathrm{\Lambda }/\sqrt{\lambda }`$. The initial energy density $`V\mathrm{\Lambda }^4`$ is dominated by oscillations of $`\sigma `$ and $`\varphi `$, while the contribution of $`\mathrm{\Psi }`$ is negligible. In this way, one can safely neglect backreaction of $`\stackrel{~}{\mathrm{\Psi }}`$ fields on the oscillation of the classical background fields $`\sigma `$ and $`\varphi `$.

The post-inflationary evolution of classical $`\stackrel{~}{\mathrm{\Psi }}`$ fields is quite interesting. The oscillation frequency of these fields is governed by the value of $`\varphi `$, and also the evolution of $`\sigma `$ leads to an additional contribution to their equation motion. On the other hand, once they become massive and can decay into leptons and higgsinos, the decay rate (26) also contributes the friction term in their equation of motion, which now read

$`\ddot{\psi }_{1,2}+(3H+\mathrm{\Gamma })\dot{\psi }_{1,2}=`$ $`\kappa ^2(2\varphi ^2+\psi _1^2+\psi _2^2)\psi _{1,2}`$ (22)
$`2\kappa \lambda \sigma \varphi \psi _{1,2},`$ (23)
Notice, that the last terms in the above equations come with an opposite sign and so the evolution of $`\psi _1`$ is different from $`\psi _2`$. In other words, in the background of the $`\sigma ,\varphi `$ fields, $`\psi _{1,2}`$ get the $`U(1)`$ invariant effective mass term $`M_+^2(\varphi )\stackrel{~}{\mathrm{\Psi }}^{}\stackrel{~}{\mathrm{\Psi }}=2\kappa ^2\varphi ^2(\psi _1^2+\psi _2^2)`$, as well as the $`U(1)`$ violating one $`M_{}^2(\sigma \varphi )(\stackrel{~}{\mathrm{\Psi }}^2+\stackrel{~}{\mathrm{\Psi }}^2)=2\lambda \kappa \sigma \varphi (\psi _1^2\psi _2^2)`$.

As we have already remarked, $`\varphi ,\sigma `$ fields make many oscillations in one Hubble time, which allows us to consider the average effect of these fields upon $`\psi _1,\psi _2`$ fields. In other words, one can replace $`\varphi ^2`$ and $`\sigma \varphi `$ by their mean values within a period of one Hubble time, for which from (19) and (20) we obtain:

$$\varphi ^2_t\frac{2\mathrm{\Lambda }^2}{\lambda },\sigma \varphi _t\frac{\sqrt{2}\mathrm{\Lambda }^2}{18\lambda }A^2(t),$$
(24)
where $`A(t)^21/t^2`$. Substituting these averages in the (22), we see that during posinflationary oscillations $`\psi _1`$ and $`\psi _2`$ have different dynamical mass terms:

$$M_{1,2}^2(t)=M_+^2\pm M_{}^2=M_\mathrm{\Psi }^2\left(1\pm \frac{A(t)^2}{25}\right)$$
(25)
This dynamical mass splitting causes their helical motion in the background of $`\sigma `$ and $`\varphi `$ and produces the lepton asymmetry in our model.

The evolution $`\stackrel{~}{\mathrm{\Psi }}`$ fields crucially depends whether the friction term in (22) is dominated by $`H`$ or $`\mathrm{\Gamma }`$. In the former case, the amplitude $`\psi `$ decreases with time as $`1/t`$, while in the latter case as $`\mathrm{exp}(\mathrm{\Gamma }t/2)`$. For the width of the decay $`\stackrel{~}{\mathrm{\Psi }}l\stackrel{~}{\phi }_2(\stackrel{~}{l}\phi _2)`$ we have:

$$\mathrm{\Gamma }_\mathrm{\Psi }=\frac{g^2}{8\pi }M_\mathrm{\Psi }\left(\frac{\kappa \eta }{\lambda }\right)^2\left(\frac{m_\nu }{0.1\mathrm{eV}}\right)\times 10^{15}\mathrm{GeV}.$$
(26)
### D Reheating

The end of inflation also marks an era of entropy creation in the Universe. The energy stored in oscillations of classical fields $`\varphi `$ and $`\sigma `$ eventually decays into relativistic particle species. In our case reheating could occur via two possible channels. Either $`\sigma `$ and/or $`\varphi `$ have some dominant channel to decay into standard particles, or, $`\varphi `$ could decay into right-handed neutrino $`\mathrm{\Psi }`$ quanta, which would subsequently decay into leptons and Higgs to produce a final thermal bath. The choice of the dominant channel shall depend on the couplings in the theory which in our case are constrained to guarantee that the reheat temperature of the Universe is not be too large. In particular, in supersymmetric theories there is an upper bound on the reheating temperature $`T_\mathrm{r}10^9`$ GeV or so . If this condition is not met, then the gravitinos are effectively produced via scattering processes in thermal bath, and their late decay products can genuinely threaten the standard nucleosynthesis unless there is some other mechanism to dilute their number density.

The reheating temperature of the Universe era can be estimated as

$$T_\mathrm{r}0.5g_{}^{1/4}\sqrt{\mathrm{\Gamma }M_\mathrm{P}}0.1\sqrt{\mathrm{\Gamma }M_\mathrm{P}},$$
(27)
where $`\mathrm{\Gamma }`$ is the decay rates of the fields $`\sigma `$ or $`\varphi `$, and $`g_{}`$ is an effective number of the particle degrees of freedom in the thermal bath. (In the supersymmetric standard model which we consider here, $`g_{}200`$).

The following remark is in order. Apart from the perturbative decay of the oscillating fields, the coherent oscillations of the inflaton could also lead to a non-thermal resonant production of particles , and in particular of gravitinos with helicity $`3/2`$ , and helicity $`1/2`$ . This would produce very stringent upper bound on the reheating temperature. However, it has been realised later on that non-thermal production of gravitinos during inflation occurs via the goldstino mode, a helicity $`\pm 1/2`$ component which can be recognised as inflatino. It eventually decays along with the inflaton to reheat the Universe, and thus non-perturbative production of gravitinos should not pose serious problems for nucleosynthesis .

One of the most important criteria for a resonant particle production is the coherent oscillations of the background fields, which is fulfilled in supersymmetric hybrid inflationary model. It has been noticed that the particle creation is quite efficient in supersymmetric hybrid models, see Ref. . In particular, quanta of the heavy RH neutrinos produced at the preheating stage by non-perturbative decay of inflaton oscillations, in their subsequent decays could produce the lepton asymmetry, provided that these decays are CP-violating .

Another interesting feature of hybrid model which we must mention here is the possibility to realise tachyonic preheating. This is because near the critical point mass squared for $`\varphi `$ field flips its sign and becomes a tachyonic mode. This violates the adiabatic vacuum condition and leads to explosive production of particles . In such a case, it is also important to consider the backreaction of the quanta on the classical fields, which actually leads to destabilising the system and the zero-mode trajectories . Eventhough, tachyonic preheating might work, it is still not very clear why there should be a resonant particle production once the backreaction of the quanta is properly taken into account.

Whatsoever be the case, we shall not worry too much upon the preheating aspects. All we assume here is that the decay of $`\varphi `$ and $`\sigma `$ is responsible for reheating the Universe. We also neglect the finite temperature effects, which could emerge if the relativistic particle species produced by the inflaton decay thermalize too early, before the inflaton decay ends up, and thus can give thermal corrections to the potential of the fields such as $`\stackrel{~}{\mathrm{\Psi }}`$ in our case. The finite temperature mass corrections can be avoided once thermalization of the Universe is delayed until the last stages of reheating.

Coming back to our situation, we have to control that the already existing couplings in the model will not produce too large reheat temperature. The dominant channel for the decay of $`\mathrm{\Phi }`$ field is governed by the last term in the superpotential (4). If $`M_\varphi >M_\mathrm{\Psi }`$, then $`\mathrm{\Phi }`$ decays into $`\mathrm{\Psi }`$ quanta which then produce the standard particles via the coupling $`g\mathrm{\Psi }l\phi _2`$. In this case we have $`\mathrm{\Gamma }_\varphi =(\kappa ^2/8\pi )M_\varphi `$ and thus $`T_\mathrm{r}0.1\sqrt{\mathrm{\Gamma }_\varphi M_\mathrm{P}}3\kappa \eta ^{1/2}\times 10^{15}`$ GeV, where we have used (9). Taking now into account the fact that in our approach the constant $`\kappa `$ should be rather large , since otherwise the decay of the orthogonal field $`\varphi `$ would not be effective in the preinflationary phase, and neither $`\eta `$ can be be very small, we obtain too a big reheat temperature.

However, if $`M_\mathrm{\Psi }>M_\varphi `$, i.e. $`\kappa >\lambda `$. situation is more interesting and the problem can be solved without any extra assumptions. Now the decay channel $`\mathrm{\Phi }\mathrm{\Psi }\mathrm{\Psi }`$ is kinematically forbidden, and the lowest order relevant operator for the decay of $`\varphi `$ into lighter particles is the $`D=6`$ one $`(\kappa g^2/M_\mathrm{\Psi }^2)\mathrm{\Phi }l^2\phi _2^2`$. Such term in the superpotential is induced by the Yukawa couplings $`g\mathrm{\Psi }l\phi _2`$ after integrating out the heavy superfields $`\mathrm{\Psi }`$. Now the decay width of the four-body decay $`\mathrm{\Gamma }(\varphi ll\phi \phi )=P\kappa ^2g^4M_\varphi ^5/M_\mathrm{\Psi }^4`$ is suppressed by a phase space factor $`P10^5`$. Finally, by taking into account (9), (10) and (11), the reheating temperature comes naturally as

$$T_\mathrm{r}\lambda \eta ^{3/2}\left(\frac{m_\nu }{0.1\mathrm{eV}}\right)\times 10^{14}\text{GeV}.$$
(28)
For instance, by taking both $`\eta `$ and $`\lambda `$ order $`10^2`$ and $`m_\nu 0.1`$ eV, we get $`T_\mathrm{r}10^9`$ GeV. This result also satisfies the thermal gravitino production bound on the reheating temperature. In addition, as we shall see in the next section, such a reheat temperature can easily produce the right amount of the baryon asymmetry.

## IV Lepton asymmetry

As we have mentioned earlier, the RH sneutrino field $`\stackrel{~}{\mathrm{\Psi }}`$ carries the global $`U(1)`$ charge, and the associated quantum number $`BL`$ (rather than $`L`$) is violated via the last term in the potential (7). This term is effective during the post-inflationary oscillations and it gives rise to a helical motion of $`\psi _{1,2}`$ in the background of $`\sigma `$ and $`\varphi `$ fields, thereby generating the lepton number in the postinflationary universe.

The $`BL`$ charge density generated during the time evolution of the RH sneutrino field is a zero component of the global $`U(1)`$ current:

$$n_{\mathrm{B}\mathrm{L}}=\frac{i}{2}(\dot{\stackrel{~}{\mathrm{\Psi }}}^{}\stackrel{~}{\mathrm{\Psi }}\stackrel{~}{\mathrm{\Psi }}^{}\dot{\stackrel{~}{\mathrm{\Psi }}})=\frac{1}{2}(\psi _1\dot{\psi }_2\dot{\psi }_1\psi _2).$$
(29)
Then from (22) we immediately obtain the equation which describes its evolution:

$$\dot{n}_{\mathrm{B}\mathrm{L}}+3Hn_{\mathrm{B}\mathrm{L}}=2\kappa \lambda \sigma (t)\varphi (t)\psi _1(t)\psi _2(t).$$
(30)
where the RH side acts as a source term which generates a net $`BL`$ asymmetry through a non-trivial motion of $`\psi _1`$ and $`\psi _2`$ fields.

Let us integrate this equation from the time moment $`t_\mathrm{i}`$, which corresponds to the end of inflation, up to a finite time interval $`t`$:

$$n_{\mathrm{B}\mathrm{L}}=\frac{2\kappa \lambda }{R_t^3}_{t_\mathrm{i}}^t𝑑t^{}R_t^{}^3\sigma \varphi _t^{}\rho ^2(t^{})\mathrm{sin}2\delta _t^{}$$
(31)
where $`R_t`$ is scale factor and we have substituted the mean value $`\sigma \varphi _t`$ given in (24). The CP-phase $`\delta _t`$ changes slowly, since for $`\lambda <\kappa `$ oscillations of the fields $`\psi _1`$ and $`\psi _2`$ have about the same oscillation frequency (see (25). In terms of the field quanta, $`n_{BL}`$ is nothing but a number density difference between the RH sneutrino $`\stackrel{~}{\mathrm{\Psi }}`$ and anti-sneutrino $`\stackrel{~}{\mathrm{\Psi }}^{}`$ states. It will be transmitted to the standard particle system via the decay of the RH sneutrinos into the ordinary leptons and Higgsinos (or sleptons and Higgses), Therefore, even if the decay rates $`\stackrel{~}{\mathrm{\Psi }}l\stackrel{~}{\phi }_2`$ and $`\stackrel{~}{\mathrm{\Psi }}^{}\overline{l}\stackrel{~}{\phi }_2`$, are exactly the same (no CP-violation in decays), we produce the different amount of $`l`$ and $`\overline{l}`$.

Once the $`BL`$ is non-zero, the net baryon number is induced via sphaleron effects which violate $`B+L`$ but conserve $`BL`$ . The sphalerons are active in a temperature range from about $`10^{12}`$ GeV down to 100 GeV. In the context of the supersymmetric model the relation between the $`B`$ and $`BL`$ is given by $`B=0.35(BL)`$ . Therefore, for obtaining the observed baryon to entropy density ratio in the range $`B=n_B/s=(0.31)10^{10}`$, as it is restricted by the primordial nucleosynthesis bounds, we need $`BL𝒪(10^{10})`$.

Now we are in grade to calculate the $`BL`$ number to entropy density ratio $`BL=n_{BL}/s`$ in the Universe. Assuming that the entropy is generated at the postinflationary reheating and there is no more entropy injection at later times, we calculate from (31) the value of $`n_{BL}`$ produced by the reheating time $`t=t_\mathrm{r}0.3g_{}^{1/2}M_\mathrm{P}/T_\mathrm{r}^2`$, and compare it to the entropy density at this time, $`s=(2\pi ^2/45)g_{}T_\mathrm{r}^3`$.

The Universe dominated by the field oscillations expands as in a matter dominated era, and so the scale factor changes as $`R_tt^{2/3}`$. On the other hand, $`\sigma \varphi _t1/t^2`$ while the CP-phase $`\delta _t`$ changes slowly, since for $`\lambda <\kappa `$ oscillations of the fields $`\psi _1`$ and $`\psi _2`$ have almost the same oscillation frequency. Therefore, the integrand in (31) goes as $`\rho ^2(t)`$, and the final answer depends whether the decay rate (26) is larger or smaller as compared to the Hubble parameter. Let us consider first the case $`\mathrm{\Gamma }<H`$. As we discussed in previous section, then we have $`\rho ^21/t^2`$. This suggests that the maximum contribution to integral (31) comes at the initial times $`1/H_\mathrm{i}`$ and hence we obtain:

$$BL\frac{\kappa \rho _\mathrm{i}^2}{H^2}\frac{T_\mathrm{r}}{M_\mathrm{P}}\mathrm{sin}2\delta _{\mathrm{eff}},$$
(32)
where $`\delta _{\mathrm{eff}}1`$ is an effective $`CP`$ violating phase which we assume to be $`𝒪(1)`$, and $`\rho _\mathrm{i}`$ is given by (21).

However, if the friction term in (22) is dominated by the decay width, i.e. $`\mathrm{\Gamma }>3H`$, then the RH sneutrino field behaves like $`\rho ^2\mathrm{exp}(\mathrm{\Gamma }t)`$ and, as a result, the magnitude of $`BL`$ will be reduced by a factor $`x=3H/\mathrm{\Gamma }`$ which can be obtained from (14) and (26): Therefore, the final result for the lepton asymmetry reads:

$$BLX\frac{T_\mathrm{r}}{M_\mathrm{P}}X\left(\frac{T_\mathrm{r}}{10^9\mathrm{GeV}}\right)\times 10^{10},$$
(33)
where the numerical factor is $`X=\mathrm{min}\{1,x\}C\kappa ^1`$, with $`C`$ being $`1`$ coefficient if the field $`\stackrel{~}{\mathrm{\Psi }}`$ gets supergravity induced mass term during inflation, or $`C=3/N_\mathrm{e}`$ otherwise.

For $`X110`$, the result (33) implies correct magnitude of the baryon asymmetry if $`T_\mathrm{r}10^9`$ GeV or so. This range for the reheating temperature can be natural in the context of our model, provided that the factor $`\lambda \eta ^{3/2}`$ in (reftr) is $`𝒪(10^5)`$. The larger $`T_\mathrm{r}`$ would contradict to the non-thermal gravitino production limit. On the other hand, if $`x>1`$, then the coefficient $`X`$ order 1 or 10 can easily emerge if $`\kappa 10^1`$ and $`C1`$ or $`1/20`$, the latter value attained to the case of 60 e-fold inflation.

## V Discussion

Up to now we have considered only one fermion generation. Let us incorporate now all three generations and discuss what happens in this case. In other words, we introduce 3 lepton species $`l_a,e_a^c`$ and $`\mathrm{\Psi }_a`$, $`a=1,2,3`$ is a generation index.

Now the last term in the superpotential (4) becomes $`\kappa _a\mathrm{\Phi }\mathrm{\Psi }^2`$. (we have taken this couplings diagonal, i.e. we work directly in a basis of the RH neutrino mass eigenstates). Let us assume for simplicity, that the neutrino Dirac terms are also diagonal in this basis: $`g_al_a\mathrm{\Psi }_a\phi _2`$. In fact, there are the fermion mass models of this type (see e.g. Ref. ), in which the neutrino mass matrix is diagonal and non-zero mixing angles in lepton sector emerge exclusively from the charged lepton mass matrix.

Let us assume that all three constants $`\kappa _{1,2,3}`$ are enough large, namely $`\kappa _a>\lambda `$, in order to evade the excessively large reheating temperature of the Universe. In this case, the hierarchy of neutrino masses goes as $`m_am_{\nu a}g_a^2/\kappa _a`$, and it should emerge from the hierarchy of the constants $`g_{1,2,3}`$.

The evolution pattern of the classic fields can be extended for the case of three RH sneutrinos in a straight forward manner. In particular, we see from (28) that the reheating temperature $`T_\mathrm{r}`$ is essentially determined by the largest neutrino mass. Recalling also that the atmospheric neutrino oscillations point to the neutrino mass in the range $`m_30.1`$ eV, we see that $`T_\mathrm{r}10^9`$ GeV can be obtained in our model provided that $`\eta 10^2`$ and $`\lambda 10^2`$, a quite natural parameter range in the hybrid model. However, the reheating temperature much smaller than this estimate is not very appealing since it would need unnaturally small value of $`\lambda `$.

On the other hand, the amount of produced $`BL`$ crucially depends on the coefficients $`x_a=3H/\mathrm{\Gamma }_a`$. In order to avoid too strong suppression of the result (33), at least one of the factors $`x_a`$ should be larger than 1 or so. Interestingly, this suppression factor is inverse proportional to neutrino mass – indeed, from (14) and (26) we see:

$$x_a=\frac{\lambda }{\kappa _a^2}\left(\frac{10^4\mathrm{eV}}{m_a}\right)$$
(34)
and thus the the largest contribution to $`BL`$ is given by the lightest neutrino mass, presumably $`\nu _1`$. Thus, the condition $`x_1>1`$ implies the upper bound on the lightest neutrino mass $`m_1<(\lambda /\kappa _1^2)\times 10^4\mathrm{eV}`$, which limit e.g. for $`\lambda <10^2`$ and $`\kappa _1>0.1`$, leads to $`m_1<10^4`$ eV. This limit can be naturally met by the by the mass of the first generation neutrino, if the neutrino mass hierarchy is about the same as that of charged leptons .

Let us conclude by summarising some interesting features of our model, which is just the simple supersymmetric hybrid model with the superpotential (4). The couplings $`\kappa _a\mathrm{\Phi }\mathrm{\Psi }_a^2`$ of the auxiliary orthogonal superfield $`\mathrm{\Phi }`$ to the RH neutrino ones $`\mathrm{\Psi }`$ puts the bridge between the inflation and particle physics sectors, thus connecting the inflation scale $`10^{15}`$ GeV to the RH neutrino mass scale needed in the context of seesaw mechanism. As a bonus, these terms can help in solving many problems of the inflationary cosmology and baryogenesis.

First of all, they allow the orthogonal field oscillations to decay enough fastly and thus can prepare the proper initial conditions for the inflation onset starting from almost arbitrary initial field configurations with the classical field values order $`M_\mathrm{P}`$.

And second, at the epoch of postinflationary field oscillations, these terms generate the dynamical lepton number breaking for the RH sneutrino fields after the end of inflation and before the end of reheating era. During de-Sitter era these fields are intrinsically massless modes and evolve very slowly due to quartic self-couplings or because of order $`H`$ mass term induced by the supergravity corrections. In either way, at the end of the inflation these fields have non-zero values order $`H`$. This is a virtue of $`\mathrm{R}`$-symmetry which actually forbids terms like $`S\mathrm{\Psi }^2`$ or $`M\mathrm{\Psi }^2`$ in the superpotential (4). After inflation they start to oscillate near origin and produce the $`BL`$ asymmetry of the Universe. This happens in an elegant way because the associated $`U(1)`$ charge is dynamically broken in the background of the oscillating inflaton fields. After the RH neutrino decay, the produced $`BL`$ number density is transfered to the standard particles, and being reprocessed by sphalerons, gives rise to the net baryon asymmetry of the Universe. In difference from the usual leptogenesis mechanisms with the RH neutrino decay , our mechanism does not require the presence of CP-violation in the lepton mixing. So, it can work in the context of predictive models  which do not contain these CP-phases but are appealing in all other respects.

In our model interesting relations emerge between the inflationary parameters, reheating temperature, $`BL`$ number density and neutrino masses. The amount of the produced $`BL`$ solely depends upon the reheat temperature $`T_\mathrm{r}`$ with some coefficient $`X`$ which incorporates the coupling and can be order 1. In this case, the correct amount of the $`BL`$ is obtained when $`T_\mathrm{r}10^9`$ GeV, close to the upper bound from the thermal gravitino production. On the other hand, the possibility of factor $`X`$ to be order 1, implies the upper limit on the lightest neutrino mass (presumably $`\nu _e`$). In general, our model is compatible also with the neutrino mass spectrum inferred from the atmospheric and solar neutrino oscillations. One could envisage, that in the context of our observation, the hybrid models  designed for explaining the reheating temperature difference between the ordinary and hidden (mirror or shadow) worlds, could also generate the non-zero baryon asymmetry in both sectors.

Concluding, our model does three jobs very neatly. First, it correctly indicates the neutrino mass range by linking the inflation scale to the RH neutrino mass scale in the context of seesaw mechanism. Second, it provides dynamically the proper initial condition for the onset of inflation in hybrid model. And finally, via the classic RH sneutrino fields, it generates the proper baryon asymmetry of the universe at the epoch of post-inflationary oscillations and reheating.

###### Acknowledgements.

We thank R. Allahverdi, M. Bastero-Gil, D. Comelli, A. Dolgov, K. Enqvist and K. Heitmann for valuable discussions. A.M. acknowledges the support of The Early Universe network HPRN-CT-2000-00152. The work of Z.B. is partially supported by the MURST research grant Astroparticle Physics.