Dark matter and low scale leptogenesis in a flavor symmetric neutrino two Higgs doublet model($\nu$2HDM)

We have studied dark matter (DM) phenomenology, neutrinoless double beta decay (NDBD) and realized low scale leptogenesis in an extension of Standard Model(SM) with three neutral fermions, a scalar doublet and a dark sector incorporating a singlet scalar and a Dirac singlet fermion. A generic model based on $A_4\times Z_4$ flavor symmetry including five flavons is used to explain both normal and inverted hierarchy mass patterns of neutrino and also to accommodate the dark matter mass. In this extension of the $\nu$2HDM, the effective neutrino mass observed in $0\nu\beta\beta$ is well within the experimental limit provided by KamLAND-ZEN. In order to validate DM within this model, we have checked relic abundance and free streaming length of the dark sector component, i.e. a Dirac singlet fermion constraining its mass in keV range. More importantly we have also realised low scale leptogenesis simultaneously within this framework and also the Dirac CP phase gets constrained with the results. Bound from LFV is also incorporated in order to constrain the Yukawa couplings. More importantly, we have analyzed the dependence of various phenomenology with decay parameter for different choice of arbitrary complex angles.


I. INTRODUCTION
The significance of physics Beyond the Standard Model(BSM) can be mainly verified by the existence of neutrino mass [1], dark datter(DM)[2, 3] and baryon asymmetry of the Universe(BAU) [4,5]. The mass of neutrinos to be of the order of sub-eV [6] got its confirmation from the neutrino oscillations [7,8]. Neutrino oscillation was first suggested by Pontecorvo in the 1960's. Whereas, its experimental discovery was made by Super-Kamiokande Observatory and Sudbury Neutrino  [9],RENO [10],T2K [11],Double-Chooz [12] have not only confirmed but also measured the neutrino oscillation parameters more accurately. In case of Type I seesaw [13,14], the Standard Model is extended with the help of three right handed neutrinos which results in the generation of neutrino masses at sub eV scale whereas the right handed neutrino(RHN) mass is related to the Grand Unified Theory(GUT) scale. Though the heavy RHN can generate leptogenesis [15] and simulataneously the small neutrino mass, a naturalness problem arises [16]. Some possible pathways to overcome the high scale leptogenesis are through resonant leptogenesis [17], ARS (Akhmedov, Rubakov and Smirnov) proposed a mechanism via neutrino oscillation wherein the baryon asymmetry of the Universe could be generated with the help of CPviolating sterile neutrino oscillations [18,19], or via Higgs decay [20,21]. Such mechanisms however demand for a degenerate RHN mass, which further is another naturalness problem. Whereas, we do have models like the ν2HDM [22] and scotogenic model [23][24][25] which accounts an explanation for these discrepencies by considering hierarchical masses of RHN and thereby producing leptogenesis at TeV scale. In this work, we have taken the ν2HDM and studied various phenomenologies persisting in it.
However, the most recent cosmological constraints on DM comes from the Planck satellite [30].
It ascertains 27% of the present Universe to constitue of DM, which accounts for five times the abundance of the baryonic matter. Depending upon the internal structure i.e. the constituent particles and symmetries, the history of thermalisation of the dark sector take different paths.
This gives rise to different types of DM on the basis of its mass range [31]. We basically have the FIMP, hot and cold type dark matter followed from the above mentioned aspect. Though the lightest of the RHN can be considered to be a sterile DM in keV scale [32,33], it faces stringent bounds from various experiments [34], thereby leaving a very small viable region for it to be a probable dark matter candidate. So, we introduce another scenario in our work by incorporating a dark sector to the ν2HDM. It includes a scalar singlet(η) and a Dirac singlet fermion(ξ), which are charged under Z 2 symmetry [35]. Interestingly, the bounds from X-ray experiment can be ignored as the stability of DM candidate ξ is maintained by the Z 2 symmetry [34]. Due to the fact that we have null results from the Direct detection [36] as well as indirect detection [37], the nature of the DM candidate is taken to be of FIMP type.
The basic idea followed in our work comprises of the realisation of the ν2HDM by the virtue of the flavor symmetries A 4 ⊗ Z 4 . Along with the constituent particles of the generic ν2HDM, we have added a singlet scalar and a Dirac singlet fermion to the model. The newly added particles incorporate the dark sector, by which the DM phenomenology can be explained. As already mentioned, the DM considered in our work is of FIMP type. We have analysed leptogenesis in our work with the consideration of a mass hierarchy of the RHN M N 1 << M N 2,3 . On finding the numerical values of the Yukawa coupling matrix, we further carry out a detailed analysis of leptogenesis as well as neutrinoless double beta decay(0νββ). Also the branching ratio of the decay (µ → eγ) is calculated in order to constrain the Yukawa couplings obtained from the model. We emphasize on a relative study due to the variation in the values of the rotational matrix angles used for parametrization of the Yukawa coupling matrix. This yields variation in decay parameter value K N 1 , which is used in determining baryon asymmetry of the Universe and also explains the various DM phenomenologies. Thus, we see how the choice of the angles has an impact on neutrino phenomenology and related cosmology carried out in this work.
We have further categorised the paper into six sections which are as follows.  In spite of being the most successful theory of particle physics, Standard model(SM) fails to explain quite a few experimental and theoritical phenomena such as smallness of neutrino mass, dark matter(DM), baryon asymmetry of the Universe(BAU), etc. In the SM, neutrino remains massless, although it can accquire a small Majorana mass through dimension five operator [38,39].
where Λ is the effective large mass scale and H is usual Higss doublet. This operator can be commonly realized in the framework of canonical seesaw [14], where SM is extended with three right handed neutrino(N i ) such that- where m D = f v is the Dirac neutrino mass considering vaccum expection value(vev) of Higgs to be < Φ 0 >= v and m N = Λf . Such kind of mechanisms does not have a direct test in the experiment beacuse of the large mass scale of m N .
ν2HDM is one such framework where Eq(2.2) can be realised naturally with mass of the right handed neutrino to be in experimental reach (order of 1 TeV) by extending the SM with three right handed neutrino and a scalar doublet(φ) having a very small vaccum expectation value (vev) [22].
We have to introduce a U (1) L global symmetry under which L Φ = 0, L φ = −1 and L N = 0, so that it forbids type-I seesaw interaction terml LΦ N . Thus, the additional scalar doublet φ couples with right handed neutrino and SM Higgs will couple with quarks and charged lepton. In this model the smallness of Dirac neutrino mass is achieved through this additional scalar doublet.
To accomodate dark matter candidate in this model, we introduce a dark sector to the generic ν2HDM(discussed above).This newly added dark sector includes a singlet scalar η and a Dirac singlet fermion ξ [31,40] which are charged under Z 2 symmetry. Considering m ξ < m η , we can say The scalar doublets Φ and φ in the model can be denoted as: Thus, the corresponding Higgs potential will be: where the lepton symmetry is explicitly and softly broken by the µ 12 term. The condition < η >= 0 should be satisfied for the unbroken Z 2 symmetry so that it donot have any interaction with the Higgs doublet. Also the mixing angle tanβ can be expressed interms of the vev of SM Higgs boson and the inert Higgs doublet, which is given by tanβ = v v . Again the minimization condition is given by the equation: which helps in expressing the vev of Higgs doublet in terms of parameter present in the Higgs potential. Radiative corrections to the term µ 2 12 is proportional to µ 2 12 itself and it is logarithimically sensitive to the cut off scale as it is the only source of lepton number violation [41]. This results in the stabilization of vev hierarchy v << v against radiative corrections. The physical Higgs boson after SSB is given by: where the mixing angle α is given as follows: (2.10) Further after neglecting O(v 2 ) and O(µ 2 12 ), the masses of the physical Higgs bosons are as follows: The effective Lagrangian for the extended ν2HDM can be expressed as: (2.12) From eq(2.12), the light neutrino mass can derived as [13]: where, m d ν = diag(m 1 , m 2 , m 3 ) is the diagonal neutrino mass matrix and U P M N S is the PMNS(Pontecorvo-Maki-Nakagawa-Sakata) matrix.
We also realize the extension of the ν2HDM with the help of discrete flavor symmetry group A 4 ⊗ Z 4 . To maintain the stability of the dark sector, a discrete Z 2 symmetry is also introduced under which the extended particles are charged. Particle content of the model with its charge assignment are given in the table(I). η and ξ fields are charged under Z 2 symmetry, considering them to be the dark sector which donot acquire any vev and thus remain invisible. We introduce five flavon fields ρ, ρ , ρ , ζ, ζ in order to break the flavor symmetry so as to generate the required mass structures.
The effective Lagrangain for Dirac mass is given by: 15) and the right handed fermions are represented by the Lagrangian: The flavon alignments considered in the model are given below: Using this alignment we can have diagonal charge lepton matrix given by: (2. 18) and the Dirac mass matrix will take the form represented as: here y 1 ,y 2 and y 3 stands for the Yukawa couplings. We further express the Dirc mass matrix in terms of the model parameters as: where, a = <φ>ω Λ y 1 , b = <φ>ω Λ y 2 and c = <φ>ω Λ y 3 are the model parameters which actually represents the Dirac masses. Because of the additional Z 4 symmetry the right handed neutrino mass matrix will be a diagonal one given as: where β 1 ,β 2 and β 3 are the Majorana couplings. Incorporating the mass matrices and Yukawa coupling matrix obtained from the model in Eq.2.13, we can finally calculate the active neutrino mass.

Bounds from Lepton flavor violating process:
It is well known that lepton flavor violating processes put significant bound on the model parameter space. The size of the LFV is controlled by the lepton number violating Yukawa couplings Y ij . The MEG collaboration has been able to set the impressive bound on muon decay Br(µ → eγ) < 4.2 × 10 −13 [43]. The branching ratio of l α → l β γ is given by [44]- where G is the Fermi constant and ∆ N i m H + ) 2 and the loop function is expressed as : The Yukawa couplings obtained from the model, i.e. found to be in the range 10 −1 − 10 −6 are consistent with the bounds of lepton flavor violating process. From Fig.(1), we see that the branching ratio value for the process µ → eγ is below the experimental bound given by MEG collaboration for both NH and IH, thereby, making the Yukawa couplings viable in explaining leptogenesis and dark matter.

III. LEPTOGENESIS
We study leptogenesis in this work, which is a consequence of the out-of-equilibrium decay of . The lepton asymmetry produced only by the decay of N 1 is converted into the baryon asymmetry of the Universe(BAU) by the electro-weak sphaleron phase transitions [45].
In case of vanilla leptogenesis, there exists an absolute lower bound on the mass of the lightest RHN to be M N 1 10 9 GeV [46,47]. However, the limit on the lightest RHN mass scale can be lowered to 10 TeV [5,25] in certain scenarios. With the consideration of a hierarchical mass spectrum m φ << M N 1 << M N 2,3 and the incorporation of the Yukawa couplings from the model, the CP-asymmetry term is given by [25]: where, Y represents the Yukawa coupling matrix and the term f (r j1 , η 1 ) is expressed as: The decay rate equation for N 1 which is given by, Meanwhile, the washout effect is checked by the decay parameter where, Γ 1 is the decay width of N 1 , H is the Hubble parameter and z = M N 1 /T with T being the temperature of the thermal bath. We can express H in terms of T and the corresponding equation is given by: In Eq.(3.5), g * stands for the effective number of relativistic degrees of freedom and M P l 1.22 × 10 19 GeV is the Planck mass.
Using the parametrization for Yukawa coupling, i.e. [48] Y = we can verify the following relation: Here, R is an orthogonal matrix satisying the condition R T R=1 which is parametrized as: where, ω 12,13,23 represents the arbitrary complex angles. Thus, we can now express the decay parameter K N 1 interms of the arbitrary angles by the following equation: which can also be verified from Eq.(3.7). In the above equation, (m diag further gives: (3.10) From Eq(3.9) and Eq(3.10), we can see that K N 1 has dependency on the arbitrary complex angles ω 12,13 but is independent of ω 23 . This relationship implies that by proper variation of these arbitrary angles, we can have a comparative analysis of the different phenomenologies associated with it.
In our work, we have choosen some benchmark values for ω 13,23 for which we have varied ω 12 .
We showcase the variations in BAU as well as DM phenomenology for two different values of ω 12 , however, keeping ω 13,23 fixed. Also, the vev of the scalar doublet φ plays a vital role in the determination of K N 1 , which further has significant impact on baryogenesis. As mentioned in various literatures [49,50], the ∆L= 0 washout processes are crucial in context of small values u. Thereby, the condition Γ 1 /M N 1 << 1 for low scale seesaw is satisfied. Thus, the Boltzmann equations for the evolution of the abundance of N 1 and N B−L are given by [46], respectively. n eq N 1 = z 2 2 K 1 (z) is the equilibrium number density of N 1 , where K i (z) is the modified Bessel function of i th type and gives the measure of the total decay rate with respect to the Hubble rate, and W 1 = Γ W Hz is the total washout rate. Again, W 1 = W 1D + W ∆L=2 [5], i.e the total washout term is the sum of the washout due to inverse decays lφ,lφ * → N 1 (W 1D = 1 4 K N 1 z 3 K 1 (z)) and the washout due to the ∆L = 2 scatterings lφ ↔lφ * , ll ↔ φ * φ * which is given by [25], (3.14) Here, g l is the internal degrees of freedom for the SM leptons, andm is the absolute neutrino mass scale, defined by:m The final B-L asymmetry n f B−L is evaluated by numerically calculating Eq.(3.11) and Eq.(3.12) before the sphaleron freeze-out. This is converted into the baryon-to-photon ratio given by [25]: In Eq.(3.16), g * = 110.75 is the effective relativistic degrees of freedom at the time when final lepton asymmetry was produced, g 0 * = 43 11 is the effective degrees of freedom at the recombination epoch and a sph = 8 23 is the sphaleron conversion factor taking two Higgs doublet into consideration. The Planck limit 2018 gives a bound on the observed BAU(n B ) to be (6.04 ± 0.08) × 10 −10 [51].
Therefore, in our work, we have chosen the free parameters such that we can generate the observed BAU for NH/IH. As mentioned earlier, we have fixed the arbitrary angles ω 13,23 bearing values 10 −3 + 10 −3 i and 10 −1 + 10 −1 i respectively. Whereas, we have chosen two benchmark values for ω 12 , i.e 10 −2 + 10 −2 i and 10 −12 + 10 −12 i which shows quite a significant change in the cosmological phenomenologies carried out in this work. In addition to this, we have taken a certain parameter space for v = 0.1 − 30 GeV to show its dependence on BAU for both NH/IH.

IV. DARK MATTER
In this work, we consider a FIMP type dark matter with the real Yukawa coupling, λ << 1 [31] between the right handed neutrino N with the dark sector, i.e. singlet fermion ξ and singlet scalar η. As of the DM type taken into account, it is known that the coupling is so weak that it never reaches thermalization even within the dark sector. We consider the singlet fermion(ξ) to be a probable dark matter candidate. Due to the feeble interaction of the DM candidate, its relic abundance is generated via the freeze-in mechanism [52]. We can obtain the expression for relic abundance by solving the Boltzmann equation: where, Y ξ and Y N 1 are the abundance of DM candidate ξ and right handed neutrino N 1 respectively, and BR ξ is the branching ratio of N 1 −→ ξη. The inverse decay process, i.e. ξη −→ N 1 is neglected and a hierarchical criteria of BR ξ ≡ BR(N 1 −→ ξη) << BR SM ≡ BR(N 1 −→ SM) 1 is considered due to the FIMP nature of ξ. Alongside, for the decay N 1 −→ ξη to obey out of equilibrium condition, it is crucial that BR ξ < 10 −2 . In our work, as we have solved the model parameters so as to find the value of K N 1 , we see that it is consistent in satisfying the bounds for relic abundance and streaming length. From Eq.4.1, an asymptotic abundance of the FIMP DM ξ can be approximated by the relation [31]: Therefore, the expression for relic abundance in agreement with the asymtotic abundance is given by: where, ρ c = 1.05371×10 −5 h 2 GeVcm −3 is the critical density of the Universe and s 0 = 2891.2cm −3 [53] is the current entropy density and h = 0.72 is the Hubble parameter. One of the frequently and significantly considered cosmological constraints that comes into play while considering light DM candidate is the free streaming limit. It provides stringent bounds on the FIMP DM mass. Due to non-trivial velocity dispersion, the free streaming of dark matter particles strike outs the matter density perturbations and consequently the structure formation on scales λ < λ F S . Subsequently, small structure formation gives the most robust constraint on the free streaming length, viz.
λ F S < O(0.1)Mpc [54]. The free streaming length is defined as the average distance covered by a particle without confronting collision [31]; where, < v ξ > is the average velocity at given time of the FIMP DM ξ, a rh and a eq are the scale factors at reheating and equilibrium respectively. The values of the cosmological parameters used are H 0 = 67.3kms −1 Mpc −1 , Ω R = 9.3 × 10 −5 and a eq = 2.9 × 10 −4 [55]. Again, the non-relativistic scale factor for FIMP DM is expressed as [40]: with g * ,0 = 3.91, g * ,rh = 106.75 and T 0 = 2.35 × 10 −4 eV. There are different decoupling and production mechanism of DM and thereby, the free streaming length is different for each of the relativistic decoupling scenarios such as hot, FIMP DM. As the free streaming length is highly dependent on the production mechanism, thus, the specific free streaming length for FIMP DM can be given by [40]: In the flavor symmetric realisation of the extended ν2HDM [40], we have explicitly studied the implication of the Yukawa coupling matrix and the rotational matrix used for parametrization.
On numerically solving the Dirac masses given by a,b and c which constitute the Yukawa coupling matrix, we eventually find the range of the coupling matrix Y. Incorporating the parameter space of the Yukawa coupling matrix in Eq.(3.6), we can deduce the active neutrino mass matrix. It has been found from the model we have worked upon that the summation of the active neutrino is consistent with the Planck limit, i.e. i=1,2,3 m i < 0.11eV [6]. Consecutively, we have also analysed the phenomenology of 0νββ [56][57][58]  eV in case of both NH/IH which is a crucial parameter required to achieve baryogenesis in TeV scale [5]. In our study, we have shown the variations that may occur due to the choice of arbitrary angles of the rotational matrix given in Eq.(3.10). Since, the dependency of BAU and dark matter phenomenology is influenced by the decay parameter K N 1 , thus we can generate desired values of K N 1 by fine tuning the arbitrary angles ω 12 and ω 13 [40]. In our case, we have kept ω 13,23 fixed at 10 −3 + 10 −3 i and 10 −2 + 10 −2 i respectively and simultaneously chosen two different values of ω 12 as 10 −1 + 10 −1 i and 10 −12 + 10 −12 i. In Fig.(2), we have shown co-relation plots between BAU and parameters such as M N 1 , K N 1 and tanβ respectively for NH and IH. From the first row of Fig.(2), we see that the observed baryogenesis is generated for M N 1 = 10 4 − 10 6 GeV for NH, however, the choice of the arbitrary angle ω 12 = 10 −1 + 10 −1 i gives large number of points satisfying the Planck limit compared to that of ω 12 = 10 −12 + 10 −12 i. For ω 12 = 10 −12 + 10 −12 i, the region M N 1 = 5 × 10 4 − 10 6 GeV has prominent points satisfying the desired BAU in case of NH. For IH, we donot observe much change in the variation of BAU w.r.t M N 1 obeying the Planck bound.
Also the variation due to the arbitrary angles are almost similar for both the heirarchies. In the second row of Fig.(2), an interesting result is seen which constraints the parameter space of K N 1 satisfying the Planck limit depending upon the choice of arbitrary angle ω 12 . The points satisfying BAU for ω 12 = 10 −12 + 10 −12 i are very scanty, whereas for ω 12 = 10 −1 + 10 −1 i it is seen that we have abundant points for the decay parameter space K N 1 = 10 − 100, thereby producing thermal leptogenesis and also serving the DM candidate to be a FIMP type for both NH/IH [31,40]. Again from third row of Fig.(2), we have obtained a clear distinction between the range of tanβ satisfying the Planck limit for BAU w.r.t the choice of ω 12 . We see that for ω 12 = 10 −1 + 10 −1 i, the range of tanβ from 10-100 is seen to produce the observed BAU, whereas for ω 12 = 10 −12 +10 −12 i the region of tanβ satisfying the BAU limit is constrained from 10-50 with very less points as compared to that for ω 12 = 10 −1 + 10 −1 i incase of NH. The results obtained for IH is almost similar to that for NH. We have obtained a crucial result from Fig.(3) which remains almost the same for both the mass heirarchies. Here, the parameter space of the Yukawa coupling element (Y † Y ) 11 which has a significant contribution in the CP asymmetry is constrained for ω 12 = 10 −12 + 10 −12 i in a small region obeying the BAU limit. Also we have obtained very few points in the parameter space satisfying the Planck bound incase of lower value of ω 12 . However, we have quite a wide range of (Y † Y ) 11 corresponding to ω 12 = 10 −1 + 10 −1 i satisfying the Planck bound for BAU, i.e.
(Y † Y ) 11 = 10 −8 − 10 −5 . Fig.(4) is a co-relation plot between the relic abundance and the BAU.   ratio in the range Br ξ = 10 −2 − 10 −7 . From Fig.(7), we see that for ω 12 = 10 −1 + 10 −1 i, we get observed relic abundance for DM mass in the range 5-100 GeV which corresponds to warm DM, whereas for ω 12 = 10 −12 + 10 −12 i, smaller mass of ξ signifying hot DM also satisfies the relic abundance limit. In a plot between relic abundance and decay parameter(K N 1 ) as shown in Fig.(8), we observe that the parameter space for K N 1 = 10 − 100 satisfies the Planck limit for relic abundance for ω 12 = 10 −1 + 10 −1 i whereas for ω 12 = 10 −1 + 10 −1 i, the same region of K N 1 hardly satisfies relic abundance constraint. However, for satifying the small structure formation constraint, K N 1 must not fall in the very weak washout region, therefore, ω 12 = 10 −12 + 10 −12 i is less favorable for both NH/IH. Also as mentioned earlier the preferable branching ratio range for warm DM satisfying the relic abundance bound is obeyed for ω 12 = 10 −1 + 10 −1 i as depicted in Fig.(9). We know that the most stringent bound on r F S comes from the small structure formation, r F S < 0.1. A relationship between r F S and K N 1 is shown in Fig.(10), where we have all the points for both the cases of ω 12 in the allowed region. As we know that K N 1 is inversely proportional to r F S , thus, very small values of K N 1 donot satisfy the small structure formation bound. An analysis on the points satisfying both relic abundance and small structure formation is been shown in Fig.(11), wherein we see that almost all the points corresponding r relic abundance bound. Also from Fig.(7), we have obtained the allowed range of warm DM mass satisying the relic abundance to be 5-100 KeV, thus we can summarise that this range of DM mass also abides the streaming length constraint. We have further shown some contour plots between the Dirac masses a, b and c w.r.t the allowed range for relic abundance of the dark matter, i.e.
Summarising the allowed parameter space of the Dirac masses from the above mentioned figures, we finally represent it in tabular form in Table.II and Table.III.

VI. CONCLUSION
Our work basically showcases a comparative analysis of different phenomenological consequences corresponding to variation in the arbitrary angles of the rotational matrix. As mentioned in the earlier sections, we have realised the extended ν2HDM with the help of flavor symmetries A 4 ⊗ Z 4 .
It is an interesting model as it can accomodate both neutrino as well as DM phenomenologies. so as to achieve TeV scale leptogenesis. We have also co-related the BAU with the values of tanβ and analyzed its variation for different ω 12 values. From Fig.(2),(3) and (4) we can conclude that the choice ω 12 = 10 −1 + 10 −1 i is more preferable in satisfying the BAU limit given by Planck data.
Also, we have obtained a constraint range of the free parameters as well as the Yukawa couplings corresponding to the two choices of the arbitrary angles. In plot (4) we obtain points which abide by the constraint from both BAU and relic abundance simultaneously which is a vital part of the analysis. We have also generated some important results in context with the decay parameter.
Depending on ω 12 values, we can generalise the transition in the relation between K N 1 and m l for both NH and IH. Considering the DM scenario, it is viable to have a warm DM candidate as it obeys the constraint coming from small structure formation and relic abundance w.r.t to the allowed parameter space of decay parameter in comparison to hot DM. Simultaneously, we can get a warm DM source when the decay parameter donot fall in the very weak washout region.
Therefore, it is seen from Fig.(5) that for ω 12 = 10 −1 + 10 −1 i we obtain the preferable range of K N 1 which further explains the warm DM. For warm DM, m ξ ∼ 10 KeV gives significant observed relic abundance, however, small DM mass correspond to hot DM which again is constrained by small structure formation. Again, a difference in the branching ratio range is obtained depending on the two different values of ω 12 . So, we can say that Br ξ 10 −5 − 10 −3 is the prefered range for ω 12 = 10 −1 + 10 −1 i obeying the Planck limit for relic abundance, whereas for ω 12 = 10 −12 + 10 −12 i, Br ξ 10 −6 − 10 −3 produces the observed relic abundance. Summarizing the above results, it can be said that for the choice of ω 12 = 10 −1 + 10 −1 i, the value of K N 1 falls in the weak washout region which further is successful in generating the desired BAU, relic abundance and also the small structure formation for both NH and IH. Also, the Yukawa couplings obtained from the model are successful in producing the branching ratio Br(µ → eγ)< 4.2 × 10 −13 as shown in Fig.(1). Thus, the model can explain the neutrino mass, leptogenesis and dark matter for the choice of the free parameters considered and the values of Dirac masses obtained.