LPT Orsay 1664
CFTP/16013
Higgs EFT for 2HDM and beyond
Hermès BéluscaMaïto, Adam Falkowski, Duarte Fontes,
Jorge C. Romão, and João P. Silva
Laboratoire de Physique Théorique, CNRS, Univ. ParisSud, Université ParisSaclay, 91405 Orsay, France
CFTP, Departamento de Física, Instituto Superior Técnico, Universidade de Lisboa,
Avenida Rovisco Pais 1, 1049 Lisboa, Portugal
We discuss the validity of the Standard Model Effective Field Theory (SM EFT) as the lowenergy effective theory for the twoHiggsdoublet Model (2HDM). Using the uptodate Higgs signal strength measurements at the LHC, one can obtain a likelihood function for the Wilson coefficients of dimension6 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.
1 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 onshell in the relevant experimental setting. Virtual effects of the heavy particles on lowenergy 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 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 modelindependence 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 be easily 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
(1.1) 
where is the SM Lagrangian, is the mass scale of BSM physics, each is an invariant operator of canonical dimension , and the parameters are called the Wilson coefficients. Terms with odd are absent assuming baryon and lepton number conservation.
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 are neglected. This corresponds to taking into account the BSM effects that scale as , 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 lowenergy phenomenology of the underlying BSM model, and 2) to what extent experimental constraints on the =6 Wilson coefficients are affected by the neglected higherdimensional 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 lowenergy 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 symmetric CPconserving twoHiggsdoublet model (2HDM). We compare the performance of the full model and its lowenergy EFT truncated at =6 to describe the Higgs signal strength measurements at the LHC. To this end, we first update the treelevel constraints on the 2HDM parameter space using the latest Higgs data from Run1 and Run2 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], 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 production and decays. If the current experimental hints of an enhanced and suppressed 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 Section 2 we review the 2HDM and its lowenergy EFT. In Section 3 we compare the direct and the EFT constraints on the parameter space imposed by the Higgs measurements. In Section 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.
2 Formalism
2.1 CPconserving 2HDM
We start by reviewing the (nonsupersymmetric) 2HDM [13, 14, 15], closely following the formalism and notation of Ref. [16]. We consider two Higgs doublets and , both transforming as under the SM gauge group. Both doublets may develop a vacuum expectation value (VEV) parametrized as , with , , and GeV. 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 symmetry under which the doublets transform as and . This symmetry is allowed to be broken only softly, that is to say, only by mass parameters in the Lagrangian. The symmetry constrains the possible form of Yukawa interactions. There are four possible classes of 2HDM, depending on how the SM fermions transform under the symmetry. They are summarized in the following table:
TypeI  TypeII  specific (TypeX)  Flipped (TypeY)  

Uptype  
Downtype  
Leptons 
It is often more convenient to work with linear combinations of and defined by the rotation
(2.1) 
It follows that , . Note that and , unlike , are not eigenstates of the symmetry. The linear combinations define the socalled Higgs basis [17], while the original doublet are referred to as the basis.
In the Higgs basis, the scalar potential takes the form:
(2.2) 
where the parameters and are all real. The symmetry is manifested by the fact that only 5 of the are independent, as they satisfy 2 relations:
(2.3)  
where . The Yukawa couplings are given by:
(2.4) 
where , and the coefficients of the Yukawa couplings are summarized in the table below:
TypeI  TypeII  TypeX  TypeY  

In the Higgs basis, the doublets can be parametrized as
(2.5) 
where and are the Goldstone bosons eaten by and , while and are the charged scalar and neutral pseudoscalar eigenstates. The two neutral scalars , are mass eigenstates, while the parameter determines their embedding in the two doublets .^{1}^{1}1The angle can be defined as the rotation angle connecting the components of the original Higgs doublets and to the mass eigenstates. In the following we will identify with the 125 GeV Higgs boson.
The equations of motion for and imply the vacuum relations:
(2.6) 
The masses of the charged scalar and the pseudoscalar are given by:
(2.7) 
The mixing angle is related to the parameters of the potential by:
(2.8) 
The masses of the neutral scalars can be written as:
(2.9)  
Finally, the couplings of the CPeven scalar , to the electroweak gauge bosons are given by
(2.10) 
and to the fermions by
(2.11) 
By convention, the sign of the couplings to and is fixed to be positive (this can always be achieved, without loss of generality, by redefining the Higgs boson field as ). On the other hand, the sign of the couplings to fermion may be positive or negative, depending on the value of and . The alignment limit is defined by , that is to say, when 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 and 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:
(2.12) 
One way to satisfy this is by making large, , which is called the decoupling limit because then , and become heavy. Another way to ensure alignment is to take small enough, . If the condition Eq. (2.12) is satisfied with then we speak of alignment without decoupling.
2.2 Lowenergy EFT
For and , Eq. (2.7) and Eq. (2.9) imply , and the spectrum below the scale is that of the SM. Consequently, we can describe Higgs production and decays at the LHC in the framework of the socalled SM EFT, where the heavy particles are integrated out, and their effects are represented by operators with canonical dimensions added to the SM. Below we discuss the Lagrangian of the lowenergy effective theory for the 2HDM, treating as the expansion parameter. We first review the known results concerning the =6 operators in the EFT with treelevel 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 =6 and tree level was discussed in Refs. [18, 11], and we will come back to it in an upcoming publication [19].
The simplest way to derive the treelevel matching is by integrating out the field and identifying with the SM Higgs doublet. The procedure is to: 1) solve the linearized equations of motion for as a function of the light fields (the scalar doublet , fermions, and gauge fields), and 2) insert the solution in the original Lagrangian. Furthermore, restricting to operators in the EFT, one can ignore all derivative terms in the equation of motion. The linearized equation of motion for with derivative terms dropped is solved as:
(2.13) 
Plugging this back, renaming , and keeping terms up to , the effective Lagrangian takes the form:
(2.14) 
The terms proportional to can be absorbed in a redefinition of the SM parameters, and they do not have observable consequences. On the other hand, the genuine =6 terms in Eq. (2.14) are in principle observable. We are interested in the impact of these =6 operators on the Higgs boson couplings probed at the LHC. Quite generally, in the SM EFT with =6 operators the CPconserving Higgs boson couplings to two SM fields can be parametrized as [9, 20]:
(2.15)  
The effect of =6 operators in Eq. (2.14) is to shift the Higgs couplings to the SM fermions and to itself:
(2.16) 
On the other hand, at treelevel and restricting to dimension6 operators in the EFT, there are no corrections to the Higgs boson interactions with gauge bosons:
(2.17) 
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 in the 2HDM expanded to linear order in , once we identify:
(2.18) 
This identification is consistent with Eq. (2.8) when .
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, fermion scattering will display a pole at the energy equals to the scalar’s mass, which cannot be captured by the 4fermion 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 lowenergy observables can still be adequately described by our EFT, even when the 2HDM has additional light scalars with mass of order . The Higgs couplings to matter are such observables, provided the 2HDM is in the alignment limit. More precisely, from the constraints on the couplings 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:
(2.19) 
This formula follows from expanding the 2HDM expressions for the Yukawa couplings of to the leading order in . Using that, Eq. (2.19) can be obtained by expanding the 2HDM Higgs couplings in and using the expression for that is also valid for alignment without decoupling. The new terms in the matching formulas are negligible in the decoupling limit , in which case they are higher order in the expansion. However they can be very important in the case of alignment without decoupling when . Such a way of extending the validity range of the EFT by adding higher order terms in the matching formula is similar to improved matching advocated in Ref. [12, 21].
3 Comparison of EFT and 2HDM descriptions of Higgs couplings
Channel  Production  Run1  ATLAS Run2  CMS Run2 
[22]  [23]  
VBF  [22]  [23]  
    
    
  [22]    
[22]  [23]  
incl.      
[22]  [24]  
VBF  [22]  [24]  
    
VBF      
    
    
    
incl.      [25]  
    
VBF      
    
    
    
VBF    [26]  [27]  
    
    
  [28]    
[29]  [30]  
incl.  [31]    
multi  cats.    [32]  [33] 
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 Run2 results. First we update the treelevel constraints on the  plane of the various 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 =6 operators. These can be subsequently recast as constraints on the 2HDM parameters using the treelevel matching in Section 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.
3.1 Update of Higgs constraints on 2HDM
We first show the constraints on various 2HDM scenarios from the LHC studies of the 125 GeV Higgs. The results are shown in Fig. 1 and Fig. 2. The constraints are the weakest in the typeI model, especially for large . The reason is that, in this case, corrections to fermionic Higgs couplings are suppressed by for a fixed . Although corrections to the Higgs couplings to and do not have that suppression, they appear only at the quadratic order in and therefore they become significant only for large . As a result, in the typeI model the 125 GeV Higgs boson can be further away from the alignment limit, with the modest bound for large enough .
The constraints on are much more stringent for the typeII, typeX, and typeY scenarios. In those cases, for a fixed , there is always a modification of some fermionic Higgs coupling that is not suppressed by . In the generic region of the parameter space the bound is  for , and even more stringent for smaller and larger . These scenarios also display a separate region of the parameter space where a large is allowed. It corresponds to the situation when corrections to the downtype quark and/or lepton Yukawa couplings flip their sign but leave the absolute values close to the SM one [37, 38, 39]. Note that 3 distinct situations can arise: when the downtype Yukawas become negative (typeY), when the lepton Yukawas become negative (typeX), and when both become negative (typeII). We refer to all these 3 cases as the “wrongsign Yukawa” region. The Higgs observables are currently weakly sensitive to the sign of the bottom and tau Yukawa, therefore these somewhat finetuned regions remain consistent with the data. Future precision tests may resolve the sign of the bottom Yukawa [40, 41, 42, 37, 38, 43, 44], but that may be challenging for the tau Yukawa.
The qualitative shape of the favored regions is the same as that obtained from Run1 Higgs data. The effect of the preliminary Run2 data is to make the constraints somewhat more stringent.
3.2 Higgs constraints on EFT
As explained in Section 2.2, at the leading order in the expansion the 2HDM induces corrections only to Higgs Yukawa and selfinteractions. The latter can be probed by nonresonant double Higgs production but, given the current level of precision, the existing limits on the Higgs cubic selfcoupling do not lead to any interesting constraints on the 2HDM parameter space. Therefore,at order , the parameters of the SM EFT relevant for the 2HDM are the three 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 Run1 and Run2 Higgs data lead to the following constraints:
(3.20) 
where the quoted uncertainties correspond to 1 , and is the correlation matrix. The central values are close to the SM point, with . These results are obtained by expanding the EFT predictions for the Higgs signal strength observables to the linear order in , and ignoring the correction of and higher. Put differently, the analysis is performed consistently at order , ignoring all effects (from =8 operators, or from the square of =6 contributions to the observables). This procedure leads to a Gaussian likelihood in the space of , in other words the corresponding function is a quadratic polynomial in . This polynomial can be unambiguously reconstructed given the central values, the 1 uncertainties, and the correlation matrix in Eq. (3.20). 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 is retained. In such a case, the likelihood is highly nonGaussian, but it can nevertheless be well approximated by a Gaussian likelihood in the parameter space region with which is preferred by the LHC Higgs data. The main qualitative consequence of using the full nonGaussian likelihood is the existence of other nearly degenerate minima (in addition to the one described by Eq. (3.20)) where some are large. For example, the full likelihood contains minima with and/or , 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.20) can be recast into a 2HDM likelihood by inserting the relation between and the 2HDM parameters. For example, in order to obtain constraints in the  plane we need to read off from Eq. (2.11)
(3.21) 
Of course this procedure cannot be in any way better than deriving the limits on and directly, as we did in Section 3.1. The purpose of this exercise is 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.20), or a more general one depending on a larger number of EFT parameters that can be subsequently projected into the 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  plane for the typeI and typeII 2HDM obtained in Section 3.1 with the ones deduced by using constraints on EFT parameters together with the EFT2HDM matching discussed earlier.
Let’s start the discussion with the typeII scenario. In this case, a recast of the Gaussian likelihood defined by Eq. (3.20) provides a very good description of the bulk of the favored region where is small. In that region, the LHC Higgs data force the deviations of the Yukawa couplings to be small, less than % of the SM value. Such small deviations can be adequately described by =6 operators of the SM EFT, and the 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 “wrongsign” 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 higherorder terms in in the likelihood function. Such a nonGaussian likelihood is capable of describing multiple minima, including the one in the “wrongsign” region where one or more are smaller than . Indeed, we can see in Fig. 3 that using the nonGaussian EFT likelihood qualitatively captures the shape of the “wrongsign” minimum, at least when is not too large. Using the nonGaussian likelihood also improves the agreement between the EFT and the direct 2HDM limits in the bulk region at small .
For typeX and typeY the results are very similar as for typeII: the EFT description captures very well the bulk of the favored parameter space with small , but it fails to capture the wrongsign Yukawa region. Again, the latter problem can be addressed by using the nonGaussian likelihood on the EFT side.
In the typeI scenario a qualitatively new issue appears. In this case the EFT provides a good approximation of the favored region for low and . However, it completely misses the relevant physics at larger and . Namely, in the typeI 2HDM the LHC Higgs data imply an upper limit on , approximately . At large , this limit is not driven by modifications of the Yukawa couplings, which are suppressed by , but rather by modifications of the Higgs couplings to and . However, these appear only at in the lowenergy EFT of the 2HDM, and are not included at all in the SM EFT truncated at =6. In other words, the typeI 2HDM at belongs to an exceptional class of BSM scenarios that are not adequately described by a SM EFT with =6 operators. Instead, in the Higgs observables, the =8 operators in the lowenergy EFT (formally ) may dominate over the ones (formally ), as the latter are suppressed by and the former are not. This is an example of selection rules in the UV theory modifying the naive power counting in the lowenergy EFT. As a consequence, the =6 EFT approach in this case misrepresents the allowed parameter space of the typeI 2HDM at large . Note that the problem is not addressed by switching from a Gaussian to a nonGaussian EFT likelihood. A more general lowenergy approach is needed to capture this scenario, for example the SM EFT truncated at the level of =8 operators, or a more phenomenological nonEFT approach.
3.3 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 , depending on the 2HDM scenario and on the value of the parameter. However, there exist further important constraints on the 2HDM. First of all, the neutral scalar and pseudoscalar, and the charged partners of the Higgs boson are targeted by direct searches in highenergy colliders. Moreover, the Higgs partners may contribute to electroweak precision observables, in particular to the , , and parameters [45], or to decays [46, 37]. Finally, the parameters of the Higgs potential should satisfy the theoretical constraints following from perturbative unitarity [47, 48, 49], and the Higgs potential should be bounded from below [50]. 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  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 GeV, and generate points with , , TeV. For the typeII scenario we demand in addition GeV, so as to satisfy the constraints [51, 52, 53, 43]. 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 nontrivial, 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 and favored by the Higgs data. This is known as the emmental effect: excluded regions in the multidimensional parameter space of the 2HDM do not show after a projection onto the twodimensional  plane.
For typeII, the bulk of the allowed parameter space with contains points with extra scalar masses ranging from very heavy ( TeV) to very light ( 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 wrongsign region, on the other hand, the extra scalars need to be relatively light, 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. Thus, while the wrongsign region is perfectly consistent with the current experimental data, the related LHC phenomenology is strictly speaking not amenable to an EFT description. We conclude that, for the typeII scenario, the SM EFT approach at is adequate in the entire parameter space allowed by the experimental data and where the additional scalars are heavy. Similar conclusions hold for the typeX and typeY models.
The situation is somewhat different for the typeI scenario. As we discussed in Section 3.2, the discrepancy between the full model and the EFT description is important, especially at large . The reason is that in this case the numerically largest effects on the Higgs boson couplings are and correspond to =8 operators, whereas the formally leading effects, captured by the =6 SM EFT, are suppressed by . This problem will always be present at large enough 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 couplings bounds is dominated by points with very light extra scalars. Thus, much like in the typeII case, most of the 2HDM parameter space where the EFT approach fails to deliver the correct bounds is anyway the one where the extra scalars do not decouple at the LHC energies.
3.4 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 about the effects of heavy new physics on lowenergy 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 =6 operators, which corresponds to retaining the effects up to in the new physics mass scale . If that is done consistently, that is the EFT predictions for the observables are expanded up then lowenergy measurements, such as the LHC Higgs signal strength observables, can be translated into a Gaussian likelihood for the =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 , the large amount of data contained in multiple Higgs analyses at the LHC can be summarized by just a handful of numbers that can be later recast to provide constraints on a large class of BSM scenarios.
The question is how much information about 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 EFT approximating the lowenergy 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:

In the wrongsign Yukawa regions of typeII, X and Y scenarios, where the relative corrections to the Yukawa bottom and/or tau couplings are large and cannot be properly described at .

For the typeI scenario at large , where the leading effects on the single Higgs production and decay come from =8 operators in the EFT, which are by default neglected when the SM EFT is truncated at =6.
One can always complicate the EFT framework such that it is capable of describing also these special cases. In particular, the wrongsign region can be captured in the EFT if one works with the =6 EFT but retains the higher order terms in =6 parameters. Indeed, the Higgs signal strength observables depend also on the squares and higher powers of the =6 EFT Wilson coefficients, which are formally or more suppressed. These are crucial to properly describe the situation when new physics contributions to observables are comparable or exceed the SM ones. As we have shown, retaining these contributions allows one to approximately reproduce the wrongsign regions in the 2HDM, at a price of introducing nonGaussian terms into the likelihood. To cover the large region of the typeI scenario the EFT Lagrangian would have to be extended to include =8 terms.
Both of these complications 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 =6 may be 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.
4 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 production mode appearing in several Higgs decay channels. On the other hand, the signal strength in the 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 Higgstop (bottom) coupling being 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 typeII and typeY models at and . 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 3parameter EFT fit in Eq. (3.20) the SM point where all is not significantly disfavored, with . The reason is that increasing the Higgstop coupling also increases the gluon fusion cross section via the 1loop top triangle diagram contribution to the amplitude. Since the measured total Higgs cross section (which is dominated by gluon fusion) agrees very well with the SM predictions, simply increasing the topHiggs couplings is not preferred by global fits. Decreasing the Higgsbottom coupling is disfavored for similar reasons. As the Higgs width is dominated by decays to bottom quarks, a smaller Higgsbottom 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 channels is outweighed by overshooting the signal strength in the precisely measured , , and final states.
The above discussion suggests a simple adhoc solution to improve the global fit in a theory with two Higgs doublets. One can arrange additional contributions to the effective Higgsgluon 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
(4.1) 
where is the gluon field strength, and is the SM strong coupling. The parameter encodes the effects of heavy colored particles beyond the 2HDM on the Higgs phenomenology. For example, integrating out a new color octet scalar of mass coupled to the Higgs sector via the interaction term , one finds . Similar extensions of the 2HDM have been considered in the past, see e.g. [54, 55, 56].
We first employ the linearized EFT approach to see whether allowing the parameter to vary freely can lead to an improvement of the Higgs fit. The Higgs boson couplings are those in Eq. (2.15) with nonzero and and the remaining coupling set to zero. With that assumption, the Run1 and Run2 Higgs data lead to the following constraints:
(4.2) 
Now the preferred values of the EFT parameters are significantly away from the SM point. Indeed, we find , which translates to the preference for BSM. We also checked that allowing for more free parameters in the EFT (e.g. ) 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 typeII and typeY scenario, once we allow for the new contributions to the Higgsgluon coupling. This is indeed the case, as shown in Fig. 5. The best fit regions move away from the SM limit where to where and are effectively generated.
For the typeII case the best fit point occurs for , , and , and has . The minimum is slightly deeper than in the 4parameter EFT fit because here we use the full (not Gaussian) likelihood function. For such low constraints from flavor physics become nontrivial and require GeV [51, 52, 53, 43]. Nevertheless, this limit does not pose consistency problems, as can be obtained with perturbative couplings in the scalar potential as long as TeV. The preferred value of requires a large contribution to the effective Higgsgluon 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 , thus the octet needs to be below the TeV scale for to remain perturbative. Note that current LHC and Tevatron data still do not exclude fairly light colored particles in a modelindependent way, see e.g. [56, 57] for a recent discussion. On the other hand, the approximate cancellation between all BSM contributions to the gluon fusion amplitude does not have a natural explanation in this model, and should be considered an accident.
For the typeY case the best fit point falls actually into the wrong sign region, at , , and has . However, it is not strongly preferred ( lower by just ) over the local minimum at , , and where all Yukawas are positive. The higher and lower at the local minimum in the typeY case are somewhat easier to accommodate than the best fit point for typeII. For example, in the scalar octet case we need , and the flavor physics bounds on are not relevant for the preferred .
For the case of typeI and typeX models we do not find any significant improvement of the fit after introducing the parameter . This is due to the fact that in these scenarios , therefore one cannot simultaneously fit the enhanced and suppressed signal.
5 Summary
In this paper we discussed the validity of the SM EFT with =6 operators as a lowenergy theory for the 2HDM. Working consistently at in the EFT expansion, the LHC Higgs signal strength measurements can be recasted 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 lowenergy EFT is known. We applied this procedure for the case of the CPconserving 2HDM, restricting to the treelevel matching. We then compared the resulting constraints on the  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 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 =6 provides a valid description of the 2HDM phenomenology, as long as the extra scalars are heavy enough such that they do not appear onshell 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 . One occurs when some SM Yukawa coupling receives corrections that are comparable to its SM value, which happens in particular in the wrongsign Yukawa regions. Another occurs for the typeI scenario at large , where the leading 2HDM effects on Higgs phenomenology are encoded in =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.
Acknowledgments
This work was supported by the Partenariats Hubert Curien programme PESSOA under project 33733UH, and by the Portuguese FCT under project 441.00 of Portugal/France cooperation program PESSOA 2015/2016. A.F is partially supported by the ERC Advanced Grant Higgs@LHC and by the European UnionÕs Horizon 2020 research and innovation programme under the Marie SklodowskaCurie grant agreements No 690575 and No 674896. The work of J.C.R. and J.P.S. is supported in part by the Portuguese Fundação para a Ciência e Tecnologia (FCT) under contract UID/FIS/00777/2013.
References
 [1] S. Weinberg, Effective Gauge Theories, Phys. Lett. B91 (1980) 51–55.
 [2] A. V. Manohar, Effective field theories, Lect. Notes Phys. 479 (1997) 311–362, [hepph/9606222].
 [3] D. B. Kaplan, Five lectures on effective field theory, 2005. nuclth/0510023.
 [4] C. Degrande, N. Greiner, W. Kilian, O. Mattelaer, H. Mebane, T. Stelzer et al., Effective Field Theory: A Modern Approach to Anomalous Couplings, Annals Phys. 335 (2013) 21–32, [1205.4231].
 [5] S. Willenbrock and C. Zhang, Effective Field Theory Beyond the Standard Model, Ann. Rev. Nucl. Part. Sci. 64 (2014) 83–100, [1401.0470].
 [6] E. Masso, An Effective Guide to Beyond the Standard Model Physics, JHEP 1410 (2014) 128, [1406.6376].
 [7] A. Pomarol, Higgs Physics, in 2014 European School of HighEnergy Physics (ESHEP 2014) Garderen, The Netherlands, June 18July 1, 2014, 2014. 1412.4410.
 [8] A. Falkowski, Effective field theory approach to LHC Higgs data, Pramana 87 (2016) 39, [1505.00046].
 [9] The LHC Higgs Cross Section Working Group collaboration, D. de Florian et al., Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector, 1610.07922.
 [10] R. Contino, A. Falkowski, F. Goertz, C. Grojean and F. Riva, On the Validity of the Effective Field Theory Approach to SM Precision Tests, JHEP 07 (2016) 144, [1604.06444].
 [11] M. Gorbahn, J. M. No and V. Sanz, Benchmarks for Higgs Effective Theory: Extended Higgs Sectors, JHEP 10 (2015) 036, [1502.07352].
 [12] J. Brehmer, A. Freitas, D. LopezVal and T. Plehn, Pushing Higgs Effective Theory to its Limits, Phys. Rev. D93 (2016) 075014, [1510.03443].
 [13] T. D. Lee, A Theory of Spontaneous T Violation, Phys. Rev. D8 (1973) 1226–1239.
 [14] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, The Higgs Hunter’s Guide, Front. Phys. 80 (2000) 1–404.
 [15] G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher and J. P. Silva, Theory and phenomenology of twoHiggsdoublet models, Phys. Rept. 516 (2012) 1–102, [1106.0034].
 [16] H. E. Haber and O. Stal, New LHC benchmarks for the conserving twoHiggsdoublet model, Eur. Phys. J. C75 (2015) 491, [1507.04281].
 [17] F. J. Botella and J. P. Silva, Jarlskog  like invariants for theories with scalars and fermions, Phys. Rev. D51 (1995) 3870–3875, [hepph/9411288].
 [18] B. Henning, X. Lu and H. Murayama, How to use the Standard Model effective field theory, JHEP 01 (2016) 023, [1412.1837].
 [19] H. BeluscaMaito, A. Falkowski, F. Duarte, J. RomãoÊ and J. Silva, in preparation, .
 [20] A. Falkowski, Higgs Basis: Proposal for an EFT basis choice for LHC HXSWG, LHCHXSWGINT2015001 (Mar, 2015) .
 [21] A. Freitas, D. LopezVal and T. Plehn, Matching Matters!, 1607.08251.
 [22] ATLAS collaboration, Combined measurements of the Higgs boson production and decay rates in and final states using collision data at 13 TeV in the ATLAS experiment, ATLASCONF2016081 (2016) .
 [23] CMS collaboration, Updated measurements of Higgs boson production in the diphoton decay channel at in pp collisions at CMS., CMSPASHIG16020 (2016) .
 [24] CMS collaboration, Measurements of properties of the Higgs boson and search for an additional resonance in the fourlepton final state at sqrt(s) = 13 TeV, CMSPASHIG16033 (2016) .
 [25] CMS collaboration, First results on Higgs to WW at , CMSPASHIG15003 (2016) .
 [26] ATLAS collaboration, Search for Higgs boson production via weak boson fusion and decaying to in association with a highenergy photon in the ATLAS detector, ATLASCONF2016063 (2016) .
 [27] CMS collaboration, VBF H to bb using the 2015 data sample, CMSPASHIG16003 (2016) .
 [28] ATLAS collaboration, Search for the Standard Model Higgs boson produced in association with a vector boson and decaying to a pair in collisions at 13 TeV using the ATLAS detector, ATLASCONF2016091 (2016) .
 [29] ATLAS collaboration, Search for the Standard Model Higgs boson produced in association with top quarks and decaying into in collisions at = 13 TeV with the ATLAS detector, ATLASCONF2016080 (2016) .
 [30] CMS collaboration, Search for production in the decay channel with pp collisions at the CMS experiment, CMSPASHIG16004 (2016) .
 [31] ATLAS collaboration, Search for Higgs bosons decaying into dimuon in collisions at = 13 TeV with the ATLAS detector, ATLASCONF2016041 (2016) .
 [32] ATLAS collaboration, Search for the Associated Production of a Higgs Boson and a Top Quark Pair in Multilepton Final States with the ATLAS Detector, ATLASCONF2016058 (2016) .
 [33] CMS collaboration, Search for associated production of Higgs bosons and top quarks in multilepton final states at , CMSPASHIG16022 (2016) .
 [34] ATLAS collaboration, G. Aad et al., Measurements of the Higgs boson production and decay rates and coupling strengths using collision data at and TeV in the ATLAS experiment, 1507.04548.
 [35] CMS collaboration, S. Chatrchyan et al., Search for a Higgs boson decaying into a Z and a photon in pp collisions at = 7 and 8 TeV, Phys.Lett. B726 (2013) 587–609, [1307.5515].
 [36] ATLAS, CMS collaboration, G. Aad et al., Measurements of the Higgs boson production and decay rates and constraints on its couplings from a combined ATLAS and CMS analysis of the LHC pp collision data at and 8 TeV, JHEP 08 (2016) 045, [1606.02266].
 [37] P. M. Ferreira, J. F. Gunion, H. E. Haber and R. Santos, Probing wrongsign Yukawa couplings at the LHC and a future linear collider, Phys. Rev. D89 (2014) 115003, [1403.4736].
 [38] D. Fontes, J. C. Romão and J. P. Silva, A reappraisal of the wrongsign coupling and the study of , Phys. Rev. D90 (2014) 015021, [1406.6080].
 [39] P. M. Ferreira, R. Guedes, J. F. Gunion, H. E. Haber, M. O. P. Sampaio and R. Santos, The CPconserving 2HDM after the 8 TeV run, in Proceedings, 22nd International Workshop on DeepInelastic Scattering and Related Subjects (DIS 2014): Warsaw, Poland, April 28May 2, 2014, 2014. 1407.4396.
 [40] I. F. Ginzburg, M. Krawczyk and P. Osland, Resolving SM like scenarios via Higgs boson production at a photon collider. 1. 2HDM versus SM, hepph/0101208.
 [41] D. Carmi, A. Falkowski, E. Kuflik and T. Volansky, Interpreting LHC Higgs Results from Natural New Physics Perspective, JHEP 07 (2012) 136, [1202.3144].
 [42] C.W. Chiang and K. Yagyu, Implications of Higgs boson search data on the twoHiggs doublet models with a softly broken symmetry, JHEP 07 (2013) 160, [1303.0168].
 [43] P. M. Ferreira, R. Guedes, M. O. P. Sampaio and R. Santos, Wrong sign and symmetric limits and nondecoupling in 2HDMs, JHEP 12 (2014) 067, [1409.6723].
 [44] T. Modak, J. C. Romão, S. Sadhukhan, J. P. Silva and R. Srivastava, Constraining wrongsign couplings with , Phys. Rev. D94 (2016) 075017, [1607.07876].
 [45] W. Grimus, L. Lavoura, O. M. Ogreid and P. Osland, The Oblique parameters in multiHiggsdoublet models, Nucl. Phys. B801 (2008) 81–96, [0802.4353].
 [46] H. E. Haber and H. E. Logan, Radiative corrections to the Z b antib vertex and constraints on extended Higgs sectors, Phys. Rev. D62 (2000) 015011, [hepph/9909335].
 [47] S. Kanemura, T. Kubota and E. Takasugi, LeeQuiggThacker bounds for Higgs boson masses in a two doublet model, Phys. Lett. B313 (1993) 155–160, [hepph/9303263].
 [48] A. G. Akeroyd, A. Arhrib and E.M. Naimi, Note on tree level unitarity in the general two Higgs doublet model, Phys. Lett. B490 (2000) 119–124, [hepph/0006035].
 [49] I. F. Ginzburg and I. P. Ivanov, Tree level unitarity constraints in the 2HDM with CP violation, hepph/0312374.
 [50] N. G. Deshpande and E. Ma, Pattern of Symmetry Breaking with Two Higgs Doublets, Phys. Rev. D18 (1978) 2574.
 [51] O. Deschamps, S. DescotesGenon, S. Monteil, V. Niess, S. T’Jampens and V. Tisserand, The Two Higgs Doublet of Type II facing flavour physics data, Phys. Rev. D82 (2010) 073012, [0907.5135].
 [52] T. Hermann, M. Misiak and M. Steinhauser, in the Two Higgs Doublet Model up to NexttoNexttoLeading Order in QCD, JHEP 11 (2012) 036, [1208.2788].
 [53] O. Eberhardt, U. Nierste and M. Wiebusch, Status of the twoHiggsdoublet model of type II, JHEP 07 (2013) 118, [1305.1649].
 [54] A. V. Manohar and M. B. Wise, Flavor changing neutral currents, an extended scalar sector, and the Higgs production rate at the CERN LHC, Phys. Rev. D74 (2006) 035009, [hepph/0606172].
 [55] G. Degrassi and P. Slavich, QCD Corrections in twoHiggsdoublet extensions of the Standard Model with Minimal Flavor Violation, Phys. Rev. D81 (2010) 075001, [1002.1071].
 [56] L. Cheng and G. Valencia, Two Higgs doublet models augmented by a scalar color octet, 1606.01298.
 [57] K. Blum, A. Efrati, C. Frugiuele and Y. Nir, Exotic colored scalars at the LHC, 1610.06582.