Common origin of modified chaotic inflation, non thermal dark matter and Dirac neutrino mass

We propose a minimal extension of the standard model of particle physics to accommodate cosmic inflation, dark matter and light neutrino masses. While the inflationary phase is obtained from a modified chaotic inflation scenario, consistent with latest cosmology data, the dark matter particle is a fermion singlet which remains out of equilibrium in the early universe. The scalar field which revives the chaotic inflation scenario by suitable modification also assists in generating tiny couplings of dark matter with its mother particle, naturally realizing the non-thermal or freeze-in type dark matter scenario. Interestingly, the same assisting scalar field also helps in realizing tiny Yukawa couplings required to generate sub-eV Dirac neutrino mass from neutrino couplings to the standard model like Higgs field. The minimality as well as providing a unified solution to all three problems keep the model predictive at experiments spanning out to all frontiers.


I. INTRODUCTION
Several experimental evidences from cosmological and astrophysical experiments strongly favour the presence of large amount of non-baryonic form of matter dubbed as Dark Matter (DM), in the present universe [1]. However, mysteries surrounding the nature and properties of DM are yet to be resolved by the experiments operating at cosmic, intensity as well as energy frontiers. It is well known that in order to accommodate DM candidate in particle physics, extension of the standard model (SM) of particle physics is required. Depending on the strength of the interactions of DM with the SM particles, different types of DM models have been proposed. Among them, the weakly interacting massive particle (WIMP) paradigm has been the most widely studied dark matter scenario for last few decades. In this framework, a DM candidate typically with electroweak scale mass and interaction rate similar to electroweak interactions can give rise to the correct DM relic abundance [2], a remarkable coincidence often referred to as the WIMP Miracle. Now, if such type of particles whose interactions are of the order of electroweak interactions really exist then we should expect their signatures in various DM direct detection experiments where the recoil energies of detector nuclei scattered by DM particles are being measured. However, the null results at direct detection experiments [3][4][5][6] have brought the several implementations of the WIMP scenario into tension. This has also resulted in an increased amount of interest in scenarios beyond the thermal WIMP paradigm where the interaction strength of DM particle with SM particles can be much lower than electroweak interaction i.e. DM may be more feebly interacting than in the thermal WIMP paradigm. One of the viable alternatives of WIMP paradigm, which may be a possible reason behind null results at various direct detection experiments, is to consider the non-thermal origin of DM [7]. For a recent review of such feebly interacting (or freeze-in) massive particle (FIMP) DM, please see [8]. In the FIMP scenario, DM candidate does not thermalise with the SM particles in the early universe due to it's feeble interaction strength and the initial abundance of DM is assumed to be zero.
At some later stage, DM can be produced non thermally from decay or annihilation of other particles thermally present in the universe.
Apart from the above mentioned DM problem which brings cosmology and particle physics close to each other in the pursuit of finding a solution, another issue faced by the standard hot big bang cosmology is related to the observed isotropy of the cosmic mi-crowave background (CMB), also known as the horizon problem. This can not be explained in the standard cosmology where the universe remains radiation dominated throughout the early stages. In order to solve this, the presence of an rapid accelerated expansion phase in the early universe, called inflation [9,10] was proposed. Originally proposed to solve the horizon, flatness and unwanted relic problem in cosmology [9,11], the inflationary paradigm is also supported by the adiabatic and scale invariant perturbations observed in the CMB [12,13]. It turns out that, one can accommodate such an early inflationary phase within different particle physics models where one or more scalar fields play the role of inflaton.
Chaotic inflation [14,15] models were one of the earliest and simplest scenarios that used power law potentials like m 2 φ 2 with a scalar field φ. It predicts specific values for inflationary parameters like the spectral index n s ∼ 0.967, tensor-to-scalar ratio r ∼ 0.133 for number of e-folds N e = 60. However Planck 2018 results [13] strongly disfavour this simple model due to its large prediction of r. In order to continue using such simple power law type inflationary scenario, one then has to modify the simplest chaotic inflation potential. There exists several proposals in literature where such attempts have been made. For example, radiative corrections to the inflationary potential [16][17][18], non-minimal coupling of the inflaton with gravity [19][20][21] or logarithimic mass correction of the inflaton [22,23] could bring down the value of r within the Planck limit as well as keeping n s within the allowed range.
Another interesting proposal exists, where the inflation sector is extended by an additional scalar field [24,25]. The mere interaction between the inflaton and the additional scalar field can revive the model successfully by reducing the magnitude of r below the Planck 2018 limit [13]. It is to be noted that the presence of linear or cubic terms of the inflaton in the Lagrangian can destroy the required flatness for having successful inflation in the simple chaotic inflationary scenarios. This can be ensured by imposing one discrete symmetry Z 2 under which φ transforms non trivially. Now, after the end of inflation, reheating of the universe is extremely essential so that the production of radiation and other matter fields can occur. This also sets the initial condition for the standard big-bang cosmology. The energy transfer of inflaton to relativistic matter fields can be realised through non perturbative preheating process or perturbative decay of inflaton. In particular, it was shown in [26] that thermalisation of the universe through instant preheating could be one simple and elegant scenario for Z 2 -odd chaotic inflation model. For a recent work on Z 2 -odd inflaton field where reheating occurs from annihilation instead of decay, please see [27,28].
Although the fundamental origins of inflation and DM could be disconnected to each other, it is very well motivating to study them in a common framework. In fact there have been proposals where a single field can play the role of inflaton as well as DM [27,28,. Motivated by such common frameworks, here we study a scenario where inflation and non-thermal DM can find a common origin. However, unlike a similar proposal [50] where non-thermal DM field acted as inflaton, here we consider a fermionic non-thermal DM which is produced dominantly from the decay of a scalar field that was present in the thermal bath of early universe. The decay is assisted by another scalar field which plays a non-trivial role in modifying the minimal chaotic inflation scenario as required by the latest cosmology data mentioned earlier. Thus, although the inflaton and DM fields are not the same, yet they are non-trivially connected by another scalar field which assists in successful inflation and DM production. The presence of DM sector also opens up the possibility of perturbative inflaton decay at the end of inflation.
In order to make our framework more minimal and predictive, we also attempt to address the origin of light neutrino masses. Non-zero neutrino mass and large leptonic mixing are well established facts by now [1] while their origin remains unknown as the SM can not explain them. Apart from neutrino oscillation experiments, cosmology experiments like Planck also constrain neutrino sector by putting an upper bound on the sum of absolute neutrino masses |m i | < 0.12 eV [2]. While the nature of light neutrinos: Dirac or Majorana, remains undetermined at oscillation experiments, popular seesaw models like [58][59][60][61], proposed to account for neutrino masses predict Majorana neutrinos. However, experiments looking for neutrinoless double beta decay (0νββ), a promising signature of Majorana neutrinos, have not yet found any positive results. Though this does not necessarily rule out the Majorana nature, yet it is motivating to study the possibility of light Dirac neutrinos. This has led to several proposals that attempt to generate tiny Dirac neutrino masses in a variety of ways . In a recent work [90], a common origin of light Dirac neutrinos and non-thermal DM was proposed where tiny couplings involved in non-thermal DM and Dirac neutrino mass had common source from higher dimensional operators. Here we extend that idea to incorporate inflation as well 1 , within a modified chaotic inflation scenario. We extend the SM by three additional scalar fields and a fermion DM field with suitable discrete symmetries in order to keep the unwanted terms away. Since light neutrinos are of Dirac-type, three right handed 1 See [29][30][31][32][33][34]40] for earlier attempts in linking neutrino and dark matter with inflation. neutrinos are present by default. While one of the scalars is the inflaton field and one is the mother particle for DM, the third one assists in inflation as well as DM production. We find that correct DM and neutrino phenomenology can be successfully reproduced in the model while the inflationary parameters predicted by the model remain allowed from Planck 2018 data.
This paper is organised as follows. In section II, we present our model and the corresponding particle spectra, including light Dirac neutrinos. In section III, we discuss the details of inflation in our model followed by the details of non-thermal fermion DM in section IV. We finally conclude in section V.

II. THE MODEL
In this section, we discuss our model, its particle content, additional symmetries and the Lagrangian of the new fields. As mentioned in the introduction, we extend the SM with two gauge singlet real scalars (φ and χ), one gauge singlet complex scalar field (η), one vector like fermion (ψ) and three right neutrinos (ν R ). All the new fields considered in the model are singlets under the SM gauge symmetry. We impose two discrete symmetries Z 4 × Z 4 in addition to the SM gauge symmetry, in order to achieve the desired terms in the Lagrangian  2 To keep our discussion minimal, we have adopted such discrete symmetries. UV completion can be achieved by suitable gauge symmetries, for example Abelian gauge extensions [93]. A recent work where B − L gauge symmetry leading to light Dirac neutrinos and a residual Z 2 × Z 2 symmetry can be found in [94].
The scalar Lagrangian as followed from the charge assignments is provided by where we identify the SM Higgs doublet as H and v, v χ are the vacuum expectation values (VEV) of the neutral component of H and χ respectively.
The relevant part of the fermionic Lagrangian consistent with the charge assignments of the fields is given by Here l L denotes the usual lepton doublet of SM and M P is the Planck mass. The bare mass term for vector like fermion ψ is the only renormalisable term involving ψ while all other terms arise only at dimension five level or higher. The discrete symmetries prevent coupling between SM lepton doublet and right handed neutrinos at tree level which would require Dirac neutrino Yukawa to be fine tuned at the level of O(10 −12 ) or even smaller.
The discrete symmetries also help in preventing the Majorana mass term of right handed neutrinos, which is needed to ensure the pure Dirac nature of light neutrinos. After χ, H acquire non-zero VEV's, the light neutrino masses arise as For v ≈ 10 2 GeV, v χ ≈ 10 9 GeV, it is possible to generate M ν ≈ 0.1 eV for Yukawa couplings y ≈ 10 −2 . The light neutrino mixing will arise from the structure of y in the flavour basis, which is not restricted by the symmetries of the model and hence it is possible to fit it with observed mixing [1] without having any consequence for DM and inflation sector as we discuss below. Some of the higher dimensional terms involving ψ can also generate a new contribution to ψ mass after one or more of the scalar fields acquire non-zero VEV.
However, since ψ already has a bare mass term it is always possible to adjust its mass at a suitable value by adjustment of different relative contributions.

III. INFLATION
Here we discuss the dynamics of inflation and its predictions in detail. The inflation is governed by the following potential, where we identify φ as the inflaton and χ is the assisting field. We also consider the coupling coefficients c 1 and λ χ to be real and positive. The global minimum of χ field is denoted by v χ which we assume to be much smaller than M P . We also ignore higher order terms of φ e.g. φ 4 by considering the associated coupling coefficients negligibly small. At the beginning of inflation the χ field acquires a negative mass squared where H Inf is the Hubble parameter during inflation. Thus the χ field is expected to be driven quickly towards its inflaton field dependent local minimum given by, where we assume v 2 χ c 1 λχ φ 2 . Around χ = χ Inf , the effective mass squared of the χ field is positive and obtained as This turns out to be bigger than the H 2 Inf (with suitable choices of c 1 and λ χ and super-Planckian φ during inflation) and hence the χ field is expected to be stabilised at χ Inf with negligible fluctuations [95][96][97][98]. Thus inflation occurs along φ field direction with the χ field stabilised at χ Inf . One can obtain the effective inflationary potential by integrating out the heavier field χ [95][96][97][98] (by replacing equation (6) into equation (5)) which is given by, where we write α = c 2 1 8λχm 2 . From here onwards, for making the analysis simple we shall work with M P = 1 unit. The parameter α determines the amount of deformation of the modified chaotic potential from the minimal chaotic inflation scenario. In left panel of figure   1, we plot the effective inflationary potential V eff Inf (normalised by m 2 ) for α = 0 (blue) and α = 0.0007 (red). As we can see, non zero α, associated with the presence of χ field, makes the inflationary potential flatter than the minimal one (α = 0). The flattening starts to occur near φ ∼ 10 in M P = 1 unit for α = 0.0007 while before that the potential merges with the original chaotic inflation potential. Although flattened, the potential becomes unbounded from below at large value of φ (e.g. after the maximum at φ ∼ 30 for α = 0.0007) as seen from left panel of figure 1. Then in order to have successful inflation we must make an important assumption that inflaton always stays below the maximum of the potential towards the flat part. However the possibility of tunnelling of the inflaton from its minimum to the unstable part of the potential still remains. In that case it is essential to confirm the metastability of the minimum of φ by calculating the corresponding tunnelling probability. We have shown a rough estimate of the decay probability in Appendix I and found it negligibly small.
Before we find the predictions of the model, let us summarize the important conditions or assumptions which we need to ensure in order to realise successful inflation. They are: (i) at the onset of inflation and afterwards the energy of the universe is dominated by φ field which implies λ χ χ 4 < 1 2 m 2 φ 2 , (ii) the χ field during inflation is massive compared to the Hubble scale i.e. m 2 χ > H 2 Inf so that we can integrate out the χ field during inflation, (iii) we keep value of χ field sub-Planckian.
Let us proceed to find out the values of the inflationary observables: spectral index (n s ) and tensor to scalar ratio (r) using the modified potential V eff Inf in equation (8). The analytic expressions for n s and r in our set up are obtained as (in M P = 1 unit) where V = ∂V ∂φ . The number of e-foldings can be determined using where φ end and φ * represent the field value at the end of inflation and the point of horizon exit respectively. The number of e-fold is connected to the inflationary parameters through the following relation [99][100][101][102], where V Inf indicates the energy scale of inflation. The energy density at the end of inflation and reheating temperature of the universe are denoted by ρ end and T rh respectively. Once these are known, N e can be easily computed. Now the spectral index and tensor to scalar ratio in slow-roll inflation model are defined as The curvature perturbation spectrum is given by The observed value of P S is found to be 2.2×10 −9 at a pivot scale k * ∼ 0.05 Mpc −1 [2]. Using be consistent with Planck TT, TE, EE+low E + lensing results [13]. However the model is ruled out if we consider Planck TT, TE, EE + low E + lensing +BICEP 2/Keck Array (BK15) + BAO data [13] which is much stringent than earlier. In addition, unboundness exists for the effective inflationary potential (equation (8)) at large φ value which gives rise to metastability issue [103]. Below we will see that addition of a mere higher dimensional term to the inflationary sector will alleviate this problem and in addition make the inflationary predictions consistent with Planck+BK15+BAO bounds as well.

A. Inflation with higher dimensional operator
In an effort to make the model consistent with Planck+BK15+BAO data we incorporate a higher dimensional operator (which is perfectly allowed from the charge assignments in table I) in the inflationary potential (equation (5)) given by, where we assume c 2 to be real and positive. We take Λ = M P as the natural cut off scale of the theory. The total inflationary scalar potential is written as V Inf T = V (φ, χ) + V Inf HO . Similar to the earlier case, here also we integrate out the heavier χ field which again receives a negative mass-squared larger than the O(H 2 Inf ) at the onset of inflation. Then we obtain the effective inflationary potential which is (in M P = 1 unit), where β 1 = λχ (in M P = 1 unit). Using this effective potential one can calculate the inflationary predictions n s and r using equations ((13)-(14)) by varying β 1 and β 2 for a fixed value of λ χ ∼ 10 −8 . The parameter m would be determined from the observed value of curvature perturbation spectrum (P S ). It is important to note that the inflationary potential in equation (17) has no maxima at large values of φ provided β 1 < 16 × β 2 with β 1 and β 2 being real and positive. This inequality is obtained using the condition V eff Inf > 0 for any arbitrary value of φ. Considering the inequality β 1 < 16 × β 2 is satisfied, the effective inflationary potential (equation (17)) due to the presence of higher dimensional operator V Inf HO in equation (16) is monotonically increasing function of the φ. Therefore, the inflaton φ can naturally roll towards the minimum φ = 0 from any arbitrary large value. This notable feature of the effective inflationary potential makes the set up more favored unlike the previous case (β 2 = 0) where the effective inflationary potential suffers from the issue of unboundeness at large φ.
We have shown a sketch of the effective inflationary potential (equation 17) in figure 2 for two benchmark points which exhibits that the effective inflationary potential (equation (17)) is much flatter than the one in minimal chaotic model with only m 2 φ 2 potential. Here the flattening starts to take place near φ ∼ 8M P . We estimate N e to be 60 corresponding to the inflationary energy scale ∼ 10 16 GeV and reheat temperature ∼ 10 14 GeV similar to the earlier case. We also take into account the limit β 1 < 16 × β 2 while scanning the parameter space. In figure 3 we show the predictions of the proposed model in presence of the higher dimensional term (equation (16) For comparison purpose, we also include the potential structure of the minimal chaotic potential (blue). high scale validity of the model also constrains the magnitude of λ χH which should be of same order as λ χ so that it does not alter the RG running of λ χ to a great extent. Another important point to note in the proposed set up is the order of magnitudes of λ φH and λ φη , the coupling coefficients of inflaton with SM Higgs and η respectively in equation (1). Following [104], it can be shown that λ φH , λ φη > O(10 −6 ) can spoil the required flatness of the chaotic inflationary potential by inducing large radiative corrections. Therefore we keep both λ φH and λ φη smaller than O(10 −6 ) in our scenario.

B. Reheating:
Once inflation ends, φ field rolls down and oscillates about its minimum φ = 0. Then the reheating process starts to take place. In our model we have various possibilities of energy transfer of inflaton to the relativistic degrees of freedom through φ 2 |H| 2 , φ 2 χ 2 , φ 2 |η| 2 interaction terms. If we consider λ φH λ φχ , λ φη , production of H fields at zero crossing of φ field during its oscillation will occur dominantly. Now H field can decay to SM gauge bosons and fermions with coupling strength determined by the SM gauge coupling constant and the Yukawa couplings respectively. Hence it is possible that the produced H field from φ may decay to SM fields before the oscillating field φ returns back to the minimum of the potential from its maximum value during oscillation. This event is known as instant preheating [26,[105][106][107][108]. In the process, the effective mass of H field grows as m 2 H = λ φH φ 2 when the field φ rolls up from the minimum of the effective potential. Then the effectively heavy Higgs field decays to SM fields at the moment when it has the greatest mass, i.e. when φ reaches its maximal value. The number density of produced H fields at zero crossing of inflaton can be obtained as [26,105] n H λ 3/4 φHφ 0 whereφ 0 is the velocity of the inflaton around φ = 0. For chaotic inflation model with quadratic potential, the amplitude of first oscillation is ∼ 0.1M P . Hence the effective mass of produced Higgs field will be around 0.1 λ φH M P which implies Higgs will be non relativistic at that moment. Then the total energy density of the Higgs field can be approximated as [26] ρ H = m H n H ∼ 10 −14 λ Now the created H fields will decay completely at the moment φ reaches its maximum during oscillation if the condition (Γ H ) −1 ∼ ∆t is obeyed, where ∆t is the required time for φ to reach its maximal value during the oscillation. Suppose the decay width of H to SM particles is provided by where δ is the coupling of Higgs with SM particles. In [26], it is shown that for δ 2 λ 1/2 φH ∼ 5 × 10 −4 , indeed the decay products of H can dominate the energy density of the universe.
This can indeed happen if λ φH 10 −8 , which holds in our model. It is to be noted that the energy dilution of the φ field can also happen through the perturbative decay φ → χηη as followed from equation (1) with the corresponding decay rate Considering the instant preheating to be the dominant process to transfer the energy of inflaton to relativistic particles, we can try to find a numerical estimate of the reheat temperature. After several oscillation of the inflaton, the complete transfer of the inflaton energy via φ → H → ff occurs. Then we can use the relation ρ T = π 2 30 g * T 4 where g * is the number of relativistic degrees of freedom in the thermal bath. At this moment we can assign a reheat temperature to the universe as where we take ρ 0 1 2 m 2 φ 2 0 with φ 0 is the initial amplitude of φ oscillation.

IV. DARK MATTER
Here we present the detailed analysis of DM phenomenology in our model. As pointed out earlier, we consider the fermion singlet field ψ as the DM candidate. From the Lagrangian of the model involving ψ upto dimension five level, it is clear that ψ does not have any decay mode and hence stable. This is also ensured by the Z 4 symmetry which remains unbroken even after χ field acquires non-zero VEV 3 . Thus, the contribution to DM is expected to come entirely from ψ as long as it is the lightest particle charged under the unbroken Z 4 symmetry. Since ψ is a gauge singlet and it does not have any renormalisable interactions with other particles of the model, it is natural that its interactions with the visible sector particles will remain out of thermal equilibrium in the early universe, a requirement to realise the freeze-in DM scenario. Before the scalar fields acquire VEV, the production of DM can occur from 2 → 2 scattering processes. However, the effective vertex of such scattering diagrams are Planck scale suppressed. This leads to a small ultra-violet (UV) freeze-in contribution to DM abundance [109]. After the scalar field χ acquires non-zero VEV, there can be an effective two body decay contributions to DM production (infrared freeze in).
This can arise from χ → ψψ and η → ψψ with the effective coupling is governed by the ratio v χ /M P which can be as small as 10 −10 for v χ ≈ 10 9 GeV. Such couplings are quite generic in FIMP scenarios where DM is produced from such two body decays of mother particles [7]. Out of these two contributions, the process η → ψψ will dominant as χ has other decay modes whose couplings are not Planck scale suppressed. For example, χ field can decay dominantly into a pair of Higgs or a pair of η. This tree level decays make the field χ to lose its abundance very quickly. On the other hand, the η field gets produced in the thermal bath in the early universe and eventually can freeze-out from the thermal plasma at a later stage. This is similar to the superWIMP scenario [110] where a metastable WIMP decays into a super-weakly interacting dark matter at late epochs. This can happen because η can decay only into ψ and the corresponding coupling is very small, thereby allowing the freeze-out to occur earlier. Thus, η can decay to ψ while it is in equilibrium and also after its thermal freeze-out. The decay width of mother particle η, where y eff = ξvχ M P is the effective coupling of ηψν R vertex. The relic abundance of both η and ψ can be found solving the following set of Boltzmann equations which can be expressed as where the comoving equilibrium number density of η is given by with T being the temperature of the thermal bath and g is the internal degrees of freedom of η. The equation (24) corresponds to the evolution of the comoving number density (Y = n s ) of mother particle (η) as a function of x = mη T where the first term on the right hand side shows the contribution from the thermal bath to the Y η whereas the second term stands for the dilution due to the decay of η to the DM particles. Similarly, the equation (25) represents the evolution of the DM particles in the universe from the non-thermal contribution coming from the decay of η. The effective relativistic degrees of freedom during thermal equilibrium and entropy degrees of freedom are denoted by the usual notations g * and g s respectively in equations (24) and (25). The σv T η in equation (24) stands for the thermally averaged total annihilation cross section of η field as provided in equation (36) of Appendix II. To estimate the relic abundance of DM, we will numerically solve the coupled differential equations (24) and (25) and use the following expression, where s 0 , Y (x=∞) ψ are the present entropy density and comoving number density respectively and ρ crit is the critical energy density of the universe. Also, here we identify m ψ as m DM .
Before performing the numerical analysis we note down the parameters which serve important roles in determining the total relic abundance of DM: {m DM , m η , y eff , λ Hη }.
In figure 4 we have shown the variation of the comoving number densities for both the fields η and ψ as a function of x (= mη T ) for a specific choice of values of the relevant parameters as mentioned in the figure. The evolution of Y η (red line) clearly shows that η was in thermal bath in the early universe and went out of equilibrium through usual freeze-out mechanism and at some later stage continues to decay to DM particle ψ. The blue line in figure 4 exhibits the evolution of Y ψ which shows that initially the ψ abundance was very small and increases gradually from the decay of the mother particle η. Figure 5 shows the relic density allowed parameter space in m DM -y eff plane and the coloured bar represents the variation of m η . Here we have varied the DM mass from 100 MeV to 100 GeV, m η from 500 GeV to 5 TeV and the y eff from 10 −12 to 10 −8 . It can be viewed that for smaller values of m DM we need larger y eff considering a fixed value of m η . As mass of DM increases, the required number density of DM i.e. Y ψ (proportional to Γ η ) has to be smaller to satisfy the correct relic abundance (see equation (27)), and hence the decrease in y eff is observed from figure   5. Also, it is seen from figure 5 that the increase in m η raises y eff in order to have correct amount of relic density. This can be understood by looking at equation (25), where larger m η causes suppression in the value of Y ψ and hence it requires comparatively large y eff .

V. CONCLUSION
We have proposed a unified framework to account for non-thermal dark matter, modified chaotic inflation and sub-eV Dirac neutrino mass by minimally extending the standard model. The minimal chaotic inflation scenario, being ruled out by Planck data, is modified by an additional scalar field to bring the predictions for inflationary parameters within allowed range. The same additional scalar field also assists in generating required tiny couplings for non-thermal DM as well as Dirac neutrino mass by virtue of Planck scale suppressed dimension five operators. We find that for suitable VEV of the additional assisting field of the order 10 9 GeV, it is natural to generate Dirac neutrino Yukawa of the order 10 −12 −10 −10 which can then generate sub eV Dirac neutrino mass by virtue of neutrino coupling to the SM Higgs. Similar couplings generated for DM coupling to its mother particle also makes the realisation of FIMP dark matter scenario natural. After showing the validity of modified chaotic inflation scenario for suitable benchmark choices of parameters, we numerically find the parameter space that can generate the correct FIMP DM abundance by scanning over DM mass, mother particle mass as well as their couplings. Future cosmology data should be able to make this model go through further scrutiny, specially in terms of the inflationary observables. Also the model can be falsified by observation of neutrinoless double beta decay which will rule out the pure Dirac nature of light neutrinos, as proposed in this model. The decay probability of the false vacuum per unit time and unit volume can be calculated by adopting semi-classical method popularly known as bounce solution. A simpler approximate analytic form of the decay probability is given by [111,112] P where φ I is the starting field value of bounce and S 4 represents the Euclidean action for the bounce configuration. It can be shown that for a λφ 4 type potential S 4 approximately turns out to be [113,114] Now the volume of the past lightcone is estimated to be ∼ e 140 M P

4
. With this, the total probability for nucleation of a bubble in present Hubble volume is [111,112] In our case inflationary potential is dominated by φ 4 term where the unboundeness starts to appear. Therefore it is legitimate to apply this simple method to calculate the tunnelling probability. Considering φ I ∼ 30 M P where the unboundness of the inflaton potential appears, and the quartic coupling coefficient of φ field as λ ∼ −m 2 α −10 −14 , P T comes out to be much smaller than unity.