Minimal renormalizable simplified dark matter model with a pseudoscalar mediator

# Minimal renormalizable simplified dark matter model with a pseudoscalar mediator

Seungwon Baek School of Physics, Korea Institute for Advanced Study, 85 Hoegiro, Seoul 02455, Korea    P. Ko School of Physics, Korea Institute for Advanced Study, 85 Hoegiro, Seoul 02455, Korea    Jinmian Li School of Physics, Korea Institute for Advanced Study, 85 Hoegiro, Seoul 02455, Korea
###### Abstract

We consider a minimal renormalizable and gauge invariant dark matter (DM) model, in which the singlet fermion DM has only axial couplings to a new pseudoscalar mediator. The mixing between the pseudoscalar mediator and the standard model (SM) Higgs boson induces the interactions between the DM and SM particles. The DM candidate in this model can provide the correct thermal relic density and evades all direct detections, while it can produce observable signals in indirect detection experiments due to its large annihilation cross section. A comparative study for DM phenomenology at the LHC is conducted for models with scalar mediators that have either scalar or pseudoscalar couplings to SM particles and the DM. We find that the three scenarios have distinguishable features in scalar decay branching ratio, DM pair production cross section as well as the signal reaches at the LHC. The LHC searches for some visible signals related to the scalar sector are also discussed.

## I Introduction

The existence of non-baryonic Dark Matter (DM) has been established only by astrophysical observations through its gravitational effects Ade et al. (2016). Since the correct abundance of DM via thermal production could be generically obtained if the DM is in the mass range of GeV and interacts with SM particles via electroweak force, the so-called Weakly-Interacting-Massive-Particle (WIMP) paradigm has been one of the most interesting scenarios for thermal DM. Given that the DM interactions with the SM particles or among themselves are unknown, effective field theory (EFT) is one viable way to simplify the study of DM phenomenology. The EFT descriptions Goodman et al. (2011, 2010); Duch et al. (2015) of DM interactions are valid only when momentum transfer is much smaller than the mass of the mediator, which is usually not true for DM productions at high energy colliders Buchmueller et al. (2014); Busoni et al. (2014a, b, c), especially since the mediator mass scale is completely unknown. Simplified DM model frameworks have been used extensively in DM searches at the LHC Abdallah et al. (2014, 2015); Abercrombie et al. (2015). Here, the DM is neutral under the Standard Model (SM) gauge group and interacting with the SM particles via the portal of a single particle Buckley et al. (2015); Harris et al. (2015, 2016); Yang (2016).

However, simplified DM models with a single mediator can often violate the SM gauge symmetry Baek et al. (2016); Kahlhoefer et al. (2016); Englert et al. (2016), thus may become invalid for describing UV-complete models 111Importance of SM gauge symmetry within the DM EFT was pointed out in Ref. Bell et al. (2015).. There are growing interests in simplified DM model that respect the gauge symmetry Bauer et al. (2016); Ko et al. (2016); Baek et al. (2016); Bell et al. (2015); Baek et al. (2013); Baker et al. (2015); Buschmann et al. (2016); Alves et al. (2015). In particular, the gauge invariant and renormalizable DM model with scalar mediators are constructed in its minimal form Baek et al. (2012) and two Higgs doublet model (2HDM) extended form Bell et al. (2016). Models with pseudoscalar mediators are more interesting, owing to the fact that stringent constraints from DM direct detection can be evaded intrinsically, while being able to explain some anomalies in DM indirect detection Ipek et al. (2014); Yaser Ayazi et al. (2016); Kim et al. (2016). The collider phenomenology of UV-complete DM models with pseudoscalar portal has been studied in Ref. Berlin et al. (2015); Fan et al. (2016); Goncalves et al. (2016); No (2016); Ghorbani and Khalkhali (2016).

In this work, a minimal renormalizable model with pseudoscalar mediator is proposed (analogy to the model in Ref. Ghorbani (2015) which focuses on the DM indirect detection signal). Compared to the models in Refs. Lopez-Honorez et al. (2012); Esch et al. (2013), the pseudoscalar mediator of this model only has an axial coupling to DM particles. We show that there is large portion of parameter space that is consistent with DM constraints while giving measurable predictions in future experiments. At the LHC, this model can be searched through signatures both with and without DM in the final state. The most remarkable DM signal is produced by recoiling the DM pair against energetic initial state radiated jet, i.e. mono-jet. We will comparatively study these signatures for models with scalar mediators that have either scalar or pseudoscalar couplings to SM particles and the DM. The pseudoscalar can also produce beyond SM (BSM) signatures without including DM. We will discuss the constraints on the signals of , and at current stage of the LHC.

## Ii Minimal renormalizable model with pseudoscalar mediator

We propose a minimal renormalizable DM model with a pseudoscalar mediator assuming DM is a SM singlet Dirac fermion that couples to a pseudoscalar which is also a SM singlet scalar with a negative parity:

 L =¯χ(i∂⋅γ−mχ−igχaγ5)χ+12∂μa∂μa−12m2aa2 −(μaa+λHaa2)(H†H−v2h2)−μ′a3!a3−λa4!a4 −λH(H†H−v2h2)2 . (1)

Note that the parity is broken by the dim-3 and terms. We remove the tadpole for and assume . This model is unique, since the mediator has a pseudoscalar coupling to the DM , and scalar couplings to the SM fields through its mixing with the SM Higgs boson (see Eq. (7) below), unlike most other renormalizable pseudoscalar mediator models based on 2HDMs and its extensions.

The term induces the mixing between the pseudoscalar and the SM Higgs boson after electroweak symmetry breaking, making two mass eigenstates and :

 H0 =hcosα+asinα , (2) A =−hsinα+acosα . (3)

So the variables , and in Eq. 1 can be expressed by physical parameters in mass eigenstate:

 λH =12v2h(m2H0cos2α+m2Asin2α) , (4) m2a =m2H0sin2α+m2Acos2α , (5) μa =sinαcosαvh(m2H0−m2A) , (6)

where is the vacuum expectation value of .

Then the interaction Lagrangian of and with the SM particles and DM will be given by

 Lint =−igχ(H0sinα+Acosα) ¯χγ5χ−(H0cosα−Asinα) ×⎡⎣∑fmfvh¯ff−2m2WvhW+μW−μ−m2ZvhZμZμ⎤⎦ (7)

The mass eigenstates of scalar fields have only scalar couplings to SM particles and have only axial couplings to DM, so we can expect that such model setup will not lead to any CP-violation effects in the SM.

On the other hand, the extended Higgs sector could affect the electroweak precision test (EWPT) Barger et al. (2008); Baak et al. (2012) by giving extra contributions to the SM gauge boson self-energy. Since the new pseudoscalar boson couples to the SM particles only through mixing with the SM Higgs doublet, constraints from the oblique parameters and the perturbative unitarity bound are exactly the same with the scalar Higgs portal case considered in Ref. Baek et al. (2012). Taking GeV, the measurements exclude the models with scalar mixing angle . Similar constraint is also obtained from the precision measurements of SM Higgs boson signal strengths at the LHC run-I Khachatryan et al. (2015); Aad et al. (2016a), which indicate  Robens and Stefaniak (2015); Cheung et al. (2015); Dupuis (2016). Moreover, if , the stringent limit from the Higgs invisible decay search Br CMS (2016a) requires .

## Iii Dark matter phenomenology

The measurements of anisotropy of the cosmic microwave background (CMB) and of the spatial distribution of galaxies find the relic density for cold non-baryonic matter to be  Ade et al. (2016). In order not to overclose the universe, the DM candidate in our model should annihilate effectively into SM particles. There are mainly three different DM annihilation mechanisms in our model framework: (1) DM mass is around the half of mass so the annihilation cross section is resonantly enhanced; (2) DM annihilate to SM gauge bosons/heavy fermions especially when is large; (3) DM mass is larger than and/or so the annihilation cross section can be enhanced by setting large scalar self-coupling.

The micrOMEGAs BÃ©langer et al. (2015) is used to calculate the observables in DM phenomenology, with the model files for Eq. (1) generated by Feynrules Alloul et al. (2014). Taking as the Higgs state with mass of 125 GeV, the model has seven free parameters:

 mA, gχ, α, mχ, λHa, μ′a, λa . (8)

In DM annihilation, varying the and can only lead to a total rescaling of the cross section, while its dependence on the is more complicate, due to the opening of new annihilation channels with increasing . Further more, as discussed in Sec. II, should be smaller than 0.4 according to the Higgs precision measurement, but not too small to guarantee sufficient signal rate at collider. So we will choose and for the discussions of this section and scan . The determines position of the pole that is due to resonant enhancement in DM annihilation. Scanning will lead to overlapped peaks in annihilation cross section thus smear out the peak structure. For clarification, we also fix GeV. The rest of parameters are scanned in the ranges listed as following.

 λHa∈±[10−3,√4π], μ′a∈[5,300] GeV, λa∈[10−3,√4π] (9)

We will adopt the exponential scan over the and in order to have more points with small , . That is we define | and perform uniform scan over between [-5.5, 1].

The relic density for models in the chosen parameter space are illustrated in Fig. 1. In the region where DM annihilating into Higgs bosons are kinematically forbidden, is the only parameter that control the relic density. The relic density becomes smaller when DM mass is approaching half of the mass. There is also a significant drop at GeV where the DM annihilating into gauge bosons are opening. When , DM can annihilate into scalar bosons through mediation. So the scalar self-couplings are important. Especially, for our parameter choice, is kinematically disfavored, the relic density is monotonically decreased with increasing ||.

The DM has been searched actively by many underground experiments through its recoiling against nuclei  Akerib et al. (2016); Tan et al. (2016). Following the notations of Ref. Kumar and Marfatia (2013), the DM-SM particles interaction can be written in terms of DM bilinear , SM bilinear as well as form factor which includes the details of the model and nuclear form factor:

 M=MχMf⋅F(s,t,u). (10)

In our model, in the limit of low momentum transfer, the DM-SM fermion scattering matrix element is

 M ∝Mχ⋅Mf=−2qi(ξ†χ^Siξχ)× [2mf(ξ†fξf)+iμmfϵijkqivj(ξ†f^Skξf)] , (11)

where is the momentum transfer, are two component spinors for nucleon and DM, is the relative velocity of the dark matter and the target nucleon, is the reduced mass of the dark matter-nucleon system. Note the and dependences are absorbed in . Eq. (11) is showing that the spin-independent (SI) DM-nucleon cross section is suppressed by the while the spin-dependent cross section is even smaller (). The results from the above semi-quantitative estimate can be seen more clearly in the full formula for the SI direct detection cross section,

 σSIχN=2πμ4m2χλ2Nv2, (12)

where

 (13)

with denoting nucleon and . Assuming the relative velocity between the DM and nuclear is given by the orbital speed of the Sun , the typical of our model is around of that in the scalar mediator model Baek et al. (2012) as also have been justified by comparing the scattering rates of and opeartors in Ref. Catena (2015). This means the DM of our model will not leave any signals in direct detection experiments.

However, the s-wave annihilation is still permitted:

 Mχ=¯χ1γ5χ2=−(E1+m1)(E2+m2)+→k2√(E1+m1)(E2+m2)ξ†χ1ξχ2 , (14)

with is the DM momentum. So the non-relativistic DM particles that concentrated at the center of galaxies may still have relatively large annihilation cross section. Thus they can be observed in final state of photons Atwood et al. (2009), positron/anti-proton Picozza et al. (2007); Kounine (2012) or neutrinos Achterberg et al. (2006).

In upper panel of Fig. 2, we plot the cross sections for all DM annihilation channels with varying , where we have chosen appropriate such that the correct relic abundance () is obtained for each point in the scanning. The exclusion bounds from the Fermi-LAT data are available only for some pure final states, e.g , , and . In order to obtain the Fermi-LAT bound to our model, especially when DM is heavy ( GeV) so that it dominantly annihilates to heavy particles (), we assume that for a given DM mass the gamma spectra of the quark and heavy particle final state have similar shape while their normalizations can be different Ackermann et al. (2014); Agrawal et al. (2015). So we can calculate the weighted total annihilation cross section by

 ⟨σv⟩′tot=⟨σv⟩totNγNγ,b¯b (15)

where the is the DM total annihilation cross section, is the number of photons within energy GeV per annihilation for a point in our model and is the corresponding number in simplified model where DM has the same mass as the point and only annihilates to . Similar methodology was also pursued in Ref. Bringmann et al. (2012). We plot the ratio in the lower panel of Fig. 2, from which we can see that the ratio is close to 1 when annihilation is dominant. However, the gauge (Higgs) boson final state can produce less (more) photons in the range GeV than the quark final state. This also leads to a double enhancement of the ratio at GeV, where multiple Higgs final state is kinematically opened. Then, the weighted total annihilation cross section can be compared to the Fermi-LAT bound on the final state directly. We can conclude that the Fermi-LAT data from dwarf galaxies can exclude the light DM mass region ( GeV) as well as the resonant region (), while all of our points are close to the bound and are expected to be discovered/excluded in the near future. It has to be noted that this limit will be weakened if our DM particle only constitutes a fraction of the total amount of DM.

## Iv LHC phenomenology

### iv.1 Invisible channel: mono-jet

In this section, we discuss the DM phenomenology at the LHC in terms of decay of scalar, production of DM and current limits from the LHC searches. To show the merit of our model setup, results are presented alongside with those of conventional theoretical frameworks for DM at collider:

 LAAint =−igχ(asinα+Acosα) ¯χγ5χ −i(acosα−Asinα)∑fmfvh¯fγ5f (16)
 LSSint =−gχ(H1sinα+H2cosα) ¯χχ−(H1cosα−H2sinα) ×⎡⎣∑fmfvh¯ff−2m2WvhW+μW−μ−m2ZvhZμZμ⎤⎦ (17)

In the following, we denote the models of Eq. (17), Eq. (16) and Eq. (7) as SS, AA and SA respectively, since they are distinguished by the scalar/axial couplings between SM particles and DM. For simplicity, in the discussion of this section and are chosen. And the DM mass is fixed to GeV to avoid SM Higgs invisible decay while we keep relatively large DM production cross section. The mass of lighter scalar (pseudoscalar) in SS (AA) scenario is chosen as GeV for comparison purpose. Then, assuming the only decay into SM particles and DM, the only parameter relevant in collider phenomenology is . This minimal decay width for (denoted by hereafter) can be written as

 Γmin(A) =Γ(A→χχ)+Γ(A→VV)+Γ(A→ff) =cos2α⋅g2χmA8π(1−4m2χm2A)i/2 +sin2α⋅Gμm3A16√2πδV ⎷1−4m2Vm2A(1−4m2Vm2A+12m4Vm4A) +sin2α⋅(mfv)23mA8π(1−4m2fm2A)j/2 , (18)

where for SA, SS, AA scenarios respectively, for AA scenario and for .

The branching ratios of are given in Fig. 3. When the is not much larger than , the factor is important. So the Br of SS scenario is smaller than that of SA scenario. As for , both scenarios give the similar branching ratios. The AA scenario always has the largest Br because of the absence of -- coupling.

The dominant DM production channel at the LHC is gluon-gluon fusion (ggF) through the top quark loop. The effective couplings for gluon-gluon-scalar/pseudoscalar after integrating the top quark are

 Lscalar =αs8πgvvτ[1+(1−τ)f(τ)]GμνGμνϕ (19) Lpseudoscalar =αs4πgvvτf(τ)Gμν~GμνA (20)

where , and

 f(τ)=⎧⎪⎨⎪⎩arcsin21√τ,τ≥1−14(log1+√1−τ1−√1−τ−iπ)2,τ<1 . (21)

However, the ggF process itself does not produce any observable signals at detectors. Extra energetic jets radiating from either initial state gluon or top quark in the loop can circumvent this issue, which raise the mono-jet signature. The leading order cross section for DM pair production in association with a jet is computed within the FeynRules/MadGraph5_aMC@NLO Alwall et al. (2014); Hirschi and Mattelaer (2015) framework, where the jet is required to have GeV. Meanwhile, the higher order corrections to the ggF cross section of Higgs production are found to be quite significant. Using the SusHi program Harlander et al. (2013), the NNLO K-factors for Higgs mass GeV are calculated to be around 2.5. So the production cross section for the DM pair associating with a jet is given by the LO cross section in MadGraph5_aMC@NLO multiplying a universal K-factor of 2.5.

The resulting cross sections for all three scenarios are presented in Fig. 4. The contributions of two propagators that mediate the DM production will interference with each other  Ko and Li (2017), leading to different degree of suppressions for different scenarios in the light region. In particular, the cross sections drop dramatically when two propagators are close in mass. Models with heavier are more interesting because of their larger production cross section. In this region, the DM productions are dominated by the on-shell production with subsequent decay. The interference effect becomes important only for GeV, where the on-shell production is kinematically suppressed to some extent. This leads to deviation in the production cross sections of SS and SA scenarios. Note the small bumps around for all scenarios are from the top quark mass effect.

The mono-jet signature has been searched by ATLAS collaboration at 13 TeV with integrated luminosity of 3.2 fb Aaboud et al. (2016). The non detection of the signal could put a constraints on our model parameters. We adopt the CheckMATE2 program Dercks et al. (2016) to calculate the LHC search constraints on our model, in which the ATLAS mono-jet search has been implemented and validated. CheckMATE2 provide the -value at the final stage of its analysis, defined as

 Rmax=maxiNmodeliNupi (22)

where and is the number of signal events of our model and number of new physics upper limit at 95% CL in the signal region , respectively.

In Fig. 5, we present the LHC search limit with signal strength () which gives the size of the cross section that is probable at current stage of the LHC. Even in the region of where the production cross section is largest, the signal rate is at least one order of magnitude below the current reach. Note that in this region, since is mostly on-shell and Br is already close to one, taking larger will not enhance the signal rate. We would expect higher luminosity of LHC to probe/exclude this region. Among three scenarios, the AA scenario has the best search sensitivity. We find that the differences are mainly originated from the production rate of mono-jet signals as shown in Fig. 4, while the kinematic distributions of final states are similar for all scenarios, i.e. similar cut efficiencies.

In a realistic model, some new decay channels of might be important, such as . This will lead to suppressed production rate of DM pair, meanwhile, the interference effect can become remarkable because of the wide width of . In Fig. 5, we also plot the signal reaches for models with ten times larger total width of than due to the opening of new decay channels. In the region with negligible interference, the signal reaches should be one order of magnitude weaker than that of models with , e.g. . The interference effect is significant when off-shell contribution is large, e.g. in the regions GeV. It shrinks the difference in signal reaches for models with narrow and broad width of , mainly because of the enhancement in production cross section. Moreover, the large interference effect can lead to distinguished signal reaches for SS and SA scenarios.

### iv.2 Visible channels

Our model also predicts BSM signals without DM in the final states. In this section, we will focus on the non-DM signals of the SA scenario as we can expect that the exclusion bounds obtained for SA scenario can be directly applied to SS scenario, since their differences only exist in DM sector. But the corresponding bounds in AA scenario could be quite different, due to different production cross section of as well as the absence of tree level couplings.

According to the Eq. (18), the heavy pseudoscalar dominantly decays into top quarks and vector bosons apart from the DM pair. The process of top quark pair production through the pseudoscalar resonance decay interferes strongly with the QCD background, leading to difficulties in its searches at hadron colliders  Dicus et al. (1994); Gori et al. (2016); Carena and Liu (2016). However, the diboson final state may still be detectable. To survey the production cross sections of visible signals in our model, we fix GeV, and varying GeV, , with the rest of parameters scanned in the range as given in Eq. (9). We note that varying and which is important in obtaining correct relic density and evading the DM indirect detections will not affect the results in the following discussions much.

For our parameter choice, the is always important when it is kinematically allowed. So the vector boson pair production cross section is suppressed by , from both production and decay. We calculate the NNLO gluon-gluon fusion production cross section at 13 TeV by using SusHi and obtain decay branching ratios of from micrOMRGAs. The results are shown in Fig. 6. At 13 TeV, the vector boson pair production cross section in our model is only around fb for GeV. The ATLAS collaboration searches the high mass diboson resonance in  ATL (2016a),  ATL (2016b) and  ATL (2016c) final states respectively with LHC run-II data. Their exclusion bounds at 95% confident level (CL) are shown in the Fig. 6 as well. It can be seen that the signal of vector boson pair production is at least two order of magnitude below the current LHC search sensitivities.

On the other hand, the production rates of scalar pairs () do not suffer from the suppression as much as those of vector boson pair, because the coupling in scalar to scalar decay is controlled by the scalar-scalar mixing and scalar self-couplings:

 λAH0H0 =−μacos3α+2(3λH−2λHa)vhcos2αsinα +2λHavhsin3α+(2μa−μ′a)cosαsin2α (23) λH0AA =−μasin3α−2(3λH−2λHa)vhsin2αcosα −2λHavhcos3α+(2μa−μ′a)sinαcos2α (24)

They can be either large or small. In the parameter space of our interest, the and can even become dominant.

When the , the pseudoscalar pair can be produced from the SM Higgs decay, which will lead to four fermion final states after . The cross section of this process can be quite large. Ref. Khachatryan et al. (2017) summarizes the recent searches for light bosons from 125 GeV Higgs decay in the final states of , , and at LHC run-I. The bounds are presented on the production cross section of each final states normalized to the SM Higgs production cross section. In our model, for GeV, the decay branching fractions of the pseudoscalar are only determined by a single parameter . So those experimental bounds for different final states can be projected to the same plane, versus , where . The projected bounds are presented by lines in different colors in Fig. 7. Further more, the precision measurements on Higgs coupling strength constrain the BSM Higgs boson decay to be Br Aad et al. (2016b) as shown by the shaded region of the same figure (it will change slightly for varying ). Finally, we plot the normalized cross section of pseudoscalar pair production of our model by dark-green dots. We can see from the Fig. 7 that the search is quite sensitive to the region where other decay modes are kinematically suppressed while searches for other final states do not have any sensitivities to our model. The bound of BSM Higgs boson decay will exclude large portion of the parameter space where the coupling is not suppressed. In the limit of small , Eq. 24 can be simplified to . We find the visible points with Br should have .

In the region , the pair can be produced through resonance decay. The cross section of production is proportional to , while the Br can be large for appropriate choice of parameters in the scalar sector. The cross section of resonant pair production from gluon-gluon fusion in our model are shown by dark green points in Fig. 8. The lines in the figure correspond to the 95% CL LHC searches constraints from  ATL (2016d),  CMS (2016b) and  CMS (2016c) channels respectively. As have been done for Fig. 7, the known decay branching ratios of have been projected out. It can been seen that the search for final state provides the best sensitivity, and the search for is better than only in the low region. For a moderate mass of the pseudoscalar GeV, some parameter points are already close to the LHC search limit. We would expect those points can be probed/excluded in the near future when larger data sample is obtained.

## V Conclusion

In this paper, we propose a minimal renormalizable and gauge invariant DM model with a pseudoscalar mediator. The singlet fermion DM has only axial couplings to the pseudoscalar, while the mixing between the pseudoscalar and SM Higgs doublet leads to the interactions of DM and SM fermions and gauge bosons. Owing to the s-wave annihilation, the DM relic density can be easily obtained and the DM indirect detection signals are remarkable. The momentum suppression in DM-nucleon scattering matrix leads to null signal in all DM direct detection experiments.

We study the most up-to-date LHC search constraints on signals of the model both with and without DM in the final state. The mono-jet signature of our model is studied comparatively with that of models with pure scalar and pure axial couplings between the mediator and SM particles/DMs. Three scenarios give different predictions on the decay branching ratio of pseudoscalar/scalar to DM and the DM pair production cross section. As a result, different mono-jet search sensitivities are obtained in different scenarios. Among them, the AA scenario has the best search sensitivity at the LHC. And the sensitivity of SA is slightly better than that of AA scenario when the inference effect between two propagators is considerable. Due to the suppression in resonant vector boson pair production, the typical production cross section of resonant vector boson pair is at least two order of magnitude below the current LHC search sensitivity. The searches for resonant scalar pairs are more promising. For light GeV, the stringent limits on the BSM Higgs boson decay branching ratio obtained from Higgs precision measurements as well as the search for light bosons from 125 GeV Higgs boson decay in final state exclude very large portion of the parameter space. As for heavy GeV, the production rate is suppressed by while the can vary freely. A much better sensitivity is obtained for this channel than that for resonant channel. Some of the parameter points are less than one order of magnitude away from the current search sensitivity, thus can be probed/excluded in the near future.

Note Added: After we submitted this paper on the arXiv.org, we came to learn that the same or similar model has been considered in Ref. Ghorbani (2015). We thank Karim Ghorbani for bringing his paper to our attention.

###### Acknowledgements.
JL would like to acknowledge Jong Soo Kim for instructions on CheckMATE. This work is supported in part by National Research Foundation of Korea (NRF) Research Grant NRF-2015R1A2A1A05001869 (SB,PK, JL), and by the NRF grant funded by the Korea government (MSIP) (No. 2009-0083526) through Korea Neutrino Research Center at Seoul National University (PK).

## References

You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters