Nonstandard Yukawa Couplings and Higgs Portal Dark Matter

We study the implications of non-standard Higgs Yukawa couplings to light quarks on Higgs-portal dark matter phenomenology. Saturating the present experimental bounds on up-quark, down-quark, or strange-quark Yukawa couplings, the predicted direct dark matter detection scattering rate can increase by up to four orders of magnitude. The effect on the dark matter annihilation cross section, on the other hand, is subleading unless the dark matter is very light -- a scenario that is already excluded by measurements of the Higgs invisible decay width. We investigate the expected size of corrections in multi-Higgs-doublet models with natural flavor conservation, the type-II two-Higgs-doublet model, the Giudice-Lebedev model of light quark masses, minimal flavor violation new physics models, Randall-Sundrum, and composite Higgs models. We find that an enhancement in the dark matter scattering rate of an order of magnitude is possible. Finally, we point out that a discovery of Higgs-portal dark matter could lead to interesting bounds on the light-quark Yukawa couplings.

The couplings to gauge bosons are found to be in agreement with the SM predictions at the O(20%) level, while the constraints on the couplings to third-generation fermions are somewhat weaker [11,12]. Much less is known experimentally about the couplings of the Higgs to the first two generations of fermions. The couplings to u, d, s, and c quarks could be as large as the SM bottom Yukawa coupling or be absent altogether [13][14][15][16]. The Higgs couplings to top and bottom quarks will be quite well known by the end of the high-luminosity LHC run. Some progress is also expected on the measurements of Higgs couplings to charm and strange quarks [13][14][15].
Large u-, d-, and s-quark Yukawa couplings, comparable in size to the b-quark Yukawa, generically require fine-tuning. A large Yukawa coupling implies a large contribution to the quark mass from the Higgs vacuum expectation value (vev). This would then need to be cancelled by a different contribution to the u-, d-, and s-quark masses, unrelated to the Higgs vev. The opposite limit, where the observed Higgs does not couple to the light quarks at all is easier to entertain. It simply requires a separate source of the light-quark masses (for an extreme example see, e.g., [17]).
Modified light-quark Yukawa couplings could, in principle, have important implications for DM phenomenology. In this article we investigate how the Higgs-portal DM predictions change if the Higgs couplings to the light quarks differ from the SM expectations. We first allow for an arbitrary flavor structure of the Higgs Yukawa couplings, only requiring that they satisfy the current experimental bounds. In Section II we derive the implications for direct DM detection, indirect DM detection and the collider searches. We show that vanishing couplings of the Higgs to light quarks only have a negligible impact on these observables.
Saturating the loose current bounds on the light-quark Yukawa couplings would, on the other hand, lead to drastically enhanced scattering cross sections on nuclei while leaving the relic density and annihilation cross sections nearly unmodified.
Clearly, an enhancement of the light-quark Yukawas by factors of O(100) or more, as allowed by current data, requires considerable fine tuning of the quark-mass terms and hence seems quite unlikely. In Section III we, therefore, explore the deviations in the Higgs Yukawa couplings for a number of beyond-the-SM scenarios and flavor models. This leads to more realistic expectations as to how large the deviations in the direct DM detection rates can be due to the poorly known Higgs couplings to the light quarks. Note that we assume the DM to be a flavor singlet and that the new flavor structure of the interactions with the visible sector is only due to the modification of the SM Higgs couplings. DM that is in a nontrivial flavor multiplet has been investigated in [18][19][20][21][22][23][24][25][26][27][28][29], while our study is closer in spirit to the work in [30][31][32] where the flavor dependence of the DM signals for flavor-singlet DM has been explored.
A somewhat surprising result of our investigation is that, if DM is discovered and turns out to be a thermal relic predominantly interacting through a Higgs portal, it could be used to constrain the light-quark Yukawa couplings. This is discussed in more detail in Section IV.
We summarize our results in Section V.

II. HIGGS PORTAL WITH NON-TRIVIAL FLAVOR STRUCTURE
We assume that DM and the SM fields are the only light degrees of freedom. The remaining new physics (NP) particles can be integrated out so that one can use an Effective Field Theory (EFT) approach. The couplings of DM to the SM are given by the Higgs-portal Above, the fermion DM can be either a Dirac or Majorana fermion (in either case we use four component notation). After electroweak symmetry breaking (EWSB) we have where v W = 246 GeV is the vacuum expectation value (vev) of the Higgs field. The above interactions therefore lead to annihilation of DM into both single Higgs, χχ → h, and double Higgs, χχ → hh, final states.
The Higgs-portal operator for fermionic DM has mass dimension five and is suppressed by the new physics scale Λ. The Higgs-portal interaction for fermionic DM can also be re-written as L χ = (g χ + ig χ )χ L χ R H † H/Λ + h.c.. Forg χ = 0 the interaction is thus both P -and CP -violating. The interaction for vector DM is most probably also due to a higherdimensional operator in the full theory. For instance, if χ µ arises from a spontaneously broken gauge symmetry in the dark sector, then g χ ∼ v 2 D /Λ 2 , where v D is the vev of the field that breaks the dark sector gauge invariance, while Λ is the mass of the mediator between DM and the Higgs.
The relevant terms, after EWSB, in the effective Lagrangian for the Higgs couplings to the SM particles are given by where the κ i are real. A sum over the SM quarks, q = u, d, s, c, b, t, and charged leptons, = e, µ, τ is implied, and we have assumed custodial symmetry. The h → γγ coupling is not relevant for DM phenomenology, since its effects are suppressed compared to the Higgs couplings to gluons. 1 The couplings are normalized such that κ q = κ = κ V = κ λ = 1 correspond to the SM. The experimental constraints on the couplings of the light quarks to the Higgs, obtained from a global fit to current data, are  [13]. Higgs couplings to the light quarks of a size comparable to the coupling to the b quark are thus still allowed. In (3) we do not allow for flavor violating Higgs couplings, since these are already tightly constrained from both Higgs decays and low-energy observables [33][34][35][36].
In the SM, keeping the Higgs on shell, the hGG coupling arises predominantly from the one-loop top-quark contribution. The κ NP g in (3) encodes only the potential NP contributions, and vanishes in the SM. In the global fits a parameter κ g is introduced that gives the total h → gg amplitude, including the SM contributions [11,12]. We have (see, e.g., [37]) Figure 1: DM annihilation channels in the Higgs-portal models.
At present, significant CP -violating Higgs couplings to fermions and gluons are still allowed experimentally (see, e.g., [38]), so we also discuss their effect on the Higgs interactions with DM: Here, theκ i are real parameters; in the SM, we haveκ i = 0. Moreover, G a µν is the gluon field-strength tensor and G a,µν = 1 2 µναβ G a αβ its dual. The normalization of the hGG term is chosen such that integrating out the top at one loop one obtainsκ g =κ t +κ N P g . Accordingly, we have Br(h → gg) ∝ κ 2 g + (3κ g /2) 2 . The CP -violating couplings of the Higgs to ZZ and W W are already well constrained, and we thus set them to zero. In the numerical analysis below, we will also assumeκ i = 0, for simplicity.
The modified Higgs couplings change the usual Higgs-portal predictions for DM annihilation rates, the relic abundance, and direct detection rates. In the following, we discuss these modifications in detail.

A. Annihilation cross sections
The dominant DM annihilation cross sections in the Higgs-portal models are χχ → bb, W + W − , ZZ, tt, and χχ → hh. The first four proceed through the s-channel Higgs exchange, while χχ → hh receives additional contributions from t-and u-channel χ exchange as well as from the four-point contact interaction (cf. Fig. 1). The χχ → bb channel is only relevant if the other channels are not kinematically allowed, i.e. for light DM masses, The χχ →bb annihilation cross section assuming SM Higgs couplings is given for scalar (S), Dirac fermion (DF ), and vector (V ) DM by where, here and below, √ s is the center-of-mass energy, r k = 4m 2 k /s, β k = √ 1 − r k is the velocity of particle k, and v rel = 2β χ is the relative velocity of the DM particles. If the Higgs coupling to the b-quarks differs from the SM value, the annihilation cross section is rescaled The annihilation cross sections σ ff to the other fermions are obtained with the obvious replacement b → f in the above expressions. Since the Higgs couplings to the light quarks are poorly constrained experimentally, the DM annihilation to two light quarks can be comparable to χχ → bb and can be important for light DM, m χ < m W .
For heavy DM, m χ > m W , the annihilation cross-sections into a pair of W or Z bosons are V = W, Z. The annihilation cross sections assuming the SM Higgs couplings to W are given by for scalar, Dirac fermion, and vector DM, respectively. The χχ → ZZ annihilation cross section is obtained by replacing W → Z, and multiplying all expressions by an extra factor of 1/2 since one has two indistinguishable particles in the final state.
The χχ → hh annihilation cross sections for scalar, Dirac fermion, and vector DM, are given by In this result we display only the leading terms in the expansion in powers of β χ . The contribution of the s-channel Higgs exchange diagram is proportional to the rescaling of the trilinear Higgs coupling, κ λ . The latter is completely unknown experimentally, at present, but can be measured to O(20%) at the end of the LHC [39,40]. Since we are mostly interested in the effects of flavor modifications we will set it to the SM value, κ λ = 1, in the numerics below. All cross sections for Majorana DM can be obtained by multiplying the corresponding Dirac DM cross sections by a factor of 4.

B. Relic abundance
The DM relic abundance Ω DM is proportional to 1/σv rel , where σ is the annihilation cross section. Assuming that the DM in our scenario accounts for all of the observed relic density, the measured value Ω DM h 2 = 0.1198(26) [41] fixes g χ for a given value of m χ . The resulting constraint in the m χ -g χ plane is denoted for the different cases by a red line in Figs. 2 to 9.
In Fig. 2 Fig. 2 denote the value of g χ for which the total Higgs decay width would be larger than   cross section to b-jets, σ(χχ → bb). For m χ m W these are the two dominant annihilation modes. Since the annihilation cross sections to b-jets and light jets are comparable, the relic abundance curve show only a small change in g χ when the b-quark threshold is reached. This should be compared with the case of the SM Yukawa couplings shown in the right panels of Figs. 2, 4, 6, and 8. In this case the χχ → jj annihilation is almost exclusively due to DM annihilating to two gluons, so that σ(χχ → jj) σ(χχ → bb), while the annihilation into two light quarks is negligible. For m χ m W and SM Yukawas, the required g χ is thus bigger by 30% − 40% than in the case of light Yukawas at their present experimental limits, and exhibits a significant jump below the b-quark threshold.

C. Invisible decay width of the Higgs
The bounds on the invisible decay width of the Higgs boson provide stringent constraints on Higgs-portal DM [7]. The partial h → χχ decay widths are given by where r χ = 4m 2 χ /M 2 h and β χ = 1 − r χ . The current best limits on the invisible branching fraction of the SM Higgs are obtained from Zh production. The CMS collaboration gives a 95% CL limit of Br(h → inv) < 0.58 for M h = 125 GeV [42] and ATLAS finds Br(h → inv) < 0.75 for M h = 125.5 GeV [43]. Note that the increased light-quark Yukawa couplings, at their presently allowed values, do not appreciably change the Higgs production cross section [13]. Their main effect is to increase the total decay width of the Higgs and thus reduce the branching ratios to the other decay modes: In Figs. 2 to 9 we denote the bound on g χ corresponding to the ATLAS upper limit on Br(h → inv) with a dashed black line and grey out the excluded region in the g χ vs. m χ plane. We see that the light DM Higgs portal, m χ m h /2 is excluded by the Higgs invisible decay width.
Vector boson fusion, gluon fusion and ttH production, with the off-shell Higgs going to two DM particles, can provide some limited sensitivity to DM masses above m h /2. A combination of the searches in the three channels at a 100 TeV collider could exclude the scalar thermal relic DM Higgs portal for DM masses in parts of the m h /2 m χ m W interval at 95% C.L. [10] (these result receives only a negligible correction if light quark Yukawa couplings are enhanced). For m χ < m h /2 the invisible Higgs decay width is, however, always the most constraining [44].

D. Indirect detection
In indirect signals of DM annihilation, the effect of changing the light-quark Yukawa couplings within the presently experimentally allowed ranges leads to at most O(1) effects.
Further, the effect is present only for DM masses below the W threshold where the dominant annihilation channel is into the bb final state. For example, Fig. 10 shows the recast of the Fermi-LAT bound from dwarf spheroidals [45] for scalar DM, following the procedure outlined below.
Photon flux measurements with γ-ray telescopes can put a strong bound on the annihilation cross-section of the DM. The strongest bound in the DM mass range of interest has been recently released by the Fermi-LAT collaboration [45] based on Pass 8 observation data of the Milky Way dwarf spheroidal satellite galaxies (dSphs). There is also a recent analysis based on the Dark Energy Survey (DES) dSph candidates using the Fermi-LAT data [46].
While this bound is competitive with the one from the known dSphs, it is still weaker on its own. One could also consider the bounds from the isotropic gamma ray background (IGRB) [47]. In our analysis, we recast the Fermi-LAT bound on the bb final state using a simple re-weighting procedure of the photon spectra which will be discussed below.
The observed differential photon flux from the annihilation of dark matter is given by where J is an astrophysical factor which depends on the distance to the source and the dark matter density profile. The factor in the brackets is the one most interesting for our purposes. It depends on the velocity-averaged cross-section and photon spectrum per DM annihilation.
The Fermi-LAT analysis gives bounds for the different final states separately while we have an admixture of final states. In order to recast the bound, we rely on the observation that for heavy DM the photon spectra from DM annihilation into quarks, gauge bosons, and the Higgs boson all peak at approximately the same photon energy and have approximately the same shape. Therefore, to extract the bound on the DM Higgs portal coupling g χ , it is sufficient to find the zeros of the polynomial where in the last term, σv fermi b is the bound from the Fermi-LAT analysis on the velocityaveraged cross-section. The photon spectra were obtained from the interpolation tables provided in the PPPC4DMID package [48]. In all cases except for the hh final state, f (g χ ) has only one zero up to a sign ambiguity. For the hh final state, however, the zero of f (g χ ) closest to the g χ corresponding to χχ → bb is the one used to rescale the Fermi-LAT bound on g χ as a function of m χ . The resulting bounds are shown in Figs. 2 to 9.

E. Direct detection
We have seen so far that most DM observables exhibit only a weak dependence on the light-quark Yukawas. This is not the case for the direct DM detection. In fact, modifying the light-quark Yukawa couplings can significantly change the predictions for DM -nucleus scattering cross sections.
The differential cross section for spin-independent DM scattering on a nucleus is given where E R is the nclear recoil energy, m A is the mass of the nucleus, is the reduced mass of the DM -nucleus system, s = (m χ +m A ) 2 is the center-of-mass energy, v rel is the DM velocity in the detector rest frame, and |M| 2 is the spin-averaged squared matrix element.
The matrix element M depends on the effective Higgs couplings to the nucleus. Since the momentum exchanges in DM scattering on nuclei are much smaller than the Higgs mass, we can calculate |M| 2 by first integrating out the Higgs and the heavy quarks (t, b, c). This gives an EFT with light quarks and gluons interacting with DM through local operators, described by the effective Lagrangians for scalar, fermion, and vector DM, respectively. The scalar current is the same in all three cases: Here, the last two terms arise from CP -violating Higgs couplings. The sums are over the light quarks q = u, d, s. The heavy quarks are integrated out and contribute only via the gluonic terms in the current. For the two corresponding dimensionless Wilson coefficients we have where the first contribution is from tree-level matching, and the remaining from one-loop matching, working in the limit of heavy quarks. This is well justified for top and bottom quarks. For scattering on heavy nuclei, e.g., on Xe or W, the maximal momentum exchanges for DM with mass above approximately 1 TeV may, however, start to become comparable to the charm-quark mass. We neglect these effects, while they may need to be included in the future if such heavy DM is discovered.
CP -violating Higgs couplings to light quarks lead to spin-dependent interactions of DM with the target nuclei. The corresponding scattering rates are suppressed relative to the spin-independent interaction rates from CP -conserving Higgs couplings. We will therefore neglect the CP -violating interactions in our numerical analysis of direct detection scattering rates; i.e., we will setκ q = 0,C q = 0 from now on.
The nucleon matrix elements of the remaining terms in the scalar current S q are conventionally parametrized by (see, e.g., [49]), In the heavy-quark limit for t, b, c the trace anomaly equation leads to the relation [49,50] f (N ) We can also define the effective Higgs coupling to nucleon as the expectation value of the scalar current, The exclusion curves from LUX [51] and Xenon100 [52], assuming a local DM density T s = 0.043 ± 0.011 [53]. The matrix elements for u and d quarks can be related to the σ πN term. A Baryon Chiral Perturbation Theory (BχPT) analysis of the πN scattering data gives σ πN = 59(7) MeV [54]. This is in agreement with BχPT fit to world lattice N f = 2 + 1 QCD data, which gives σ πN = 52(3)(8) MeV [55]. Including both ∆(1232) and finite spacing parametrization in the fit shifts the central value to σ πN = 44MeV. To be conservative we use σ πN = (50 ± 15)MeV, which gives f T d = (3.8 ± 1.1) × 10 −2 , using the expressions in [56]. For the effective Higgs coupling to nucleons this gives We use the results in [57] to relate the nuclear matrix elements to actual scattering rates on nuclei via nuclear form factors.
We  Fig. 9). This means that if thermal relic DM is discovered, it would immediately place an upper bound on κ u , κ d , κ s , assuming Higgs-portal mediation (unless in the case of fermion DM that has purely pseudoscalar couplings). We comment in more detail on that observation in Section IV.
Since the DM -nucleus scattering cross section is the only DM observable that exhibits a rather pronounced dependence on the values of the light-quark Yukawas, we study this dependence in more detail. In Fig. 11 (left) we show how the direct detection bounds on g χ are affected by changes in the values of the light-quark Yukawas. We plot the ratio where g max χ (κ q ) is the upper bound on g χ obtained from direct detection experiments for a given value of κ q , with q = u, d, s. Hence, g max χ (1) is the bound obtained assuming SM Yukawa couplings. Its value depends on m χ , on whether DM is a scalar, fermion, or vector, and on the experiment that measured the bounds. Similarly, also g max χ (κ q ) depends on m χ , the spin of DM, and the experiment; however, all these dependences cancel in the ratio ξ gχ .
The ratio ξ gχ thus only depends on κ q and on which target material was used to derive the direct detection bounds. In Fig. 11 we show ξ gχ for a Xenon target, varying in turn κ u (dark red line), κ d (light red) and κ s (blue), while keeping all other parameters fixed to their SM values. We set the hadronic matrix element f (N ) T q to their present central values, anticipating that in the future their uncertainties will be further reduced. In Fig. 11 (right) we show a closely related quantity -the ratio of the scattering cross sections with varied κ u,d,s and the scattering cross section with SM Yukawa couplings, σ d.d. /σ SM d.d. . Fig. 11 illustrates clearly that the difference between the bounds where the light quark Yukawa couplings are small or vanish completely, and the bounds where all the couplings are SM-like, is very small, O(5%). Saturating the present experimental bounds on κ u or κ d , the allowed value of g χ could lie two orders of magnitude below what one obtains for the case of SM Yukawa couplings. Such large values for the light-quark Yukawas are not very likely to be realized in a concrete model, as we will discuss in the next section. However, it is very interesting to observe that even a moderate increase of the values of the light-quark Yukawa couplings to only a few times their SM values can have a significant effect on the direct detection bounds, enhancing the scattering cross sections by up to a factor of ten.  [58,59], in the MSSM at tree level, the Giudice-Lebedev model of quark masses (GL) [60], in NP models with minimal flavor violation (MFV) [61], in Randall-Sundrum models (RS) [62], and in models with a composite Higgs, realized as a pseudo-Nambu-Goldstone boson (pNGB) [63][64][65][66]. For completeness, we include both the flavor-conserving and flavor-violating Yukawa interactions, and allow for CP violation. The Higgs couplings to quarks are thus described by where a sum over the SM quark fields is understood. The first two terms are flavor diagonal, with the first term CP conserving and the second term CP violating, and coincide with the definitions in eqs. (3) and (5), respectively. The terms in square brackets are flavor violating, with the real (imaginary) part of the coefficient CP conserving (violating). In the SM we have κ q = 1, whileκ q = κ qq =κ qq = 0. The flavor-violating couplings in the above set of NP models are collected in Tables III and IV. These tables complement the analyses in [67][68][69] (see also [70], where implications of a negative top-quark Yukawa were explored).

A. Dimension-Six Operators with Minimal Flavor Violation
We start our discussion by considering dimension-six operators arising from integrating out NP at a high scale Λ. In addition, we assume that the flavor breaking in the NP sector is only due to the SM Yukawas, i.e. that NP satisfies the Minimal Flavor Violation (MFV) hypothesis [61,[71][72][73][74][75][76]. Integrating out the new physics states gives for the Higgs couplings Model κ ij /κ t , are the same as the CP -conserving ones, apart from substituting "Im" for "Re" in the "MFV" row.
to quarks where Λ is the scale of new physics and H c = iσ 2 H * . We identify the NP scales in the upand down-quark sectors for simplicity. There are also modifications of quark kinetic terms through dimension-six derivative operators. These can be absorbed in (35)

using equations
Model are the same as the CP -conserving ones, apart from substituting "Im" for "Re" in the "MFV" row.
of motion [77]. The quark mass matrices and Yukawa couplings after EWSB are thus Because Y q and Y q appear in two different combinations in M q and y q , the two, in general, cannot be made diagonal in the same basis.
In MFV the coefficients of the dimension-six operators can be expanded in terms of Y u,d , with a q , b q , c q ∼ O(1). Working to first order in dimension-six operator insertions we can thus write for the Yukawa couplings, in the mass eigenbases for up and down quarks respectively, Here y u,d SM are the diagonal matrices of the SM Yukawa couplings, while V is the Cabibbo-Kobayashi-Maskawa (CKM) matrix. In general, the coefficients a q , b q , c q are complex so that CP -violating Higgs couplings arise at O(v 2 W /Λ 2 ). Flavor-violating Higgs couplings arise first from the terms proportional to c u,d in the Yukawa expansion and are thus suppressed by the corresponding CKM matrix elements. In Tables I-IV we  Under the assumptions above, the mass and Yukawa terms can be diagonalized in the same basis, so that there is no flavor violation and no CP violation in the Yukawa interactions: We obtain a universal shift in all up-quark Yukawa couplings, and a different universal shift in all down-quark Yukawa couplings, given by Since the shifts are universal over generations, the precise measurements of the Higgs couplings to top and bottom quarks will also determine the Higgs couplings to light quarks.
Both κ t and κ b are expected to be known with O(5%) precision after the end of the highluminosity LHC run [78,79]. The uncertainties in the DM direct detection rates due to uncertainties in the Yukawa couplings will thus be negligible, assuming NFC. Note that the Higgs portal with an additional SM singlet mixing with the Higgs is also described by the above modifications of fermion couplings, with a completely universal shift κ i = cos θ, where θ is the singlet-Higgs mixing angle [80].
Our analysis of modified Higgs-portal DM phenomenology given in Section II applies in the somewhat special limit where the DM only couples to the lightest mass-eigenstate h.
For instance, for scalar DM the general Higgs portal is If the hermitian matrix of couplings g χ,ij is such that it has h as the only eigenstate with nonzero eigenvalue, then our analysis in Section II applies unchanged. In general, however, all the expressions in Section II get corrected by terms of order 1/m 2 H i due to exchanges of heavy Higgs bosons with masses m H i . If DM is heavy, m χ > m H i , the presence of heavy Higgs bosons would also open new annihilation channels.

C. Type-II Two-Higgs-Doublet Model
The MSSM tree-level Higgs potential and the couplings to quarks are the same as in the type-II two-Higgs-doublet model (2HDM), see, e.g., [81]. This is an example of a 2HDM with natural flavor conservation in which v u = sin β v W , v d = cos β v W . The mixing of h u,d into the Higgs mass-eigenstates h and H is given by where h is the observed SM-like Higgs. Thus while the flavor-violating and/or CP -violating Yukawas are zero. In the decoupling limit (β − α = π/2) the heavy Higgs bosons become infinitely heavy, while the Yukawa couplings tend toward their SM value, κ i = 1. The global fits to Higgs data in type-II 2HDM already constrain β − α π/2 [37,82,83] so that in this case the corrections to Higgs-portal DM phenomenology due to non-standard Higgs Yukawa couplings are small.
As in the case of MHDM, the DM phenomenology of Section II remains unchanged only in the limit where the DM couples to the light Higgs h but not to the heavy Higgs H. In the general case, our analysis gets corrections that are relatively suppressed by O(m 2 h /m 2 H ). If we are not too far away from the decoupling limit these corrections can be neglected, while in parts of the parameter space, where cancellation can occur, the extra contributions are numerically important [8].

D. Higgs-dependent Yukawa Couplings
In the model of quark masses introduced by Giudice and Ledebev (GL) [60] the Higgsquark interactions are written in terms of effective operators They can be thought of as arising from integrating out heavy mediators at a large mass scale M . In this model the light quarks couple to the Higgs only through operators with mass dimension higher than four, i.e., for light quarks we have n u,d ij = 0. The values of the integers n u,d ij , and of the coefficients c u,d ij that take values of order unity, are chosen such that the hierarchies of the observed quark masses and mixing angles are explained, after EWSB, in terms of the expansion parameter ≡ v 2 W /M 2 ≈ 1/60. Thus, the Yukawa couplings are of the form After mass diagonalization the SM Yukawas are diagonal in the same basis as the quark masses, (y u,d ij ) SM ∝ δ ij m u,d i , while the y u,d ij are not diagonal in the same basis 3 . Using the ansatz n u,d ij = a i + b u,d j with a = (1, 1, 0), b d = (2, 1, 1), and b u = (2, 0, 0), this gives the deviations in the Yukawa couplings collected in Tables I-IV in the row denoted by "GL".
Since the couplings to the bottom quark is enhanced by a large factor, κ b 3, the simplest version of the GL model is already excluded by the Higgs data on h → W W , h → ZZ and h → γγ decays.
We therefore modify the initial GL proposal and assume that we have two Higgs doublets in (44), H u that only gives masses to up-type quarks and H d that only gives masses to down-type quarks. The correct mass and CKM angle hierarchy is obtained by using b d = (1, 0, 0) in the ansatz for n d ij , and leaving a and b u unchanged. This gives satisfactory Higgs phenomenology at present as long as κ b = sin α/ cos β 1 up to O(20%). In this limit also κ t = cos α/ sin β 1.
The coupling of DM to gluons (25) C g 4κ t +κ b , so that C g ∼ (5/3)C SM g , andC g ∼ O(C g ). Taking κ b 1, this means that the effective Higgs coupling to nucleons, governing the direct DM detections rates, gets enhanced compared to the SM Higgs Yukawa couplings by Here most of the enhancement over the SM comes from enhanced κ u and κ d , which is also the reason for enlarged isospin breaking (the difference between f (p) S and f (n) S ). As a result of larger couplings to light quarks the spin-independent DM scattering cross section can thus be enhanced by an order of magnitude in the GL2 model of light-quark masses.

E. Randall-Sundrum models
The Randall-Sundrum (RS) warped extra-dimensional models with the SM fields propagating in the bulk provide a solution to the hierarchy problem and simultaneously explain the hierarchy of the SM fermion masses without large hierarchies in the initial five-dimensional (5D) Lagrangian [62,84,85]. The fermion zero modes are either localized toward the UV brane (for lighter fermions) or toward the IR brane (the top, the left-handed b quark and potentially the right-handed c quark) [86,87]. The Higgs field and the Higgs vev are both localized toward the IR brane. Integrating out the Kaluza-Klein (KK) modes and working in the limit of a brane-localized Higgs, the SM quark mass matrices are given, to leading order in v 2 W /m 2 KK , by [88] (see also [89][90][91][92][93][94][95][96][97]) Here, m KK is the KK mass scale. The F q,u,d are diagonal 3 × 3 matrices of fermion wavefunctions for the left-handed electroweak quark doublets and the right-handed electroweak up and down quark singlets, respectively, evaluated at the IR brane. Assuming flavor anarchy, the 5D Yukawa matrices for up and down quarks, Y 5D 1,2 , are general 3 × 3 complex matrices with O(1) entries. For a Higgs field propagating in the bulk, 5D gauge invariance [88].
At leading order in v 2 W /m 2 KK the Higgs Yukawas are aligned with the quark masses, i.e., The misalignment arises from dimension-six operators that are generated by tree-level KK quark exchanges. They give whereȲ is a typical value of the dimensionless 5D Yukawa coupling and is in numerical analyses typically taken to be belowȲ 4 (see, e.g., [92]). The Higgs mediated FCNCs are thus suppressed by the same zero-mode wave-function overlaps that also suppress the quark masses, giving rise to the RS GIM mechanism [98][99][100].
Using that the CKM matrix elements are given by V ij ∼ F q i /F q j for i < j, Eq. (49) can be rewritten as This yields the κ i collected in Tables I-IV. In the numerical analysis of ref. [88] the diagonal values κ i were typically found to be smaller than one, with deviations in κ t up to 30%, κ b up to 15%, in κ s,c up to ∼ 5%, and in κ u,d of 1% (these estimates were obtained fixing the mass of the first KK gluon excitation to 3.7 TeV, above the present ATLAS bound [101] Finally, we investigate the possibility that the Higgs is a pseudo-Goldstone boson arising from the spontaneous breaking of a global symmetry in a strongly coupled sector [63][64][65][66].
We assume that the SM fermions couple linearly to composite operators O L,R [102], where i, j are flavor indices. This is the 4D dual of the fermion mass generation in 5D RS models. The Higgs couples to the composite sector with a typical coupling y * . The SM masses and Yukawa couplings then arise from expanding the two-point functions of the O L,R operators in powers of the Higgs field [103], giving rise to four-and higher-dimensional Higgs operators, such as in (35).
The new ingredient, related to the pNGB nature of the Higgs, is that the shift symmetry dictates the form of the higher-dimensional operators. The flavor structure and the composite Higgs coset structure completely factorize if the SM fields couple to only one composite operator. The general decomposition of Higgs couplings then becomes [103] (see also [104,105]) and similarly for the down quarks. Here P (h/f ) = a 0 + a 2 (h/f ) 2 + . . . is an analytic function whose form is fixed by the structure of the spontaneous breaking and the embedding of the SM fields in the global symmetry of the strongly coupled sector, while f is the equivalent of the pion decay constant and is of order v W . Since the flavor structure of the coefficients of the dimension-four and dimension-six operators is the same, they can be diagonalized simultaneously. All corrections to the quark Yukawa couplings from this effect are therefore strictly diagonal, and we have For example, for the models based on the breaking of SO(5) to SO(4), the diagonal Yukawa couplings can be written as [106] where n, m are positive integers. The MCHM4 model corresponds to m = n = 0, while MCHM5 is given by m = 0, n = 1.
The flavor-violating contributions to the quark Yukawa couplings then arise only from corrections to the quark kinetic terms. That is, they are related to dimension-six operators of the form [103]q These operators arise from the exchange of composite vector resonances with typical mass M * ∼ Λ. After using the equations of motion they contribute to the misalignment between the fermion masses and the corresponding Yukawa couplings. The NDA estimates for these corrections are, neglecting relative O(1) contributions in the sum [15,103,107], and similarly for the down quarks. If the strong sector is CP violating, thenκ u,d ij ∼ κ u,d ij . The exchange of composite vector resonances contributes also to the flavor diagonal Yukawa couplings, shifting the estimate (53) by (note the different normalizations of κ q and κ qq in (34)) top, are expected to be composite, explaining the large top mass (i.e., λ q L,3 ∼ λ u R,3 ∼ 1). In the anarchic flavor scenario, one expects the remaining quarks to be mostly elementary (so the remaining λ i 1). However, if there is some underlying flavor alignment, it is also possible that the light quarks are composite. This is most easily achieved in the right-handed sector [105,108,109]. Taking all right-handed up-type quarks fully composite, and assuming that this results in a shift ∆κ u ∼ ∆κ c ∼ ∆κ t ∼ 1, this would lead to an increase in the effective Higgs coupling to nucleons, f S , of about 50%, and an increase in the DM-nucleon scattering rate of about 100%.

IV. CONSTRAINING THE LIGHT-QUARK YUKAWA COUPLINGS
If DM is a thermal relic interacting with ordinary matter predominantly via SM Higgs exchange, direct detection scattering rates immediately give information about the lightquark Yukawa couplings once the coupling of the DM particle to the Higgs particle is fixed.
In fact, DM scattering in direct detection searches would be one of the very few possible probes of the light-quark Yukawa couplings. The interactions of the Higgs boson with u, d, or s quarks give rise to flavor-conserving neutral currents. Off-shell Higgs contributions in processes with only SM external particles always compete with other, much larger flavorconserving neutral currents induced by gluon, photon, or Z exchange. This leaves us with two options: either to consider on-shell Higgs decays [13][14][15], or to use new probes, such as DM scattering in direct detection experiments.
In principle, there is enough information to make a closed argument. Suppose that indirect DM searches yield a positive DM annihilation signal for m χ > m h /2. At the end of the high-luminosity LHC run, the Higgs couplings to W , Z, t, and b will be precisely determined. Assuming that DM is a thermal relic interacting only through the Higgs portal, this fixes the value of g χ since the annihilation cross section for m χ > m h /2 is otherwise almost completely controlled by the Higgs couplings to W , Z, and t. In principle, a consistency check that the DM is really interacting through a Higgs portal could be provided, for a very limited range of DM masses m χ m h /2, by a 100-TeV hadron collider [10].
After the discovery of DM, the direct detection searches would immediately imply an upper bound on the light-quark Yukawa couplings. As an illustration, consider the excess in γ-ray emmission in the recently discovered dwarf spheroidal galaxy Reticulum 2 [46].
Let's take the bold step of interpreting this signal as originating from DM annihilating into bb pairs (see Ref. [110] for details). Assuming the Dirac-fermion DM scenario with purely CP-violating couplings, we obtain a 1σ region in the m χ -g χ plane that is not yet excluded by direct detection constraints, denoted by the orange lines in Fig. 12. Note that part of this region is consistent with DM furnishing the dominant component of the observed relic density while at the same time not being excluded by the invisible Higgs decay width.
Concentrating on the overlap region, m χ ∼ 75 GeV, a comparison with the ratios shown in   Figure 12: The γ-ray excess in the recently discovered dwarf spheroidal galaxy Reticulum 2, interpreted as a signal of DM annihilating into bb pairs, is shown as the black 1σ contour (see Ref. [110] for details). The orange lines show the 95% CL exclusion limits at the 14-TeV LHC (solid line) and a prospective 100-TeV hadron collider (dashed line), obtained by rescaling the bounds given in Ref. [10]. The remaining color coding is the same as in Fig. 2. See text for more details.
These estimates could potentially be loosened by uncertainties in the DM velocity profile and the local DM density. On the other hand, if DM is discovered in direct detection the relative size of the light-quark Yukawas could be probed by comparing scattering rates on different target materials.
An additional cross check of our scenario could be provided by searches for DM production at hadron colliders. In Fig. 12 we denote the 95% CL exclusion limits, assuming 3000/fb of data, at the 14-TeV LHC by a solid orange line and at a prospective 100-TeV hadron collider by a dashed orange line. These curves have been obtained by converting the bounds in Ref. [10] to the case of Dirac DM using FeynRules [111] and MadGraph5 [112]. We see that, while the LHC will be sensitive to a part of the interesting region in parameter space, the scenario of m χ = 75 GeV DM can be excluded only at a 100-TeV collider.

V. CONCLUSIONS
Not much is known experimentally about the couplings of the Higgs to light quarks. It is entirely possible that the Higgs couples only to the third generation of fermions. Experimen-tally equally viable is the possibility that the light-quark Yukawas are significantly enhanced, up to O(50) for κ s and up to O(10 3 ) for κ u and κ d . Such extremely large enhancements are not natural from a model-building point of view as they require a large fine tuning of the light-quark masses, but at present cannot be excluded experimentally.
Modified Yukawa couplings to light quarks could have implications for DM searches. In this paper we focused on Higgs-portal DM. We considered constraints on scalar, vector, and fermionic Higgs-portal models of DM from relic density, direct and indirect detection, and the invisible Higgs width. A central result of our analysis is that, for phenomenologically The effects in MFV models, multi-Higgs models with natural flavor conservation, and the type-II two-Higgs-doublet model (i.e., the tree-level Higgs sector of the MSSM), on the other hand, are expected to be much smaller.
Finally, we point out that a discovery of Higgs-portal DM in indirect searches would immediately imply an upper bound on the light-quark Yukawa couplings due to the upper bounds in direct DM searches.