New LHC Benchmarks for the CP-conserving Two-Higgs-Doublet Model

We introduce a strategy to study the parameter space of the general, CP-conserving, two-Higgs-doublet Model (2HDM) with a softly broken Z_2-symmetry by means of a new"hybrid"basis. In this basis the input parameters are the measured values of the mass of the observed Standard Model (SM)-like Higgs boson and its coupling strength to vector boson pairs, the mass of the second CP-even Higgs boson, the ratio of neutral Higgs vacuum expectation values, and three additional dimensionless parameters. Using the hybrid basis, we present numerical scans of the 2HDM parameter space where we survey available parameter regions and analyze model constraints. From these results, we define a number of benchmark scenarios that capture different aspects of non-standard Higgs phenomenology that are of interest for future LHC Higgs searches.


Introduction
The observation in 2012 [1,2] of a new boson with a mass close to 125 GeV [3] has been widely viewed as the discovery of the long sought after Higgs boson [4][5][6]. Despite the limited data set, the detailed analyses of the Run 1 LHC data by the ATLAS and CMS Collaborations have confirmed that the observed signal strengths (cross section times branching ratio) of the Higgs boson candidate relative to that expected of the Standard Model (SM) Higgs boson are consistent with SM predictions to within the accuracy of the current measurements [7,8]. Nevertheless, the current precision is at best 20% in the bosonic channels and considerably less accurate in the fermionic channels. Thus, there is still considerable room for new physics beyond the SM to contribute to the properties of the already discovered scalar state.
The Standard Model posits that the dynamics of one complex hypercharge-one, weak doublet of scalar fields is solely responsible for electroweak symmetry breaking. Three of the four degrees of freedom of this scalar doublet provide the longitudinal modes of the massive W ± and Z gauge bosons. The remaining scalar degree of freedom is the SM Higgs boson. But, why should the scalar sector of the SM be of minimal form? The spin-1/2 quark and lepton degrees of freedom of the SM appear in three generations. The origin of the non-trivial flavor structure of the SM is presently unknown. By analogy, one might expect a replication in the scalar sector as well-a non-minimal Higgs sector consisting of multiple doublets. Adding additional doublets yields new phenomena in the scalar sector-charged Higgs scalars, neutral Higgs scalars of opposite CP quantum numbers (if CP is conserved by the scalar potential) or neutral Higgs of indefinite CP (if CP is violated by the scalar potential). The current Higgs data do not rule out the existence of an extended Higgs sector which, as we will discuss in some detail, could be realized at mass scales close to (or in the extreme case even below) 125 GeV. Thus, apart from any other theoretical motivation, it behooves us to devote a dedicated program at the LHC to search for evidence of a non-minimal Higgs structure.
Of course, theoretical arguments have been also advanced in support of a non-minimal Higgs sector. Perhaps the most persuasive is based on the fact that the SM is unnatural [9,10]-namely, it is difficult to understand how the scale of electroweak symmetry breaking (v = 246 GeV) arises in a more fundamental framework that includes gravity and its associated Planck scale, M PL ≃ 10 19 GeV. In the context of the SM coupled to gravity, electroweak symmetry is achieved only by fine-tuning the squared-mass parameter of the scalar potential to an accuracy of 34 decimal places. New physics beyond the SM that attempts to address this issue introduces new phenomena that enters at or near the TeV scale. Many such approaches invoke non-minimal Higgs sectors. The most well studied example of this kind is the minimal supersymmetric extension of the SM (MSSM), which employs a two-Higgs doublet scalar sector with quartic terms in the scalar potential that respect supersymmetry (SUSY) [11][12][13][14].
Since we do not presently know the precise nature of the new physics beyond the SM, it is prudent to be open minded about the possible structure of the non-minimal Higgs sector. However, the observed value of the electroweak parameter ρ = m 2 W /(m 2 Z cos θ W ) ≃ 1 strongly favors non-minimal Higgs sectors comprised only of singlet and doublet scalar fields [15,16]. 1 Among theories in this category, we find the two-Higgs doublet model 2 (2HDM) particularly attractive as the minimal extension that yields new scalar phenomena, including new charged states (H ± ) and neutral states with different (or mixed) CP properties. Direct searches for the additional 2HDM Higgs states based on run-I LHC data have appeared both from ATLAS [21] and CMS [22,23]. Interpretations of the recently discovered Higgs boson at 125 GeV in the context of the 2HDM and implications for future LHC searches have been presented in [24][25][26][27][28][29][30][31][32][33][34].
In this paper, we develop a set of benchmark scenarios for LHC Higgs searches that capture different aspects of 2HDM phenomenology. Our scenarios are devised taking into account that the observed Higgs boson already possess properties that are close to those expected in the SM. The experimental data already provide important constraints on the 2HDM framework. For example, the absence of observable flavor changing neutral currents (FCNCs) mediated by tree-level exchange of neutral Higgs bosons requires that the structure of the Higgs-fermion Yukawa couplings must be of a special form [35,36]. The simplest way to eliminate tree-level Higgs-mediated FCNCs is to impose a particular discrete Z 2 symmetry [cf. Table 1 in Section 2.4] that is exactly respected by the dimensionfour quartic terms of the Higgs Lagrangian [37,38]. 3 We make one further assumption in our analysis by imposing CP symmetry on the scalar potential. This assumption is not required based on data (even though a pure CP-odd state has been ruled out as explanation for the 125 GeV state [39]). The restriction to CP-symmetry provides an additional simplification to the analysis. In particular, the neutral scalar spectrum consists of two CP-even states h and H with m h < m H , and the CP-odd state A. Indeed, the methods developed in this paper are rather easily generalized to a 2HDM with a CP-violating scalar potential or CP-violating vacuum. This possibility will be addressed elsewhere.
In Section 2, we provide the theoretical background relevant to the 2HDM. We introduce the Higgs basis [40,41], which is especially useful in our analysis, as it provides a very clear way of parametrizing the terms that yield deviations of the properties of the observed Higgs boson from those of the SM. Ultimately, we construct a "hybrid basis" of parameters, which we will employ in 2HDM parameter scans. These parameters include the two CP-even neutral Higgs boson masses (m h and m H , where one of these masses is identified with the mass of the observed Higgs boson), the parameter c β−α which parametrizes the deviation of the lightest CP-even neutral Higgs boson h from that of the SM, the ratio of neutral Higgs vacuum expectation values, tan β, and three dimensionless quartic couplings of the Higgs basis. Two of these three dimensionless couplings are related in a simple way to the masses of the charged Higgs and neutral CP-odd Higgs boson of the 2HDM.
In Section 3, we present our numerical scans of the 2HDM parameter space using the hybrid basis of parameters. Based on these scans, which are performed using the code 2HDMC [42,43], we develop seven different benchmark scenarios: (1) a SM-like h, with c β−α small but non-zero as allowed by current data to yield interesting phenomenology of the heavier H; (2) a "flipped" scenario with a SM-like H (which implies that |c β−α | is near 1); (3) overlapping CP-even and CP-odd scalars with masses around 125 GeV; (4) a SM-like h and a heavy non-SM-Higgs mass spectrum with short cascade decays H → ZA or H → W ± H ∓ ; (5) a SM-like h and a heavy non-SM Higgs mass spectrum with long cascade decays H ± → W ± A → W ± ZH or A → W ± H ∓ → W ± W ∓ H; (6) SM-like hV V and htt couplings (V = W ± and Z), but with opposite sign hbb and hτ + τ − couplings as compared with the SM; and (7) an MSSM-like scenario, in which the scalar potential parameters (with one exception) are fixed by the corresponding tree-level MSSM relations. Finally, in Section 4, we present our conclusions. The proposed benchmark scenarios are summarized in a set of tables presented in Appendix A. The two-Higgs doublet model (2HDM) consists of two complex SU(2) doublet, hypercharge-one fields, Φ 1 and Φ 2 , and an SU(2)×U(1)-invariant scalar potential, where the parameters m 2 12 , λ 5 , λ 6 and λ 7 are potentially complex. All other scalar potential parameters are manifestly real. We assume that the parameters are chosen such that the minimum of the scalar potential spontaneously breaks the SU(2)×U(1) electroweak gauge symmetry to U(1) EM [44,45]. That is, at the minimum of the scalar potential, the neutral components of the complex doublet scalar fields acquire vacuum expectation values (vevs), where by convention v 1 and v 2 are real and non-negative. The combination v 2 ≡ v 2 1 + v 2 2 ≃ (246 GeV) 2 is fixed by its relation to the Fermi constant and the W boson mass, v 2 = 1/( √ 2G F ) = 4 m 2 W /g 2 . In the most general 2HDM, the fields Φ 1 and Φ 2 are indistinguishable. Thus, it is always possible to define two orthonormal linear combinations of the two scalar doublet fields without modifying any prediction of the model. Performing such a redefinition of fields (henceforth called a change of basis of the scalar doublet fields) leads to a new scalar potential with the same form as Eq. (1) but with modified coefficients. In this paper we shall focus on the case where the scalar potential and the vacuum state are CP-conserving, leaving the more general case for future work.
The scalar potential is explicitly CP-conserving if and only if there exists a basis choice for the scalar fields in which m 2 12 , λ 5 , λ 6 and λ 7 are simultaneously real. The general conditions that guarantee the existence of such a basis (called a real basis) were developed in [46]. Henceforth, we shall assume that all scalar potential parameters given in Eq. (1) are real. However, it is still possible that the vacuum spontaneously breaks CP. In particular, spontaneous CP-violation takes place if and only if the scalar potential is explicitly CP-conserving, but no real basis exists in which the scalar vacuum expectation values are simultaneously real. Sufficient conditions for a CP-conserving vacuum are easily obtained. The minimization of the scalar potential fixes v 1 ≡ v cos β, v 2 ≡ v sin β (where 0 ≤ β ≤ 1 2 π) and the relative phase of the two vevs, ξ ≡ ξ 2 − ξ 1 through the equations, where s β ≡ sin β and c β ≡ cos β and all scalar potential parameters are real by assumption. If sin ξ = 0, then Eq. (5) can be used to obtain Moreover, this value of ξ corresponds to the minimum [maximum] of the scalar potential if λ 5 > 0 [λ 5 < 0]. Thus, it follows (e.g., see Appendix B of ref. [47]) that no scalar potential minimum occurs for sin ξ = 0 if Thus if Eq. (7) is satisfied, then the scalar potential is minimized for sin ξ = 0 (i.e., ξ = nπ for integer n). One can then perform a hypercharge gauge transformation, Φ i → e −iξ 1 Φ i , followed by the field redefinition, Φ 2 → (−1) n Φ 2 , which yields real non-negative vevs Φ 0 i = v i / √ 2 (for i = 1, 2). If Eq. (7) is not satisfied, then the scalar potential is minimized for sin ξ = 0. Nevertheless, one still must check that there exists no basis transformation of the scalar fields such that the resulting vevs are real. This can be accomplished without explicitly considering all possible basis choices for the scalar fields by making use of the so-called Higgs basis. In Section 2.2, we define the Higgs basis and explain how it can be used to verify that a 2HDM with an explicitly CP-conserving scalar potential also has a CP-conserving vacuum.

The Higgs basis
It is convenient to define new Higgs doublet fields, It follows that H 0 1 = v/ √ 2 and H 0 2 = 0. This is the Higgs basis [40,41], which is uniquely defined up to an overall rephasing, H 2 → e iχ H 2 .
In the Higgs basis, the scalar potential takes the same form as Eq. (1) but with new coefficients, where Y 1 , Y 2 and Z 1 , . . . , Z 4 are real and uniquely defined, whereas Y 3 , Z 5 , Z 6 and Z 7 are complex and transform under the rephasing of H 2 , For an explicitly CP-conserving scalar potential where all the coefficients given in Eq. (1) are real, the real coefficients of the scalar potential in the Higgs basis are given by where s 2β ≡ sin 2β, c 2β ≡ cos 2β, etc. The potentially complex coefficients of the scalar potential in the Higgs basis are given by One can check that by using Eqs. (3)-(5), one recovers the expected scalar potential minimization conditions in the Higgs basis, Note that if sin ξ = 0 then Y 3 , Z 5 , Z 6 and Z 7 are all real and the scalar potential and the vacuum are CP-conserving. More generally, if all the coefficients of the Higgs basis are real for some choice of χ [cf. Eq. (10)], then it follows that any basis related to this real Higgs basis by a real orthogonal transformation of the two doublet fields is also a real basis with real vevs. Thus, it follows that the scalar potential and vacuum are CP-conserving if and only if [41,46,48,49] Im(Z * 5 Z 2 6 ) = Im(Z * 5 Z 2 7 ) = Im(Z * 6 Z 7 ) = 0 .
Thus, we shall parametrize the general CP-invariant 2HDM scalar potential by its real Higgs basis form given in Eq. (9), where all the Y i and Z i are real. The scalar potential minimum conditions given in Eq. (21) fix Y 1 and Y 3 in terms of Z 1 and Z 6 , respectively. Since H 0 1 = v/ √ 2 and H 0 2 = 0, it follows that Y 2 , Z 1 , . . . , Z 5 and the product Z 6 Z 7 are uniquely defined, whereas Z 6 and Z 7 separately change sign under the only possible transformation among real Higgs bases, H 2 → −H 2 . That is, Z 6 and Z 7 are pseudo-invariant quantities with respect to the real Higgs basis transformation H 2 → −H 2 . However, the relative sign of Z 6 and Z 7 is meaningful.
Physical observables must be invariant with respect to any possible Higgs basis transformation. In the most general 2HDM, tan β is also unphysical [41,50] since there is no physical significance to an arbitrary real basis of scalar fields (apart from the Higgs basis and the basis defined in terms of neutral Higgs mass eigenstates). However, it is often the case that the form of the Higgs-fermion Yukawa couplings pick out a special scalar field basis, in which case tan β is promoted to a physical parameter.
It is convenient to employ the real Higgs basis to evaluate the spectrum of Higgs masses in the CP-conserving model. The physical charged Higgs boson is the charged component of the Higgs-basis doublet H 2 , and its mass is given by The three physical neutral Higgs boson mass-eigenstates are determined by diagonalizing a 3 × 3 real symmetric squared-mass matrix that is defined in the Higgs basis [40,50] We immediately identify the CP-odd Higgs boson A = √ 2 Im H 0 2 with squared mass, Note that the real Higgs mass-eigenstate field A is defined up to an overall sign change, which corresponds to the freedom to redefine H 2 → −H 2 . The upper 2 × 2 matrix block given in Eq. (25) is the CP-even Higgs squared-mass matrix, where c β−α ≡ cos(β − α) and s β−α ≡ sin(β − α) are defined in terms of the mixing angle α that diagonalizes the CP-even Higgs squared-mass matrix when expressed in the original basis of scalar fields, {Φ 1 , Φ 2 }. The real Higgs mass-eigenstate fields H and h are defined up to an overall sign change. This implies that β − α is defined modulo π. In particular, It follows that the product s β−α c β−α is a pseudo-invariant quantity with respect to the real Higgs basis transformation H 2 → −H 2 . The squared masses of h and H are then given by, The following identity therefore holds, Hence, diagonalizing M 2 H yields the following expressions: where s β−α ≡ sin(β − α) and c β−α ≡ cos(β − α). Note that Eq. (35) imply that Indeed, Eq. (37) is invariant with respect to the real Higgs basis transformation H 2 → −H 2 , so the sign of the quantity Z 6 s β−α c β−α is physically meaningful. Using the fact that β − α is defined modulo π, we shall establish a convention where In this convention, s β−α is non-negative, the sign of the field h is fixed and c β−α is pseudo-invariant with respect to the real Higgs basis transformation H 2 → −H 2 . One can then derive expressions for c β−α and s β−α from Eqs. (34) and (35), where the signs of the corresponding quantities are fixed by Eqs. (37) and (38): Note that we have used Eq. (33) to obtain the second form for c β−α and s β−α in Eqs. (39) and (40), respectively.

2.3
The CP-conserving 2HDM scalar potential with a softly-broken discrete Z 2 symmetry When we introduce the Higgs-fermion Yukawa couplings in Section 2.4, there will be some motivation to restrict the parameter freedom of the most general scalar potential given in Eq. (1) by requiring the invariance of the scalar potential under the discrete Z 2 symmetry Φ 1 → +Φ 1 and Φ 2 → −Φ 2 . Imposing this discrete symmetry implies that m 2 12 = λ 6 = λ 7 = 0 in Eq. (1). In this case, λ 5 is the only potentially complex scalar potential parameter, which can be rendered real by an appropriate rephasing of Φ 1 . It then follows from Eqs. (6) and (7) that the scalar potential is minimized for sin ξ = 0 if λ 5 < 0 and cos ξ = 0 if λ 5 > 0. In the latter case, Φ 0 2 / Φ 0 1 = ±i tan β. However, a redefinition of Φ 1 → ∓iΦ 0 1 yields real vevs while λ 5 → −λ 5 . Thus, the Z 2 -invariant scalar potential and the vacuum are CP-invariant.
One can relax the discrete symmetry by allowing for m 2 12 = 0 in Eq. (1) which softly breaks the Z 2 symmetry. The quartic terms in Eq. (1) still respect the Z 2 symmetry, so that λ 6 = λ 7 = 0. However, the scalar potential is now CP-violating unless Im(λ * 5 [(m 2 12 ] 2 ) = 0. In what follows, we shall assume that the latter condition is satisfied. In this case, one can rephase Φ 1 such that m 2 12 and λ 5 are both real. That is, the CP-conserving scalar potential of interest is given by where all scalar potential parameters are real. We denote this basis of scalar fields as the Z 2 -basis.
, then the vevs are also real in light of Eq. (7), which implies that the vacuum is also CP-conserving. Otherwise, there is the potential for spontaneous CP-violation. To guarantee that the vacuum is CP-conserving, we check the conditions given by Eq. (22). Setting λ 6 = λ 7 = 0 in Eq. (23) yields The other two conditions given in Eq. (22) simplify considerably when λ 6 = λ 7 = 0, where cos ξ = m 2 12 /(λv 1 v 2 ) if λ 5 > |m 2 12 |/(v 1 v 2 ) and sin ξ = 0 otherwise. Although the vacuum is CP-violating for generic values of the softly broken Z 2 -invariant 2HDM scalar potential when sin ξ = 0, special cases can arise in which sin ξ = 0 and yet the vacuum is CPconserving. We have already encountered one such example when m 2 12 = λ 6 = λ 7 = 0 and λ 5 > 0, in which case the scalar potential is minimized for ξ = 1 2 π. More generally, consider the case of m 2 12 = 0, λ 1 = λ 2 , λ 5 = |λ 1 − λ 3 − λ 4 | and λ 6 = λ 7 = 0. In this case, Eqs. (42)-(44) all vanish, which implies that a real Higgs basis exists. Indeed, for this particular example, Eqs. (18)-(20) yield for λ 5 = ±(λ 1 − λ 3 − λ 4 ). Thus, one can rephase the Higgs basis field H 2 to render Z 5 , Z 6 and Z 7 real. Henceforth, we assume that a Higgs basis has been chosen such that all the Z i are real and sin ξ = 0. The most general parametrization of the CP-conserving 2HDM scalar potential and vacuum is easily specified in the Higgs basis. To determine the constraints on the Higgs basis parameters that guarantee the existence of a softly-broken Z 2 -invariant, CP-conserving scalar potential, one simply sets λ 6 = λ 7 = sin ξ = 0 in Eqs. (11)- (20). It is also convenient to redefine tan β as the ratio of vevs. That is, henceforth we shall define tan β ≡ Φ 0 2 / Φ 0 1 = e iξ v 2 /v 1 = ±v 2 /v 1 for ξ = 0 and ξ = π, respectively. In particular, the Z i are now given by where λ 345 ≡ λ 3 + λ 4 + λ 5 . Since there are five nonzero λ i and seven nonzero Z i , there must be two relations. One can check that the following two identities are satisfied: where Eliminating tan 2β yields the following relation among the Z i , Indeed, in a general CP-conserving 2HDM, if Eq. (54) is satisfied then there must exist a basis in which λ 6 = λ 7 = 0. The corresponding value of tan β (which specifies the basis) can be determined from either Eq. (51) or (52). At this stage, the parameter tan β is a pseudo-invariant, since its sign can be flipped by redefining H 2 → −H 2 [e.g., Φ 2 → −Φ 2 and ξ 2 → ξ 2 + π in Eq. (8)]. It is convenient to work in a convention where the ratio of vevs is non-negative, in which case we can take In this convention, the signs of the pseudo-invariants Z 6 and Z 7 are now fixed, since we can no longer flip their signs by redefining H 2 → −H 2 . This implies that the relative sign of s β−α and c β−α is fixed by Eq. (37), in which case the sign of c β−α is determined by combining Eqs. (37) and (38). Finally, in the CP-conserving, softly-broken Z 2 -invariant 2HDM, it is convenient to introduce the squared-mass parameter, One can express m 2 in terms of Y 2 , Z 1 and Z 6 , Then, combining Eq. (57) with Eqs. (26), (35) and (36) yields

Special forms for the Higgs-fermion Yukawa couplings
We next turn to the Higgs-fermion Yukawa couplings. One starts out initially with a Lagrangian expressed in terms of the scalar doublet fields Φ i (i = 1, 2) and the interaction-eigenstate quark fields. After electroweak symmetry breaking, one can identify the 3 × 3 quark mass matrices. By redefining the left and right-handed quark and lepton fields appropriately, the quark and charged lepton mass matrices are transformed into diagonal form, where the diagonal elements are real and non-negative. The resulting Higgs-fermion interaction Lagrangian in terms of the quark and lepton mass-eigenstate fields, U = (u, c, t), D = (d, s, b), N = (ν e , ν µ , ν τ ), and E = (e, µ, τ ), is given by where K is the CKM quark mixing matrix, h U,D,L are 3 × 3 Yukawa coupling matrices and there is an implicit sum over a = 1, 2. The diagonal quark and charged lepton mass matrices are given by However, the couplings of the neutral Higgs bosons to the fermions are not flavor-diagonal. Thus, Eq. (59) would yield large tree-level Higgs-mediated flavor changing neutral currents (FCNCs) which is in conflict with observed data.
In a general extended Higgs model, tree-level Higgs mediated FCNCs are absent if for some choice of basis of the scalar fields, at most one Higgs multiplet is responsible for providing mass for quarks or leptons of a given electric charge, as first pointed out by Glashow, Weinberg and Pascos (GWP) [35,36]. This GWP condition can be imposed by a symmetry principle, which guarantees that the absence of tree-level Higgs-mediated FCNCs is natural. By an appropriate choice of symmetry transformation laws for the fermions and the Higgs scalars, the resulting Higgs-fermion Yukawa interactions take on the required form in a specific basis of scalar fields. The symmetry also restricts the form of the Higgs scalar potential in the same basis. These considerations were first applied in the 2HDM in Refs. [37] and [38].
The GWP condition can be implemented in four different ways [51-54]: 1. Type-I Yukawa couplings: The four types of Yukawa couplings can be implemented by a discrete symmetry as shown in Table 1.  Table 2: Neutral Higgs-fermion couplings in the 2HDM subject to the Z 2 symmetries given in Table 1. The couplings listed above are normalized relative to the SM Higgs couplings to U U , DD, and EE.
The neutral Higgs Yukawa couplings (relative to the corresponding couplings of the SM Higgs boson) are conveniently summarized in Table 2 for the four possible implementations of the GWP condition. The imposition of the discrete symmetry also restricts the form of the Higgs scalar potential given in Eq. (1) by setting m 2 12 = λ 6 = λ 7 = 0. As discussed in Section 2.3, the condition m 2 12 = 0 can be relaxed. In the case of a softly-broken Z 2 -invariant 2HDM with m 2 12 = 0, Higgs-mediated FCNCs are still absent at tree-level, although they are generated at one-loop order. Nevertheless, the size of these FCNCs may be phenomenologically acceptable, depending on the region of the 2HDM parameter space. This motivates us to focus on the Higgs scalar potential of the form given in Eq. (41). Note that the parameter tan β ≡ Φ 0 2 / Φ 0 1 is defined in terms of the vevs with respect to the Z 2 -basis of scalar fields, where the discrete Z 2 symmetry of the Higgs-fermion Yukawa interactions is manifest (i.e., where three of the six Higgs-fermion Yukawa matrices vanish). Indeed as previously advertised, the parameter tan β has been promoted to a physical parameter of the theory.

Parametrizing the softly-broken Z 2 -symmetric CP-conserving 2HDM
The scalar potential in the Z 2 -basis [Eq. (41)] is governed by three squared-mass parameters and five dimensionless quartic coupling parameters. Minimizing the scalar potential, we re-express m 2 11 and m 2 22 in terms of v = 246 GeV and tan β. Excluding v (which determines the W and Z masses), we are left with seven real parameters: m 2 12 , tan β and the λ i (i = 1, 2, . . . , 5). In the Higgs basis, the counting is also straightforward: after imposing the scalar potential minimum conditions the relevant real parameters (excluding v) are Y 2 and Z i (i = 1, 2, . . . , 7). Imposing one relation, Eq. (54), to guarantee the existence of a basis where λ 6 = λ 7 = 0, we are again left with seven real parameters.
A more physical choice of parameters would consist of α, β and the four scalar masses, m h , m H , m A and m H ± . This leaves one additional parameter, which is usually chosen to be m 2 12 or λ 5 (cf. Appendix D of [47]). In many of the previous studies of the 2HDM parameter space [24-26, 29-31, 33, 34], scans were performed over the parameters m h , m H , m A , m H ± , m 2 12 , α and β. Acceptable points in the scan must satisfy unitarity and perturbativity constraints [55,56]. However, for random choices of m 2 12 and masses of the three non-SM-like Higgs bosons, one finds that the unitarity and perturbativity constraints are often violated. This is easily understood by examining the decoupling limit where Y 2 ≫ v. In this limit, m H ≫ v in which case c β−α → 0 [cf. Eq. (39)] and the properties of h approach that of the SM Higgs boson. Moreover, eqs. (24), (26), (36) and (56) yield squared-mass differences, The conditions of unitarity and perturbativity constrain these squared-mass differences to be of O(cv 2 ), where c < ∼ 10. To maximize the efficiency of costly scans over the 2HDM parameter space (e.g. in a global fit), one should define the input parameters of the model to include the mass of the observed SM-like Higgs boson and at most one additional Higgs mass. The other two Higgs masses and m 2 are determined in terms of the quartic Higgs self-coupling parameters. As long as the magnitudes of the quartic parameters are constrained to be less than O(10), the scan over the other Higgs masses and m 2 will not violate the unitarity and perturbativity constraints.

The hybrid basis of parameters
In light of the discovery of a SM-like Higgs boson, we have information on two of the parameters of the 2HDM. Since the observed boson is not CP-odd, we shall identify it with either h or H. Given that the hV V and HV V couplings (V = W ± or Z) relative to the SM Higgs boson h SM are given by the values of the SM-like Higgs boson mass (m h or m H ) and c β−α are already known with some precision. We therefore propose a "hybrid" strategy for specifying the input parameters for the softly-broken Z 2 -invariant 2HDM. The masses of the physical CP-even Higgs scalars are given directly as input parameters together with c β−α , which determines the phenomenologically important couplings of the CP-even scalars to the W ± and Z bosons [cf. Eq. (61)], and tan β which specifies the basis of scalar fields where the discrete symmetry of the Higgs-fermion Yukawa interactions is manifest. In addition, we specify the real Higgs basis self-coupling coefficients Z 4 , Z 5 , and Z 7 as input parameters. We henceforth designate the input parameter set {m h , m H , c β−α , tan β, Z 4 , Z 5 , Z 7 } as the hybrid basis of parameters as indicated in Table 3. In the convention adopted in Eq. (38), we have 4  Table 3: 2HDM input parameters in the hybrid basis. By convention, we take tan β to be non-negative. For a more detailed description of their phenomenological relevance, see the text. 4 Different convention choices for β − α appear in the 2HDM literature. For example, the convention 0 ≤ c β−α ≤ 1 and −1 ≤ s β−α ≤ 1 is employed internally by the 2HDMC code [42,43]. Another common choice in the literature is to take − 1 2 π ≤ α ≤ 1 2 π. Indeed, one is always free to define either α or β − α modulo π. It is a simple matter to translate among the various conventions. With these seven input parameters, one may compute the real Higgs basis self-coupling coefficients Z 1 , Z 2 , Z 3 and Z 6 as well as the CP-odd Higgs and charged Higgs scalar masses m A and m H ± as follows. First, we employ Eqs. (34) and (35), In order to ensure that unitarity and perturbativity constraints are not violated, one should choose the input parameters m H and c β−α such that m 2 H c β−α ∼ O(cv 2 ), where c < ∼ O(10) in order to avoid a potentially large value of Z 6 . Next, we use Eq. (51) to compute Z 2 , There is a danger that Z 2 may become too large if β is near 0 or 1 2 π. This is not a true singularity since in the formal limit of tan 2β = 0, it follows that Z 6 = Z 7 = 0 [cf. Eqs. (46) and (47)] corresponding to the inert 2HDM limit discussed briefly at the end of this Section. In practice, we restrict our scans over regions of tan β such that Z 2 is never too large.
The parameter Z 3 is determined from Eq. (52), For this quantity, values of β near 0, 1 4 π or 1 2 π appear to be problematical, yielding potentially large values of Z 3 . The case of β near 0 or 1 2 π has already been mentioned below Eq. (65). The case of β = 1 4 π is not a singular limit, as it corresponds to Z 1 = Z 2 and Z 6 = Z 7 with Z 345 an independent parameter [cf. Eqs. (46)- (50)]. In particular, for β = 1 4 π, we must replace Z 7 with Z 3 as an input parameter in the hybrid basis. We can sidestep this special parameter regime by avoiding values of tan β too close to 1.
In practice, taking the parameters Z 4 and Z 5 as input parameters is not very intuitive. Clearly a more physical approach is to take m A and m H ± as input parameter and then compute Z 4 and Z 5 using Eqs. (67) and (68). Then these parameters can be employed in evaluating Z 3 via Eq. (66). The only danger with such an approach is that a poor choice of m A and m H ± will yield values of Z 4 and Z 5 that violate unitarity and perturbativity constraints, as discussed in Section 2.5. In practice, we first employ our hybrid basis to perform our parameter scans. Once we have identified suitable parameter regimes, we will then take a more physical approach by employing m A and m H ± to define the corresponding benchmark scenarios.
The parameter Z 7 is also not particularly intuitive, so one could advocate taking m 2 as the input parameter and then compute Z 7 using Eq. (69). However, m 2 is not a physical squared-mass parameter, so there is no particular advantage for adopting this strategy. A poor choice of m 2 would yield a value of Z 7 that violates unitarity and perturbativity constraints. Indeed, the CP-conserving, softly-broken Z 2 -invariant 2HDM always requires one extra parameter beyond the four physical Higgs masses, α and β. There is an advantage to choosing this parameter to be dimensionless for the reasons noted above. Another possible choice for this parameter that is often found in the literature is λ 5 . Using Eqs. (48)-(50), one can check that so one can always trade in λ 5 for Z 7 and vice versa. The case of m 2 = 0 is special and corresponds to a 2HDM with an exact Z 2 symmetry. In this case, Z 7 is determined by Eq. (69) in terms of m h , m H and c β−α . It is noteworthy that one cannot take m H arbitrarily large in the Z 2 -invariant 2HDM without violating unitarity and perturbativity bounds. 5 Thus, there is no consistent decoupling limit where the heavy Higgs states H, A and H ± decouple from the low-energy theory. In contrast, in the softly-broken Z 2 -invariant 2HDM with m 2 = 0, the decoupling limit exists in which H, A and H ± are all very heavy with squared-masses of O(m 2 ) [cf. Eq. (57)], while keeping all Higgs self-coupling parameters bounded.
Finally, the inert 2HDM [57][58][59], corresponding to Z 6 = Z 7 = 0 and Type-I Yukawa couplings [20], must be treated separately. In the inert 2HDM, we have either β = 0 or β = 1 2 π and either c β−α = 0 (in which case h is identified as the SM-like Higgs boson) or s β−α = 0 (in which case H is identified as the SM-like Higgs boson). In both cases, the SM-like Higgs boson is denoted by h SM and the other inert neutral Higgs bosons will be denoted by H I and A I . The states H I and A I are relatively CP-odd, although (despite the notation) the individual CP-quantum numbers of these two states are not well defined. By convention, we shall define H I to be the heavier of the two neutral inert Higgs scalars. Then, using Eqs. (24) and (25), the charged and neutral inert Higgs masses are given by The relevant input parameter set in this case is Note that parameters Z 2 and Z 3 , which govern the self-couplings of the inert scalar doublet, must be provided separately since Eqs. (65) and (66) are not applicable in this case. We shall not consider the inert 2HDM further here, but phenomenological constraints and benchmark scenarios for this model have been recently discussed in [60].

Numerical analysis
To capture the 2HDM phenomenology that is interesting for upcoming LHC Higgs searches, we have performed numerical scans over the softly broken Z 2 -symmetric 2HDM parameter space using the hybrid basis of parameters for inputs for the scans. 6 Based on these results, we identify a set of useful benchmark scenarios that we will now describe in more detail. Most of the scenarios are presented in the form of 2-dimensional benchmark planes, whereas some of them are in the form of lines where only one parameter is varied. As a general rule, these scenarios could be made into benchmarks of higher dimensions by promoting additional parameters with fixed values to vary within certain ranges. Below we discuss a few examples of how this could be done.
For the numerical evaluations we use the code 2HDMC [42,43] (v. 1.7.0), where the hybrid basis of parameters has been implemented according to Section 2.6. Constraints on the quartic couplings from (absolute) vacuum stability and S-matrix unitarity are evaluated at the input scale. For the case with Z 2 symmetry (as is implicit already in the definition of the hybrid basis), the condition that the Higgs potential is positive definite is equivalent to the well-known relations [57], For the unitarity constraints, we impose an upper limit which corresponds to saturation of the unitarity bound by the tree-level contribution. Using the formulation of [56], this is equivalent to the constraint |Λ| < 16 π on individual matrix elements. As an example, the regions in (c β−α , tan β ) compatible with theoretical constraints from perturbative unitarity and positivity of the potential are shown in Branching ratios and LHC cross sections for Higgs production are evaluated using 2HDMC and SusHi [61] in accordance with the recommendations of [62]. We recommend that numerical values used in an experimental analysis are also calculated with HIGLU [63] and HDECAY [64], since such a comparison can provide a first estimate of the theoretical uncertainties in the treatment of (missing) higher-order corrections. These programs provide cross section predictions for the heavier 2HDM Higgs bosons, H and A, in the dominant gluon fusion and bb associated production modes. One of the benchmark scenarios discussed below (Scenario E) allows for the production of heavy charged Higgs bosons. The cross section for pp → tbH − /tbH + depends only on m H ± and tan β. A detailed numerical analysis of this process has recently been presented in [65], with numbers that are applicable to our Scenario E.
Constraints from direct Higgs searches at LEP, the Tevatron, and the LHC are evaluated using HiggsBounds [66][67][68] (v. 4.2.0), which selects for each parameter point the most sensitive exclusion limit (at 95% C.L.). The compatibility of the 2HDM with the observed 125 GeV Higgs signal is calculated in terms of a χ 2 value taking into account the full LHC run-I data using HiggsSignals [69] (v. 1.3.0). To determine the viable parameter regions in a benchmark scenario with two free parameters, we demand compatibility with the best-fit point, which usually is in very good agreement with the SM, within 2 σ (∆χ 2 = χ 2 − χ 2 min < 6.18). The best-fit (minimal χ 2 value) is reevaluated for each benchmark scenario. We make no quantitative statements about the relative degree of compatibility between different scenarios and the data, but it is checked explicitly that the best fit points obtained in all cases lie very close to the SM predictions. Fig. 2 demonstrates how these constraints work on the 2HDM parameter space under the assumption that the observed LHC signal corresponds to the lightest 2HDM Higgs boson, h. Demanding that the Higgs rates are within 2 σ of the measurements restricts c β−α to be close to the SM limit (c β−α ∼ 0) in the case of Type-II Yukawa couplings. For Type-I couplings, the deviation from the SM limit can be somewhat larger. The shape of the allowed region arises from the dependence of the production and decay rates on the mixing angles, most importantly the total width which in the Type-II model is dominated by the hbb coupling.
One additional set of constraints to keep in mind when designing viable benchmark scenarios are the electroweak precision tests, where in particular the oblique T parameter [70] can receive sizable contributions due to a large mass splitting between the non-SM-like charged and neutral Higgs states [71][72][73][74][75]. In particular, where α ≡ e 2 /(4π) with the electromagnetic coupling e and the weak couplings g and g ′ defined in the MS scheme evaluated at m Z . The function F is defined by with F(m 2 , m 2 ) = 0. Typically, the O(g ′ 2 ) contribution to T is numerically small. Hence, to ensure that the 2HDM contribution to T is within the current bound, T ≤ 0.2 at 95% C.L. [76], one must be close to the custodial limit. In particular, in the decoupling limit of the 2HDM, Eq. (73) reduces to [75] αT from which one can see that the 2HDM contribution to the T parameter will be sufficiently small if either m 2 Consequently, the choice Z 4 = ±Z 5 is special in the sense that the corresponding Higgs sector contributions to the T parameter vanish in the decoupling limit. Several of the benchmark scenarios that we propose below satisfy this property. Finally, the general 2HDM (without supersymmetry) can also be constrained by various low-energy (flavor physics) processes (see, e.g., [77]). Since our main objective is to define scenarios capturing interesting LHC phenomenology, we will not be concerned with the details of these constraints, nor will they be explicitly applied in our numerical results. In all the following, it should be noted that there exist generic lower bound on the charged Higgs boson mass in the Type-II 2HDM, m H ± 480 GeV at 95% C.L., from measurements of the BR(B → X s γ) [78].

Scenario A (non-alignment)
Our first benchmark scenario has the "normal" interpretation of the 125 GeV signal as the lightest CP-even Higgs boson, h, with SM-like properties. The h is SM-like in the so-called alignment limit 7 where c β−α → 0, in which case the hV V coupling approaches the corresponding SM value. On the other hand, to allow for some interesting phenomenology of the heavier CP even state, H, we define the scenario with a non-alignment (c β−α = 0) as allowed by the present constraints [cf. Fig. 2]. The scenario focuses on searches for the heavier CP-even state, H, in SM final states (including the H → hh decay). The remaining two Higgs bosons, A and H ± (which are kept mass-degenerate), are decoupled to a sufficient degree to create a small hierarchy For m H > 150 GeV, this can be achieved by setting Z 4 = Z 5 = −2, which leads to values of m H ± satisfying the b → sγ constraint for Type-II models. The value of c β−α is fixed close to the maximum allowed by the LHC Higgs constraints, c β−α = 0.1 for Type-I and c β−α = 0.01 for Type-II couplings. Consequently, we keep m H and tan β as free parameters. As shown in Fig. 3, these parameter choices  We now look at the predictions for the production and decay modes of the heavy CP-even Higgs in this scenario. Fig. 4 shows contours of the cross section for the dominant gluon fusion production mode, gg → H versus the free parameters (m H , tan β ). As can be seen from this figure, the cross section is maximized as low values of tan β due to the suppression of the Htt coupling by 1/ tan β . This holds for both Type-I and Type-II couplings due to the universality of the couplings to up-type quarks. In the Type-II case, we also see an enhancement of the cross section at high values of tan β where the bottom loops become dominant, although in this scenario the parameter space regions with a significant enhancement is already disfavored by LHC constraints. In the Type-I case, the decoupling property of the H couplings to fermions sets in at high tan β , which can lead to arbitrarily small values of the cross section. This scenario can therefore provide a useful benchmark to gauge and interpret experimental progress in probing smaller and smaller values of σ × BR in particular channels. In the Type-II model, at high tan β , the cross section for associated production with b quarks, bb → H, can become important since it scales as tan 2 β, whereas for Type-I couplings this process is suppressed.
Turning now to the decay modes, the case of Type-I Yukawa couplings is shown in Fig. 5 for tan β = 1.5 and 7. The Type-I couplings gives a very distinct pattern of dominant decay modes, where the vector bosons dominate below 250 GeV. Above this value, the cascade H → hh opens up and quickly becomes dominant, while there are still substantial contribution from the decays to W W and ZZ. At low tan β this changes above the tt threshold, but for higher values of tan β this mode (like all fermionic decays) is suppressed. Scenario A with Type-I Yukawa couplings is a suitable benchmark to replace the SM as the standard parametrization of 2ℓ and 4ℓ searches at higher masses, where the rates can basically be adjusted to anything below the present limit. Note that arbitrarily high values of m H are not allowed in this scenario for all tan β values, owing to the restrictive unitarity and stability constraints when c β−α is non-zero [cf. Fig. 1].  With Type-II couplings, the restriction on m H from theoretical constraints is much less severe, since the LHC data already forces the scenario to be defined with c β−α much closer to the alignment limit. This also has the effect of suppressing bosonic decay modes, which can be seen in Fig. 6 where fermionic modes (bb, τ τ and tt) are clearly dominant over the full mass range. This scenario can be used for combinations of searchers utilizing different fermionic decay modes. It is also useful to gauge performance of e.g. τ τ searches in a setting where only a single CP-even Higgs boson is present at a given mass (in contrast to the MSSM case where typically H and A are nearly mass-degenerate). It also avoids other complications associated with SUSY interpretations, such as the non-holomorphic ∆ b corrections to the bottom Yukawa coupling and/or the presence of SUSY decay modes [82,83].

Scenario B (low-m H )
Scenario B corresponds to a "flipped" 2HDM benchmark scenario. In this scenario both h and H are light, but it is the heavier of the two which has m H = 125 GeV and is SM-like. Since m h < m H , the lighter of the two CP-even Higgs bosons must have strongly suppressed couplings to vector bosons to be compatible with direct search limits which forces s β−α → 0 [cf. Eq. (61)]. 8 Similarly to Scenario A, the hybrid basis quartic parameters are chosen to decouple the other two Higgs states, A and H ± , to a sufficient degree not to affect the phenomenology (here: Z 4 = Z 5 = −5). This is to ensure compatibility with the strong indirect constraints (from flavor physics) on a light charged Higgs boson with mass of order m H , in particular for Type-II Yukawa couplings. In the parameter space region 90 < m h < 120 GeV LHC constraints (from h → bb, τ τ ) apply, which leads to an upper limit on tan β. This limit depends only weakly on c β−α in the limit c β−α ≪ 1. We focus here on the benchmark line with tan β = 1.5 as a function of m h , although it is straightforward to generalize this analysis to a benchmark plane by varying tan β. The colors indicate statistical compatibility with the 125 GeV signal at 1 σ (green), 2 σ (yellow) and 3 σ (blue). The gray (dashed) region is excluded at 95% C.L. by constraints from direct searches at LEP and the LHC.
Since this scenario implements a completely different interpretation of the 125 GeV signal, the assumptions for the parameter space fit presented in Fig. 2 are not fulfilled. In Fig. 7 we therefore show the corresponding results with the assumption of m H ≃ 125 GeV for a region of parameter space near the SM limit, |c β−α | → 1. From this figure it can also be seen that the region where H → hh is open is very hard to reconcile with LHC measurements, we will therefore not consider this possibility any further.
As suitable benchmark scenarios we select the case of exact alignment, c β−α = 1, with either Type-I or Type-II Yukawa couplings. As a third scenario, we choose a non-aligned value of c β−α = 0.9, which however is only relevant for Type-I couplings. The dominant decay modes of the light, non-SM-like, Higgs bosons are similar for all the three different scenarios, with BR(h → bb) ∼ 75-80% and BR(h → τ + τ − ) ∼ 8% the two most interesting from the phenomenological point of view. 8 In light of Eq. (40), the alignment limit of s β−α → 0 is achieved in the limit of Z6 → 0, in which case the heavier of the two CP-even Higgs scalars is approximately aligned with the Higgs basis field, H ≃ √ 2 Re H 0 1 − v, and mH ≃ Z1v 2 [cf. footnote 7]. In contrast to Scenario A, the decoupling limit can never be realized in Scenario B since at least one of the Higgs states is lighter than the SM-like Higgs boson.

Scenario C (CP-overlap)
In this work we have restricted ourselves to benchmarks for a 2HDM Higgs sector with CP-conservation, while a detailed analysis of the hybrid basis extended to the CP-violating case is postponed to future work. Nevertheless, we now consider a scenario where overlapping CP-odd and CP-even Higgs bosons simultaneously have mass close to 125 GeV [84]. Since the CP-odd Higgs boson does not couple to vector bosons at tree level, there are surprisingly few channels where it is possible to distinguish this scenario from the case with a single light Higgs, h. The most important channel where the CP-odd contribution to the total rate could reach O(1) is through gluon (bb) fusion, followed by the decay h/A → τ + τ − . We shall therefore analyze this process in more detail.
Requiring that m A = m h (= 125 GeV) in the hybrid basis of parameters, Eq. (67) yields Requiring in addition that m H ± = m H , Eq. (68) implies that The remaining quartic parameter, Z 7 , does not enter in the mass determination. To maintain a SMlike h, we focus on the alignment limit and set c β−α = 0. In this case, Eqs. (77) and (78) reduce to We fix Z 7 , which has only minor impact on the phenomenology, such that m 2 12 = m 2 A s β c β , or equivalently λ 5 = 0. Setting m H = 300 GeV, this leaves tan β as the only remaining free parameter in Scenario C. Other choices for m H (or even varying its value continuously) would lead to a benchmark plane generalization of Scenario C. However, varying m H has no impact on the properties of the overlapping CP states (h and A).
Scenario C can be considered with both Type-I and Type-II Yukawa couplings. In Fig. 8, we show the production cross sections in the hadronic modes gg → h/A and bb → h/A (and their sum) for Type-I couplings. In this case the CP-even production dominates, except at low values of tan β 1.6. Since there is no tan β enhancement of the bottom Yukawa, the bb-induced production remains small and can be neglected. Combining the (total) production cross sections (gg + bb) with the branching ratios BR(h/A → τ τ ), we form the inclusive quantity for the τ τ rate relative to the SM. Predictions for R τ τ in Scenario C are shown in Fig. 9. The figure shows that the total rate (green curve) approaches the SM value (R τ τ = 1) in the limit of high tan β due to the decoupling property for the A contribution with Type-I couplings. In the opposite limit, the CPodd contribution becomes increasingly important and dominates for tan β 1.6. In effect, the total predicted τ τ rate significantly exceeds the SM prediction. Experimental constraints limit the maximal total rate, with the currently measured rates being µ ATLAS τ τ = 2.1 +0. 9 −0.8 [85] and µ CMS τ τ = 0.34±1.09 [86]. 9  To determine the fractions of the total rate which can be due to the CP-even and CP-odd 2HDM components, we show in the right panel of Fig. 9 their relative contributions R i τ τ /R tot τ τ . Due to the decrease in the cross section, the CP-odd component can be seen to decrease monotonously with increasing tan β. This scenario therefore provides a consistent model to parametrize an arbitrarily small CP-odd admixture in the 125 GeV signal.
With Type-II Yukawa couplings, the production cross sections for the CP-odd state, A, have a different dependence on tan β . The effect of this change is shown in Fig. 10. As can be seen from this figure the CP-odd contribution now becomes much more important at high tan β , both due to an increase in gluon fusion which has a minimum around tan β ∼ 4 and due to a contribution from the bb → A process with a cross section increasing as tan 2 β. In the case of Type-II couplings the  The total τ τ rate (adding gg and bb production modes), relative to the SM, from h (long dashes), A (short dashes) and their sum (green, solid) in Scenario C with Type-II Yukawa couplings. Right: the respective fractions of the inclusive τ τ rate resulting from h (long dashes) and A (short dashes). total cross section does not approach the SM value for any value of tan β . This also has consequences for the total τ τ rate, as can be seen in Fig. 11. In this scenario the minimum rate, obtained for tan β ∼ 3.5, is R τ τ ≃ 1.5. While 50% above the SM expectation, this is not experimentally excluded. Interestingly, the similar magnitude of the CP-even and CP-odd contributions leads to fractions that are nearly equal over a large range in tan β (Fig. 11, right). This scenario is therefore particularly suitable as a "best-case" test of the possibility to distinguish CP-properties of the 125 GeV Higgs boson in the τ τ channel.

Scenario D (Short cascade)
This scenario is constructed with a SM-like h by fixing c β−α to be zero (exact alignment). The mass hierarchy can be modified allow for either one-or both-of the decay modes H → ZA or H → W ± H ∓ to be open, thus resulting in a "small cascade" of Higgs-to-Higgs decays. 10 These decays can typically be made dominant in the mass window 250 GeV < m H < 350 GeV (below tt threshold). A recent search for H → ZA by the CMS Collaboration [23] already places interesting constraints on the parameter space of this type of scenario. Other modes that can be potentially of simultaneous interest is H → hh and H → AA (when A is very light).
We present realizations of Scenario D for all the interesting cases below. We define these scenarios with a single free parameter by fixing tan β = 2, but this assumption could easily be relaxed. For simplicity, we choose two of the three non-SM-like Higgs masses to be equal. In the hybrid basis of parameters, these mass degeneracies can be implemented by an appropriate choice of Z 4 and Z 5 as follows [cf. Eqs. (67) and (68)]: In addition to the short cascade that we are investigating, additional "exotic" decay modes may be accessible, such as e.g. A → W ± H ∓ (when H ∓ is light). However, with the possibility of nondegeneracy between the heavy Higgs bosons (H and A in this case) there is in principle no guarantee for these other exotic decay modes to co-exist at any appreciable rate. We therefore choose to discuss only the H decays here, but there could be additional modes of similar type to exploit.
Starting with the case of only a low m A , this can be realized in the hybrid basis of parameters by choosing Z 4 = −1, Z 5 = 1. For m H close to 250 GeV, the decay H → AA can be open, with a rate that can be adjusted by varying Z 7 . Here we choose Z 7 = −1, which satisfies stability requirements in the whole mass range. 11 Fig. 12 shows the branching ratios of the heavy CP-even Higgs boson, H, into the interesting final states. As this figures shows, there is nearly no difference in the decay pattern for the two Yukawa types (although, of course, the production cross sections are different). One thing that can be noteworthy is that the cascade decay, in this case H → ZA maintains an appreciable branching ratio also beyond the top threshold, with percent-level rates up to m H ∼ 380 GeV. In Type-II models it is possible to suppress H → tt further by going to higher tan β , but we find for the chosen set of parameters that this requires adjusting Z 7 to maintain positivity of the potential.
An alternative mass hierarchy with m H ± < m A = m H can be arranged by setting Z 4 = 2 and Z 5 = 0 while keeping the remaining parameters fixed. This generates the new possible decay modes H → W ± H ∓ (and H → H + H − for very low m H ± ). Note that these low m H ± might be in a region which is disallowed by flavor constraints. The decay branching ratios in this scenario are shown in Fig. 13. Similarly to the case with light A, the cascade decay H → W ± H ∓ is dominant over the range 250 < m H < 350 GeV, and with a sizable tail towards higher masses. In this case percent-level branching ratios are obtained for m H ≃ 420 GeV. Note that the widths for the leading-order predictions for the decays involving one Higgs boson and one gauge bosons, e.g. H → ZA, are proportional to gauge couplings (and therefore relatively fixed), whereas the modes involving triple-Higgs couplings, such as H → AA, are proportional to scalar couplings appearing in the Higgs potential. For example, 12 The strength of such modes is therefore more sensitive to the specific choice of the 2HDM parameters.

Scenario E (Long cascade)
The short cascade (Scenario D) was defined with a one-step decay involving two (or more) Higgs bosons. In Scenario E, we extend this to a "long cascade" where two-step decays involving all three non-SM-like Higgs bosons are possible, assuming that H is always the lighter of the three. Assuming first (E1) that H ± is heavier than both A and H (realised in the hybrid basis by choosing Z 4 = −6, Z 5 = −2 (with Z 7 always zero in this scenario), can give rise to a "long" cascade and a second complementary direct decay that is always present [88], The other hierarchy, with m A > m H ± > m H is achieved by setting Z 4 = 1, Z 5 = −3 (E2). This leads to long cascades which are inverted compared to Eq. (85) [87]: Interesting signatures in this scenario include a heavy Higgs boson H in the τ + τ − or bb channels augmented with multiple leptons and/or additional jets. As in scenario D, the recent search for A → ZH [23] places constraints on the parameter space. It should be interesting to analyze the interplay between this mode and other available signatures in this scenario in more detail. In Table 4 we present the branching ratios for the various cascade decay modes in the two different incarnations of Scenario E for input masses m H = 200 GeV and m H = 300 GeV with Type-I Yukawa couplings. Note that for this low value of tan β there is no significant change in any of these branching ratios in a Type-II setting. The main difference between the two types are instead in the decay of H at the end of the cascade, where the final states with down-type fermions (b, τ ) are more favored with Type-II couplings.  Table 4: Mass spectrum and Type-I Higgs branching ratios to bosonic decay modes in Scenario E.
There are several things of generic interest that can be noted from Table 4. The heaviest (parent) Higgs boson has a rather high probability, 70-85%, to decay into a lighter Higgs and a vector boson. However, most of this decay goes directly into the lighter of the two states (H) due to the larger available phase space. The aggregated branching ratio for a "long" cascade is thus suppressed, reaching typical values at the (few) percent-level (up to 5-6% in some cases). Nevertheless, to be able to study these types of final states would be an intriguing possibility, that would contain a lot of information about the structure of an extended Higgs sector.

Scenario F (Flipped Yukawa)
The flipped Yukawa scenario is characterized by SM-like couplings for the light Higgs, h, except for the couplings to down-type fermions which has a change of sign relative to the SM. This scenario is realized with Type-II Yukawa couplings for values of (c β−α , tan β ) solving the equation [89] g hdd g SM This unconventional solution is what generates the second branch of 2 σ-allowed parameter space in Fig. 2 (right), with valid solution for not too small values of tan β. We define Scenario F by simply fixing the value of c β−α from the chosen tan β using this relation, while the remaining parameters are given the same values as in Scenario A: Z 4 = Z 5 = −2, Z 7 = 0. The (small) effects on the h phenomenology induced by the opposite sign bottom Yukawa coupling have been studied in detail in [89]. However, the comparably large allowed values of c β−α also open possibilities for the decay modes of the heavier Higgs boson beyond what can be achieved near the alignment limit as considered, for example, in Scenario A. To illustrate this point, we show in Fig. 15 the decay branching ratios of H as function of m H and tan β . As this figure shows, the decays into vector bosons can be sizable over a large part of the remaining parameter space, clearly beyond what is possible in the alignment case. This observation holds even in the region above the top threshold, which makes the Flipped Yukawa scenario a suitable benchmark to replace the SM in interpretations of future heavy Higgs searches in the dilepton and 4ℓ final states over a larger mass range.

Scenario G (MSSM-like)
This scenario is inspired by the Higgs potential of the MSSM. The tree-level MSSM is defined by the following values for the quartic couplings (in the "SUSY-basis" where supersymmetry is manifest), Following Eq. (56), the remaining parameter is given by m 2 12 = 1 2 m 2 A sin 2β in terms of the more usual MSSM inputs: the CP-odd Higgs mass, m A , and tan β. In the tree-level MSSM, the predicted mass of the lightest Higgs boson, m h < m Z , is not compatible with the LHC measurements. However, the CP-even mass matrix receives corrections beyond leading order. To a good approximation, the leading radiative corrections can be parametrized as an additional contribution to λ 2 [90], corresponding to the shift λ 2 → λ 2 + δ, while (sub-leading) MSSM contributions to the other quartic couplings are neglected. This is reminiscent of the approach pursued in [91]. Since the leading radiative corrections in the MSSM can most easily be specified in terms of the SUSY-basis, we use this instead of hybrid basis input for Scenario G. Working in this approximation, we define our MSSM-inspired 2HDM scenario by three parameters: m h , m A and tan β. Using this input, together with the MSSM relations for the λ i (i = 2) given by Eq. (88), we solve for the value of δ necessary to reproduce the desired value of m h for the chosen values of m A and tan β. In practice, this is done with 2HDMC using an iterative procedure. 13 To satisfy the constraints of unitarity and perturbativity, we impose the additional condition λ 2 < 4π.
Although this scenario is inspired by the MSSM, it is (like all our other benchmarks) completely defined within the CP-conserving, softly-broken Z 2 -symmetric 2HDM. There is therefore no principal restriction to the Type-II structure for the Yukawa couplings, although we choose to stick to this familiar pattern here to keep the connection to the MSSM tree-level structure. Using a fixed value m h = 125 GeV, we scan over the remaining parameters m A and tan β to determine the viable regions of parameter space. The results are shown in Fig. 16, where green color indicates regions compatible with m h = 125 GeV (and λ 2 < 4π). The shaded (gray) regions show the excluded regions at 95% C.L. from direct Higgs searches. In particular the limit H/A → τ τ plays a very important role to constrain this scenario both for high and low values of tan β for values m A < 2 m t . Constraints from heavy Higgs searches benefit in this scenario from the near-degeneracy of the heavy Higgs bosons H and A.
The resulting phenomenology of this scenario is very similar to the MSSM Higgs sector in the absence of additional low-energy degrees of freedom, which in the MSSM could provide additional decay channels for the heavy non-SM-like Higgs bosons. This can be seen, for example, from the region allowed by the LHC Higgs measurements in Fig. 16. In the absence of large radiative corrections affecting the CP-even Higgs mixing [81], the LHC Higgs data forces the MSSM into the decoupling limit. Hence, we find m A > 360 GeV at 95% C.L. almost independently of tan β. A related analysis in the MSSM context at low values of tan β with a heavy supersymmetric spectrum that makes use of the framework described above has recently been presented in [92].

Additional scenarios for consideration
There are a number of additional scenarios worthy of consideration that we have not included in this work. For example, one can consider the inert 2HDM with Type-I Yukawa couplings, which provides a plausible dark matter candidate [59,[93][94][95][96]. As noted at the end of Section 2.6, phenomenological constraints and benchmark scenarios for this model have been recently discussed in [60].
It is also possible to consider simplified models of 2HDM phenomena, where h is SM-like and one of the non-SM-like Higgs states among H, A and H ± is significantly lighter than the two heaviest states (keeping in mind that the split spectrum must be consistent with constraints due to the T parameter). In this case, one can focus on the phenomenology of a single non-SM-like scalar. Scenarios A and B already provide examples of this type, in which A and H ± are significantly heavier than h and H. But other scenarios could be considered that would feature the production and decay of A or H ± alone. In these latter scenarios, it would be appropriate to slightly modify the hybrid basis, which specifies two scalar masses as input parameters. When focusing on either A or H ± , it would be more appropriate to specify the SM-like Higgs mass and one other scalar mass [m A or m H ± ] as the two input scalar masses. Of course, one could then use Eqs. (67) and (68) to determine the masses of the remaining non-SM-like Higgs scalars.

Conclusions
We have introduced a new "hybrid" basis to define input parameters in the general, CP-conserving, two-Higgs-doublet Model with a softly broken Z 2 symmetry in a way that naturally accommodates constraints from perturbativity and unitarity. In this basis, the input parameters are the two CP-even Higgs masses, m h and m H , three quartic couplings defined in the Higgs basis, Z 4 , Z 5 , and Z 7 , and the mixing angles cos(β − α) and tan β which are very relevant for the phenomenology. In this hybrid approach, tan β can be interpreted as specifying the basis where the soft Z 2 -breaking is manifest.
The hybrid basis of parameters has been implemented as an input option for the public computer code 2HDMC. Using this setup, we have performed a scan and a detailed numerical analysis of the parameter space constraints and predictions of cross sections and branching ratios for different 2HDM configurations. Using the results of this analysis, we have defined a set of 2HDM benchmark scenarios that we deem relevant for the design and interpretation of Higgs searches beyond the SM at LHC run-II and beyond. The scenarios are defined to cover different aspects of the 2HDM phenomenology, including interpretation of the 125 GeV Higgs signal as either the light or the heavy CP-even Higgs boson. We also devise a new scenario with overlapping CP-even and CP-odd Higgses at 125 GeV, which could potentially be very useful for the interpretation of CP-studies in the τ τ final state, as well as a scenario where the Yukawa coupling of the Higgs to down-type fermions has the "flipped" sign (relative to the coupling to vector bosons). Another two sets of scenarios provide Higgs cascade decays as a prominent feature, giving rise to completely non-standard final states that may contain additional leptons and/or jets. Finally, we define the 2HDM equivalent of the decoupled MSSM, where the Higgs spectrum follows MSSM relations with the dominant radiative corrections to the Higgs masses completely captured in a shift of the quartic coupling λ 2 appearing in the Higgs potential. Results presented for this scenario should be immediately familiar to anyone following MSSM Higgs searches at the LHC, which could make it useful for comparison and communication of results.
We hope that the presented scenarios will be useful and serve to inspire new discoveries!

A. Parameter values for 2HDM benchmark lines and planes
For convenience, this appendix contains a list of the input parameters that can be used to realize the 2HDM benchmark scenarios described in this paper.