Minimal Signatures of Naturalness

We study the naturalness problem using a model independent bottom up approach by considering models where only the interaction terms needed to cancel the Higgs quadratic divergences are present. If quadratic divergences are canceled by terms linear in the Higgs field, then the collider phenomenology is well covered by current electroweakino and fourth generation searches. If quadratic divergences are canceled by terms bilinear in the Higgs field, then the signatures are highly dependent on the quantum numbers of the new particles. Precision Higgs measurements can reveal the presence of new particles with either vevs or Standard Model charges. If the new particles are scalar dark matter candidates, their direct and indirect detection signatures will be highly correlated and within the reach of XENON100 and Fermi. Observation at one of these experiments would imply observation at the other one. Observable LHC decay channels can also arise if the new particles mix with lighter states. This decay channel involves only the Higgs boson and not the gauge bosons. Observation of such decays would give evidence that the new particle is tied to the naturalness problem.


Introduction
The LHC has recently discovered a Higgs like boson with mass 125 GeV [1,2]. In the Standard Model (SM), the Higgs boson receives quadratically divergent contributions to its mass from radiative corrections. The biggest contributions to these quadratic divergences come from top quark and gauge boson loops and are shown in Fig 1. Such radiative corrections will depend on the new physics scale, being proportional to the square of the cutoff scale. If new physics appears only at the GUT scale, these radiative corrections to the bare Higgs mass will be of the order of 10 32 GeV. The observed Higgs mass would then result from a miracle cancellation between the bare Higgs mass and the radiative corrections, requiring a large amount of fine tuning of UV parameters. Various models such as supersymmetry, extra dimensions [3,4], Little Higgs [5,6], Twin Higgs [7] and folded SUSY [8] have all been proposed to solve this naturalness issue, but the existence of other possible solutions should not be ignored. The current solutions give rise to very diverse signatures at colliders and dark matter detectors. However, most of these signatures are uncorrelated with the cancelation of the Higgs quadratic divergence.
There have been bottom up approaches to naturalness in the context of SUSY [9,10,11,12,13] or Little Higgs [14,15,16]. However, as mentioned before, the phenomenology of these models is not necessarily tied to the cancelation of the Higgs quadratic divergences. Applying Occam's razor to models of naturalness, the simplest solution of the naturalness problem would cancel quadratic divergences without any additional interactions or matter content 1 .
Thus in this paper we focus on the signatures of interactions necessary for the cancelation of quadratic divergences. These signatures vanish in the limit where the quadratic divergences are not canceled. This approach leads us to consider the only two terms in the Lagrangian which contribute to the 1-loop quadratic divergences of the Higgs mass, where the quantum numbers of the fields ψ i are only constrained by gauge invariance. Ideally, these couplings would be measured at a collider and shown to cancel the quadratic divergences induced by the SM fields. λ 1 can be directly observed through the production and decay of ψ i . However, even proving the existence of λ 2 is difficult, since it is not directly related to the production or decay of the ψ i . Direct observation at a low luminosity LHC is not feasible. Depending on the quantum numbers of ψ i , this term may be the only renormalizable interaction, preventing any other means of detection. 1 Anthropics may favor using additional quartic interactions to fine tune the Higgs mass rather than using the bare mass. As generic particles have no other interactions with the SM, the resulting theory cancels quadratic divergences with no additional interactions.  This article studies the signatures of models where either λ 1 or λ 2 are used to cancel quadratic divergences. Models which use λ 1 to cancel quadratic divergences have signatures that are similar to MSSM electroweak phenomenology [17] or fourth generation models [18,19]. Models which use λ 2 to cancel quadratic divergences can have interesting correlated dark matter direct and indirect signatures. The λ 2 term can also lead to up to 10% modifications of the Higgs production cross section or of its decay width to gauge bosons.
Predicting other, more visible, collider signatures coming from the λ 2 term requires introducing additional terms in the Lagrangian. Generically, these terms would lead to signatures which are not directly related to the quartic term of Eq. 1.1. Introducing mass mixing with a lighter state is the only way of getting an observable decay channel for ψ i without making additional assumptions on the UV symmetry enforcing naturalness. If ψ i mixes with a lighter field, then λ 2 induces a decay of ψ i to this field and a Higgs boson.
This article is organized as follows. Sec. 2 analyzes the signatures of simplified models which use a coupling linear in the Higgs field to cancel quadratic divergences. Sec. 3 studies simplified models which use a coupling bilinear in the Higgs field to cancel quadratic divergences. Sec. 4 studies signatures that result from mass mixing and the bilinear term in Eq. 1.1. We conclude in Sec. 5.

Simplified models with couplings linear in H
We consider minimal extensions of the SM where there is only the term L ⊃ λHψ 1 ψ 2 (2.1) in the Lagrangian. In order for this coupling to contribute to the Higgs quadratic divergence at 1-loop, we need ψ 1 and ψ 2 to be fermions. For a hard UV cutoff, the 1-loop quadratically divergent contribution to the Higgs mass is always negative and is − |λ| 2 Λ 2 8π 2 . Thus if this term in the Lagrangian is to cancel quadratic divergences, it must cancel the positive divergences due to the gauge bosons. Due to gauge invariance, either ψ 1 or ψ 2 has to be charged under SU (2). The minimal models considered in this section are labeled by whether ψ 2 is a SM particle or not. Other non-minimal signatures can be obtained by combining these models.
If both ψ 1 and ψ 2 are non-SM particles, they are pair produced through intermediate gauge bosons and then decay through gauge bosons and Higgses. In general, the associated decay pattern closely resembles the SUSY signatures involving the gauginos [17]. If the lighter of the ψ is colored or charged, then we can have the pair production of R-hadrons [22,23,24], CHAMPS [25,27] or electroweakino-like decays ending in R-hadrons or CHAMPS. The LHC signatures of this scenario are decays that involve gauge and/or Higgs bosons and either missing E T , CHAMPS or R-hadrons. The similarity of these signatures with electroweakino decays is not an accident as these are the couplings that SUSY theories use to cancel the gauge quadratic divergences.
If only ψ 1 is a new particle, it can be both pair produced and singly produced. Its pair production cross section is determined by its gauge quantum numbers while its single production cross section is determined by what standard model particle ψ 1 decays into. e.g. if ψ 1 is charged electroweakly and decays into leptons, then it is both pair and singly produced through intermediate W/Zs. Since canceling the gauge quadratic divergences requires gauge and Yukawa couplings to be of the same order, both single and pair production matrix elements are comparable. Because ψ 1 is heavier than its Standard Model partner, we will see single production before pair production. For ψ 1 which are colored, the pair production through gluons may be seen before the single production through W/Zs. The cross section for single production depends on λ in Eq. 2.1 so that its magnitude can be used to determine which gauge divergence is canceled. The decay of ψ 1 is the same as in other fourth generation models [18,19] or models with new vector-like fermions such as Little Higgs models [5,6].
If ψ 1 is lighter than ψ 2 , then the SM particle ψ 2 can have new three body decay modes. The cases of ψ 1 being colored or having an electroweak charge are both ruled out by the R-hadron [22,23,24] and LEP bounds [26] respectively. If ψ 1 is a SM singlet, then it is neutrino like and after the Higgs obtains a vev, Eq. 2.1 generates a large mass term for the neutrinos. Thus we find that without making additional assumptions, this scenario is also excluded.

Simplified models with couplings bilinear in H
In this section, we consider couplings of the form where the Higgs vev v will be defined as v = 246 GeV all throughout this article. For a hard UV cutoff, the quadratic divergence can be either positive or negative depending on the sign of λ. Thus we have the option of canceling any quadratic divergence induced by either gauge bosons or fermions.
ψ is either a scalar or a vector-like fermion but does not have to be charged under the SM. If ψ is lighter than half the mass of the Higgs, the second term of Eq. 3.2 opens a new decay channel for the Higgs. If ψ is a scalar, it can get a vev and mix with the Higgs after EWSB. If ψ is charged under the SM, then it modifies the Higgs couplings to the gauge bosons at 1-loop. Finally, if ψ is a SM singlet, then it is a dark matter candidate and can be accessible at the current dark matter detection experiments. All these possibilities give rise to various signatures, summarized in Tab. 2.
If the new particles are lighter than half the Higgs mass, then they can give rise to Higgs alternate decay modes. Light stable standard model singlets or neutral components of standard model triplets will lead to invisible Higgs decay modes. If these particles cancel top or gauge quadratic divergences then the decay widths into these new particles would be orders of magnitude above the decay width into bottoms. Current LHC searches rule out decays widths this large [34]. Canceling other quadratic divergences yields modifications too small to observe.
Light electrically charged of colored particles have been ruled out by LEP [26] and R-hadron searches respectively. Decays of new light SU (2) charged particles to standard model states would lead to model dependent modifications of the Higgs branching ratios with respect to their standard model values. If the final states associated to h → ψψ do not significantly contribute to any of the main Higgs search channels, all the detection rates in these channels will be uniformely suppressed. In this last case, current LHC Higgs coupling measurements rule out the possibility of cancelling the top or gauge quadratic divergences [48, 29].

ψ = 0
This subsection considers the case where the new particle which cancels quadratic divergences is a scalar that obtains a vev. If this scalar is an SU (2) doublet, then the model is a two Higgs doublet model that satisfies the Veltman conditions [47]. If the new particle is an SM singlet, the minimal Lagrangian is The production of φ m occurs only through mixing with the Higgs and is proportional to sin 2 α ≤ 14%. If φ m is lighter than twice the Higgs mass, it acts like a heavy SM Higgs boson with uniformly suppressed couplings; the signal strength is suppressed by sin 2 α in all channels. If φ m is heavier than twice the Higgs mass, φ m can decay into 2 Higgses. By the goldstone boson equivalence theorem, the branching ratios into Ws/Zs/hs is 2:1:1 in the large mass limit. The best limit on heavy SM Higgs arises from H → ZZ → llll [49, 50]. The bounds on sin 2 α are generally of about 20% and can be as low as 10% for very specific masses. Thus we find that precision Higgs physics is the best method to constrain this class of models.

Scenario
Signature m ψ < m h 2 Higgs decays to missing E T , CHAMPS, R-hadrons ψ = 0 suppressed Higgs couplings and a heavy Higgs-like scalar ψ charged under SM O(10%) changes to H → γγ/gg ψ is dark matter Correlated direct and indirect detection signatures Table 2: The scenarios and collider signatures that results from minimal models of naturalness which cancel the quadratic divergences through the term shown in Eq. 3.1. The quantum numbers and spin of ψ are allowed to vary when not specified.

ψ is charged under the SM
If ψ has SM quantum numbers, the three point couplings in Eq. 3.2 can lead the new particles running in loops to contribute to the Higgs decays to gauge bosons. Measuring the deviation of these decay rates from their SM values gives information as to which quadratic divergence could be cancelled. As an example, consider a singlet fermion which cancels the top quadratic divergence and has charge 1. The relevant terms in the Lagrangian are The modification of the Higgs decay rate to two photons is given by the low energy theorem [35,36] to be If this uncolored top partner had a mass of 600 GeV, the diphoton decay width of the Higgs will be larger than in the standard model by 10.5%. The current ATLAS and CMS 95% confidence limits on the Higgs to diphoton rate allow excesses of up to 90% for ATLAS and 20% for CMS [30,31]. Although 10% modifications of the Higgs to diphoton rate should be within reach of the future LHC precision Higgs measurements, determining the Higgs couplings to other gauge bosons with a similar precision would take much longer.

ψ is a dark matter candidate
If ψ is a dark matter candidate, Eq. 3.2 shows that after EWSB, Higgs exchange with nuclei and annihilation into H → bb → π 0 s giving many photons gives rise to direct and indirect detection signatures respectively. The dark matter particle's coupling to the Higgs boson is For fermions/complex scalars, we have where N c is the color factor (N c = 3) and N s , N f are the multiplicity of scalar and fermionic dark matter candidates. Dark matter has been considered to interact using this term in other contexts [39,40,41]. For simplicity, we will be assuming a single complex dark matter field. Direct detection constrains the spin independent cross section where a = 4 if X is a real scalar, a = 1 if X is a majorana fermion or a complex scalar and a = 1/4 if X is a dirac fermion [42]. We used the values f T u = f T d = 0.025 and f T s = 0.0532 as suggested by lattice studies [43]. The results are shown in Fig. 2 and show that sub-TeV singlet dark matter particles canceling the top quadratic divergences are excluded by XENON100. Continuum photons constrain the annihilation rate of dark matter into bottoms/Ws/Zs. Fermi provided the bound on dark matter annihilation into a pair of bottoms [33] while, in this particular model, dark matter annihilates into a pair of Higgses. Higgses dominantly decay into bottoms and Ws which in turn give many pions and hence many photons. Using pythia, it was found that the number and distributions of the photons from dark matter annihilating into bottoms and Higgses were within 20% of each other. Therefore, we compare the dark matter annihilation rate into a pair of Higgses to the bound on dark matter annihilation into a pair of bottoms. Fig. 2 gives a rough estimate on what the exclusion limits should be. microOMEGAs was used to calculate the annihilation cross sections. For fermions, the vanishing of the cross section for annihilation into Higgs bosons can be understood because the J = 0 initial state of two fermions has CP = −1 while the final state of two identical scalar particles cannot have spin 0 and CP = −1. Thus the annihilation amplitude must vanish in the v = 0 limit.
If the dark matter sector solves the Higgs naturalness problem, it leads to unique predictions. Fermionic dark matter is visible in direct detection experiments but not in indirect detection experiments. Scalar dark matter has a much more distinct signature, which is that the ratio of direct to indirect detection cross sections is in the large dark matter mass limit. If scalar dark matter cancels the top quadratic divergences, it annihilates too efficiently so obtaining the correct relic abundance requires non-thermal cosmology. Fermionic dark matter annihilates less efficiently as its cross sections are velocity suppressed so it can potentially give the correct thermal relic abundance. Unfortunately, the parameter points which give the correct abundance are ruled out by direct detection constraints. Relic abundances were calculated using micrOMEGAs [46].
As can be seen from Fig. 2, sub-TeV singlet scalar dark matter cannot cancel the top quadratic divergences; if it cancels gauge quadratic divergences, then its detection could be just around the corner for both direct and indirect detection experiments. Simultaneous detection in both experiments with cross sections obeying Eq. 3.14 would be strong evidence that dark matter is a scalar involved in the cancelation of quadratic divergences and has a non-thermal production mechanism.
Direct evidence of new particles canceling the Higgs quadratic divergences can be obtained through dark matter direct and indirect detection experiments or through LHC precision Higgs measurements.  Figure 2: The abundance, direct and indirect detection bounds on dark matter particles which cancels quadratic divergences. The indirect detection bound is the Fermi bound on XX → bb [33], while the annihilation channel is into two Higgses. Using pythia, it was seen that the number and distributions of photons from a pair of Higgses and a pair of bottoms was within 20% of each other. Thus the plot can be used to find the approximate indirect detection bound. The direct detection bound is taken from XENON100 [32] and the relic abundance is taken from [44]. Bounds are placed assuming X makes up all of the dark matter. The blue line signals a cancelation of the top quadratic divergence while a red line is the gauge quadratic divergence. The black lines are the current bounds. Solid lines are complex scalar dark matter while dashed lines are dirac fermion dark matter.
we examine how the current LHC searches bound the simple case of top quark, light quark ner. Seeing a signal in a multiple higgs search would be strong evidence for a new particle uadratic divergences. Details about Monte Carlo generation and exclusions are in App. A. nce of an explicit model to study, it is not clear which of the three di↵erent three body ortant. In the example of a little higgs theory at large mass, the goldstone boson equivas a good approximation. While this theorem does not hold at the lower masses of interest, that it does to obtain a handle on the various decay modes. Thus we consider the case rix elements for (T ! W + W t), (T ! ZZt) and (T ! hht) are 2:1:1 respectively 1 . th will then vary from this ratio due to phase space factors.
any di↵erent searches that all have di↵erent reaches. As the dominant decay involves bosons and higgses, we expect that the most important searches are those involving s and b jets. nce of an explicit model to study, it is not clear which of the three di↵erent three body ortant. In the example of a little higgs theory at large mass, the goldstone boson equivas a good approximation. While this theorem does not hold at the lower masses of interest, that it does to obtain a handle on the various decay modes. Thus we consider the case rix elements for (T ! W + W t), (T ! ZZt) and (T ! hht) are 2:1:1 respectively 1 . th will then vary from this ratio due to phase space factors.
any di↵erent searches that all have di↵erent reaches. 7 divergences should be studied. The particular clas quartic coupling responsible for canceling quadrati particles. These models are unique because discove have something to do with the quadratic divergenc m mix the simplified models of Sec. 3. Finally, we conclude in Sec. 5.

Simplified models with couplings linear in H
We consider simplified models where there is the term L H¯ 1 2 in the Lagrangian. In order for this coupling to contribute to the Higgs quadratic divergence we need 1 and 2 to be fermions. This quadraticly divergent contribution to the Higgs mass negative and is | | 2 ⇤ 2 8⇡ 2 . Thus if the new particles are to cancel the quadratic divergences, cancel the positive divergences due to the gauge bosons.
If both 1 and 2 are both new particles and the only terms in the Lagrangian are those cancel quadratic divergences, then we notice that the lightest new particle is stable. There is a which all of the new particles are odd while all of the SM particles are even. These particles t through gauge bosons and higgs particles yielding a cascade that is similar to the electroweakin from supersymmetric theories. If the lightest is uncharged under the SM, then it makes a matter candidate and the cascade is exactly like an electroweakino cascade. If is colored, th than missing energy, the cascade will end in a R-hadron, though this option is cosmologically d 1 with the Higgs to cancel quadratic divergences. Sec. 3 outlines simplified models which use a bilinear coupling with the Higgs to cancel quadratic divergences. Sec. 4 shows what are the current bounds on the simplified models of Sec. 3. Finally, we conclude in Sec. 5.

Simplified models with couplings linear in H
We consider simplified models where there is the term in the Lagrangian. In order for this coupling to contribute to the Higgs quadratic divergence at 1-loop, we need 1 and 2 to be fermions. This quadraticly divergent contribution to the Higgs mass is always negative and is | | 2 ⇤ 2 8⇡ 2 . Thus if the new particles are to cancel the quadratic divergences, they must cancel the positive divergences due to the gauge bosons.
If both 1 and 2 are both new particles and the only terms in the Lagrangian are those needed to cancel quadratic divergences, then we notice that the lightest new particle is stable. There is a Z 2 under which all of the new particles are odd while all of the SM particles are even. These particles then decay through gauge bosons and higgs particles yielding a cascade that is similar to the electroweakino cascades from supersymmetric theories. If the lightest is uncharged under the SM, then it makes a good dark matter candidate and the cascade is exactly like an electroweakino cascade. If is colored, then rather than missing energy, the cascade will end in a R-hadron, though this option is cosmologically disfavored.
1 is organized as following. Sec. 2 outlines simplified models which use a linear coupling to cancel quadratic divergences. Sec. 3 outlines simplified models which use a bilinear he Higgs to cancel quadratic divergences. Sec. 4 shows what are the current bounds on odels of Sec. 3. Finally, we conclude in Sec. 5. ed models with couplings linear in H plified models where there is the term an. In order for this coupling to contribute to the Higgs quadratic divergence at 1-loop, 2 to be fermions. This quadraticly divergent contribution to the Higgs mass is always Thus if the new particles are to cancel the quadratic divergences, they must ive divergences due to the gauge bosons. nd 2 are both new particles and the only terms in the Lagrangian are those needed to c divergences, then we notice that the lightest new particle is stable. There is a Z 2 under new particles are odd while all of the SM particles are even. These particles then decay osons and higgs particles yielding a cascade that is similar to the electroweakino cascades etric theories. If the lightest is uncharged under the SM, then it makes a good dark te and the cascade is exactly like an electroweakino cascade. If is colored, then rather ergy, the cascade will end in a R-hadron, though this option is cosmologically disfavored. 1 1 Figure 3: If there is mass mixing, the term responsible for canceling quadratic divergences can give rise to decay channels. This signature is the unique indirect signature that one is led to in a bottom up approach to finding naturalness.
Probing naturalness through precision Higgs measurements will take many years. Some of the fields ψ in Eq. 3.1 have also none of the characteristics mentioned in Tab. 2. Such fields would have no visible signatures in current experiments if only Eq. 3.1 is present. While as shown in [45], future experiments such as the ILC can probe naturalness, the LHC has limited reach. In order to obtain visible signatures at the LHC, one is led to tying the cancelation of quadratic divergences to the decays of new particles.

LHC signatures with minimal assumptions
In any theory where the Higgs quadratic divergences cancel, the bilinear term is related to the Yukawa couplings by a symmetry. If a specific model is assumed, then the structure associated to this symmetry can be tested to obtain indirect evidence for the cancelation of quadratic divergences. e.g. in SUSY the gauginos have gauge strength trilinear interactions or Higgsinos having SM Yukawa strength trilinear interactions. Instead, motivated by our bottom up approach, we are lead to look for IR effects independent of the UV symmetry which can shed light on the bilinear term. Additional interaction terms in the Lagrangian typically induce phenomenology of their own, independent of whether the quadratic divergences are canceled or not. However there is a unique effect which can combine with quadratic divergences to yield detectable signals, mass mixing. As shown in Fig. 3, combining mass mixing with the bilinear interaction term given in Eq. 4.1 yields a decay channel. Observation of this new decay channel gives indirect evidence that the new particle observed is involved with the Higgs quadratic divergences.
While non-renormalizable terms are necessarily UV dependent, renormalizable terms are not. e.g. the soft IR breaking of SUSY is expressed in terms of renormalizable soft terms whereas the UV supersymmetry in manifest in the marginal terms 2 . We can then use the renormalizable mass mixing term, that is potentially unconstrained by the UV symmetry, in combination with the bilinear term to give an observable decay channel. It should be noted however that one could simply have started with a term ψ 1 ψ 2 HH † in the Lagrangian which emulates this effect. As such, finding this decay channel provides only indirect evidence for the cancelation of quadratic divergences rather than direct evidence. ψ 1 and ψ 2 in Fig. 3 could both be new particles or one could be a SM particle. They decay in cascades involving only Higgses and not the associated Ws or Zs that one expects from the familiar trilinear Yukawa interactions. If the ψ 1,2 are charged under the SM gauge group, then they can be pair produced via gauge bosons. They decay to a Higgs and either a SM particle, a CHAMP or a R-hadron. If the ψ 1,2 are not charged under the standard model, then additional assumptions are required for how they might be produced at the LHC.
If ψ 2 is a SM field, then ψ 1 has the same quantum numbers as a Standard model particle and we have a vector like fourth generation. As an explicit example, assume that ψ 1 has the same quantum numbers as a right handed up-type quark. Writing down all interactions up to dimension 5, we have where i and j are flavor indices. ψ 1 has tree level mixing with the quarks. Working in the small vev limit and diagonalizing the mass terms in the small vev limit yields The mass eigenstate of ψ 1 is called U . λ U U is related to the structure of the quadratic divergences and is not directly measurable. Due to mixing effects it is related to λ U u . In the case where λ 3 = 0 or λ 4 = 0, as in Little Higgs models, we notice that λ U u and λ U U are directly related by mixing angles. Observation of the decay channel U → uHH † gives strong indirect evidence that λ U U = 0 and that the particle U may be involved in canceling the quadratic divergences. This evidence is indirect because it is always possible that λ 3 and λ 4 are chosen exactly such that λ U u = 0 while λ U U = 0. Generically, this miraculous cancelation does not occur.
The new particle U has decay channels resulting from λ U u and λ i U . If λ i U > λ i U u v m U , the decays from λ i U are dominant, and the situation is simply that listed in Sec. 2 where a term linear in H dominates the phenomenology. Thus, we will focus on the case where λ U u dominates over λ i U . This limit can be realized in Little Higgs theories as shown in App. B.
Aside from the desire to categorize all models solving the naturalness problem, there are two other reasons why one might expect that the signatures associated with λ U u dominate. Generically, we expect O(1) mixing with no ad hoc cancelations so that we have λ U ∼ λ SM 1 and λ U U ∼ λ U u ∼ λ 2 SM 2 2m U , where λ SM 1 and λ SM 2 are the Yukawa couplings associated with the SM particle ψ mixes with and the quadratic divergences cancelled by ψ respectively. In the case where ψ is mixing with a light quark, but is canceling the top quadratic divergence, we get that λ i U is suppressed by the light quark Yukawa so that λ i U u dominates.
Another reason why λ i U terms may be suppressed is that λ i U plays a role in flavor physics and contributes to flavor changing neutral currents (FCNC). In fourth generation and Little Higgs models, the typical assumption is that there is some ad hoc texture such that FCNCs are all avoided. Analogously, we can assume that λ i U is small due to flavor physics so that FCNCs are not an issue. λ i U small also implies that the STU [38] parameters are small. λ i U u in Eq. 4.3 can give four decay channels. The decay channels are The three body decays can be the dominant or subdominant decay mode of U . Their relative branching ratios are shown in Fig. 4. Unlike two body decays resulting from a Yukawa term, two body decays resulting from Eq. 4.1 involve only the Higgs and not the gauge bosons. As shown in App. A, this unique structure is important for maintaining the goldstone boson equivalence theorem. At high m U , the goldstone boson equivalence theorem implies that there will be only three body decays into pairs of Ws, Zs and hs with a ratio of 2:1:1. This structure is specific to the interaction term in Eq. 4.1. Observation of three body decays into a pair of Higgs bosons and another particle is evidence that a new particle couples to the Higgs with a bilinear coupling.
Although the previous example illustrates the case of a new particle mixing only with an up-type quark, mixing with any other SM particle leads to similar conclusions, as shown in Table 4.4. For masses that the LHC can probe, the two body decays will be seen first. The new particles will be pair produced through intermediate gauge bosons.
If there are two new particles which mix with each other rather than with the SM, the associated collider signatures will be cascade decays through where the final state H † H can stand for h, W + W − , Z 0 Z 0 , and hh. ψ 2 can result in missing E T , a CHAMP, or a R-hadron depending on its quantum numbers. In supersymmetric models, none of the new Quantum numbers SM particle LHC final states (SU (3), SU (2)) U (1) spin that is mixed with  Fig. 4. For low masses, the signature of canceling a quadratic divergence is the observation of a decay to a Higgs but not to W or Z bosons.
particles have the same quantum numbers as each other. After electroweak symmetry breaking left and right handed sparticles can mix, but the effects of these mixings are too small to observe.

Collider bounds on top quark, light quark and lepton partners
This subsection studies the sensitivity of the current ATLAS searches to the decays of quark and lepton partners of the form Pair produced quark or lepton partners give rise to collider signatures involving two Higgs bosons and two Standard Model particles. Due to the wide variety of possible decay channels for a 125 GeV Higgs boson, there are many final states, often involving one or more lepton and/or b-jet. Missing E T and leptons are needed to reduce the QCD background so that the identity of the Standard Model particle produced in association with the Higgs determines which search will be the most sensitive to the signal studied.
Each of the processes shown in Eq 4.9 is studied separately and only signals corresponding to pair produced quark and lepton partners are investigated. In what follows, hadronically decaying τ leptons are treated as jets. In the simulated samples, the leptonic partners L decay to electrons and Higgs bosons and the light quark partners decay to up quarks and Higgs bosons. Details about the event generation and search validation are given in App. C.  8%  33%  34%  3 leptons  6%  13%  0%  4 leptons  3%  3%  0%   Table 4: Probabilities (in %) of getting a given number of leptons and b-quarks from decays of pair produced electron or muon partners. Detector effects and particle identification efficiencies are not included.

Top quark partner
requiring exactly one lepton and four or more b-jets. The associated mass reach is about 800 GeV if T → th is the only available decay channel.

Electron partner
Decays of pair produced electron partners to an electron and a Higgs boson have particularly characteristic final states with a large number of leptons and b-jets, as shown in have extremely low background rates and are sensitive to processes like the one studied. In addition to the four lepton requirement, the effective mass m ef f is required to be bigger than 600 GeV. m ef f is defined as As can be seen from the distributions in Fig 5, the efficiencies of these cuts are very high even for lepton partner masses as low as m L = 500 GeV. For m L ≥ 500 GeV, the 95% confidence limit on the lepton partner production cross section has very low sensitivity to m L and is O(40) fb. For lower masses, signal regions with 50 and 75 GeV missing E T cuts start being more competitive and the corresponding exclusion bounds are shown in Fig. 6. The electron partner pair production cross section has been computed at leading order (LO) using MadGraph 4.5.1. As shown in Fig. 6, getting reasonable sensitivity to L → l+H processes requires either going to higher energy and luminosity or designing a dedicated search. Three lepton searches [53] and searches requiring b-jets produced in association with two leptons [54, 55] could also have sensitivity to the final states shown in Table 4. Most of the final states with three leptons are characterized by an opposite sign same flavor lepton pair (OSSF) with large invariant mass. The corresponding search region in the ATLAS 3 lepton search also requires m T > 110 GeV, with m T being the transverse mass of the lepton not belonging to the OSSF pair and the missing E T . Since in most of the final states, the third lepton and the missing E T come from the decay of a W or a τ , this transverse mass requirement suppresses most of the signal. In the searches for two leptons produced in association with b-jets, the two leading leptons are required either to have the same sign or to come from the decay of an on-shell Z. The two leptons directly produced through the decay of the two lepton partners do not satisfy these requirements.   Figure 6: 95% confidence limits on the pair production cross section of an electron partner which decays to a Higgs and an electron with a branching ratio of 1 (blue). The black line shows the LO pair production cross section of the electron partner

Up quark partner
Unlike top quark partners, no existing LHC search studies the decays of up quark partners to Higgs bosons and jets. The current exclusion bounds on quark partner masses are obtained by assuming a non-zero branching ratio to W bosons plus jets [56]. The associated cuts heavily depend on the decay topology and cannot be transposed to the case studied here. Another search, specific to heavy bottom partners has been performed [57], but it requires final states with exactly two leptons of the same sign, which can be obtained only through B → W t.
Decays of pair produced light quark partners lead to the production of two Higgs bosons and two jets. Higgs searches, and especially searches looking for Higgs bosons produced through vector boson fusion (VBF), seem like they may be suited to this kind of signal. However, searches looking for Higgs bosons produced in association with two jets through VBF require a large rapidity gap. As can be seen from Fig. 7, this separation requirement proves particularly harmful for our signal. In all of the Higgs searches except the H → γγ search, signal regions where the dominant single Higgs production mode is not VBF either require at most one jet or have exclusion bounds that are too loose at 125 GeV.
The ATLAS h → γγ search [58] sets reasonable exclusion bounds on the signal strength at 125 GeV and most of its signal regions do not have specific requirements on jets. This search is then the Higgs search providing the best sensitivity to our signal. A 95% confidence bound on the UŪ production cross section can be derived by requiring that we do not produce more photons than the single SM Higgs does. Since the current exclusion bound on the signal strength for a 125 GeV Higgs is µ ∼ 1.7 and the single SM Higgs production cross section is about 20 pb, it is not surprising that our exclusion bounds on the UŪ production cross section are of the order of 10 pb. Since no tight kinematic cuts are applied in the event preselection, these bounds are roughly independent of m U . The associated mass reach is about 300 GeV.
Decays of the Higgs bosons through taus or vector bosons can lead to signals with 3 and 4 leptons. Unfortunately, due to the low branching ratio and the limited lepton identification efficiencies, the bounds found using the ATLAS 3 and 4 lepton searches are looser than the ones obtained using the h → γγ search.

Conclusion
This article studies model independent signatures of solutions to the naturalness problem. The Higgs 1-loop quadratic divergences can be cancelled by new particles interacting with the Higgs through either a Yukawa term or a quartic term. The first case is being probed by existing electroweakino or fourth generation searches. In the second case, if the new particle is a dark matter candidate, direct evidence of the existence of a bilinear quartic term could be around the corner for dark matter direct detection experiments. For scalar dark matter, a measurable correlation between direct and indirect detection signals is also predicted. More generally, if the particles canceling the quadratic divergences are charged under the SM gauge group, they could lead to modifications of the Higgs couplings to gauge bosons of up to 10%. Alternatively, if the new particles obtain a vev, the Higgs couplings are suppressed.
In the case of a single new particle, we find that introducing mass mixing with a SM particle is the sole model independent method for obtaining an observable collider signature linked with quadratic divergences. This mass mixing gives rise to two and three-body decay modes involving Higgs bosons. The two body decays arising from the quartic bilinear term dominate at sub-TeV masses and only involve Higgs bosons and not gauge bosons.
Particles mixing with the top quark and decaying to a top quark and a Higgs are already constrained by the LHC to be heavier than 800 GeV. Observing particles decaying to leptons or other quarks would require either going to higher energy and luminosity or designing dedicated searches. At the high luminosity LHC, new searches could also isolate the three body decays of the heavy new particles to a pair of Higgses or gauge bosons and another SM particle. A thorough study of the associated couplings and of possible correlations between the two and three-body decay rates would shed light on whether the decaying particle could cancel the Higgs quadratic divergences.
rrent LHC searches bound the simple case of top quark, light quark a multiple higgs search would be strong evidence for a new particle etails about Monte Carlo generation and exclusions are in App. A. el to study, it is not clear which of the three di↵erent three body of a little higgs theory at large mass, the goldstone boson equiva-. While this theorem does not hold at the lower masses of interest, a handle on the various decay modes. Thus we consider the case W + W t), (T ! ZZt) and (T ! hht) are 2:1:1 respectively 1 . this ratio due to phase space factors. that all have di↵erent reaches. As the dominant decay involves we expect that the most important searches are those involving ein, we have included the bounds from the following searches, ??. . ??. e the tightest bounds are ttom up study of naturalness at the LHC. As it is not clear how is canceled, all low energy e↵ective models which cancel quadratic articular class of models studied in this article are those where the ling quadratic divergences is also the discovery channel for the new cause discovery of the new particles immediately indicate that they tic divergences.
lo Generation 5. Showered in Pythia. Details about lepton isolation, b tagging we get limits.
sponsible for setting a smaller v for the top sector, we expect that the dominant ple of a little higgs theory at large mass, the goldstone boson equ tion. While this theorem does not hold at the lower masses of inter ain a handle on the various decay modes. Thus we consider the c ! W + W t), (T ! ZZt) and (T ! hht) are 2:1:1 respective m this ratio due to phase space factors. ches that all have di↵erent reaches. As the dominant decay invo s, we expect that the most important searches are those involv is vein, we have included the bounds from the following searches, Figs. ??. lace the tightest bounds are bottom up study of naturalness at the LHC. As it is not clear h ggs is canceled, all low energy e↵ective models which cancel quadr e particular class of models studied in this article are those where nceling quadratic divergences is also the discovery channel for the because discovery of the new particles immediately indicate that t dratic divergences.
arlo Generation ph 5. Showered in Pythia. Details about lepton isolation, b tagg id we get limits.
re responsible for setting a smaller v for the top sector, we expect that the domi udy, it is not clear which of the three di↵erent three body tle higgs theory at large mass, the goldstone boson equivae this theorem does not hold at the lower masses of interest, le on the various decay modes. Thus we consider the case t), (T ! ZZt) and (T ! hht) are 2:1:1 respectively 1 . io due to phase space factors. ll have di↵erent reaches. As the dominant decay involves ect that the most important searches are those involving have included the bounds from the following searches, ??.
ghtest bounds are p study of naturalness at the LHC. As it is not clear how eled, all low energy e↵ective models which cancel quadratic r class of models studied in this article are those where the dratic divergences is also the discovery channel for the new iscovery of the new particles immediately indicate that they ergences. neration wered in Pythia. Details about lepton isolation, b tagging limits.
e for setting a smaller v for the top sector, we expect that the dominant 7 HC searches bound the simple case of top quark, light quark iple higgs search would be strong evidence for a new particle about Monte Carlo generation and exclusions are in App. A. udy, it is not clear which of the three di↵erent three body tle higgs theory at large mass, the goldstone boson equivae this theorem does not hold at the lower masses of interest, le on the various decay modes. Thus we consider the case t), (T ! ZZt) and (T ! hht) are 2:1:1 respectively 1 . io due to phase space factors. ll have di↵erent reaches. As the dominant decay involves ect that the most important searches are those involving have included the bounds from the following searches, ??.
ghtest bounds are p study of naturalness at the LHC. As it is not clear how eled, all low energy e↵ective models which cancel quadratic r class of models studied in this article are those where the dratic divergences is also the discovery channel for the new iscovery of the new particles immediately indicate that they ergences. neration wered in Pythia. Details about lepton isolation, b tagging limits.
e for setting a smaller v for the top sector, we expect that the dominant the limit of large m T where the goldstone boson equivalenc dition to the three body decays, there is the two body deca v/m T relative to the three body decay. The two three body decays that will be compared are T ! th agrams are shown in Fig. ??
De ta i l s o f M o n te Ca rl o Ge n e ra ti o n vents were generated in Madgraph 5. Showered in Pythia. d our detector simulator. How did we get limits.
ners ion we examine how the current LHC searches bound the simple case of top q partner. Seeing a signal in a multiple higgs search would be strong evidence els quadratic divergences. Details about Monte Carlo generation and exclusio absence of an explicit model to study, it is not clear which of the three di↵ important. In the example of a little higgs theory at large mass, the goldst em is a good approximation. While this theorem does not hold at the lower m ume that it does to obtain a handle on the various decay modes. Thus we matrix elements for (T ! W + W t), (T ! ZZt) and (T ! hht) are 2 width will then vary from this ratio due to phase space factors. are many di↵erent searches that all have di↵erent reaches. As the domina auge bosons and higgses, we expect that the most important searches ar ptons and b jets. In this vein, we have included the bounds from the follow obtained are shown in Figs. ??. tly the searches which place the tightest bounds are clusion icle, we have initiated a bottom up study of naturalness at the LHC. As it tic divergence to the Higgs is canceled, all low energy e↵ective models which s should be studied. The particular class of models studied in this article ar pling responsible for canceling quadratic divergences is also the discovery ch hese models are unique because discovery of the new particles immediately hing to do with the quadratic divergences.
ails of Monte Carlo Generation e generated in Madgraph 5. Showered in Pythia. Details about lepton iso tector simulator. How did we get limits.
se where extra dimensions are responsible for setting a smaller v for the top sector, we expe important. In the example of a little higgs theory at large mass, m is a good approximation. While this theorem does not hold at me that it does to obtain a handle on the various decay modes atrix elements for (T ! W + W t), (T ! ZZt) and (T ! width will then vary from this ratio due to phase space factors. re many di↵erent searches that all have di↵erent reaches. As uge bosons and higgses, we expect that the most important tons and b jets. In this vein, we have included the bounds from obtained are shown in Figs. ??. ly the searches which place the tightest bounds are clusion cle, we have initiated a bottom up study of naturalness at the ic divergence to the Higgs is canceled, all low energy e↵ective m should be studied. The particular class of models studied in th pling responsible for canceling quadratic divergences is also the d hese models are unique because discovery of the new particles im hing to do with the quadratic divergences.
ails of Monte Carlo Generation e generated in Madgraph 5. Showered in Pythia. Details abou ector simulator. How did we get limits.
ing a signal in a multiple higgs search would be strong evidence for a new p divergences. Details about Monte Carlo generation and exclusions are in A n explicit model to study, it is not clear which of the three di↵erent three In the example of a little higgs theory at large mass, the goldstone boson e approximation. While this theorem does not hold at the lower masses of in does to obtain a handle on the various decay modes. Thus we consider th ents for (T ! W + W t), (T ! ZZt) and (T ! hht) are 2:1:1 respect hen vary from this ratio due to phase space factors. erent searches that all have di↵erent reaches. As the dominant decay in and higgses, we expect that the most important searches are those inv jets. In this vein, we have included the bounds from the following searche shown in Figs. ??. hes which place the tightest bounds are initiated a bottom up study of naturalness at the LHC. As it is not clea e to the Higgs is canceled, all low energy e↵ective models which cancel qua tudied. The particular class of models studied in this article are those whe sible for canceling quadratic divergences is also the discovery channel for th are unique because discovery of the new particles immediately indicate tha ith the quadratic divergences.
onte Carlo Generation in Madgraph 5. Showered in Pythia. Details about lepton isolation, b ta ator. How did we get limits.
dimensions are responsible for setting a smaller v for the top sector, we expect that the do ) In the absence of an explicit model to study, i decays are important. In the example of a little hi lence theorem is a good approximation. While this we will assume that it does to obtain a handle on where the matrix elements for (T ! W + W t), The decay width will then vary from this ratio du There are many di↵erent searches that all hav multiple gauge bosons and higgses, we expect th multiple leptons and b jets. In this vein, we have The limits obtained are shown in Figs. ??.
Currently the searches which place the tightest

Conclusion
In this article, we have initiated a bottom up stu the quadratic divergence to the Higgs is canceled, divergences should be studied. The particular clas quartic coupling responsible for canceling quadrati particles. These models are unique because discove have something to do with the quadratic divergenc

A Details of Monte Carlo Gener
Events were generated in Madgraph 5. Showered and our detector simulator. How did we get limits partners In this section we examine how the current LHC search and lepton partner. Seeing a signal in a multiple higgs which cancels quadratic divergences. Details about Mo In the absence of an explicit model to study, it is decays are important. In the example of a little higgs lence theorem is a good approximation. While this the we will assume that it does to obtain a handle on the where the matrix elements for (T ! W + W t), (T The decay width will then vary from this ratio due to There are many di↵erent searches that all have d multiple gauge bosons and higgses, we expect that multiple leptons and b jets. In this vein, we have incl The limits obtained are shown in Figs. ??.
Currently the searches which place the tightest bou

Conclusion
In this article, we have initiated a bottom up study o the quadratic divergence to the Higgs is canceled, all l divergences should be studied. The particular class of quartic coupling responsible for canceling quadratic div particles. These models are unique because discovery o have something to do with the quadratic divergences.

A Details of Monte Carlo Generati
Events were generated in Madgraph 5. Showered in P and our detector simulator. How did we get limits. 1 In the case where extra dimensions are responsible for setting decay is instead (T ! hht) 7 in the limit of large m T where the goldstone boson equiva addition to the three body decays, there is the two body of v/m T relative to the three body decay.
The two three body decays that will be compared are T diagrams are shown in Fig. ??

A Goldstone Boson Equivalence Theorem
In this section, we clarify how the goldstone boson equivalence theorem work for quartic interactions. In particular, we will study the system in the limit of large m T where the goldstone boson equivalence theorem holds. In addition to the three body decays, there is the two body decay T → th which is suppressed by a factor of v/m T relative to the three body decay.
The two three body decays that will be considered are T → thh and T → tW + W − . The two Feynmann diagrams are shown in Fig. 8. Schematically the matrix elements are . The last term results from the polarization vectors of the massive gauge bosons which scale in the large energy limit as k µ m W . The usual expectation is that putting the Higgs on-shell would be the dominant contribution to Eq. A.3, i.e. the two body decay would be dominant. In the large m T limit, this intuition is not valid because the large momentum flowing in the propagator is cancelled by the large k µ in the polarization vectors of the gauge bosons. Schematically, scalar and fermion final states contribute k 0, 1 2 µ respectively so do not overcome the off-shell propagator suppression.
Thus in the large m T limit, decays proceeding through an off-shell intermediate Higgs dominate. This calculation shows that even though the two body decay of the Higgs is suppressed, its existence is still important for maintaining the goldstone boson equivalence theorem.

B Little Higgs, Quadratic Divergences, and Three Body Decays
In Sec 4, we considered bilinear couplings to the Higgs where the Yukawa term was suppressed. In this section, we show how this limit can be reached in a toy Little Higgs model. Consider a toy Little Higgs models describing the spontaneous breaking of a SU (3) down to an SU (2) by a fundamental vev f [14]. The cut off scale at which new physics must enter is 4πf . The breaking gives 5 pseudo-goldstone bosons, 4 of which can be made into the Higgs field and the other will be ignored as it plays no role in subsequent analysis. The goldstones are combined into a non-linear sigma field Σ.
A pair of colored weyl fermions, u and u c , are added to cancel the top quadratic divergence by collective symmetry. One of the weyl fermions combines with Q 3 to form a SU (3) triplet χ = (Q 3 , u ). The lagrangian is then The leading 1/f terms in the lagrangian are then , Γ(T → Zt) and Γ(T → ht) respectively. The dashed blue/red/black lines show Γ(T → W + W − t), Γ(T → ZZt) and Γ(T → hht) respectively. The goldstone equivalence theorem was used to calculate the three body decays. In the case of a top partner, it is hard for the 3 body decay to be parametrically larger than the 2 body decay, but for a bottom partner it happens quite readily.
To further illustrate this model, consider top and bottom partners with various λ 2 . As before, we work in the limit where v f . We use the goldstone boson equivalence theorem to obtain the three body decay rates so that the ratio of their matrix elements squared are thus 2 : 1 : 1. Since the three body phase space is especially sensitive to the mass of the decay products, it is only at high masses that the decay widths approach this ratio as well.
In Fig. 9, we see that one can interpolate between making the 3 body decay dominant and the 2 body decay dominant by varying λ 2 . For a top quark partner, it is difficult to make the 3 body decay parametrically dominant while for a bottom quark partner, the three body decay is expected to be the dominant decay mode.

B.2 v ∼ f
In the previous subsection, the critical assumption of v f was made. The assumption is invalidated at low masses because if m T ∼ 500 GeV and λ 2 ∼ 10 then f ∼ 50 GeV. The reason why v = 0 leads to such strong constraints is that HH † u c 3 u always leads to a two body decay when H obtains a vev. One can estimate what the maximal relative branching ratios are as a function of m T . Assuming massless decay products, one obtains the relationship This ratio shows that for large λ 2 , the three body decay to two Higgses dominates over the two body decay around 8 TeV. The exact Little Higgs computation confirms this expectation. Fig. 10 shows the differential cross section as a function of m T . Except for rather large masses, it is impossible to make a three body decay dominate over all of the two body decays. In the large λ 2 limit, the decay to bW and tZ are suppressed so that the only competing two body decay channel is T → th. This suppression is an indication that what is causing that particular two body decay channel is in fact the quartic coupling rather than the Yukawa interaction. This effect is simply what was discussed in Sec. 4.

C Details of Monte Carlo Generation
Events were generated using MadGraph 4.5.1 [59] and showered with Pythia 6.4 [61]. After showering, all hadron-level events are passed to the PGS 4 detector simulation, which parameterizes the detector response. The detector parameters used are those of the default ATLAS PGS card except for jets, which are clustered using the anti-k t algorithm with a radius of 0.4. Heavy flavor jets are tagged using the p T dependent efficiencies found in [62]. In order to validate the ATLAS three and four lepton searches, light leptons identified by PGS are tagged using the p T and η dependent tight electron and muon identification efficiencies of respectively [63] and [64]. Hadronic taus are treated as jets.
For most SUSY searches, ATLAS gives the cut flow associated to a few benchmark models, which can then be used to validate our results. The results of this validation process for the ATLAS multilepton searches studied in Sec 4.1 are shown in Table 5     . The benchmark model tested here is χ ± 1 χ 0 2 → W ± Zχ 0 1 χ 0 1 . For each step of the event selection, the number of events obtained using our simulated samples and our code is given on the left and the number of events given by ATLAS is shown on the right. 20