Higgs EFT for 2HDM and beyond

We discuss the validity of the Standard Model Effective Field Theory (SM EFT) as the low-energy effective theory for the two-Higgs-doublet Model (2HDM). Using the up-to-date Higgs signal strength measurements at the LHC, one can obtain a likelihood function for the Wilson coefficients of dimension-6 operators in the EFT Lagrangian. Given the matching between the 2HDM and the EFT, the constraints on the Wilson coefficients can be translated into constraints on the parameters of the 2HDM Lagrangian. We discuss under which conditions such a procedure correctly reproduces the true limits on the 2HDM. Finally, we employ the SM EFT to identify the pattern of the Higgs boson couplings that are needed to improve the fit to the current Higgs data. To this end, one needs, simultaneously, to increase the top Yukawa coupling, decrease the bottom Yukawa coupling, and induce a new contact interaction of the Higgs boson with gluons. We comment on how these modifications can be realized in the 2HDM extended by new colored particles.


Introduction
Effective field theories (EFTs) allow one to describe the lowenergy dynamics of a wide class of quantum theories [1][2][3]. The idea is to keep only the subset of light degrees of freedom, while discarding the heavy ones that cannot be produced on-shell in the relevant experimental setting. Virtual effects of the heavy particles on low-energy observables are represented by an infinite series of operators constructed out of the light fields.
In the context of the LHC experiments, the light degrees of freedom are those of the Standard Model (SM), and the heavy a e-mail: hermes.belusca@th.u-psud.fr b e-mail: afalkows017@gmail.com c e-mail: duartefontes@ist.utl.pt d e-mail: jorge.romao@tecnico.ulisboa.pt e e-mail: jpsilva@cftp.ist.utl.pt ones correspond to hypothetical new particles. The lowenergy effective description of such a framework is called the SM EFT; see e.g. [4][5][6][7][8][9] for reviews. The SM EFT allows for a unified description of many possible signals of physics beyond the SM (BSM), assuming the new particles are too heavy to be directly produced. This model-independence is a great asset, given we currently have little clue about the more complete theory underlying the SM. Another strength of this approach is that constraints on the EFT parameters can easily be translated into constraints on masses and couplings in specific BSM constructions. Thus, once experimental results are interpreted in the EFT language, there is no need to reinterpret them in the context of every possible model out there.
A less appealing feature of EFTs is that the Lagrangian contains an infinite number of interaction terms and parameters, in contrast to renormalizable theories. In the SM EFT, these terms are organized in an expansion L eff = L SM + i c (6)  In practice, the series in Eq. (1.1) must be truncated, such that one works with a finite set of parameters. In most applications of the SM EFT, terms with D ≥ 8 are neglected. This corresponds to taking into account the BSM effects that scale as O(m 2 W / 2 ), and neglecting those suppressed by higher powers of . It is important to discuss the validity of such a procedure for a given experimental setting [10]. More precisely, the questions are: (1) whether the truncated EFT gives a faithful description of the low-energy phenomenology of the underlying BSM model, and (2) to what extent experimental constraints on the D = 6 Wilson coefficients are affected by the neglected higher-dimensional operators. Generically, in the context of LHC Higgs studies the truncation is justified if is much larger than the electroweak scale. But, to address the validity issue more quantitatively and identify exceptional situations, it is useful to turn to concrete models and compare the description on physical observables in the full BSM theory with that in the corresponding low-energy EFT. Such an exercise provides valuable lessons about the validity range and limitations of the SM EFT.
In this paper we perform that exercise for the Z 2symmetric CP-conserving two-Higgs-doublet model (2HDM). We compare the performance of the full model and its low-energy EFT truncated at D = 6 to describe the Higgs signal strength measurements at the LHC. To this end, we first update the tree-level constraints on the 2HDM parameter space using the latest Higgs data from Run-1 and Run-2 of the LHC. We use the same data to derive leading-order constraints on the parameters of the SM EFT. Given the matching between the EFT and the 2HDM parameters [11][12][13], the EFT constraints can be subsequently recast as constraints on the parameter space of the 2HDM. By comparing the direct and the EFT approaches, we identify the validity range of the EFT framework where it provides an adequate description of the impact of 2HDM particles on the LHC Higgs data.
We also remark that neither the SM nor the 2HDM provides a very good fit to the Higgs data, mostly due to some tension with the measured rate of the tth production and h → bb decays. If the current experimental hints of an enhanced tth and suppressed h → bb are confirmed by the future LHC data, the 2HDM alone will not be enough to explain these. Here the EFT approach proves to be very useful in suggesting extensions of the 2HDM that better fit the current Higgs data. In particular, we show that a good fit requires simultaneous modifications of the EFT parameters controlling the top and bottom Yukawa couplings and the contact interaction of the Higgs boson with gluons. We show how these modifications can be realized in the 2HDM extended by new colored particles coupled to the Higgs. This paper is organized as follows. In Sect. 2 we review the 2HDM and its low-energy EFT. In Sect. 3 we compare the direct and the EFT constraints on the parameter space imposed by the Higgs measurements. In Sect. 4 we discuss how to improve the fit to the LHC Higgs data by extending the 2HDM with new colored states coupled to the Higgs.

CP-conserving 2HDM
We start by reviewing the (non-supersymmetric) 2HDM [14][15][16], closely following the formalism and notation of Ref. [17]. We consider two Higgs doublets 1 and 2 , both transforming as (1, 2) 1/2 under the SM gauge group. Both doublets may develop a vacuum expectation value (VEV) parametrized as We assume that all parameters in the scalar potential are real, which implies the Higgs sector preserves the CP symmetry at the leading order.
Furthermore, we assume that the Lagrangian is invariant under a discrete Z 2 symmetry, under which the doublets transform as 1 → + 1 and 2 → − 2 . This symmetry is allowed to be broken only softly, that is to say, only by mass parameters in the Lagrangian. The Z 2 symmetry constrains the possible form of Yukawa interactions. There are four possible classes of 2HDM, depending on how the SM fermions transform under the Z 2 symmetry. They are summarized in the following It is often more convenient to work with linear combinations of 1 and 2 defined by the rotation , H 0 2 = 0. Note that H 1 and H 2 , unlike i , are not eigenstates of the Z 2 symmetry. The linear combinations H i define the so-called Higgs basis [18], while the original doublet i are referred to as the Z 2 basis.
In the Higgs basis, the scalar potential takes the form where the parameters Y i and Z i are all real. The Z 2 symmetry is manifested by the fact that only 5 of the Z i are independent, as they satisfy 2 relations: , where Z 345 ≡ Z 3 + Z 4 + Z 5 . The Yukawa couplings are given by (2.4) whereH i = iσ 2 H * i , and the coefficients of the H 2 Yukawa couplings are summarized in the table below: In the Higgs basis, the doublets can be parametrized as where G ± and G z are the Goldstone bosons eaten by W ± and Z , while H ± and A are the charged scalar and neutral pseudo-scalar eigenstates. The two neutral scalars h, H 0 are mass eigenstates, while the parameter c β−α ≡ cos(β − α) determines their embedding in the two doublets H i . 1 In the following we will identify h with the 125 GeV Higgs boson. The equations of motion for H 1 and H 2 imply the vacuum relations The masses of the charged scalar and the pseudo-scalar are given by The mixing angle is related to the parameters of the potential by The masses of the neutral scalars can be written as 1 The angle α can be defined as the rotation angle connecting the components of the original Higgs doublets 1 and 2 to the mass eigenstates. (2.9) Finally, the couplings of the CP-even scalar, h, to the electroweak gauge bosons are given by (2.10) and to the fermions by By convention, the sign of the h couplings to W W and Z Z is fixed to be positive (this can always be achieved, without loss of generality, by redefining the Higgs boson field as h → −h). On the other hand, the sign of the h couplings to a fermion may be positive or negative, depending on the value of c β−α and tan β. The alignment limit is defined by c β−α → 0, that is to say, when h has SM couplings. There is a strong evidence, both from Higgs and from electroweak precision measurements, that the couplings of the 125 GeV boson to W and Z bosons are very close to those predicted by the SM. Therefore the 2HDM has to be near the alignment limit to be phenomenologically viable. From Eq. (2.8), the condition for alignment is One way to satisfy this is by making Y 2 large, Y 2 v 2 , which is called the decoupling limit because then H 0 , A and H + become heavy. Another way to ensure alignment is to take |Z 6 | small enough, |Z 6 | 1. If the condition Eq. (2.12) is satisfied with Y 2 v 2 then we speak of alignment without decoupling.

Low-energy EFT
Eqs. (2.7) and (2.9) imply m A ∼ m H + ∼ m H 0 ∼ , and the spectrum below the scale v is that of the SM. Consequently, we can describe Higgs production and decays at the LHC in the framework of the so-called SM EFT, where the heavy particles are integrated out, and their effects are represented by operators with canonical dimensions D > 4 added to the SM. Below we discuss the Lagrangian of the low-energy effective theory for the 2HDM, treating 1/ as the expansion parameter. We first review the known results concerning the D = 6 operators in the EFT with tree-level matching [11,12]. This is enough for the purpose of this paper, in which the main focus is the accuracy of the EFT to describe the current LHC Higgs measurements. Matching beyond D = 6 and tree level was discussed in Refs. [11,13,19,20], and we will come back to it in an upcoming publication [21].
The simplest way to derive the tree-level matching is by integrating out the field H 2 and identifying H 1 with the SM Higgs doublet. The procedure is to: (1) solve the linearized equations of motion for H 2 as a function of the light fields (the scalar doublet H 1 , fermions, and gauge fields), and (2) insert the solution in the original Lagrangian. Furthermore, restricting our study to D = 6 operators in the EFT, one can ignore all derivative terms in the H 2 equation of motion. The linearized equation of motion for H 2 with derivative terms dropped is solved as Plugging this back, renaming H 1 → H , and keeping terms up to 1/ 2 , the effective Lagrangian takes the form 14) The terms proportional to Y 3 can be absorbed in a redefinition of the SM parameters, and they do not have observable consequences. On the other hand, the genuine D = 6 terms in Eq. (2.14) are in principle observable. We are interested in the impact of these D = 6 operators on the Higgs boson couplings probed at the LHC. Quite generally, in the SM EFT with D = 6 operators the CP-conserving Higgs boson couplings to two SM fields can be parametrized as [9,22] The effect of D = 6 operators in Eq. (2.14) is to shift the Higgs couplings to the SM fermions and to itself: On the other hand, at tree level and restricting our study to dimension-6 operators in the EFT, there are no corrections to the Higgs boson interactions with gauge bosons: One can check that the couplings of the Higgs in the effective theory described by the Lagrangian in Eq. (2.14) are the same as the couplings of h in the 2HDM expanded to linear order in c β−α , once we identify: This identification is consistent with Eq. (2.8) when Y 2 v 2 . We also comment on the interesting case of alignment without decoupling. Our EFT is supposed to be a good description of the 2HDM in the decoupling limit where all the additional scalars are heavy. In general, the EFT will not work when one or more scalars are light, even in the alignment limit. Indeed, if one of the new Higgs scalars are light, 2 → 2 fermion scattering will display a pole at the energy equal to the scalar's mass, which cannot be captured by the 4-fermion operators in Eq. (2.14). Similarly, double Higgs production will have a pole at the new mass (if the other neutral scalar is light), which again cannot be described by the operators Eq. (2.14).
However, it is possible that certain low-energy observables can still be adequately described by our EFT, even when the 2HDM has additional light scalars with mass of order m h . The Higgs couplings to matter are such observables, provided the 2HDM is in the alignment limit. More precisely, from the constraints on the couplings δy f one can correctly infer constraints on the parameters of the 2HDM in the limit of alignment without decoupling. However, to this end, the mapping between the parameters of the EFT and the 2HDM has to be modified: instead of Eq. (2.16), we have to use the following map: This formula follows from expanding the 2HDM expressions for the Yukawa couplings of h to the leading order in c β−α . Using this result, Eq. (2.19) can be obtained by expanding the 2HDM Higgs couplings in c β−α and using the expression for c β−α , which is also valid for alignment without decoupling.
The new terms in the matching formulas are negligible in the decoupling limit Y 2 v 2 , in which case they are higher order in the v 2 / 2 expansion. However, they can be very important in the case of alignment without decoupling when Y 2 v 2 . Such a way of extending the validity range of the  EFT by adding higher-order terms in the matching formula is similar to v-improved matching advocated in Refs. [12,23].

Comparison of EFT and 2HDM descriptions of Higgs couplings
In this section we discuss constraints from the Higgs signal strength measurements at the LHC. To this end, we use the results summarized in Table 1, which also include preliminary Run-2 results. First we update the tree-level constraints on the c β−α -tan β plane of the various Z 2 -symmetric versions of the 2HDM. The same LHC data can also be used to derive leading-order constraints on the parameters of the SM EFT with D = 6 operators. These can be subsequently recast as constraints on the 2HDM parameters using the treelevel matching in Sect. 2.2. As long as the extra scalars of the 2HDM are heavy, we expect that the EFT should give an adequate description of the Higgs physics, and then the constraints should be the same regardless whether we obtain them directly or via the EFT. The goal of this section is to validate this expectation and quantify the validity range of the EFT for the 2HDM. Finally, we will also compare the results obtained by the above analyses with more sophisticated parameter scans of the 2HDM, that take into account the limits from precision measurements, unitarity, and boundedness of the Higgs potential.

Update of Higgs constraints on 2HDM
We first show the constraints on various 2HDM scenarios from the LHC studies of the 125 GeV Higgs. 2 The results are shown in Figs. 1 and 2. The constraints are the weakest in the type-I model, especially for large tan β. The reason is that, in this case, corrections to fermionic Higgs couplings are suppressed by tan β for a fixed c β−α . Although corrections to the Higgs couplings to W and Z do not have that suppression, they appear only at the quadratic order in c β−α and therefore they become significant only for large |c β−α |. As a result, in the type-I model the 125 GeV Higgs boson can be further away from the alignment limit, with the modest bound |c β−α | 0.4 for large enough tan β.
The constraints on c β−α are much more stringent for the type-II, type-X, and type-Y scenarios. In those cases, for a fixed c β−α , there is always a modification of some fermionic Higgs coupling that is not suppressed by tan β. In the generic region of the parameter space the bound is |c β−α | 0.1-0.2 for tan β ∼ 1, and even more stringent for smaller and larger tan β. These scenarios also display a separate region Footnote 2 continued c β−α -tan β would be obtained if we instead used experimental constraints on the κ coupling modifiers. of the parameter space where a large c β−α is allowed. It corresponds to the situation when corrections to the down-type quark and/or lepton Yukawa couplings flip their sign but leave the absolute values close to the SM one [41][42][43]. Note that 3 distinct situations can arise: when the down-type Yukawas become negative (type-Y), when the lepton Yukawas become negative (type-X), and when both become negative (type-II). We refer to all these 3 cases as the "wrong-sign Yukawa" region. The Higgs observables are currently weakly sensitive to the sign of the bottom and tau Yukawa, therefore these somewhat fine-tuned regions remain consistent with the data. Future precision tests may resolve the sign of the bottom Yukawa [41,42,[44][45][46][47][48], but that may be challenging for the tau Yukawa.
The qualitative shape of the favored regions is the same as that obtained from Run-1 Higgs data. The effect of the preliminary Run-2 data is to make the constraints somewhat more stringent.

Higgs constraints on EFT
As explained in Sect. 2.2, at the leading order in the 1/ 2 expansion the 2HDM induces corrections only to Higgs Yukawa-and self-interactions. The latter can be probed by non-resonant double Higgs production but, given the current level of precision, the existing limits on the Higgs cubic self-coupling do not lead to any interesting constraints on the 2HDM parameter space. Therefore, at order 1/ 2 , the parameters of the SM EFT relevant for the 2HDM are the three δy f characterizing corrections to the SM Higgs Yukawa couplings; see Eq. (2.15). As an intermediate step in connecting the 2HDM to the SM EFT, we can derive the constraints on these 3 EFT parameters. We find that the Run-1 and Run-2 Higgs data lead to the following constraints:  1). Incidentally, the constraints change very little (by less than 20%) if the EFT predictions are not expanded to a linear level, but instead the full nonlinear dependence on δy f is retained. In such a case, the likelihood is highly non-Gaussian, but it can nevertheless be well approximated by a Gaussian likelihood in the parameter space region with δy f 1, which is preferred by the LHC Higgs data. The main qualitative consequence of using the full non-Gaussian likelihood is the existence of other nearly degenerate minima (in addition to the one described by Eq. (3.1)) where some δy f are large. For example, the full likelihood contains minima with δy d ∼ −2 and/or δy e ∼ −2, as the Higgs data are currently very weakly sensitive to the sign of the bottom and tau Yukawa.
The EFT likelihood defined by Eq. (3.1) can be recast into a 2HDM likelihood by inserting the relation between δy f and the 2HDM parameters. For example, in order to obtain constraints in the c β−α -tan β plane we need to read off from Eq. (2.11) Of course this procedure cannot be in any way better than deriving the limits on c β−α and tan β directly, as we did in Sect. 3.1. The purpose of this exercise is to investigate how useful EFT is as a tool to constrain various BSM scenarios. The idea is that the LHC experiments present the EFT likelihood like the one in Eq. (3.1), or a more general one depending on a larger number of EFT parameters that can be subsequently projected into the δy f subspace. That likelihood function can be recast to quickly obtain constraints on a host of BSM models. Our exercise is a case study for the validity of the EFT approach to LHC Higgs data, which allows one to understand limitations of the EFT and avoid possible pitfalls.
The results of our exercise are shown in Fig. 3. We compare the favored regions in the c β−α -tan β plane for the type-I and type-II 2HDM obtained in Sect. 3.1 with the ones deduced by using constraints on EFT parameters together with the EFT-2HDM matching discussed earlier.
Let us start the discussion with the type-II scenario. In this case, a recast of the Gaussian likelihood defined by Eq. (3.1) provides a very good description of the bulk of the favored region where c β−α is small. In that region, the LHC Higgs data force the deviations of the Yukawa couplings to be small, less than ∼30% of the SM value. Such small deviations can be adequately described by D = 6 operators of the SM EFT, and the O(1/ 4 ) contributions to Higgs observables can be safely neglected. As the experimental precision increases, and assuming that no large deviations from the SM are reported, this conclusion will only be strengthened, and the agreement between the EFT and the complete description will further improve. On the other hand, we can see that the Gaussian EFT likelihood completely misses the existence of the "wrong-sign" Yukawa region. This is inevitable, as a Gaussian likelihood has only one minimum, and therefore it cannot capture other degenerate minima in the parameter space. The situation can be improved by complicating the description on the EFT side, and instead including all higher-order terms in δy f in the likelihood function. Such a non-Gaussian likelihood is capable of describing multiple minima, including the one in the "wrong-sign" region where one or more δy f are smaller than −1. Indeed, we can see in Fig. 3 that using the non-Gaussian EFT likelihood qualitatively captures the shape of the "wrong-sign" minimum, at least when c β−α is not too large. Using the non-Gaussian likelihood also improves the agreement between the EFT and the direct 2HDM limits in the bulk region at small c β−α .
For type-X and type-Y the results are very similar to those for type-II: the EFT description captures very well the bulk of the favored parameter space with small c β−α , but it fails to capture the wrong-sign Yukawa region. Again, the latter Constraints from LHC Higgs data on the parameter space of the type-I (left) and type-II (right) 2HDM. We show the 68% (darker green) and 95% (lighter green) CL region in the c β−α -tan β plane. We compare it with the 95% CL region obtained by recasting the EFT limits in Eq. (3.1) which are derived from a Gaussian likelihood function (black dashed). We also show the boundary of the 95% CL region obtained by recasting the limits derived from a non-Gaussian EFT likelihood where higher-order terms in δy f are kept (red dotted) problem can be addressed by using the non-Gaussian likelihood on the EFT side.
In the type-I scenario a qualitatively new issue appears. In this case the EFT provides a good approximation of the favored region for low c β−α and tan β 2. However, it completely misses the relevant physics at larger c β−α and tan β. Namely, in the type-I 2HDM the LHC Higgs data imply an upper limit on |c β−α |, approximately |c β−α | 0.4. At large tan β, this limit is not driven by modifications of the Yukawa couplings, which are suppressed by tan β, but rather by modifications of the Higgs couplings to W W and Z Z. However, these appear only at O(1/ 4 ) in the low-energy EFT of the 2HDM, and they are not included at all in the SM EFT truncated at D = 6. In other words, the type-I 2HDM at tan β 1 belongs to an exceptional class of BSM scenarios that are not adequately described by a SM EFT with D = 6 operators. Instead, in the Higgs observables, the D = 8 operators in the low-energy EFT (formally O( −4 )) may dominate over the D = 6 ones (formally O( −2 )), as the latter are suppressed by tan β and the former are not. This is an example of selection rules in the UV theory modifying the naive power counting in the low-energy EFT. As a consequence, the D = 6 EFT approach in this case misrepresents the allowed parameter space of the type-I 2HDM at large tan β. Note that the problem is not addressed by switching from a Gaussian to a non-Gaussian EFT likelihood. A more general low-energy approach is needed to capture this scenario, for example the SM EFT truncated at the level of D = 8 operators, or a more phenomenological non-EFT approach.

Parameter scans
So far we have limited ourselves to studying the constraints on the 2HDM resulting from the LHC studies of the 125 GeV Higgs boson. These select an interval(s) for the allowed values of the Higgs mixing angle c β−α , depending on the 2HDM scenario and on the value of the tan β parameter. However, there exist further important constraints on the 2HDM. First of all, the neutral scalar and pseudo-scalar, and the charged partners of the Higgs boson are targeted by direct searches in high-energy colliders. Moreover, the Higgs partners may contribute to electroweak precision observables, in particular to the S, T , and U parameters [49], or to Z → bb decays [41,50]. Finally, the parameters of the Higgs potential should satisfy the theoretical constraints following from perturbative unitarity [51][52][53], and the Higgs potential should be bounded from below [54]. In the EFT approach, all of this information is not used. In particular, the heavy Higgses are integrated out from the spectrum. The natural question then is whether the region in the c β−α -tan β plane selected by the 125 GeV Higgs data can be realized in the full 2HDM given the existing constraints.
To address this question, we have performed scans of the 2HDM parameter space. Our results are shown in Fig. 4. We set m h = 125 GeV, and generate points with −π/2 ≤ α ≤ π/2, 0 ≤ tan β ≤ 10, 125 GeV ≤ m H 0 ,A,H + ≤ 2 TeV. For the type-II scenario we demand in addition m H + ≥ 480 GeV, so as to satisfy the b → sγ constraints [47,[55][56][57][58]. In that range, we search for parameter points that pass the 125 GeV Higgs constraints evaluated earlier, as well as the unitarity, boundedness, and electroweak precision constraints. The latter constraints turn out to be non-trivial, in the sense that they eliminate a fraction of points that would pass the Higgs constraints alone. However, in our simulations they seem not to eliminate any particular value of c β−α and tan β favored by the Higgs data. This is known as the emmental effect: excluded regions in the multi-dimensional parameter space of the 2HDM do not show after a projection onto the twodimensional c β−α -tan β plane.
For type-II, the bulk of the allowed parameter space with |c β−α | 1 contains points with extra scalar masses ranging from very heavy ( 1 TeV) to very light (∼125 GeV), corresponding to alignment with and without decoupling, respectively. That region is also recovered (to a good approximation) by recasting the Gaussian EFT analysis of Higgs data into constraints on 2HDM parameters. Our scan shows that this entire region can be realized in the 2HDM with all extra scalars decoupled at the LHC energies. In such a case, the heavy states are not accessible directly, and their only observable effect is to modify the couplings of the 125 GeV Higgs.
In the wrong-sign region, on the other hand, the extra scalars need to be relatively light, m H 0 ,A,H + 700 GeV. This implies the heavy states are not decoupled at the LHC energies, and can be relevant for direct searches and resonant double Higgs production analyses. The main theoretical constraints in this region come from the unitarity bounds, because large values of some scalar couplings in the Higgs potential are required, λ i 5. Our unitarity bounds on λ i are those obtained at LO by Refs. [51][52][53]. Even stronger conclusions were obtained in Ref. [59] which takes into account unitarity constraints on NLO (one-loop) partial wave amplitudes. These, together with the bounds on the charged Higgs mass, disfavor at 95% CL the entire wrong-sign region in the type-II and type-Y scenarios, although still allow for a wrong-sign solution of the lepton couplings in the type-X scenario. In this paper we only impose the LO unitarity bounds, since a full one-loop renormalized calculation of theoretical bounds, relations between λ i and masses, mixing angles, decay rates, and production cross-sections is still missing to our knowledge.
In any case, whether or not the wrong-sign region is consistent with the latest experimental and theoretical constraints, the related LHC phenomenology is strictly speaking not amenable to an EFT description. We conclude that, for the type-II scenario, the SM EFT approach at O( −2 ) is adequate in the entire parameter space allowed by the experimental data and where the additional scalars are heavy. Similar conclusions hold for the type-X and type-Y scenarios.
The situation is somewhat different for the type-I scenario. As we discussed in Sect. 3.2, the discrepancy between the full model and the EFT description is important, especially at large tan β. The reason is that in this case the numerically largest effects on the Higgs boson couplings are O( −4 ) and correspond to D = 8 operators, whereas the formally leading O( −2 ) effects, captured by the D = 6 SM EFT, are suppressed by tan β. This problem will always be present at large enough tan β even when precision of Higgs measurements is improved significantly. However, the scan in Fig. 4 shows that the parameter space where the two descriptions disagree about the Higgs coupling bounds is dominated by points with very light extra scalars. Thus, much like in the type-II case, most of the 2HDM parameter space where the EFT approach fails to yield the correct bounds is anyway the one where the extra scalars do not decouple at the LHC energies.

Discussion
Working with the SM EFT one always needs to make a compromise between generality and simplicity. In principle, the full EFT Lagrangian contains all information as regards the effects of heavy new physics on low-energy observables, but that information is encoded in an infinite number of parameters. The usual approach is to truncate the EFT expansion at the level of D = 6 operators, which corresponds to retaining the effects up to O( −2 ) in the new physics mass scale . If that is done consistently, that is the EFT predictions for the observables are expanded up to O( −2 ) then lowenergy measurements, such as the LHC Higgs signal strength observables, can be translated into a Gaussian likelihood for the D = 6 EFT parameters. This allows for a very concise presentation of results, as a Gaussian likelihood is fully specified by the central values, 1 σ uncertainties, and correlation matrix of the parameters. Thus, using the SM EFT at O( −2 ), the large amount of data contained in multiple Higgs analyses at the LHC can be summarized by just a handful of numbers that can later be recast to provide constraints on a large class of BSM scenarios.
The question is how much information as regards the UV physics is lost due to these dramatic simplifications. This can be addressed quantitatively by comparing the performance of complete UV models and the D = 6 EFT approximating the low-energy physics of those models. Our case study of 2HDM scenarios and their EFTs leads to a few interesting conclusions. First of all, the Gaussian likelihood provides a very good approximation of the new physics effects in the bulk of the allowed parameter space. We, however, identified the exceptional situations where this is not the case: 1. In the wrong-sign Yukawa regions of type-II, -X and -Y scenarios, where the relative corrections to the Yukawa bottom and/or tau couplings are large and cannot be properly described at O( −2 ). 2. For the type-I scenario at large tan β, where the leading effects on the single Higgs production and decay come from D = 8 operators in the EFT, which are by default neglected when the SM EFT is truncated at D = 6.
One can always complicate the EFT framework such that it is capable of describing also these special cases. The wrong-sign region can be captured in the SM EFT with only dimension-6 operators if, when calculating the Higgs signal strength observables, one retains in the final expressions the terms which are non-linear in D = 6 Wilson coefficients. In this approach the likelihood function contains terms which are higher than quadratic in the parameters of interests, therefore it is non-Gaussian. Generally speaking, in the SM EFT with dimension-6 operators one controls the expansion only up to the order 1/ 2 , and thus keeping O( −4 ) and higherorder terms may not be consistent. However, as discussed in detail in Ref. [10], keeping higher-order terms can be justified in some circumstances without invalidating the EFT expansion. In particular, these higher-order terms are crucial to properly describe the situation when new physics contributions to observables are comparable or exceed the SM ones. We have shown here that retaining these terms allows one to approximately reproduce the wrong-sign regions in the 2HDM, at the price of introducing non-Gaussian terms into the likelihood.
To cover the large tan β region of the type-I scenario the EFT Lagrangian would have to be extended to include D = 8 terms. One can understand why and when these terms can be important by looking back at Eqs. (2.10) and (2.11). In the full theory, the terms proportional to s β−α come from the rotation from the Higgs basis H 1 , which carries all the VEV v, into the lightest scalar in the mass basis, h. Using these expressions in the EFT, by expanding 1 − c 2 β−α to higher orders in c β−α and making the identification c β−α ↔ − Z 6 v 2 2 , we get a perfect match to the full theory. Granted, there will be other D > 6 terms in a general EFT, but these are the ones needed for the matching with the 2DHM.
One interesting point is how such terms come about when integrating out H 2 . The matrix taking the Higgs basis (H 1 , H 2 ) into the mass basis (h, H ), is orthogonal with angle β−α. But, by integrating out H 2 and discarding H , the "truncated rotation" of H 1 into h by s β−α violates unitarity, turning the derivative terms into non-canonical form. This is what happens when matching to the EFT, where the problem is fixed at each order in v/ by rescaling the field h. The procedure introduces new D > 6 terms which are exactly those which would appear by expanding Eqs. (2.10) and (2.11) to higher order in c β−α . Thus, in the 2HDM case, the validity or lack thereof of the D = 6 EFT is determined by c β−α . From Fig. 3 we conclude that the dimension-8 operators become non-negligible for |c β−α | 0.3.
Accommodating the wrong-sign regions in the 2HDM and the large tan β region of the type-I scenario would make it more challenging to perform EFT analyses at the LHC and present their results. In our opinion, the Gaussian approach with the EFT Lagrangian truncated at D = 6 may be more productive in the context of Higgs signal strength observables. This simple approach is sufficient in generic situations, while the special cases described above can be treated separately. Indeed, our parameter scans show that the special cases are always associated with the extra scalars being not much heavier than the 125 GeV Higgs boson, and therefore they should be probed directly using the complete model description and without passing through the EFT.

Beyond 2HDM
The LHC measurements of the Higgs signal strength summarized in Table 1 show some tension with the SM predictions. On the one hand, there is an excess in the tth production mode appearing in several Higgs decay channels. On the other hand, the signal strength in the h → bb decay channel is suppressed for several production modes. Assuming for a moment this is not merely a statistical fluctuation, the data may point to the Higgs-top (-bottom) coupling being 30% larger (smaller) than in the SM. Within the 2HDM, it is straightforward to arrange the Higgs couplings to top quarks to be enhanced, and the Higgs coupling to bottom quarks to be simultaneously suppressed. This happens in the type-II and type-Y models at tan β ∼ 1 and c β−α > 0. However, these regions of the parameter space are not favored by the global fits showed in Fig. 1 or in Fig. 2. More generally, in the 3-parameter EFT fit in Eq. (3.1) the SM point where all δy f = 0 is not significantly disfavored, with χ 2 SM −χ 2 min ≈ 2. The reason is that increasing the Higgs-top coupling also increases the gluon fusion cross section via the 1-loop top triangle diagram contribution to the gg → h amplitude. Since the measured total Higgs cross section (which is dominated by gluon fusion) agrees very well with the SM predictions, simply increasing the top-Higgs couplings is not preferred by global fits. Decreasing the Higgs-bottom coupling is disfavored for similar reasons. As the Higgs width is dominated by decays to bottom quarks, a smaller Higgs-bottom coupling increases the Higgs branching fractions (and thus the signal strength) into other final states. In a global fit, the gain from fitting better the suppressed h → bb channels is outweighed by overshooting the signal strength in the precisely measured W W , Z Z, and γ γ final states.
The above discussion suggests a simple ad hoc solution to improve the global fit in a theory with two Higgs doublets. One can arrange additional contributions to the effective Higgs-gluon coupling beyond those from integrating the top quark and other SM fermions. If the sign of these contributions was opposite to that induced by the top, the new physics could cancel the effect of the increased Yukawa in the gluon fusion Higgs production cross section. We can parametrize these new contributions by adding a new term in the 2HDM Lagrangian where G a μν is the gluon field strength, and g s is the SM strong coupling. The parameter c gg encodes the effects of heavy colored particles beyond the 2HDM on the Higgs phenomenol-ogy. For example, integrating out a new color octet scalar S a of mass m S coupled to the Higgs sector via the interaction term +λ S |H 1 | 2 S a S a , one finds c gg = − . Similar extensions of the 2HDM have been considered in the past; see e.g. [60][61][62].
We first employ the linearized EFT approach to see whether allowing the parameter c gg to vary freely can lead to an improvement of the Higgs fit. The Higgs boson couplings are those in Eq. (2.15) with non-zero δy f and c gg and the remaining coupling set to zero. With that assumption, the Run-1 and Run-2 Higgs data lead to the following constraints: Now the preferred values of the EFT parameters are significantly away from the SM point. Indeed, we find χ 2 SM −χ 2 min ≈ 11, which translates to the 2.3σ preference for BSM. We also checked that allowing for more free parameters in the EFT (e.g. c γ γ ) does not lead to further significant improvement of the fit.
The EFT results in Eq. (4.2) suggest that the Higgs fit can be improved also in the context of the type-II and type-Y scenario, once we allow for the new contributions to the Higgs-gluon coupling. This is indeed the case, as shown in Fig. 5. The best fit regions move away from the SM limit where c β−α = 0 to c β−α > 0 where δy u > 0 and δy d < 0 are effectively generated.
For the type-II case the best fit point occurs for c β−α ≈ 0.29, tan β ≈ 1.0, and c gg ≈ −3.9 × 10 −3 , and has χ 2 SM − χ 2 min ≈ 15. The minimum is slightly deeper than in the 4-parameter EFT fit because here we use the full (not Gaussian) likelihood function. For such low tan β constraints from flavor physics become non-trivial and require m H + 480 GeV [47,[55][56][57]. Nevertheless, this limit does not pose consistency problems, as c β−α ∼ 0.3 can be obtained with perturbative couplings in the scalar potential as long as m H + 1.5 TeV. The preferred value of c gg requires a large contribution to the effective Higgs-gluon coupling from new particles, approximately one half (in magnitude) that of the top quark in the SM. In the example with a scalar octet we need λ S ≈ 11 × (m S /TeV) 2 , thus the octet needs to be below the TeV scale for λ S to remain perturbative. Note that current LHC and Tevatron data still do not exclude fairly light colored particles in a model-independent way; see e.g. [62,63] for a recent discussion. On the other hand, the For the type-Y case the best fit point falls actually into the wrong-sign region, at c β−α ≈ 0.35, tan β ≈ 4.3, c gg ≈ −2.7 × 10 −3 and has χ 2 SM − χ 2 min ≈ 15. However, it is not strongly preferred (χ 2 lower by just 0.3) over the local minimum at c β−α ≈ 0.16, tan β ≈ 1.9, and c gg ≈ −2.6 × 10 −3 where all Yukawas are positive. The higher tan β and lower c β−α at the local minimum in the type-Y case are somewhat easier to accommodate than the best fit point for type-II. For example, in the scalar octet case we need λ S ≈ 7×(m S /TeV) 2 , and the flavor physics bounds on m H + are not relevant for the preferred tan β.
For the case of type-I and type-X models we do not find any significant improvement of the fit after introducing the parameter c gg . This is due to the fact that in these scenarios δy u = δy d , therefore one cannot simultaneously fit the enhanced tth and suppressed h → bb signal.

Summary
In this paper we discussed the validity of the SM EFT with D = 6 operators as a low-energy theory for the 2HDM. Working consistently at O( −2 ) in the EFT expansion, the LHC Higgs signal strength measurements can be recast into a Gaussian likelihood for the EFT Wilson coefficients. That likelihood can then be used to place constraints on the parameter space of various extensions of the SM, once the matching between the BSM model and its low-energy EFT is known. We applied this procedure for the case of the CP-conserving 2HDM, restricting our study to the tree-level matching. We then compared the resulting constraints on the c β−α -tan β plane with those derived directly without passing through the EFT. We find that, in the bulk of the allowed parameter space of the 2HDM where c β−α is small, the Gaussian likelihood approximates very well the effects of the new scalars on the Higgs phenomenology. In those regions, the SM EFT truncated at D = 6 provides a valid description of the 2HDM phenomenology, as long as the extra scalars are heavy enough such that they do not appear on-shell in LHC Higgs observables.
However, we also identified the situations where our EFT procedure miscalculates the impact of the 2HDM on Higgs physics, even when m h . One occurs when some SM Yukawa coupling receives corrections that are comparable to its SM value, which happens in particular in the wrong-sign Yukawa regions. Another occurs for the type-I scenario at large tan β, where the leading 2HDM effects on Higgs phenomenology are encoded in D = 8 operators of the lowenergy EFT. These two exceptions are important to keep in mind when EFT results are interpreted as constraints on BSM, as they are representative of a wider class of models. It is possible to generalize the EFT approach such that it becomes adequate also in the above situations, but that would come at the price of a greater complexity of the analysis and a less transparent presentation.
We also applied the EFT approach to investigate what deformations of the SM Higgs couplings are needed to improve the fit to the Higgs data. According to Eq. (4.2), this requires simultaneously: (1) increasing the top Yukawa coupling, (2) decreasing the bottom Yukawa coupling, and (3) inducing the contact interaction of the Higgs boson with gluons. We discussed how these modifications can be realized in the 2HDM extended by new colored particles coupled to the Higgs. Future analyses of the LHC data from the 2016 run will tell whether the current small tension between the measurements and the SM predictions is just due to a statistical fluctuation, or due to new physics contributions to the Higgs boson couplings.