AdS/QCD approach to the scaleinvariant extension of the standard model with a strongly interacting hidden sector
Abstract
In this paper, we revisit a scaleinvariant extension of the standard model (SM) with a strongly interacting hidden sector within AdS/QCD approach. Using the AdS/QCD, we reduce the number of input parameters to three, i.e. hidden pion decay constant, hidden pion mass and that is defined as the ratio of the vacuum expectation values (VEV) of the singlet scalar field and the SM Higgs boson. As a result, our model has sharp predictability. We perform the phenomenological analysis of the hidden pions which is one of the dark matter (DM) candidates in this model. With various theoretical and experimental constraints we search for the allowed parameter space and find that both resonance and nonresonance solutions are possible. Some typical correlations among various observables such as thermal relic density of hidden pions, Higgs boson signal strengths and DMnucleon cross section are investigated. We provide some benchmark points for experimental tests.
pacs:
Valid PACS appear hereI Introduction
Although the SMlike Higgs boson has been discovered at the Large Hadron Collider (LHC) Aad:2012tfa (); Chatrchyan:2012xdj (), there are still a number of questions that call for physics beyond the SM (BSM): (i) the origin of the mass of Higgs particle or the origin of weak scale, (ii) the nature of nonbaryonic dark matter (DM), (iii) the origin of neutrino masses and mixing, (iv) matterantimatter asymmetry of the universe, to name a few.
The first question is often phrased as hierarchy problem, that addresses why the electroweak (EW) scale GeV is much smaller than the Planck scale. One of the nice ways to understand this is through quantum dimensional transmutation, which explains why the proton mass is much suppressed compared with the Planck mass in Quantum Chromodynamics (QCD) Hill:2002ap (); Wilczek:2005ez (). Technicolor (TC) provides an answer in this way, but the naive version of it is strongly disfavored by the electroweak precision test (EWPT) Peskin:1990zt ().
Since the observation of W. Bardeen Bardeen:1995kv (), softly broken scale invariance has been considered as a possible solution for the hierarchy problem. If the model is scale invariant at classical level, no dimensionful parameters are allowed and the scale symmetry is broken only logarithmically through scaleanomaly. Many authors have studies this type of models where the EW symmetry is dynamically broken via dimensional transmutation in the hidden sector with new confining strong interactions Hur:2007uz (); Ko:2008ug (); Ko:2009zz (); Ko:2010rj (); Hur:2011sv (); Heikinheimo:2013fta (); Heikinheimo:2013xua (); Holthausen:2013ota (); Jung:2014zga (); Salvio:2014soa (); Kubo:2014ova (); Kubo:2014ida (); Schwaller:2015tja (); Ametani:2015jla (); Kubo:2015joa (); Haba:2015qbz (), or ColemanWeinberg mechanism Meissner:2006zh (); Foot:2007as (); Meissner:2007xv (); Hambye:2007vf (); Foot:2007iy (); Iso:2009ss (); Iso:2009nw (); Holthausen:2009uc (); AlexanderNunneley:2010nw (); Ishiwata:2011aa (); Lee:2012jn (); Okada:2012sg (); Iso:2012jn (); Gherghetta:2012gb (); Das:2011wm (); Carone:2013wla (); Khoze:2013oga (); Farzinnia:2013pga (); Antipin:2013exa (); Hashimoto:2014ela (); Hill:2014mqa (); Guo:2014bha (); Radovcic:2014rea (); Binjonaid:2014oga (); Allison:2014zya (); Farzinnia:2014xia (); Pelaggi:2014wba (); Farzinnia:2014yqa (); Foot:2014ifa (); Benic:2014aga (); Guo:2015lxa (); Oda:2015gna (); Fuyuto:2015vna (); Endo:2015nba (); Plascencia:2015xwa (); Hashino:2015nxa (); Karam:2015jta (); Ahriche:2015loa (); Wang:2015cda (); Haba:2015nwl (); Ghorbani:2015xvz (); Helmboldt:2016mpi (); Jinno:2016knw (); Ahriche:2016cio (); Ahriche:2016ixu (); Das:2016zue (); Khoze:2016zfi ().
Some of the present authors have proposed a scaleinvariant extensions of the SM with a strongly interacting hidden sector, namely hidden QCD models Hur:2007uz (); Ko:2008ug (); Ko:2009zz (); Ko:2010rj (); Hur:2011sv (); Hur:2011sv (). At the classical level, the scale invariance is imposed so that all dimensionful parameters are forbidden in the classical Lagrangian. The Higgs mass term arises at quantum level through the dimensional transmutation driven by asymptotically free gauge theories in the hidden sector. Hidden sector couples to the Higgs through singlet scalar field only and there are stable or longlived particles (lightest hidden mesons and hidden baryons) that can make good DM candidates. In those works, hidden QCD sector was studied in the chiral effective Lagrangian approach and nonperturbative parameters were estimated by naive dimensional analysis. Then the same model was analyzed in the Nambu–JonaLasinio (NJL) approach in Ref. Holthausen:2013ota (); Kubo:2014ida ().
In this paper, we consider the same model using another approximation method, the AdS/QCD DaRold:2005zs (); DaRold:2005vr (), in order to analyze nonperturbative strong dynamics in the hidden QCD models. First we reformulate the hidden QCD sector in terms of the linear sigma model, in which the sigma and pi mesons are effective degrees of freedom. We consider a linear sigma model coupled with a scaleinvariant Higgssinglet sector and analyzed the masses and mixing between the SM Higgs boson, a singlet scalar messenger and the sigma meson. In the AdS/QCD we successfully reduce the number of free parameters by matching the mass spectra of the lightest scalar, vector and axial vector mesons.
Next, we apply this model to dark matter phenomenology. In our model, since the hidden quarks do not couple to any U(1) gauge fields, the hidden pions cannot decay through the U(1) anomaly and are found to be stable lightest particles coupling weakly with the SM fields. Hence the hidden pions become candidates of the weakly interacting massive particle (WIMP) DM. With the free parameters reduced by the AdS/QCD, we identify the parameter space that satisfies the recent observations. Then we study the distinctive features of the allowed parameter region and also some typical correlations among various observables. We address on the possible signatures of the model that can be further scrutinized in the future experiments such as LHC RunII, ILC and so forth.
This paper is organized as follows. In Sec. II, we revisit the original hidden QCD models Hur:2007uz (); Hur:2011sv () by reformulating the hidden QCD sector with the linear sigma model. Then in Sec. III, we apply the idea of the AdS/QCD to the linear sigma model described in Sec. II. In Sec. IV, numerical results on the Higgs and the dark matter phenomenologies are presented. Then Sec. V is devoted to summary and discussions.
Ii The Model
Scaleinvariant extension of the SM with a strongly interacting hidden sector contains the SM fields plus a singlet scalar and a scaleinvariant hidden QCD sector Ko:2008ug (); Ko:2009zz (); Ko:2010rj (); Hur:2011sv (). The corresponding Lagrangian is given by
(1)  
where is the mass parameter of the SM Higgs boson. We have replaced all the mass parameters (the Higgs boson mass, the RH neutrino masses and the current quark masses of hiddensector quarks) by real singlet scalar operators or following the idea of classical scale invariance. is the field strength of the hidden QCD with gauge symmetry. The SM singlet scalar couples to the hiddensector quarks through the Yukawa interaction. Since there are no dimensionful parameters in the Lagrangian, this system is scaleinvariant at the classical level. At quantum level and lowenergy scale, the hiddenQCD quarks can condensate. Such condensates induce a linear term in . Then the potential of can be tilted and can develop a VEV. The VEV of the singlet scalar generates a Higgs boson mass term , as well as the RH neutrino masses and the current quark masses of the hiddensector quarks. Thus all the mass scales in this model are generated by , which is a result of nonperturbative dynamics in the strongly interacting hidden sector. For this to happen, we assume that so that nonzero triggers the electroweak symmetry breaking.
Hereafter we consider the case in which and , for which we can use the known results from the hadronic system with and mesons. Then and for simplicity we assume the hidden quarks have isospin symmetry . In such a case the lowenergy effective theory of the hidden QCD is described by the pi meson triplets and the sigma meson. It would be written in the form of a linear sigma model
(2) 
where and represents sigma and pi meson fields.
We parameterize the VEVs and fluctuations of scalars as
(3) 
To minimize the potential energy
(4) 
three minimization conditions () should be satisfied. These conditions reduce the number of free parameters. Furthermore, two parameters , are traded with the pion mass and a sigma meson mass parameter . Hence the scalar mass matrix takes the form of
(5) 
with
(6)  
(7)  
(8)  
(9) 
where and . Since a offdiagonal part satisfies , the Higgssinglet mixing can be large when is small. and are traded with pion mass and sigma meson mass by
(10) 
The couplings , are given by
(11) 
where . Since one of the physical scalar should be the Higgs boson with mass , is one of the eigenvalues of . Thus a condition ( is a unit matrix) yields
(12) 
where we have parameterized .
The mixing matrix is defined as
(13) 
where is the SMlike Higgs boson with GeV and are extra scalar particles with .
At this stage we have four free parameters: , , and (or ). To reduce the number of free parameters, in particular, to relate the with , we use a holographic treatment of the hidden QCD.
Iii AdS/QCD analysis
In the AdS/QCD DaRold:2005zs (); DaRold:2005vr (), the hidden QCD sector is described by gauge theory on space with metric
(14) 
where and is the curvature radius of . breaks the conformal symmetry in the infrared (IR) regime, while one can take to be arbitrary small, . The nonperturbative breaking of chiral symmetry is regarded as the spontaneous breaking of symmetry by the VEV of the bulk scalar which is a bidoublet . The parity transformation corresponds to the exchange .
Then the 5D bulk Lagrangian is given by
(15) 
where , , and the bulk mass parameter is chosen so as to relate the bulk scalar field with the dimensionthree operator . The profile of the VEV is obtained by solving the zeromode equation of motion. We have
(16) 
where and can be written in terms of the value of at boundaries
(17) 
from boundary conditions
(18) 
Here nonzero corresponds to the spontaneous breaking of the chiral symmetry in the IR, while the boundary condition at corresponds to the explicit breaking of chiral symmetry. The boundary condition at is induced by the scalar potential localized on boundary. The boundary interaction is
(19) 
After the symmetry breaking , vector and axialvector gauge bosons , are expanded into KaluzaKlein (KK) modes
(20) 
where , correspond to hadronic vector and axialvector currents, respectively.
The bulk scalar is decomposed into and gauge fixing conditions are
(21) 
and have the following KK expansions:
(22) 
where and are scalar and pseudoscalar hadronic states, respectively. In particular corresponds to the pion.
In the QCD the two point correlators for the scalars and pseudoscalars are defined as
(23) 
where the twopoint correlator can be obtained from the generating function according to
(24) 
where and are the scalar and pseudoscalar external sources coupled to QCD:
(25) 
According to the AdS/CFT correspondence, the generating function is obtained by integrating bulk fields restricted to a given Ultraviolet(UV)boundary value which play the role of the external sources coupled to QCD. For the 5D scalar field we have
(26) 
( is a constant which will be determined in the matching of correlation in UV as DaRold:2005vr ()) or explicitly,
(27) 
Since the quark masses are given by , the singlet scalar fluctuation can be related to the scalar source term as
(28) 
This correspondence can be used to obtain the couplings of to the meson states.
This AdS/QCD model has five relevant free parameters: , , , and . is traded with the pion mass by the GellMann–Oaks–Renner relation
(29) 
or where . is related with the betafunction of the QCD and we fix
(30) 
where we consider the case . is related with the mass of the first KK state of which corresponds to the rho meson mass by . The value of can be fixed by adjusting the mass of the first KK of and the first KK vector meson mass
(31) 
with and , which yields DaRold:2005zs ()
(32) 
The pion decay constant is written in terms of , , as
(33) 
when . Here is the gamma function. With and one has
(34) 
which well agree with the experimental value PDG () when . Hence for and we have . The hidden rho meson mass and hidden axial vector meson mass is estimated as
(35) 
Now we fix the value of . In the original paper the author estimated and identified the lightest scalar meson as . In the present study, we regard sigma meson as lightest scalar resonance state, . In the AdS/QCD, since the wave functions of , vanish at
(36) 
there are no direct interactions between meson states with source term. In the AdS/QCD, sourcepionpion interactions are given by
(37) 
() where , and is the mass, decay constant and the coupling. These terms arise due to the mixing. In particular, the first term of r.h.s. of eq. (37) is induced by coupling through the source mixing. We assume that , and that the mixing is given by
(38) 
Together with the GellMann–Oaks–Renner relation eq. (29) we obtain
(39) 
In the AdS/QCD,
(40) 
and is obtained for , . and is obtained by formulas summarized in the Appendix. We find numerically that eq. (39) is satisfied when
(41) 
(See fig. 1),
and hence we obtain a relation (See fig. 2)
(42) 
When we take a scale normalized by so that we have , we obtain , (fig. 2).
This result well agree with the experimental bound PDG (), if we identify as the sigma meson.
Iv Hidden Pion Phenomenology
As a result of the previous sections, we now have three free parameters in this model: and . For numerical analysis, we scan the threedimensional parameter space . Note that the SMlike Higgs boson with mass GeV is termed as , and extra scalar particles as and with . We have considered several theoretical and experimental constraints.
Let us first consider the theoretical constraints. From the stability of the potential, the dimensionless couplings should satisfy the relation,
(43) 
which are translated into the following relation with the minimization conditions
(44) 
Since the strongly interacting hidden sector are now treated as linear sigma model, the symmetry breaking condition in the sector would constrain the value of to be positive. With the help of eq. (10), the condition reads
(45) 
We also adopt the perturbativity bound on with our definition of the Lagrangian:
(46) 
Experimental constraints considered in the analysis are listed in the following:

Signal strength for the SM Higgs boson Khachatryan:2014jba (); Aad:2015gba (),
(47) 
Bounds for extra scalar particles from the LEP Barate:2003sz () and the LHC CMS:2012bea (); TheATLAScollaboration:2013zha (); ATLAS:2014aga ().

Relic density from Planck satellite Ade:2015xua ().
(48) 
Neutrino signals through the DM capture by the Sun, mostly from SuperKamiokande for upward muon flux Tanaka:2011uf ().

FermiLAT 6year results for DM annihilation Ackermann:2015zua ().

Higgs invisible width from the LHC Aad:2014iia (); Chatrchyan:2014tja (),
(49) 
Direct detection bound, mostly from LUX Akerib:2013tjd (), SuperCDMS Agnese:2014aze (), and CRESSTII Angloher:2014myn ().
We apply bounds with these experimental constraints except for the relic density, for which we use the measured value as an upper bound. This is because there could be additional contributions from hidden baryons to DM thermal relic density, which we do not include in this paper. We vary the up to 2 TeV and the ranges for and are fixed with eq. (44) and eq. (45). We use micrOMEGAs Belanger:2013oya () for evaluating DMrelated observables.
The result of scanning is depicted in fig. 3. Here we see that there is definite lower bound for , around 100 GeV. Also is bounded from below, , mainly because of the perturbativity of since it can be written in the form
(50) 
So if is too small, will have very large value, above the perturbativity bound. The island on the leftmost side are the solution points where , i.e. the SMlike Higgs resonances. Other points include light scalar resonanaces with , heavy scalar resonances with and nonresonance solutions with . The nonresonance solutions favor relatively small . For example, if the hidden pion mass is away from both resonance regions more than , i.e. , then is constrained to be smaller than about . No points survive if the hidden pion mass is far away from the resonance regions, more than about . This is because when is small, the offdiagonal term of the mass matrix presented in eq. (5), , is enhanced compared with large case so that the mixing between the SMlike Higgs boson and singlet scalar fields are enhanced.
These features for the solution points can be more easily checked with another forms of plots. In fig. 4, the solution points are shown in plane. The thin branch in the left plot is the collection of solution points where , light scalar resonances. Note that there is no solution points when , i.e. all extra scalar particles are heavier than the SMlike Higgs boson. The right plot, where the points are shown in plane, also includes the thin branch that is corresponding to the solution points with heavy scalar resonances. The shape of the plot is almost same as fig. 3, since a relation generally holds in this model. As a result, has a definite lower bound as does. Both plots also include the nonresonance cases, of which DM mass are far from both resonances and other parameters are tuned to satisfy the all constraints.
As a distinctive observable, we show the deviation of the triple Higgs coupling from the SM prediction in fig. 5. The triple Higgs coupling in the model can reach of the SM prediction for relatively small . If we take larger values for and , it approaches to the SM one and cannot be a distinctive observable. Especially when is larger than 1 TeV, the triple Higgs coupling is very close to the SM prediction and the deviation cannot be detected.
Let us consider another observables. To be more specific, we separate the cases as the SMlike Higgs resonance and the other cases. fig. 6 shows the correlations between the Higgs signal strength and DMnucleon cross section in Higgs resonance case. Deviation of from 1 is generated by the mixing angle in eq. (13). In this case, the SMlike Higgs boson can decay to a pair of DM’s and these decay modes contribute to the Higgs invisible decay width. As one can see in the left plot, the invisible decay width of the SMlike Higgs boson increases when the signal strength decreases and vice versa, just like DMnucleon cross section. We can understand this by the fact that the deviation of the signal strength is determined entirely by nonzero mixing angles among the SMlike Higgs and other extra scalar particles. The more they are mixed, the more wide the channel between the visible and hidden sector is open. We show in the right plot the same correlation with relic density contours. The variation of the relic density is caused by the small variation of and as long as is close to . In both plots we apply bound for the signal strength such as .
The same correlation is shown in fig. 7 for the other cases than the SMlike Higgs resonance. As mentioned before, the cases include light and heavy scalar resonances and nonresonance solutions. In this plot, the color contour represents the mass of the DM. Unlike the Higgs resonance case, relatively large values of the signal strength are favored, with generically larger values of DMnucleon cross section. Note that there is an upper limit for . This bound is originated from the mixing, that should be different from zero for avoiding the overclosure of the universe by the DM.
fig. 8 is showing the correlation between the signal strength and triple Higgs coupling normalized with the SM prediction. The left plot is for the SMlike Higgs resonance solutions and the right one for other cases. Both cases are predicting sharp linear correlations, but with different slopes. Two solution points are disjoint with each other, so the measurements of them can be used for the identification of the scenarios, though their values are small.
Finally, some benchmark points are collected in Table 1. We classify the points as (A) SM Higgs resonance solutions, (B) light scalar resonances, (C) heavy scalar resonances, and (D) nonresonance solutions. The solution points with large relic densities (close to upper bound, ) are labelled as (I) and (II,III) correspond to the cases of small relic densities. Note that the closer we take the DM mass to the exact resonance, (or ), the smaller the relic density becomes.
type  
AI  692.19  55.27  0.7881  346.281  3460.96  0.1220  0.754  0.206  0.932  
AII  1524.33  62.49  1.648  415.13  7621.67  0.0027  0.981  0.994  
BI  1244.17  1000.08  4.288  2033.54  6224.10  0.1088  0.986  0  1.000  
BII  1607.03  999.63  5.638  2000.65  8036.53  0.0101  0.986  0  1.000  
BIII  775.72  199.22  2.761  396.81  3878.60  0.0061  0.984  0  0.998  
CI  227.49  549.85  1.479  570.38  1153.54  0.1193  0.985  0  0.988  
CII  387.01  999.03  6.867  378.17  1938.82  0.0681  0.985  0  1.000  
DI  185.66  319.49  1.337  323.97  930.67  0.1192  0.964  0  0.965  
DII  208.73  231.86  0.906  405.30  1044.98  0.0101  0.967  0  0.966 
V Summary and Conclusion
In this paper, we have analyzed the scaleinvariant extension of the SM with vectorlike confining gauge theory in the hidden sector by using the AdS/QCD proposed in Refs. DaRold:2005zs (); DaRold:2005vr (). The model contains the singlet scalar field that connects the confining hidden sector and the scaleinvariant SM sector. Hidden sector fermions develop nonzero chiral condensates and generate the linear term in the potential of the singlet scalar field . As a result, the singlet scalar field develops a nonzero VEV and it provides the tachyonic mass term for the SM Higgs field. Therefore the origin of the EWSB in the SM sector lies in the new strong dynamics in the hidden sector.
We have used the AdS/QCD approach to describe nonperturbative dynamics of the hidden QCD sector. By the AdS/QCD, strongly interacting gauge theory in the hidden sector with twoflavors can be described by gauge theory on . The spectrum of the mesonic states then can be calculated up to overall scale by considering the twopoint correlators. We first fixed the values of the AdS/QCD parameters that reproduce the known spectra of the mesons by identifying first KK mode of the vector state as rho meson. We applied the results to the hidden QCD. In this case, hidden rho meson mass is be treated as overall scale of the hidden QCD. By this, we successfully found out the relation between hidden sigma meson mass and hidden pion decay constant etc. As a result, we reduced the number of free parameters of the model to three, i.e. and .
The hidden pions can be the DM candidates since the hidden sector flavor symmetry becomes an accidental symmetry of hidden sector strong interaction. We have analyzed these “hidden pion” properties as the DM. Many results of the DM search experiments were considered. In addition to the SMlike Higgs boson, we have two extra neutral scalar fields in the model. Those extended scalar sectors are constrained by the LHC data, for example Higgs signal strengths and nonobservation of another scalar particles, etc. By scanning the threedimensional parameter space , we found that the nonresonance solutions are also possible in addition to the resonance solutions. We also considered various correlations among the experimental observables. For example, there is the correlation between the Higgs signal strength and DMnucleon cross section, and also between and the triple SMlike Higgs coupling. Especially for the latter, we found that their values and correlations behave differently depending on whether hidden pions have the SMlike Higgs resonance or not. Though the Higgs signal strength has been measured quite precisely and seems to be consistent with the SM prediction, there is still room for the physics beyond the SM as discussed in this paper. If the Higgs signal strength is measured more precisely, according to the sharp correlations we found, we can give peculiar predictions on the DM properties and others such as triple Higgs coupling etc. This could be seen in the benchmark points we presented at the end of the analysis.
Let us comment on the future prospects. Our model contains two extra neutral scalar bosons that mix with the SMlike Higgs boson. Mass spectra of those two scalar bosons are constrained by the uptodate experimental results on the Higgs signal strengths in such a way that both of them are heavier than the 125 GeV SMlike Higgs boson. Besides the resonance solutions by the SMlike Higgs, extra light and heavy scalar particles, the nonresonance solutions are also possible for moderate values of hidden pion mass and decay constant. In that case, the mass of the light extra scalar particle will be around a few hundred GeV, which can be accessible at the LHC RunII. The model also predicts the values of other observables such as relic density, DMnucleon cross section and triple Higgs coupling and so forth. Especially, the Higgs signal strength will be sharply determined. The more detailed study on the collider phenomenologies, for example, the pair production of the SMlike Higgs boson, could be possible. In addition, more complete studies with the hidden baryons, another DM candidates, can be pursued with the AdS/QCD.
Acknowledgements.
This work is supported in part by National Research Foundation of Korea (NRF) Research Grant NRF2015R1A2A1A05001869 (HH, DWJ, PK), NRF2015R1D1A1A01059141, NRF2015R1A2A1A15054533 (DWJ) and by the NRF grant funded by the Korea government (MSIP) (No. 20090083526) through Korea Neutrino Research Center at Seoul National University (PK).Appendix A AdS/QCD formulas
Scalar meson mass and decay constant are obtained DaRold:2005vr () by comparing the scalar correlator
(51)  
(52)  
(53) 
(where is the Euclidean momentum) with the correlator in Large QCD
(54) 
The masses of scalar resonances are determined by finding the poles, or
(55) 
and corresponding residues gives the scalar decay constants
(56) 
References
 (1) G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716 (2012) 1 [arXiv:1207.7214 [hepex]].
 (2) S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716 (2012) 30 [arXiv:1207.7235 [hepex]].
 (3) C. T. Hill and E. H. Simmons, Phys. Rept. 381 (2003) 235 [Phys. Rept. 390 (2004) 553] [hepph/0203079].
 (4) F. Wilczek, Int. J. Mod. Phys. A 21 (2006) 2011 [physics/0511067].
 (5) M. E. Peskin and T. Takeuchi, Phys. Rev. Lett. 65 (1990) 964.
 (6) W. A. Bardeen, FERMILABCONF95391T, C950827.3.
 (7) T. Hur, D. W. Jung, P. Ko and J. Y. Lee, Phys. Lett. B 696 (2011) 262 [arXiv:0709.1218 [hepph]].
 (8) P. Ko, Int. J. Mod. Phys. A 23, 3348 (2008) doi:10.1142/S0217751X08042109 [arXiv:0801.4284 [hepph]].
 (9) P. Ko, AIP Conf. Proc. 1178, 37 (2009). doi:10.1063/1.3264554
 (10) P. Ko, PoS ICHEP 2010, 436 (2010) [arXiv:1012.0103 [hepph]].
 (11) T. Hur and P. Ko, Phys. Rev. Lett. 106 (2011) 141802 [arXiv:1103.2571 [hepph]].
 (12) M. Heikinheimo, A. Racioppi, M. Raidal, C. Spethmann and K. Tuominen, Mod. Phys. Lett. A 29, 1450077 (2014) [arXiv:1304.7006 [hepph]].
 (13) M. Heikinheimo, A. Racioppi, M. Raidal, C. Spethmann and K. Tuominen, Nucl. Phys. B 876, 201 (2013) [arXiv:1305.4182 [hepph]].
 (14) M. Holthausen, J. Kubo, K. S. Lim and M. Lindner, JHEP 1312, 076 (2013) [arXiv:1310.4423 [hepph]].
 (15) D. W. Jung and P. Ko, Phys. Lett. B 732, 364 (2014) [arXiv:1401.5586 [hepph]].
 (16) A. Salvio and A. Strumia, JHEP 1406, 080 (2014) [arXiv:1403.4226 [hepph]].
 (17) J. Kubo, K. S. Lim and M. Lindner, Phys. Rev. Lett. 113, 091604 (2014) [arXiv:1403.4262 [hepph]].
 (18) J. Kubo, K. S. Lim and M. Lindner, JHEP 1409, 016 (2014) [arXiv:1405.1052 [hepph]].
 (19) P. Schwaller, Phys. Rev. Lett. 115, no. 18, 181101 (2015) [arXiv:1504.07263 [hepph]].
 (20) Y. Ametani, M. Aoki, H. Goto and J. Kubo, Phys. Rev. D 91, no. 11, 115007 (2015) doi:10.1103/PhysRevD.91.115007 [arXiv:1505.00128 [hepph]].
 (21) J. Kubo and M. Yamada, PTEP 2015, no. 9, 093B01 (2015) [arXiv:1506.06460 [hepph]].
 (22) N. Haba, H. Ishida, N. Kitazawa and Y. Yamaguchi, Phys. Lett. B 755, 439 (2016) doi:10.1016/j.physletb.2016.02.052 [arXiv:1512.05061 [hepph]].
 (23) K. A. Meissner and H. Nicolai, Phys. Lett. B 648 (2007) 312 [hepth/0612165].
 (24) R. Foot, A. Kobakhidze and R. R. Volkas, Phys. Lett. B 655 (2007) 156 [arXiv:0704.1165 [hepph]].
 (25) K. A. Meissner and H. Nicolai, Phys. Lett. B 660 (2008) 260 [arXiv:0710.2840 [hepth]].
 (26) T. Hambye and M. H. G. Tytgat, Phys. Lett. B 659 (2008) 651 [arXiv:0707.0633 [hepph]].
 (27) R. Foot, A. Kobakhidze, K. L. McDonald and R. R. Volkas, Phys. Rev. D 77 (2008) 035006 [arXiv:0709.2750 [hepph]].
 (28) S. Iso, N. Okada and Y. Orikasa, Phys. Lett. B 676 (2009) 81 [arXiv:0902.4050 [hepph]].
 (29) S. Iso, N. Okada and Y. Orikasa, Phys. Rev. D 80, 115007 (2009) [arXiv:0909.0128 [hepph]].
 (30) M. Holthausen, M. Lindner and M. A. Schmidt, Phys. Rev. D 82 (2010) 055002 [arXiv:0911.0710 [hepph]].
 (31) L. AlexanderNunneley and A. Pilaftsis, JHEP 1009, 021 (2010) [arXiv:1006.5916 [hepph]].
 (32) K. Ishiwata, Phys. Lett. B 710, 134 (2012) [arXiv:1112.2696 [hepph]].
 (33) J. S. Lee and A. Pilaftsis, Phys. Rev. D 86, 035004 (2012) [arXiv:1201.4891 [hepph]].
 (34) N. Okada and Y. Orikasa, Phys. Rev. D 85, 115006 (2012) [arXiv:1202.1405 [hepph]].
 (35) S. Iso and Y. Orikasa, PTEP 2013, 023B08 (2013) [arXiv:1210.2848 [hepph]].
 (36) T. Gherghetta, B. von Harling, A. D. Medina and M. A. Schmidt, JHEP 1302, 032 (2013) [arXiv:1212.5243 [hepph]].
 (37) M. Das and S. Mohanty, Int. J. Mod. Phys. A 28, 1350094 (2013) [arXiv:1111.0799 [hepph]].
 (38) C. D. Carone and R. Ramos, Phys. Rev. D 88, 055020 (2013) [arXiv:1307.8428 [hepph]].
 (39) V. V. Khoze and G. Ro, JHEP 1310, 075 (2013) [arXiv:1307.3764 [hepph]].
 (40) A. Farzinnia, H. J. He and J. Ren, Phys. Lett. B 727, 141 (2013) [arXiv:1308.0295 [hepph]].
 (41) O. Antipin, M. Mojaza and F. Sannino, Phys. Rev. D 89, no. 8, 085015 (2014) [arXiv:1310.0957 [hepph]].
 (42) M. Hashimoto, S. Iso and Y. Orikasa, Phys. Rev. D 89, no. 5, 056010 (2014) [arXiv:1401.5944 [hepph]].
 (43) C. T. Hill, Phys. Rev. D 89, no. 7, 073003 (2014) [arXiv:1401.4185 [hepph]].
 (44) J. Guo and Z. Kang, Nucl. Phys. B 898, 415 (2015) [arXiv:1401.5609 [hepph]].
 (45) S. Benic and B. Radovcic, Phys. Lett. B 732, 91 (2014) [arXiv:1401.8183 [hepph]].
 (46) M. Y. Binjonaid and S. F. King, Phys. Rev. D 90, no. 5, 055020 (2014) [Phys. Rev. D 90, no. 7, 079903 (2014)] [arXiv:1403.2088 [hepph]].
 (47) K. Allison, C. T. Hill and G. G. Ross, Phys. Lett. B 738, 191 (2014) [arXiv:1404.6268 [hepph]].
 (48) A. Farzinnia and J. Ren, Phys. Rev. D 90, no. 1, 015019 (2014) [arXiv:1405.0498 [hepph]].
 (49) G. M. Pelaggi, Nucl. Phys. B 893, 443 (2015) [arXiv:1406.4104 [hepph]].
 (50) A. Farzinnia and J. Ren, Phys. Rev. D 90, no. 7, 075012 (2014) [arXiv:1408.3533 [hepph]].
 (51) R. Foot, A. Kobakhidze and A. SpencerSmith, Phys. Lett. B 747, 169 (2015) [arXiv:1409.4915 [hepph]].
 (52) S. Benic and B. Radovcic, JHEP 1501, 143 (2015) [arXiv:1409.5776 [hepph]].
 (53) J. Guo, Z. Kang, P. Ko and Y. Orikasa, Phys. Rev. D 91, no. 11, 115017 (2015) [arXiv:1502.00508 [hepph]].
 (54) S. Oda, N. Okada and D. s. Takahashi, Phys. Rev. D 92, no. 1, 015026 (2015) [arXiv:1504.06291 [hepph]].
 (55) K. Fuyuto and E. Senaha, Phys. Lett. B 747, 152 (2015).
 (56) K. Endo and K. Ishiwata, Phys. Lett. B 749, 583 (2015) [arXiv:1507.01739 [hepph]].
 (57) A. D. Plascencia, JHEP 1509, 026 (2015) [arXiv:1507.04996 [hepph]].
 (58) K. Hashino, S. Kanemura and Y. Orikasa, arXiv:1508.03245 [hepph].
 (59) A. Karam and K. Tamvakis, Phys. Rev. D 92, no. 7, 075010 (2015) [arXiv:1508.03031 [hepph]].
 (60) A. Ahriche, K. L. McDonald and S. Nasri, arXiv:1508.02607 [hepph].
 (61) Z. W. Wang, T. G. Steele, T. Hanif and R. B. Mann, arXiv:1510.04321 [hepph].
 (62) N. Haba, H. Ishida, R. Takahashi and Y. Yamaguchi, arXiv:1511.02107 [hepph].
 (63) K. Ghorbani and H. Ghorbani, JHEP 1604, 024 (2016) doi:10.1007/JHEP04(2016)024 [arXiv:1511.08432 [hepph]].
 (64) A. J. Helmboldt, P. Humbert, M. Lindner and J. Smirnov, arXiv:1603.03603 [hepph].
 (65) R. Jinno and M. Takimoto, arXiv:1604.05035 [hepph].
 (66) A. Ahriche, K. L. McDonald and S. Nasri, arXiv:1604.05569 [hepph].
 (67) A. Ahriche, A. Manning, K. L. McDonald and S. Nasri, arXiv:1604.05995 [hepph].
 (68) A. Das, S. Oda, N. Okada and D. s. Takahashi, arXiv:1605.01157 [hepph].
 (69) V. V. Khoze and A. D. Plascencia, arXiv:1605.06834 [hepph].
 (70) L. Da Rold and A. Pomarol, Nucl. Phys. B 721 (2005) 79 [hepph/0501218].
 (71) L. Da Rold and A. Pomarol, JHEP 0601, 157 (2006) [hepph/0510268].
 (72) K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38, 090001 (2014).
 (73) V. Khachatryan et al. [CMS Collaboration], Eur. Phys. J. C 75 (2015) 5, 212 [arXiv:1412.8662 [hepex]].
 (74) G. Aad et al. [ATLAS Collaboration], arXiv:1507.04548 [hepex].
 (75) R. Barate et al. [LEP Working Group for Higgs boson searches and ALEPH and DELPHI and L3 and OPAL Collaborations], Phys. Lett. B 565 (2003) 61 [hepex/0306033].
 (76) CMS Collaboration [CMS Collaboration], CMSPASHIG13027.
 (77) The ATLAS collaboration [ATLAS Collaboration], ATLASCONF2013067, ATLASCOMCONF2013082.
 (78) G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 92 (2015) 1, 012006 [arXiv:1412.2641 [hepex]].
 (79) P. A. R. Ade et al. [Planck Collaboration], arXiv:1502.01589 [astroph.CO].
 (80) T. Tanaka et al. [SuperKamiokande Collaboration], Astrophys. J. 742 (2011) 78 [arXiv:1108.3384 [astroph.HE]].
 (81) M. Ackermann et al. [FermiLAT Collaboration], arXiv:1503.02641 [astroph.HE].
 (82) G. Aad et al. [ATLAS Collaboration], Phys. Rev. Lett. 112 (2014) 201802 [arXiv:1402.3244 [hepex]].
 (83) S. Chatrchyan et al. [CMS Collaboration], Eur. Phys. J. C 74 (2014) 2980 [arXiv:1404.1344 [hepex]].
 (84) D. S. Akerib et al. [LUX Collaboration], Phys. Rev. Lett. 112 (2014) 091303 [arXiv:1310.8214 [astroph.CO]].
 (85) R. Agnese et al. [SuperCDMS Collaboration], Phys. Rev. Lett. 112 (2014) 24, 241302 [arXiv:1402.7137 [hepex]].
 (86) G. Angloher et al. [CRESSTII Collaboration], Eur. Phys. J. C 74 (2014) 12, 3184 [arXiv:1407.3146 [astroph.CO]].
 (87) G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 185 (2014) 960 [arXiv:1305.0237 [hepph]].