The Impact of Two-Loop Effects on the Scenario of MSSM Higgs Alignment without Decoupling

In multi-Higgs models, the properties of one neutral scalar state approximate those of the Standard Model (SM) Higgs boson in a limit where the corresponding scalar field is roughly aligned in field space with the scalar doublet vacuum expectation value. In a scenario of alignment without decoupling, a SM-like Higgs boson can be accompanied by additional scalar states whose masses are of a similar order of magnitude. In the Minimal Supersymmetric Standard Model (MSSM), alignment without decoupling can be achieved due to an accidental cancellation of tree-level and radiative loop-level effects. In this paper we assess the impact of the leading two-loop O(alpha_s h_t^2) corrections on the Higgs alignment condition in the MSSM. These corrections are sizable and important in the relevant regions of parameter space and furthermore give rise to solutions of the alignment condition that are not present in the approximate one-loop description. We provide a comprehensive numerical comparison of the alignment condition obtained in the approximate one-loop and two-loop approximations, and discuss its implications for phenomenologically viable regions of the MSSM parameter space.


Introduction
Since the initial discovery of a new scalar particle with mass of about 125 GeV [1,2], detailed studies of the data from Run 1 and 2 of the Large Hadron Collider (LHC) at CERN are beginning to establish the phenomenological profile of what appears to be the Higgs boson associated with electroweak symmetry breaking. Indeed, the measurements of Higgs production cross sections times decay branching ratios into a variety of final states appear to be consistent with the Higgs boson predicted by the Standard Model (SM) [3]. One can now say with some confidence that a "SM-like" Higgs boson has been discovered. Nevertheless, the limited precision of the current Higgs data from the LHC still allows for deviations from SM expectations. If such deviations were to be confirmed, new physics beyond the SM would be required.
Deviations from the SM Higgs behavior can be accommodated by introducing additional Higgs scalars to the electroweak model. Typically, the SM Higgs sector is extended by adding additional electroweak scalar doublets and/or singlets in order to avoid deviations of the approximate relation between the W and Z boson mass, M W M Z cos θ W , where θ W is the weak mixing angle. However, the existence of a SM-like Higgs boson already imposes significant constraints on any extended Higgs sector. We can always define a neutral Higgs field that points in the direction of the scalar doublet vacuum expectation value (vev) in field space. The tree-level couplings of such a scalar field to the SM gauge bosons and fermions are precisely those of the SM Higgs boson. However in general, this aligned scalar field is not a mass eigenstate field, since it will mix with other neutral scalar fields of the extended Higgs sector. Thus, the current Higgs data is consistent with an extended Higgs sector only if the observed scalar particle with mass 125 GeV is approximately aligned in field space with the doublet vev. This so-called alignment limit [4][5][6][7][8] is either the result of some symmetry of the scalar sector [9,10], or it is the result of some special choice of the scalar sector parameters.
An example of the latter is the decoupling regime of the extended Higgs sector [4,11]. The scalar potential typically contains a number of mass parameters. One of those mass parameters is fixed by the doublet scalar vev, which must be set to v = 246 GeV to explain the observed value of the Fermi constant G F . If other scalar sector mass parameters are characterized by a scale M that is significantly larger than v, then one of the neutral scalar mass eigenstates will be of O(v), whereas all other scalar masses will be of order M v. One can then integrate out the heavy scalar states below the mass scale M . The resulting effective scalar theory will be that of the SM with a single Higgs doublet, which will yield one neutral Higgs boson state whose couplings are approximately those of the SM Higgs boson. Of course, in such a scenario, additional scalar states would be quite heavy and may be difficult to discover at the LHC.
One can also achieve alignment independently of the masses of the non-SM like Higgs bosons. Generically, the aligned scalar field (which possesses the couplings of the SM Higgs boson) is not a mass eigenstate. However, if the parameters of the scalar sector (either accidentally or due to a symmetry) yield suppressed mixing between the aligned scalar field and the other neutral scalar interaction eigenstates, then approximate alignment is realized. In any multi-Higgs doublet model, an exact alignment condition can be specified, in which the aligned scalar field is a mass eigenstate (and thus its mixing with all other scalar eigenstate fields vanishes). Hence, if this alignment condition is approximately fulfilled, it is possible to have a SM-like Higgs boson along with additional scalar states with masses that are not significantly larger than the electroweak scale and thus more amenable to discovery in future LHC runs. We denote this latter scenario as alignment without decoupling [4][5][6][7][8][12][13][14].
Extended Higgs sectors in isolation suffer from the same problem as the SM Higgs sector, namely there is no natural explanation for the origin of the electroweak scale. There have been numerous attempts in the literature to devise models of new physics beyond the SM (BSM) that can provide a natural explanation of the electroweak scale, either via new dynamics or a new symmetry. All such approaches invoke new fundamental degrees of freedom, and many models of BSM physics incorporate enlarged scalar sectors. One of the most well studied models of this type is the minimal supersymmetric extension of the SM (MSSM) [15][16][17][18], which requires a second Higgs doublet in order to avoid anomalies associated with the supersymmetric fermionic partners of the SM Higgs doublet. In light of the fact that no supersymmetric particles have yet been discovered, it follows that the scale of supersymmetry (SUSY)-breaking, M S , must lie somewhat above the electroweak scale. This already leads to some tension with the requirements of a natural explanation of the electroweak scale (sometimes called the little hierarchy problem [19][20][21][22]). Nevertheless, if supersymmetric particles are ultimately discovered at the LHC, it would provide a significant amelioration of the large hierarchy problem associated with the fact that the electroweak scale is 17 orders of magnitude smaller than the Planck scale.
Numerous searches for supersymmetric particles at the LHC (as well as at previous lower energy colliders such as LEP and Tevatron) provide important constraints on the allowed MSSM parameter space [23,24], with additional constraints from considerations of virtual supersymmetric particle contributions to SM processes (see, e.g., Ref. [25] for a review). Finally, due to the enlarged Higgs sector of the MSSM, the properties of the observed Higgs boson and the absence of evidence for additional Higgs scalars yield additional constraints. In particular, given that the observed Higgs boson appears to be SM-like, it follows that the Higgs sector of the MSSM must be close to the alignment limit. In the MSSM, the scale of the non-SM-like Higgs boson is governed by a SUSY-breaking mass parameter. Although this mass parameter is logically distinct from the mass parameter M S that governs the mass scale of the heavy supersymmetric particles, one might expect these two parameters to be of a similar order of magnitude. If this is the case, then the approximate alignment limit of the MSSM Higgs sector is a result of the decoupling of heavy Higgs states. On the other hand, one may wonder whether the approximate alignment limit of the MSSM Higgs sector can be achieved outside of the decoupling limit, in which case one might expect the possibility that additional non-SM like Higgs scalars could soon be discovered in future LHC running.
The possibility of alignment without decoupling has been analyzed in detail in Refs. [4][5][6][7][8][12][13][14]. 1 More recently, the connection of Higgs alignment without decoupling in the MSSM with dark matter has been investigated in Ref. [27]. In Ref. [28], a parameter scan of the phenomenological MSSM (pMSSM) with eight parameters was performed, taking into account the experimental Higgs boson results from Run I of the LHC and further low-energy observables. One of the central questions considered in Ref. [28] was whether parameter regimes with approximate Higgs alignment without decoupling are still allowed in light of the current LHC data. Two separate cases were considered in which either the lighter or the heavier of the two CP-even neutral Higgs bosons of the MSSM is identified with the observed Higgs boson of mass 125 GeV. In the first case, we identified allowed regions of the MSSM parameter space in which the non-SM-like Higgs bosons could be as light as 200 GeV. In the second case, we demonstrated that the heavy CP-even Higgs boson is still a viable candidate to explain the Higgs signal -albeit only in a highly constrained parameter region. Both cases correspond to parameter regimes of approximate alignment without decoupling.
In the MSSM, alignment without decoupling arises due to an approximate accidental cancellation between tree-level and loop-level effects. Given the current precision of the Higgs data, we concluded in Ref. [28] that this region of approximate cancellation, while accidental in nature, does not require an extreme fine-tuning of the MSSM parameters. Indeed, such regions must appear in any comprehensive scan of the MSSM parameter space. In Ref. [28], we showed that the result of our numerical scans could be understood using simple analytical expressions in which the leading one-loop and two-loop radiative corrections to the MSSM Higgs sector are included. In this paper, we provide a detailed treatment of this analytic approximation and demonstrate the importance of the leading two-loop radiative effects in determining the allowed parameter regions for approximate alignment without decoupling.
The remainder of this paper is structured as follows. In Section 2, we review the alignment limit at tree-level in the context of the general CP-conserving two Higgs doublet model 2HDM). Both the decoupling limit and the limit of alignment without decoupling are discussed. We can apply these results to the MSSM by treating the MSSM Higgs sector as an effective non-supersymmetric 2HDM at tree-level, obtained by integrating out heavier supersymmetric particles. The effects of the SUSY-breaking lead to corrections that are logarithmic in the supersymmetry breaking scale, M S , as well as finite threshold corrections that can be of O(1). The leading one-loop corrections to the exact alignment condition are treated in Section 3. However, it is known that the two-loop corrections to the MSSM Higgs sector can be phenomenologically relevant. Employing a procedure first introduced in Ref. [29] and later extended in Ref. [30], the one-loop results of Section 3 are modified to obtain the leading two-loop corrections to the exact alignment condition in Section 4. In Section 5, a numerical comparison of the impact of the corresponding leading one-loop and two-loop corrections is given. In addition, we discuss the M S values required to achieve a SM-like Higgs boson mass of 125 GeV, and give a criterion on the CP-odd Higgs mass, M A , that determines whether the lighter or the heavier CP-even Higgs boson is aligned with the SM Higgs vev. In the latter scenario in which the heavier of the two CP-even Higgs bosons is identified with the observed Higgs scalar at 125 GeV, a new decay mode H → hh is possible if m H > 2m h . We discuss the magnitude of the relevant triple Higgs coupling and the resulting branching fraction for this decay in Section 6. Finally, we present our conclusions and outlook in Section 7.

The alignment limit in the two Higgs doublet model
In light of the LHC Higgs data, which strongly suggests that the properties of the observed Higgs boson are SM-like [3], we seek to explore the region of the MSSM parameter space that yields a SM-like Higgs boson. Since the Higgs sector of the MSSM is a constrained CP-conserving 2HDM, we first review the limit of the 2HDM that yields a SM-like Higgs boson. In a multi Higgs doublet model, a SM-like Higgs boson arises in the alignment limit, in which one of the neutral Higgs mass eigenstates is approximately aligned with the direction of the Higgs vacuum expectation value (vev) in field space.
The 2HDM contains two hypercharge-one weak SU(2) L doublet scalar fields, Φ 1 and Φ 2 . By an appropriate rephasing of these two fields, one can choose their vevs, Φ 0 to be real and non-negative. In this convention, It is convenient to introduce the following linear combinations of Higgs doublet fields, Re H 0 1 −v were a Higgs mass eigenstate, then its tree-level couplings to SM particles would be precisely those of the SM Higgs boson. This would correspond to the exact alignment limit. To achieve a SM-like Higgs boson, it is sufficient for one of the neutral Higgs mass eigenstates to be approximately given by with a corresponding squared-mass Z 1 v 2 . The observed Higgs mass implies that Z 1 0.26.
The CP-even neutral Higgs squared-mass matrix given by Eq. (6) is controlled by two independent mass scales, v 246 GeV and Y 2 , where the latter enters via the parameter M 2 A [cf. Eq. (5)]. In addition, the scalar potential parameters Z 1 , Z 5 and Z 6 are typically of O(1) or less (in the MSSM, they are of order the square of a gauge coupling). Consequently, a SM-like neutral Higgs boson can arise in two different ways: This corresponds to the so-called decoupling limit, where h is SM-like and In particular, one can achieve alignment without decoupling if |Z 6 | 1, independently of the value of the non-SM-like Higgs states H, A and H ± . Indeed, if the heavier of the two neutral CP-even Higgs states is SM-like, then one must have |Z 6 | 1 in a non-decoupling parameter regime.
After diagonalizing the CP-even neutral Higgs squared-mass matrix, one obtains the CPeven Higgs mass eigenstates h and H (where m h < m H ), 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, Since the SM-like Higgs field must be approximately We can now apply the above results to the MSSM Higgs sector. In the usual treatment of the MSSM, one introduces two Higgs doublets, H U and H D of hypercharge Y = +1 and Y = −1, respectively. 2 To make contact with the notation of the 2HDM presented above, we can relate these fields to the hypercharge Y = +1 scalar fields, where 12 = − 21 = 1 and 11 = 22 = 0, and there is an implicit sum over the repeated SU(2) L index j = 1, 2. The tree-level quartic couplings Z i can be expressed in terms of the electroweak SU(2) L and U(1) Y gauge couplings g and g , respectively, where c 2β ≡ cos 2β and s 2β ≡ sin 2β. We have already noted that the squared-mass of the SM-like Higgs boson is approximately given by Z 1 v 2 , which is equal to M 2 Z c 2β at tree-level in the MSSM, and thus incompatible with the observed Higgs mass of 125 GeV. Moreover, if the existence of a SM-like Higgs boson is due to alignment without decoupling, then the relation Z 6 = 0 must be approximately fulfilled, which implies that sin 4β = 0 (i.e. β = 0, 1 4 π or 1 2 π). Of course, the extreme values of β = 0 or β = 1 2 π are not phenomenologically realistic, whereas β = 1 4 π would yield a massless CP-even Higgs boson at tree-level. In order to achieve a realistic MSSM Higgs sector, radiative corrections must be incorporated [34][35][36] (see, e.g., Refs. [25,[37][38][39] for reviews). It is well-known that the observed Higgs mass of 125 GeV is compatible with a radiatively-corrected Higgs sector in certain regions of the MSSM parameter space [40]. Moreover, a SM-like Higgs state is easily achieved in the decoupling limit where M 2 A v 2 , where h is identified as the observed Higgs boson. In this paper, we focus on the alternative scenario in which a SM-like Higgs boson is a consequence of approximate alignment without decoupling. When loop corrections are taken into account, the possibility of alignment without decoupling must be reconsidered.

Alignment without decoupling at the one-loop level
In the MSSM, exact alignment via Z 6 = 0 can only happen through an accidental cancellation of the tree-level terms with contributions arising at the one-loop level (or higher). In this case the Higgs alignment is independent of the values of M 2 A , Z 1 and Z 5 . The leading one-loop contributions to Z 1 , Z 5 and Z 6 proportional to h 2 t m 2 t , where m t is the top quark mass and is the top quark Yukawa coupling, have been obtained in Ref. [12] in the limit M Z , M A M S (using results from Ref. [41]): where s β ≡ sin β, s 2β ≡ sin 2β, c 2β ≡ cos 2β, M S ≡ √ mt 1 mt 2 denotes the SUSY-breaking mass scale that governs the top-squark (stop) sector, given by the geometric mean of the light and heavy stop masses, and 3 The approximate expression for Z 6 v 2 given in Eq. (13) depends only on the unknown parameters µ, A t , tan β and M S . Exact alignment arises if Z 6 = 0. Note that Z 6 = 0 is trivially satisfied if β = 0 or 1 2 π (corresponding to the vanishing of either v 1 or v 2 ). However, this choice of parameters is not relevant for phenomenology as it leads to a massless b quark or t quark, respectively, at tree-level. 4 Henceforth, we assume that tan β is finite and non-zero; by convention, we take tan β to be positive. Regarding the other parameters, µ, A t and M S , we generously allow for rather large parameter values in this work. However, one should keep in mind that parameter points with |µ/M S | and |A t /M S | larger than about 3 are often severely restricted by vacuum (meta-)stability requirements, in particular the absence of a color and/or electric charge-breaking global minimum of the full MSSM scalar potential [43][44][45][46][47][48][49] (for recent analyses, see also Refs. [50][51][52].) Furthermore, for values of |X t /M S | > ∼ 3, the theoretically predicted loop-corrected Higgs squared-mass decreases rapidly from its maximal value (which at one-loop is achieved at X t /M S √ 6), and is ultimately driven to negative values. In our numerical analysis, we consider only |A t /M S | values up to about 3, and highlight regions of the parameter space that exhibit |X t /M S | ≥ 3 where our analysis is untrustworthy.
We simplify the analysis by solving Eq. (11) for ln(M 2 S /m 2 t ) and inserting the result back into Eq. (13). The resulting expression for Z 6 now depends on Z 1 , tan β, and the dimensionless ratios Using Eq. (14) to rewrite the final expression in terms of A t and µ, we obtain Setting Z 6 = 0, we can identify Z 1 v 2 with the mass of the observed (SM-like) Higgs boson (which may be either h or H depending on whether s β−α is close to 1 or 0, respectively). We can then numerically solve for tan β for given values of A t and µ. Indeed, t β ≡ tan β is the solution to a seventh order polynomial equation, 3 The elements of the top-squark squared-mass matrix are governed by the supersymmetric higgsino mass parameter µ and the soft-SUSY-breaking trilinear H 0 U t L t * R coupling A t . For simplicity, we ignore potential CP-violating effects by taking µ and A t to be real parameters in this work. 4 A potential loophole to this last remark arises in models, dubbed "uplifted supersymmetry", in which down-type fermion masses are absent at tree-level but are generated radiatively by loop-induced couplings to the up-type Higgs doublet, H U . Further details are described in Ref. [42]. A seventh order polynomial has either one, three, five or seven real roots. In light of the comments below Eq. (14), we are only interested in real positive solutions of Eq. (17); i.e., we exclude the possibility of t β = 0. Moreover, we can interpret the negative t β solutions at the point ( µ , A t ) as corresponding to positive t β solutions at the point (− µ , A t ). 5 Finally, the solution to Eq. (17) is invariant under the simultaneous inversion of µ → − µ and A t → − A t (keeping the sign of tan β fixed). The latter is a consequence of the symmetry properties of the approximate one-loop expressions for Z 1 , Z 5 and Z 6 . It therefore follows that a negative t β solution at the point ( µ , A t ) corresponds to a positive t β solution at the point ( µ , − A t ).
In the left [right] panel of Fig. 1 we show the number of real [positive] solutions to the above polynomial, Eq. (17), in the ( µ , A t ) plane. We observe that there is one real root of Eq. (17) for | µ| < ∼ 5 to 8 (depending on the value of A t = 0). For larger values of | µ| in the ( µ , A t ) plane, there are three real roots. The transition between these two regions occurs when two of the three roots coalesce (yielding a degenerate real root) and then move off the real axis to form a complex conjugate pair. In the quadrants with µ A t > 0 ( µ A t < 0) with large | µ|, these two real roots are always positive (negative), whereas the sign of the third root, which also exists at smaller | µ|, depends on the value of , it is of the opposite sign. To see how the roots evolve in a continuous manner in the ( µ, A t ) parameter plane, consider a path in the left panel of Fig. 1 that begins at A t ∼ 3 and µ ∼ −9. For these values, Eq. (17) possesses three negative roots and no positive roots. Keeping µ fixed and reducing A t , one of the negative roots decreases without bound until it reaches −∞ at A t ∼ √ 6. Taking A t below √ 6, the root switches over to +∞, and then steadily decreases. When A t crosses from positive to negative values, all positive and negative roots interchange. Consequently, when we enter the quadrant where A t is negative, we now have two positive roots and one negative root. This happens because one of the negative roots goes to −∞ as A t approaches zero from above, and then switches over to +∞ after crossing A t = 0. Finally, the remaining negative root goes to −∞ as A t approaches − √ 6 from above, and then switches over to +∞. For values of A t < − √ 6, there is now a third positive root. For smaller values of | µ|, Eq. (17) possesses only one real root, since the two other roots that were real at larger values of | µ| are now complex, as noted above. Finally, the two panels of Fig. 1 are symmetric under µ → − µ and A t → − A t , reflecting the symmetry of Eq. (17). We shall discuss the tan β values of all these roots in greater detail in Section 5, when we present the results of our numerical analysis.
It is instructive to obtain an approximate analytic expression for the value of the largest real root. Assuming µ A t tan β 1 the following approximate alignment condition, first written in Ref. [12], is obtained, where t < 2 and | µ| is sufficiently large such that the numerator of Eq. (18) is negative. Keeping in mind that Eq. (18) was derived under the assumption that µ A t tan β 1, we have observed in our numerical evaluation in Section 5 that the largest of the three roots of Eq. (17) always satisfies the stated conditions above. Another consequence of Eq. (18) is that by increasing the value of | µ A t | (in the region where 2 < A 2 t < 6), it is possible to lower the tan β value at which alignment occurs.
If | A t | 1, then Eq. (18) is no longer a good approximation. Returning to Eq. (17), we set A t = 0 and again assume that tan β 1. We can then solve approximately for tan β, For example, in the parameter regime where A t 0 and | µ| 1, we obtain tan β | µ|/ √ 6. Once the value of t β corresponding to exact alignment is known at a specific point in the ( µ , A t ) plane, we can use Eq. (11) to determine the value of the SUSY mass scale, M S , such that Z 1 v 2 = (125 GeV) 2 is the observed Higgs squared mass. We shall explore the numerical values of M S in Section 5 for each of the physical solutions of the alignment condition.
The question of whether the light or the heavy CP-even Higgs boson possesses SM-like Higgs couplings in the alignment without decoupling regime depends on the relative size of Z 1 v 2 and A . Combining Eqs. (12) and (13), it follows that in the limit of exact alignment where Z 6 = 0, we can identify Z 1 v 2 as the squared mass of the observed SM-like Higgs boson and We define a critical value of M 2 A , where Z 1 v 2 = (125 GeV) 2 and Z 5 v 2 is given by Eq. (20). Note further that the squared-mass of the non-SM-like CP-even Higgs boson in the exact alignment limit, M 2 A + Z 5 v 2 , must be positive, which implies that the minimum value possible for the squared-mass of the CP-odd Higgs boson is That is, if Z 5 is negative, then the minimal allowed value of M 2 A is non-zero and positive. If we compute Z 5 from Eq. (20) using the value of tan β obtained by setting A given in Eq. (22)], then H can be identified as the SM-like Higgs boson with M H 125 GeV. We shall exhibit numerical results for M A,c in Section 5 for each of the realistic tan β solutions.
Finally, we note that using the same one-loop approximations employed in this section, the leading contribution to the squared-mass splitting of the charged Higgs boson and CP-odd Higgs boson is given by [41], In particular, in the parameter regime in which H is identified as the SM-like Higgs boson, there is an upper bound on the charged Higgs mass obtained by inserting M A = M A,c into Eq. (23). In this case, collider and flavor constraints relevant to such a light charged Higgs boson can significantly reduce the allowed MSSM parameter space [28,53].
In Ref. [30], it was shown that the dominant part of these two-loop corrections can be obtained from the corresponding one-loop formulae with the following very simple two step prescription. First, we replace where m t (m t ) 163.6 GeV is the MS top quark mass [54], and the running top quark mass in the one-loop approximation is given by In our numerical analysis, we take α s = α s (m t (m t )) 0.1088 [54]. Second, when m 4 t multiplies the threshold corrections (i.e., the one-loop terms proportional to X t and Y t ), then we make the replacement, where Note that the running top-quark mass evaluated at M S includes a threshold correction at the SUSY-breaking scale that is proportional to X t . Here, we only keep the leading contribution to the threshold correction under the assumption that m t M S (a more precise formula can be found in Appendix B of Ref. [30]). The above two step prescription can now be applied to Eqs. (11)- (13), which yields a more accurate expression for the radiatively corrected Higgs mass and the condition for exact alignment without decoupling.
In applying the prescription outlined above, we formally work to O(α s ) while dropping terms of O(α 2 s ) and higher. For example, The end results are the following approximate expressions for Z 1 , Z 5 and Z 6 that incorporate the leading two-loop O(α s m 2 t h 2 t ) effects, where we have defined, and In the above equations, m t ≡ m t (m t ) is the MS top quark mass. Note that the approximate loop-corrected formulae for Z 1 , Z 5 and Z 6 are no longer invariant under X t → −X t , Y t → −Y t (or equivalently A t → −A t , µ → −µ) due to the asymmetry introduced by Eq. (27) at O(α s ).
We can now derive analogous expressions to Eqs. (17) and (20) that incorporate the leading two-loop effects at O(α s m 2 t h 2 t ). First, we note that Eq. (29) yields where B 1 is to be determined. Inserting Eq. (34) into Eq. (29), the O(1) terms cancel exactly.
Keeping only terms of O(α s ), we end up with the following expression for B 1 , Now we substitute Eq. (35) back into Eq. (34) to obtain Finally, we insert Eq. (36) into Eq. (31) and set Z 6 = 0 to obtain, That is, t β ≡ tan β is the solution to a 11th order polynomial equation, As previously noted, solutions to this equation for negative tan β at a point in the ( µ , A t ) plane can be reinterpreted as positive tan β solutions at the point (− µ , A t ).
In order to obtain two-loop improved versions of M 2 A,c and M 2 A,m [cf. Eqs. (21) and (22)], we need to impose the alignment limit condition, Z 6 = 0, on the two-loop expression for Z 5 given by Eq. (30). Our strategy is similar to the one employed above in deriving Eq. (37). First, we derive another expression for L based on Eq. (30); the steps leading to Eq. (36) are modified by the following substitutions, The end result is, Finally, we insert Eq. (40) into Eq. (31) and set Z 6 = 0 to obtain, Solving for Z 5 , and again expanding out in α s and dropping terms of O(α 2 s ) and higher, which yields the O(α s ) correction to Eq. (20). One can now define the two-loop improved versions of M 2 A,c and M 2 A,m via Eqs. (21) and (22). Likewise, the two-loop improved formula for the charged Higgs mass is obtained by replacing m t in Eq. (23) by m t (M S ) according to Eq. (27). The end result is In the left [right] panel of Fig. 2 we show the number of real [positive] solutions to the polynomial given in Eq. (38), corresponding to the two-loop condition for alignment without decoupling which determines tan β as a function of µ and A t . Compared to the one-loop results of Section 3, there are a few notable changes, which we now discuss. First, in our scan of the ( µ , A t ) plane, we have observed numerically that there are three real roots of Eq. (38) for | µ| < ∼ 8-10 (depending on the value of A t ), whereas for larger values of | µ|, a region opens up in which there are five real roots. As previously discussed, the transition between these two regions occurs when two of the real roots in the large | µ| regime coalesce (yielding a degenerate real root) and then move off the real axis to form a complex conjugate pair as the value of | µ| is reduced. Comparing with the roots of Eq. (17), we see that two new roots have come into play. We have analyzed these two roots and find that one is positive and one is negative. However, the positive root always corresponds to a value of |X t | > 3M S , which lies outside our region of interest. Henceforth, we simply discard this possibility. What remains then are at most two real roots at a given point in the ( µ , A t ) plane that can be identified as the two-loop corrected versions of the corresponding one-loop results obtained earlier.
We can now see the effects of including the leading O(h 2 t α s ) corrections. The regions where positive solutions to Eq. (38) exist, shown in the right panels of Fig. 2 (excluding the positive solution corresponding to |X t | > 3M S as noted above), have shrunk considerably in the two quadrants where µ A t > 0, as compared to the corresponding positive solutions to Eq. (17) shown in the right panel of Fig. 1. In contrast, in the two quadrants where µ A t < 0, the respective sizes of the regions where positive solutions to Eq. (17) and Eq. (38) exist are comparable.
One new feature of the two-loop approximation not yet emphasized is that we must now carefully define the input parameters µ and A t . In the formulae presented in this section, we interpret these parameters as MS parameters. However, it is often more convenient to reexpress these parameters in terms of on-shell parameters. In Ref. [30], the following expression was obtained for the on-shell squark mixing parameter X OS t in terms of the MS squark mixing parameter X t , where only the leading O(α s ) corrections are kept, Since the on-shell and MS versions of µ are equal at this level of approximation, we also have The approximations employed in the section capture some of the most important radiative corrections relevant for analyzing the alignment limit of the MSSM. However, it is important to appreciate what has been left out. The analysis of this section ultimately corresponds to a renormalization of cos(β − α), which governs the couplings of the Higgs boson in the effective 2HDM theory below the SUSY-breaking scale and its departure from the alignment limit. However, radiative corrections also contribute other effects that modify Higgs production crosssections and branching ratios. It is well-known that for M A M S , the effective low-energy theory below the scale M S is a general two Higgs doublet model with the most general Higgsfermion Yukawa couplings. These include the so-called wrong-Higgs couplings of the MSSM [55], which ultimately are responsible for the ∆ b and ∆ τ corrections that can significantly modify the coupling of the Higgs boson to bottom quarks and tau leptons. 6 In addition, integrating out heavy SUSY particles at the scale M S can generate higher dimensional operators that can also modify Higgs production cross-sections and branching ratios [57]. None of these effects are accounted for in the analysis presented in this section.
In this section we present the numerical results for the physical (i.e. real positive) tan β solutions of the alignment condition, and, in particular, compare the results obtained in the one-loop and two-loop approximations given in Sections 3 and 4, respectively. Moreover, we shall discuss for each of these solutions their implications for the correlated parameters, i.e. the SUSY-breaking mass scale, M S , and the critical M A value, M A,c , which determines whether the light or the heavy CP-even Higgs boson is the one aligned with the SM Higgs vev in field space.
As shown in Figs. 1 and 2, there may be more than one value of tan β corresponding to exact alignment for a given µ and A t . In the left and right panels of Previous works on Higgs alignment without decoupling in the MSSM [8,12,28] have largely focused on the tan β solution displayed in the two bottom panels of Fig. 3. This value of tan β, which can appear already at moderately large | µ| values, has also been employed in the definition of MSSM benchmark scenarios with Higgs alignment without decoupling [12,58]. However, this solution is associated with a large trilinear scalar coupling, A t , and thus part of the parameter space exhibiting this solution may yield a color or electric charge breaking vacuum and/or feature an unreliable theoretical prediction of the Higgs mass. In order to highlight this, we overlay the region where |X t |/M S ≥ 3 with a blue shading in Fig. 3. Since for the relevant parameter space, µ A t tan β 1, this tan β solution is approximated by Eq. (18), as first employed in Ref. [12]. Alignment without decoupling at moderately small values of tan β 10, as suggested by constraints from LHC H/A → τ + τ − searches [59,60], can be found for | µ| ∼ 2-3 and | A t | ∼ 3. Comparing the numerical values of this tan β solution obtained in the approximate one-loop and two-loop descriptions, we observe that the improved description at the two-loop level yields rather small corrections, which slightly increase the value of tan β. Furthermore, note the small asymmetry between the two sectors ( µ > 0 and A t > 0 vs. µ < 0 and A t < 0) introduced by the finite threshold correction proportional to X t entering at the two-loop level.
The smallest of the tan β solutions, displayed in the top panel of Fig. 3, was only briefly mentioned in Ref. [8,12], but was subject to detailed discussions in Ref. [28] in the context of scenarios where the observed SM-like Higgs boson was interpreted in terms of the heavy CP-even Higgs boson. In fact, such scenarios were found viable in this parameter region at | µ| ∼ 6-8, partly because for such large µ values, large ∆ b corrections suppress the light charged Higgs contribution to the rare flavor physics decays B → X s γ. In the top panel of Fig. 3, we observe that this tan β solution extends over all four sectors of the ( µ, A t ) parameter space, however, with the restriction that for | µ| 5 (7), the parameters µ and A t have to be of opposite sign in the one-loop (two-loop) description. In the latter case, and as long as | µ|| A t | tan β 1, the tan β alignment solution derived at the one-loop level [top left panel of Fig. 3] is approximately described by Eq. (18). The impact of the two-loop improved calculation on the numerical values of this solution is significant and again shifts the tan β values towards larger values. Whereas alignment without decoupling at moderately small values of tan β 10 is achieved in the oneloop description for µ 2.2, A t ∼ −1.3 and for µ < ∼ − 2.2, A t ∼ 1.3 (due to the µ → − µ, A t → − A t symmetry), the two-loop description pushes these results to higher absolute values of µ > ∼ 4.2 and µ < ∼ − 3.4, respectively. Even lower tan β values 5 can be obtained by allowing even larger µ values, as can be seen in Fig. 3. However, with further increasing µ 1, this turns over into the approximate behavior tan β | µ|/ √ 6, found in the limit µ 1 and small A t of Eq. (19), and thus tan β starts to increase with µ. At such large µ values alignment solutions are also found in the parameter regions with µ and A t having the same sign, which feature small values of tan β 5 (for the µ range considered here).
The remaining tan β solution, displayed in the middle panels in Fig. 3, has not been discussed previously in the literature (except for some brief comments in our previous work [28]). It occurs only in the regions where µ and A t have the same sign, and only for very large | µ| 5 (7 − 8) in the one-loop (two-loop) description. The tan β value of this solution is small at large | A t |, and approaches +∞ as | A t | → 0. This solution is only found in regions of the parameter space that also exhibit the solution shown in the top panels of Fig. 3, and its tan β values are always larger. Therefore, and because of the very large µ values required especially after taking into account the two-loop corrections, this alignment solution is phenomenologically not relevant.
In our numerical scans in the ( µ , A t ) plane, the size of the SUSY-breaking mass scale, M S , varies as required by the condition of exact alignment (Z 6 = 0) such that the SM-like Higgs mass is fixed to its observed value of 125 GeV. That is, given the value of t β for exact alignment at a point in the ( µ , A t ) plane, one can use Eq. (11) [Eq. (29)] in the one-loop [twoloop] approximation to determine the value of M S such that Z 1 v 2 = (125 GeV) 2 . These M S values are exhibited in the three rows of Fig. 4, which are in one-to-one correspondence with the three rows of Fig. 3, i.e. each row shows a different solution of the alignment condition, and on the left [right] we show the one-loop [two-loop] result. We define maximal mixing in the top squark sector to correspond to the value of X t /M S that maximizes the value of Z 1 v 2 given by Eq. (11) [Eq. (29)] in the one-loop [two-loop] approximation, prior to fixing the Higgs mass at its observed value of 125 GeV. In the one-loop approximation, maximal mixing occurs at X t /M S = √ 6, and the maximal value of the Higgs mass corresponds to maximal mixing with tan β 1. Turning this around, if we fix the value of the Higgs mass to be its observed value of 125 GeV, then the minimal value of M S occurs at maximal mixing with tan β 1. This can be seen in the top panels of Fig. 4, where the minimal M S contour is located in the region of A t = √ 6 and µ A t < 0. In this region, tan β 1 so that A t X t /M S , corresponding to the region of maximal mixing at large tan β. Maximal mixing at large tan β is also evident in the bottom panels of Fig. 4. At smaller values of | A t |, it is still possible to reach maximal mixing at large values of | µ|, albeit with smaller values of tan β shown in Fig. 3. In contrast, maximal mixing is never reached in the middle panels of Fig. 4, as the corresponding X t /M S values are closer to the minimal mixing value of X t = 0. In this case, the smaller values of M S are associated with the larger values of tan β, which occur when | A t | → 0. Under the assumption that the heavy CP-even Higgs boson H is identified with the observed Higgs boson at 125 GeV, the M A,c values can easily be translated into upper bounds on the charged Higgs boson mass, M H + , according to Eqs. (23) and (43) in the leading oneand two-loop description, respectively. Here one should keep in mind that the leading radiative corrections are negative and proportional to µ 2 . Consequently, at large µ, these radiative corrections can substantially decrease the M H + prediction with respect to its tree-level prediction, M tree H + = (M 2 A + M 2 W ) 1/2 . Collider and flavor constraints on such scenarios arising from a light charged Higgs boson have been extensively discussed in Ref. [28] (see also Ref. [53] for a similar analysis in the framework of the 2HDM). We close this section with a few comments on how these results compare with the numerical fit results found in Ref. [28]. In the global fit, Ref. [28] identified two distinct parameter regions with phenomenologically viable points near the limit of alignment without decoupling. These regions resemble the parameter regions that are exhibited in the top and bottom panels in Figs. 3-5. Specifically, under the assumption that h is identified as the SM-like Higgs boson at 125 GeV, the preferred points with low M A (i.e. in the non-decoupling regime) found in Ref. [28] were located near the alignment solution displayed in the bottom panels. The main reason for this, however, is the restriction |µ|/M S ≤ 3 imposed in the fit for the light Higgs interpretation, which essentially excludes the other possible alignment solutions identified in this work. In contrast, assuming that H is identified as the SM-like Higgs boson at 125 GeV, Ref. [28] found viable points only near the parameter regions displayed in the top panels of Figs. 3-5. Here, the restriction |µ|/M S ≤ 3 was not imposed, and the main reason for this observation was a coupling suppression of the charged Higgs contribution to the branching fraction of the B meson decay B → X s γ, thus yielding phenomenologically acceptable values despite the presence of a very light charged Higgs boson. Lastly, a word of caution is in order: the numerical results displayed in this work are based on the exact alignment limit, whereas in the global fit studies of the MSSM parameter space in the non-decoupling regime, the parameter points only need to be near the alignment limit in order to be phenomenologically viable. In particular, Ref. [28] quantified the maximal values of |Z 6 |/Z 1 for the parameter points allowed at the 2σ level in the light Higgs (with low M A ) and heavy Higgs interpretation, resulting in ∼ 0.3 and ∼ 0.2, respectively. Such values indicate that these parameter regions are not yet parametrically finetuned, and non-negligible deviations from the alignment limit are still allowed by the current data. 7 6 SM-like Higgs branching ratios in the alignment limit In the exact alignment limit, the tree-level couplings of the SM-like Higgs boson are precisely those of the Higgs boson of the SM. Nevertheless, in the case of alignment without decoupling, deviations from SM Higgs boson properties can arise because the effective theory at the electroweak scale contains additional fields beyond the fields of the SM. In this work, we have employed the framework of the MSSM under the assumption that the SUSY-breaking scale M S M Z , M H ± . Thus the effective electroweak theory at energy scales below M S is the 2HDM. Moreover, SUSY-breaking effects can generate so-called wrong-Higgs couplings with coefficients that in some cases are tan β-enhanced (see footnote 6). Thus, we are led to consider the 2HDM with the most general Higgs-fermion Yukawa interactions as the effective theory below M S . In particular, the masses of the additional scalar states are assumed to be of the same order as the scale of electroweak symmetry breaking. In the exact alignment limit, deviations of the SM-like Higgs boson branching ratios from the corresponding SM predictions can arise due to two possible effects: (i) new loop-contributions due to the exchange of non-SM Higgs scalars that modify partial decay rates, and (ii) new decay channels in which the SM-like Higgs boson decays into a pair of lighter scalars, if kinematically allowed.
If new tree-level Higgs decays are present, these will typically yield the dominant contributions to the deviations of the Higgs branching ratios from their SM values. In particular, we expect that any additional deviations that arise from the exchange of non-SM Higgs scalars (which compete with the SM loop corrections) would result only in small shifts of the Higgs decay rates away from their corresponding SM predictions, and will be difficult to isolate experimentally. In contrast, consider the loop-induced Higgs couplings to γγ and Zγ, which have no tree-level counterpart. In this case, new loop corrections due to charged Higgs exchange can compete with the corresponding SM loop contributions, since by assumption M H ± does not differ appreciably from the mass of the SM-like Higgs boson [14]. In practice, due to the domination of the W -loop contribution to the loop-induced Higgs couplings to γγ and Zγ relative to the fermion and scalar loop contributions, the shift in the loop-induced Higgs couplings from their SM values due to the contribution of charged Higgs exchange will typically be small.
The most significant deviation from SM Higgs branching ratios in the alignment limit with- out decoupling arises if new decay channels are present in which the SM-like Higgs boson decays into a pair of lighter scalars. In Ref. [28], we demonstrated that regions of the MSSM parameter space in which the heavier of the CP-even scalars, H, is SM-like and m h < m H /2 are still allowed after taking into account the experimental constraints from SUSY particle searches and the measurement of Higgs boson properties at the LHC. In such a scenario, the decay mode H → hh is kinematically allowed, which has an impact on the predicted SM Higgs branching ratios. At tree-level, the Hhh coupling survives in the exact alignment limit where s β−α = 0. Indeed, when expressed in terms of the coefficients of the scalar potential in the Higgs basis, the tree-level Hhh coupling is given by [4,13,14] g where Thus, in the alignment limit, g Hhh → −vZ 345 . In the one-loop corrected MSSM in the limit of M Z , M A M S , we make use of the results of Ref. [41] to obtain, 8 Including the approximate leading O(α s m 2 t h 2 t ) corrections we obtain where X 34 ≡ 1 4 (x t + y t ) 2 − 1 12 x 2 t y 2 t and X 5 ≡ x t y t (1 − 1 12 x t y t ) and the other relevant quantities have been defined in Eq. (32).
In the left and right panel of Fig. 6 we show the contours of Z 345 for the first alignment solution [shown in the top panels in Fig. 3] derived in the one-loop and two-loop description, respectively. For most of the parameter space Z 345 is negative, with large negative values found for large | µ| values. A small parameter region with small | A t | 0.4 and | µ| between 2.5 and 5 in the one-loop [5 and 7 in the two-loop] description exhibits small positive Z 345 values, shown by the dark green color in Fig. 6. At the boundary between the light and dark green region the coupling Z 345 vanishes. Thus, in these regions the decay H → hh becomes coupling-suppressed and the branching fraction can even become zero, irrespective of the available phase-space. This feature has been numerically observed in Ref. [28] and in particular exploited in Ref. [27] in the definition of low mass light Higgs boson benchmark scenarios for dark matter studies.
The other two alignment solutions [middle and bottom panels in Fig. 3] do not exhibit phenomenologically relevant parameter regions where the coupling Z 345 vanishes, and thus we do not exhibit them here. 9 For the second alignment solution [shown in the middle panels in Fig. 3], the Z 345 values are positive, whereas in the third alignment solution [shown in the bottom panels in Fig. 3], the Z 345 values are negative, with larger magnitudes found at larger values of | µ|.
In order to further illustrate the feature of a vanishing branching fraction BR(H → hh) in the first alignment solution, as mentioned above, we show BR(H → hh) for two choices of the light Higgs mass, M h = 10 GeV and 60 GeV, in the left and right panel of Fig. 7. Here, we exhibit the alignment solution in the two-loop description, and calculate the branching ratio, where Γ SM tot = 4.1 MeV is the SM Higgs boson total decay width [62], and 10 By comparing Fig. 7 with Fig. 6 one can clearly observe that the branching fraction vanishes in the region where Z 345 = 0. Nevertheless, while such an "accidental" parameter constellation leading to a vanishing coupling Z 345 can occur in this alignment solution, it should be noted that generically the Hhh coupling does not vanish in the alignment limit in the scenario where the heavy CP-even Higgs boson is identified with the observed Higgs boson at 125 GeV. Thus, precision measurements of the properties of the observed Higgs boson can be employed to further constrained the heavy SM-like Higgs scenario. 9 In the second solution [middle panels in Fig. 3], one can achieve a vanishing Z 345 in a small parameter strip in the second solution [middle panels in Fig. 3], which we have discarded as phenomenologically irrelevant. 10 Employing Eq. (49) for Z 345 in Eq. (51) incorporates some of the leading one and two-loop corrections to the H → hh decay rate. A more complete one-loop computation can be found in Ref. [63].

Conclusions and Outlook
Given the current precision of Higgs boson measurements at the LHC, the observed state with mass 125 GeV is consistent with a Higgs boson that possesses the spin, CP quantum number and coupling properties predicted by the SM. In an extended Higgs sector of a BSM theory, a scalar mass-eigenstate would possess the properties of the SM Higgs boson if it is aligned in field space with the scalar vacuum expectation value responsible for electroweak symmetry breaking. This defines the so-called alignment limit. Higgs alignment can be achieved due to the decoupling of the non-SM Higgs states, under the assumption that these scalars are considerably heavier than the SM-like Higgs boson, or via the suppression of the mixing between the aligned scalar state and the other scalar states of the Higgs sector. In the latter scenario, the possibility of alignment without decoupling arises if the masses of the non-SM-like Higgs states are of the same order of magnitude as the mass of the SM-like Higgs boson.
In the MSSM, the simplest way to achieve approximate Higgs alignment is in the decoupling limit in which M A M h and h is identified as the observed SM-like Higgs boson. In this paper, we have addressed the possibility of achieving approximate Higgs alignment without decoupling, which can arise due to an accidental cancellation of tree-level and loop-level effects, independently of the mass scale of the other Higgs states. Under the assumption that the SUSYbreaking scale, M S , is significantly larger than the masses of the non-SM-like Higgs scalars, the properties of the Higgs sector are well described by an effective two-Higgs doublet extension of the SM (2HDM). In previous work, alignment without decoupling was achieved due to the cancellation of tree-level and one-loop contributions to the effective 2HDM Lagrangian. 11 The present work goes beyond previous studies and provides a detailed analysis of the leading twoloop corrections and their impact on the parameter regions that exhibit an exact realization of the alignment without decoupling scenario. In particular, we assessed the leading radiative corrections proportional to the strong coupling constant α s , which first enters at the two-loop level, by employing an approximation scheme developed in Refs. [29,30]. This scheme employs an optimal choice of the renormalization scale for the running top quark mass, which captures the leading logs at the two-loop level as well as a significant part of the two-loop corrections proportional to the stop mixing parameter X t .
Taking the observed Higgs boson mass of 125 GeV as an additional constraint, the alignment condition (i.e., the equation corresponding to exact alignment, independent of the value of M A ) can only be fulfilled for a specific value of tan β and M S that depends on the location in the ( µ ≡ µ/M S , A t ≡ A t /M S ) parameter plane. We discussed all physical solutions of the alignment condition both at the one-loop and two-loop level. We found that up to three physical solutions exist simultaneously, out of which at most two appear to be phenomenologically relevant. Comparing the one-and two-loop approximations we found some significant differences in the number of physical solutions in the ( µ, A t ) parameter plane. Nevertheless, the gross qualitative features of the one-loop solutions are maintained in the two-loop improved results.
We presented a detailed numerical comparison of the resulting tan β and M S values obtained in the one-loop and two-loop approximation in the exact alignment limit. We found that the two-loop corrections are sizable and lead to significant changes of the phenomenology. In particular, the tan β values are corrected towards larger values, and the SUSY mass scale M S is corrected towards smaller values, with respect to the corresponding values obtained in the one-loop approximation. Because tan β is a parameter that significantly influences the collider phenomenology of the non-SM Higgs bosons (in particular the CP-odd Higgs boson A), the twoloop corrections to the alignment condition cannot be neglected in a detailed phenomenological study of the viability of the alignment without decoupling scenario in the MSSM [28]. We found that the SUSY-breaking mass scale M S varies in the ( µ, A t ) plane from below 500 GeV up to values in the multi-TeV range.
We furthermore defined a critical mass of the CP-odd Higgs boson, M A,c . For parameter points with M A below this value, the heavy CP-even Higgs boson plays the role of the SM-like Higgs boson at 125 GeV, whereas parameter points with M A > M A,c feature a SM-like light CP-even Higgs boson. We exhibited numerical results for M A,c for all viable solutions to the alignment condition in the one-loop and two-loop approximations. Again, we noted a significant impact of the two-loop corrections, which in general lead to a substantial downward shift of M A,c , thus narrowing the parameter space that exhibits a SM-like heavy Higgs boson H.
In the heavy Higgs interpretation, i.e. the scenario where the heavy CP-even Higgs boson plays the role of the SM-like Higgs boson at 125 GeV, a new decay mode H → hh is possible if m H > 2m h . We discussed the magnitude of the relevant triple Higgs coupling and the resulting branching fraction BR(H → hh) for two choices of the light Higgs mass, M h = 10 and 60 GeV, in the ( µ, A t ) plane. We find that generically the relevant coupling is unsuppressed in the limit of alignment without decoupling, thus leading to a value of BR(H → hh) that is in conflict with the LHC Higgs data. However, in one of the solutions to the alignment condition, the responsible triple Higgs coupling (accidentally) vanishes in certain regions of the parameter space. These regions are found at small | A t | 0.4 and large | µ| values around 3 to 5 (5 to 7) in the one-loop (two-loop) description. Parameter points exhibiting this accidental suppression of BR(H → hh) in the heavy Higgs interpretation have previously been observed numerically in Refs. [27,28].
In the effective 2HDM Lagrangian, exact alignment corresponds to setting the effective Higgs basis parameter Z 6 to zero. Given that exact Higgs alignment in the MSSM is achieved by an accidental cancellation between tree-level and loop-level contributions to Z 6 , the astute reader may object that this scenario is of no interest as it represents a set of measure zero of the MSSM parameter space. To address this concern, we first note that the present Higgs data implies that the observed state at 125 GeV is consistent with that of the SM Higgs boson with an accuracy that is roughly 20-30%. Consequently, as long as the parameters of the MSSM Higgs sector yield a result close to the alignment limit, such MSSM parameter regions are presently not ruled out by the Higgs data. Indeed, in Ref. [28], a detailed numerical scan of the parameter space of a phenomenological MSSM governed by eight parameters revealed the existence of regions in which approximate Higgs alignment without decoupling is satisfied. In the preferred region, some points with values of Z 6 as large as |Z 6 /Z 1 | ∼ 0.3 were within two standard deviations of the best fit point. Thus, the present Higgs data does not require excessive fine-tuning of the MSSM parameters to achieve approximate Higgs alignment without decoupling.
The analysis of the exact alignment limit given in this paper provides an understanding of the regions of the MSSM parameters where Higgs alignment without decoupling can occur, and the impact of including or neglecting the leading two-loop effects. Our analytic approximations include the leading effects proportional to the fourth power of the top quark Yukawa coupling and include leading logarithmic terms (sensitive to the mass scale of SUSY-breaking, M S , arising from the top squark sector), and the leading threshold effects at M S due to top squark mixing. However, subdominant effects proportional to the square of the top quark Yukawa coupling, the bottom quark and tau lepton Yukawa couplings, and the electroweak gauge couplings have been neglected, as well as non-leading logarithmic terms and non-leading threshold effects due to bottom squark mixing. It is straightforward to include such effects analytically (see, e.g. Ref [29]). Additional corrections not treated in this work include effects arising from higher dimensional operators (ultimately arising from integrating out the heavy SUSY sector) as well as genuine electroweak radiative corrections to the low-energy effective 2HDM. Nevertheless, the impact of including all such corrections, while modifying some of the precise details of the cancellation between tree-level and loop-level contributions in achieving exact Higgs alignment, will not change the overall qualitative understanding of the MSSM parameter regime that yields the approximate alignment limit without decoupling.
Further experimental Higgs studies at the LHC will improve the precision of the properties of the 125 GeV Higgs boson, while further constraining or discovering the existence of new scalar states of the extended Higgs sector. Both endeavors will be critical for providing a more fundamental understanding as to why the observed 125 GeV scalar resembles the SM Higgs boson.