Muon $g-2$ in the Aligned Two Higgs Doublet Model

We study the Two-Higgs-Doublet Model with the aligned Yukawa sector (A2HDM) in light of the observed excess measured in the muon anomalous magnetic moment. We take into account the existing theoretical and experimental constraints with up-to-date values and demonstrate that a phenomenologically interesting region of parameter space exists. With a detailed parameter scan, we show a much larger region of viable parameter space in this model beyond the limiting case Type X 2HDM as obtained before. It features the existence of light scalar states with masses $3$ GeV $\lesssim m_H^{} \lesssim 50$ GeV, or $\ 10$ GeV $\lesssim m_A^{} \lesssim 130$ GeV, with enhanced couplings to tau leptons. The charged Higgs boson is typically heavier, with $200$ GeV $ \lesssim m^{}_{H^+} \lesssim 630$ GeV. The surviving parameter space is forced into the CP-conserving limit by EDM constraints. Some Standard Model observables may be significantly modified, including a possible new decay mode of the SM-like Higgs boson to four taus. We comment on future measurements and direct searches for those effects at the LHC as tests of the model.


I. INTRODUCTION
After the discovery of a Higgs boson with mass around 125 GeV [1][2][3][4], the ATLAS and CMS collaborations have continued to accumulate data and measured the detailed properties of this new particle. Within the attainable accuracy, they appear consistent with the elementary Higgs boson of the Standard Model (SM). This milestone discovery strongly encourages the exploration for new physics, in particular the search for a richer structure of the Higgs sector beyond the SM. Although no clear indication exists yet for new physics beyond the SM, this search will be actively pursued at the LHC and elsewhere.
The apparent discrepancy with a 3 − 4 sigma deviation between the theoretical [5][6][7][8] and measured [9] values of the muon anomalous magnetic moment (g − 2) µ is a long-standing puzzle which may point to new physics [10][11][12][13][14]. The possibility that the observed excess in (g − 2) µ can be explained by an extended Higgs sector has been raised by several authors [15][16][17][18][19][20][21]. Extended Higgs sectors arise in a variety of new physics scenarios, including supersymmetry, Grand Unified Theories, dark matter models, flavor models, and others. A feature common to many models, and one of the simplest extensions, is the addition of a second Higgs doublet. We thus will focus on a general formulation of the two-Higgs-doublet-model (2HDM) [22][23][24][25][26][27][28] in this work in the hope to reach some general conclusions.
2HDMs have been extensively studied in the past, most often in the context of global flavor symmetries and with the assumption of no new sources of CP violation. In this paper we focus on an interesting broader formulation, the Aligned Two Higgs Doublet Model (A2HDM) [29]. We find that the A2HDM can account for the experimental value of (g − 2) µ [9]. However, it can only account for such an excess in a very restricted range of the otherwise-allowed parameter space. In particular, mixing between the new states and the SM Higgs boson must be very small, one of the new neutral states must be quite light, the new states will couple strongly to taus, and CP-violating effects must be negligible.
The paper is organized as follows: We describe the A2HDM model for our study in Sec. II, and then present the current experimental bounds on the model parameters in Sec. III.
Accounting for these constraints, we perform a detailed multiple-dimensional parameter scan in light of the (g − 2) µ excess, and show our main results in Sec. IV. We discuss the physical implication of these results on SM observables in Sec. V and comment on direct searches for the new Higgs states at the LHC in Sec. VI. In Sec. VII, we draw our conclusions.
An appendix is provided to summarize our (g − 2) µ calculation in the A2HDM model.

II. BRIEF REVIEW OF A2HDM
For two SU(2) L Higgs doublets with hyper-charge one, we can redefine a theoretical basis by rotating the doublets according to a global U (2) transformation [30]. We choose to work in the Higgs basis in which only one doublet, H 1 , contains the SM electroweak symmetry breaking vacuum expectation value (vev) [31]: The scalar potential can then be expressed using Higgs basis fields as follows: where Y 1 , Y 2 and Z 1,2,3,4 are real-valued and invariant in this basis, whereas Y 3 and Z 5,6,7 are complex "pseudo-invariants". The pseudo-invariants can be re-phased by the U (1) transformation H 2 → e iχ H 2 . The SM-like minimization condition for the scalar potential requires The neutral Higgs boson mass-eigenstates can be determined by diagonalizing a 3 × 3 squared-mass matrix given by 23 ) −Im(Z 6 e −iθ 23 ) Re(Z 6 e −iθ 23 ) A 2 /v 2 + Re(Z 5 e −2iθ 23 ) − 1 2 Im(Z 5 e −2iθ 23 ) −Im(Z 6 e −iθ 23 ) − 1 2 Im(Z 5 e −2iθ 23 ) where A 2 ≡ Y 2 + 1 2 [Z 3 + Z 4 − Re(Z 5 e −2iθ 23 )]v 2 . The diagonalizing matrix R, can be written as where c ij = cos θ ij and s ij = sin θ ij . In general θ 23 appears in the diagonalizing matrix R for the matrix M in some particular pseudo-invariant basis. Under the rephasing H 2 → e iχ H 2 , however, θ 23 → θ 23 − χ [32]. Thus R and M as written are invariant and we may use the remaining freedom in our choice of basis to eliminate θ 23 , which has no physical significance.
Here, we choose a convention where − 1 2 π ≤ θ 12 , θ 13 < 1 2 π. We will identify one of the mass eigenstates as the observed SM-like Higgs near 125 GeV. The charged Higgs mass is given by The Higgs couplings to the fermions can be written as where K is the CKM matrix and κ i = The A2HDM is defined by the relation [29] where A i is a potentially complex number. Thus the Yukawa couplings of the two Higgs doublets are aligned 1 in that they are proportional to each other and can be simultaneously diagonalized, which ensures that flavor changing neutral currents (FCNCs) do not arise at tree level. The A2HDM is thus an example of Minimal Flavor Violation [34]. FCNC effects will develop from loop diagrams involving the charged Higgs, which always depend on powers of the CKM matrix. The often studied discrete Z 2 2HDMs can be understood as special cases of this model which correspond with the relations among the A i parameters summarized in Table I. In these cases, the ratio of Yukawa couplings allows us to rotate to a basis where the symmetry is manifest and only one doublet couples to each fermion type 1 Note that this alignment is different from the "alignment limit" discussed in certain literature which corresponds with small mixing between H 1 and H 2 [33]. (up, down and lepton, respectively). In such cases, |A u | is equivalent to cot β ≡ v 1 /v 2 , the ratio of vevs in that basis.
In the 2HDM, CP is potentially broken in the scalar potential and vacuum as well as in the neutral Higgs Yukawa interactions. CP is preserved if the following terms are real [32] where i = u, d and l. In our chosen basis this is equivalent to setting Z 5 , Z 6 , Z 7 , A u , A d , A l to be real and it implies that s 13 = 0.
Details of the general CP-violating A2HDM and calculations relevant to our scan can be found in a forthcoming publication Ref. [35].

III. CONSTRAINTS ON A2HDM
There exist a number of relevant theoretical and experimental constraints on the A2HDM.
For theoretical considerations, we require that the quartic coupling parameters Z i satisfy the partial-wave unitarity bounds, which are taken from the Appendix of [32]. We also require that the SM-like electroweak symmetry breaking vev be at a local minimum, and that the potential be positive at large field values. For convenience, we set the parameter Z 2 = 4π, corresponding to the maximum value allowed by partial-wave unitarity. Z 2 does not enter any of our phenomenological considerations so we choose a value that contributes maximally to the stability of the potential.
On the experimental side, we take the following results into account:

A. Precision EW Data
We directly impose the LEP bound on the charged Higgs mass [36] m H ± ≥ 80 GeV.
For other searches applicable to exotic Higgs, we make use of the program HiggsBounds [37], which checks for exclusion of potential signals at the 95% level incorporating a large number of searches from the LEP, Tevatron, LHC and other experiments. We also check for compatibility with the Peskin-Takeuchi S, T and U parameters [38] at the Z pole.
A further constraint of interest for the large |A l | values we allow comes from measurements of the Zτ τ coupling measured at the Z pole. This is characterized by effective couplings g L τ and g R τ which have been determined to high precision. Experimentally they are found to be [39]: SM best fit predictions based on precision data yield [36] g L τ (SM ) = −0.26919 ± 0.0002, g R τ (SM ) = 0.23274 ± 0.0002.
Our calculation for the effects of the new Higgs states is based on Ref. [40].

B. B Meson Rare Decays
Typically, some of the strongest experimental constraints on the 2HDM come from flavor physics [41][42][43]. Although highly suppressed in the SM, FCNCs are experimentally observed and provide stringent constraints on new physics. The A2HDM guarantees that tree-level FCNCs vanish, but they will appear at the loop level. We make use of the NLO prediction of BR(B → X s γ) in the 2HDM available from the SusyBSG code provided by the authors of Ref. [44]. The numerical estimation of BR(B → X s γ) has been performed up to NNLO for the SM [45]. The most updated SM prediction is The most recent average of experimental data on b → sγ rate gives [36] BR(b → sγ)| exp = (3.43 ± 0.22) × 10 −4 .
NNLO corrections for the 2HDM are discussed in Ref. [46].
A second rare decay of interest is B s → µ + µ − . This decay is of high theoretical significance and may receive contributions from both charged and neutral Higgs bosons in the A2HDM via box and penguin diagrams. Details of the calculation can be found in Ref. [47].
The SM prediction works out to be Measurements of this decay have been made at both LHCb [48] and CMS [49] and found to be consistent with the SM expectation. The combined result can be expressed in a ratio for comparison with the Standard Model:

C. Lepton Universality
Many additional flavor phenomena can potentially be affected by new physics in the A2HDM. A large selection of relevant processes is analyzed in Ref. [29] in the A2HDM under the assumption that the charged scalar H + provides the dominant effects. These include leptonic tau decays, leptonic decays of heavy mesons such as B → τ ν and D (s) → µ(τ )ν, measurements of Z → bb, B 0 −B 0 mixing and K 0 −K 0 mixing. These can be translated into constraints at 95% confidence on the charged Higgs mass and fermion couplings as follows: The most relevant of these limits for our purposes is the first, which comes from tests of lepton flavor universality in tau decays. In particular, it derives from the ratio of tau decay widths Γ(τ → µνν)/Γ(τ → eνν), where Γ indicates the partial width normalized to its SM value. This ratio is equivalent to the ratio of fitted effective coupling parameters (g µ /g e ) 2 given by the HFAG collaboration [50]. In Ref. [19], the authors make use of the additional fitted value (g τ /g e ) 2 which can be translated into the ratio Γ(τ → µνν)/Γ(µ → eνν). The most recent HFAG fits to purely leptonic decays are g µ g e = 0.0018 ± 0.0014, g τ g e = 0.0029 ± 0.0015, g τ g µ = 0.0011 ± 0.0015.
Neglecting the electron Yukawa coupling, only τ → µν is altered by the tree-level charged Higgs graph.Within the measured bounds of the Michel parameters which characterize the decay distribution, the only allowable effects from the tree-level charged Higgs will tend to decrease Γ(τ → µνν). Thus the bound depends strongly on the fact that the current fit to (g τ /g e ) 2 is almost a 2σ excess over the Standard Model while the charged Higgs gives small negative corrections.
In light of this, Ref. [19] includes the leading one-loop effects from the new Higgs states which tend to decrease the coupling of W τ ν. The one-loop effect applies equally to both tau decays. Hence, the ratio g µ /g e is only affected by the tree-level Higgs graph, while g τ /g e should include both tree-level and leading one-loop effects. They exclude parameters which exceed the 95% confidence limits on the combined fit, accounting for correlations in the data. We adopt their calculation of the one-loop effects, which is consistent with Ref. [51].
In principle we may work with the ratio g τ /g µ instead of g τ /g e since only two of the ratios in Eq. (21) are independent. The ratio g τ /g µ is equivalent to Γ(τ → eνν)/Γ(µ → eνν) and recieves a negative correction from the one-loop graphs but not from the tree-level charged Higgs. Hence, the strength of the bound comes from the fact that both g τ /g µ and g µ /g e are high in the purely leptonic fits. However, HFAG also provides fits to g τ /g µ from Γ(τ → hν)/Γ(h → µν) where h = K, π, which are found to be low compared to the SM expectation. These ratios should also be affected by the one loop corrections but not the tree level term since we neglect the Higgs coupling to light quarks. The combined fit is then reported as [50] ( We use this number to constrain the one loop corrections at 95% confidence and g µ /g e to constrain the tree-level charged Higgs correction. This gives a weaker exclusion than the purely leptonic fits.

D. Heavy Quarkonium Decay
As will be seen, our results allow for rather light new neutral scalars which may have a non-negligible coupling to heavy quarks. For scalars with masses below 10 GeV we consider bounds from measurements of rare Υ decays. In particular, the branching fraction for constrain the coupling of new light scalars to b quarks [52]. Based on Refs. [53,54] we impose the following limits:  Table II, where the superscripts of the signal strength modifiers denote the production channels of the Higgs boson.
We perform a global χ 2 fit of the model predictions to the observed Higgs signal strengthŝ µ i . We assume each channel listed can be treated as an independent measurement, giving us 18 degrees of freedom from the LHC measurements as listed in Table II. Note that some analyses in the literature give separate fits to associated (V H) and vector-boson fusion (VBF) production. In a given model, they are not independent. In our treatment we use the combined fit, taking into account the fact that these two processes are determined by the scalar mixing parameter θ 12 for both W s and Zs. Similarly, some analyses have performed a fit to associated top quark (ttH) production, but at present these have large error bars and we choose to use the combined fit with gluon-gluon fusion (ggF ) when appropriate since the latter is dominated by the top quark loop. We require a predicted χ 2 value consistent with experiment at the 95% level.
In the framework of the A2HDM, the SM-like Higgs boson could decay to two lighter states when kinematically accessible, which then subsequently decay to four taus. An estimate of exclusion limits for such decays based on recast 3-lepton searches was made in Ref. [67]. They found that If the new particles are too light to decay to taus then we expect them to decay largely to muons, which is bounded by experiment to be less than ∼ 1 fb [70]. We include these limits as a constraint on our results.
The SM value of the muon anomalous magnetic moment has been calculated up to fiveloop contributions in QED [5] and two-loop in weak interactions [6][7][8]. Among the diagrams, the hadronic contributions are particularly challenging to reliably estimate [13]. The current In the SM, the Higgs contribution to a µ is suppressed by a factor of m 2 µ /m 2 h compared to the electroweak contributions [7,8]. However, the Higgs sector contributions to a µ can be considerably enhanced in a 2HDM. The significance of the a µ constraint for 2HDMs was emphasized in Refs. [15,73,74]. Representative Feynman diagrams in the A2HDM contributing to (g − 2) µ are depicted in Fig. 1. New Higgs bosons may contribute to (g − 2) µ at leading order at the one-loop level. However, for a Higgs boson mass larger than ∼ 10 GeV, dominant Higgs contributions to (g − 2) µ come from the two-loop Barr-Zee diagram with a fermion in the loop [75]. It is also known that in the CP-conserving 2HDM of Type X, a light pseudo-scalar, together with a large tan β value can explain the measured ∆a µ via such diagrams [15][16][17][18][19]. Recently, the importance of additional contributions arising from charged Higgs bosons in the A2HDM has been emphasized in Ref. [76]. We include these in our predictions.

IV. RESULTS OF THE NUMERICAL SCAN
We perform a scan over the free parameters of the Higgs potential and Yukawa couplings.
As mentioned above, without constraining our results Z 2 is fixed at 4π. We randomly generate points with a flat distribution over the range for each parameter allowed by unitarity and consistent with positive eigenvalues for the masses, specifically 2 , Without loss of generality we then choose |Z 6 | so as to guarantee that one mass sits at 125 GeV. We record generated points which pass all the constraints discussed above within 2σ of the experimental values, including the requirement that the χ 2 fit to LHC Higgs data is consistent within 95% bounds. We generate billions of points, most of which are discarded due to the stringent experimental constraints. We interpret our solutions only as representative since there may still be corners of parameter space not being fully sampled. if |A d | 40. Generally then, one must rely on an enhanced coupling to leptons with large 2 Note these conditions are necessary but not sufficient and we make additional checks to ensure that a generated point satisfies our requirements.
A l , which increases the contribution from all loops due to the muon coupling and doubly enhances the graph with a tau loop. This implies that light new Higgs states will typically decay to taus and any mixing of the 125 GeV Higgs with the scalar from H 2 will lead to significant deviations from an SM-like coupling to taus, which is strongly constrained by current LHC results. Hence, (g − 2) µ forces us into the large A l region with small mixing angle θ 12 , which in turn requires that Z 1 v 2 is close to (125 GeV) 2 and |Z 6 | is small. A light Higgs with large A l also favors small A d so as not to violate various flavor bounds.
Note that the effect of scalars versus pseudo-scalars on (g − 2) µ is not as simple in the A2HDM as in the more restricted models which have been considered in the literature [16][17][18][19]. At the one-loop level it is generically true that the pseudo-scalar gives a negative In the upper left panel of Fig. 2 we show the ranges for input parameters which pass our constraints. As discussed, Z 1 is largely fixed, varying within about 15% of the observed Higgs mass value m 2 h /v 2 = 0.26, and Z 6 is also small, 0.5, corresponding to the low-mixing limit. Z 3 is positive and can be as large as the unitarity bound at 8π. Z 4 is negative and Low • In the first region, which occupies most of the allowed parameter space, the pseudo-scalar GeV.
In all regions, the upper bound on the mass of the heavier Higgs is determined by the allowed mass splitting between H and A. This splitting is largely controlled by |Z 5 |, which is constrained by unitarity to be less than 2π. At the upper end of our allowed range for |A l | the pseudo-scalar mass can approach 130 GeV with a light charged Higgs helping to enhance (g − 2) µ . Two light neutral particles which decay to taus and couple to quarks strongly enough to account for the magnetic moment are largely excluded. Thus when A is light, the scalar H must usually be significantly heavier, such that it can decay instead to AA, to AZ, or to top quarks. This implies that typically the heavier neutral state is above 200 GeV. This bound also follows from the approximate degeneracy between the heavier neutral state and the charged Higgs. The latter is preferably heavier than the top quark to avoid constraints from top decays.
If |A l | is large enough to generate ∆a µ through the tau loop alone, then in principle the quark couplings can be arbitrarily small, making the tau search from direct production irrelevant. This can be taken as the relevant limit for the Type X model. In this case, however, without the sizeable contribution of top loops to (g − 2) µ , the pseudo-scalar must be lighter or |A l | even larger than the typical values shown in our plots. This parameter space is highly constrained by the measurements of the Zττ coupling and lepton universality.
When tau loops dominate the scalar terms give a negative contribution to ∆a µ which must be suppressed relative to the pseudo-scalar, so a large mass splitting is still required.
The second region, with very light scalars and intermediate mass pseudo-scalars, corresponds with relatively mild enhancement of the lepton coupling, i.e. |A l | between 5 and 35.
These scalars are light enough, 3 − 20 GeV, that their positive one-loop contributions to (g − 2) µ , which arise from the 1-loop graphs in Figure 1 can become larger than the two-loop terms which are dominant in the first region.
The third region represents the case where two-loop terms dominate, with a scalar mass in between the previous two cases. This is qualitatively somewhat similar to the first region with the roles of A and H reversed. However, this possibility requires a larger mass splitting than in the light A scenario for two reasons. One is that the new light scalar would need SM-like coupling to the top to generate a sufficient positive contribution, and thus the limits We compare |A l | and |A u | in the lower left panel in Fig. 2. For the light pseudo-scalars, with m A ∼ 100 GeV, |A l | and |A u | are anti-correlated. This is because, for fixed masses, high values of |A l | require a suppressed |A u | to avoid the upper R Bs→µ + µ − bound, while lower values of |A l | require |A u | ∼ 1 to give a sufficiently large contribution to ∆a µ . As m A is lowered, smaller values of |A u | are favored to avoid bounds from Higgs searches with the low mass accounting for the enhanced (g − 2) µ instead. For sufficiently small |A u | however, tau loops dominate and constraints from direct production, B s → µ + µ − and b → sγ become unimportant. In the lower right panel of Fig. 2 we plot |A u | against |A d |. There is an upper bound which is driven by the constraint from b → sγ. For very light scalars, couplings become constrained by searches for rare Υ decays to γµ + µ − and γτ + τ − . Since light scalars will decay almost entirely to these leptons, this puts a strong bound on their coupling to b-quarks, requiring |A d | < 1. One can also see from these plots that the Type-X model limit, is only a small fraction of the viable A2HDM parameter space.

B. (g − 2) µ in the CP-violating A2HDM
New sources of CP violation, which are generically possible in the A2HDM, will give rise to new electric dipole moment (EDM) contributions for the electron and the neutron.
Leading contributions to these dipoles are generated by Barr-Zee-type graphs which are directly analogous to those which are required to explain ∆a µ . Unless we have a CPpreserving symmetry which ensures that the EDM contributions cancel, they will typically be too large to agree with the experimental limit. In the CP violating case with EDM constraints imposed, the largest ∆a µ generated in our general scan is an order of magnitude too small to account for the lower limit of the measured value [35]. Nonetheless, although the CP-conserving limit most naturally allows one to generate a large ∆a µ in the A2HDM, it is interesting to consider how far one can move away from this limit and still account for the experimental data. 5 To investigate this we allow the parameters φ u , φ d , φ l , 2φ 5 , φ 6 and φ 7 to deviate from a common phase φ 0 where each δ i is randomly distributed within a restricted range. For |δ i | < 2 × 10 −4 there is essentially no constraint from EDMs. The distribution of masses and parameter ranges is qualitatively the same as in the exact CP-conserving case. If we increase the allowed deviations to |δ i | < 2 × 10 −3 , then the electron EDM becomes pertinent. Roughly two thirds of model points which would otherwise be allowed are excluded by the electron EDM [35].
At this level the neutron EDM predictions are still uniformly below the measured values.
For |δ i | < 2 × 10 −2 , approximately 98% of the sample after other cuts is excluded by the electron EDM, and a few percent of the remaining points are now excluded by the neutron EDM.
Surviving points can still be found across the larger range of individual δ i , which suggests that sufficiently small EDMs can be accommodated by cancellations of several terms even for relatively large CP-violating phases. However, these canceling conditions become increasingly rarefied as the CP-violating phases grow. Without some principle which enforces such cancellations, we regard them as highly fine-tuned. Thus, while an explanation for the observed (g − 2) µ in the A2HDM does not strictly rule out significant new sources of CP-violation, they appear highly unnatural in this model. We will leave further discussions to a later publication [35].

A. Deviations of the SM-like Higgs couplings
The solutions to the (g − 2) µ excess in the A2HDM rely on the relatively light new neutral states, which in turn modify certain SM observables. We first consider modifications to the SM Higgs properties at the observed 125 GeV resonance. The range of predicted coupling strengths κ i , defined by κ i ≡ g i /g SM i is shown in Fig. 3 compared with the χ 2 fit.
As expected, deviations are driven by the large values of |A l |, which are apparent in the potentially large modifications of κ τ . One can see that κ τ may even vanish without violating our bound on the combined fit to data, although this is disfavored by the τ τ channel itself.
The fit to LHC data directly limits these excursions, which in turn requires that the mixing angle θ 12 is small. This mixing determines κ V , which is thus constrained to be very close to 1. We find κ V > 0.998, except for the case of very light scalars, which may accomodate κ V > 0.982. In either case, this is much more restrictive than the current or even expected limits from direct fits to κ V . This is equivalent to a mixing angle between the two scalar states of sin θ 12 < 0.06 (0.19 for very light scalars.) Since mixing is constrained to be so small and since the enhancement of other Yukawa couplings cannot be as large as |A l |, we don't see large deviations in the other tree level couplings. In general κ b can vary by about 10% relative to the Standard Model. The effective coupling to γγ and to γZ is typically enhanced by a few percent, although the former can be as large as 40% high and the latter as large as 10% high. In the case of light scalars only 1 − 5% deviations are expected for the γγ and γZ couplings. The coupling to gluons is predicted to be within a few percent of the Standard Model value with ∼ 10% deviations possible with light scalars.
With these relatively mild deviations, the variance in χ 2 values strongly tracks the fit to τ τ data. The best fit value is χ 2 8, well within the 1σ expectation for the number of degrees of freedom (18). For comparison the Standard Model value in our fit is χ 2

13.
The total width of the SM-like Higgs boson can differ in our results in the range 3.5 − 7 MeV. The upper range comes from decays to pseudo-scalar pairs. As discussed above, we limit these to a 25% branching fraction based on three-lepton searches for h 1 > 30 GeV and lighter masses are constrained by published results. The χ 2 fit by itself would rule out values higher than ∼ 50%.   Fig. 6. The E989 experiment at Fermilab is expected to improve experimental errors on the (g − 2) µ measurement by a factor of four [77].
In the upper-right panel of Fig. 6, we display the predicted values for δg L,R τ . The magnitudes of the left and right-handed couplings are strongly correlated with δg L τ > 0 and δg R τ < 0. This is opposite to the small shifts from Standard Model predictions found experimentally. In particular, since g R τ is high by roughly one standard deviation, it provides an important constraint on our model. In the lower panels of Fig. 6, we display the predicted coupling shifts as a function of |A l |. One can see that the measured value of g R τ becomes a significant bound for |A l | > 50 and excludes a significant fraction of otherwise viable points.
Since δg L τ is highly correlated with g R τ and the experimental fit for the latter is closer to the SM, δg L τ does not impose any relevant additional constraint. must typically be massive enough to decay to the lighter Higgs state or to tops in order to limit the decay to taus and to suppress negative contributions to (g − 2) µ , which is strongly constrained by LHC searches.

VI. SEARCHES FOR NEW HIGGS STATES
In Fig. 7 we present the leading branching fractions for the decay of the heavier neutral For all cases we find that the heavier new neutral state may decay to tops with at most 30% branching fraction and to taus with at most 50%, as seen in the two middle panels. In general the maximum branching fraction to taus or to h 1 h 1 declines at higher masses, while the fraction for decay to Zh 1 (H or A) increases. Decays to bb are relatively rare since only a mild enhancement over SM-like couplings is allowed. As seen in the two lower panels in Fig. 7, the charged Higgs can have a large branching fraction to τ ν if |A l | is large which tends to dominate at lower m H + with light pseudo-scalars. At higher masses, and generally in the case of light scalars, the charged Higgs will decay primarily to W h 1 . After accounting for these two decays any remaining fraction will be largely H + → tb, which may reach about 20%.
In Fig. 8 The gluon-fusion channel will continue to be a leading search mode at the LHC. However, given the wide range of the possible κ g values as seen in Fig. 8, more quantitative evaluation would be needed to draw a conclusion for the observability of those states at the LHC. The b-associated channel can be mildly enhanced when |A d | > 1, but the t-associated channel will be suppressed due to |A u | < 1.
A second production mechanism is through the electroweak pair production channels qq → W ± → H ± A, H ± H, H ± h, and qq → Z → AH, Ah.
The charged current channel W ± → H ± A is via the pure SU(2) gauge interaction independent of other model parameters, and the process Z → AH is also close to the full EW strength due to the small h − H mixing. By the same token, production of h via these modes is quite small. These channels have the advantage of being present even if the quark couplings are small. It is conceivable to search for Z * → AH → 4τ /4τ Z/6τ , and similarly for W + * → H + h 1 → 3τ /4τ W + could be pursued. W + * → H + h 2 is also possible, although typically it will decay to a more complicated final state. Estimates of the sensitivity of such searches in the Type X model were recently published in Ref. [79]. The Z * channel can also be exploited at a lepton collider, particularly at the ILC which, with an energy of 500 GeV or 1 TeV could kinematically access a large part of the interesting parameter space.

VII. SUMMARY AND CONCLUSIONS
In this work we have explored the viable parameter space of the CP-conserving A2HDM, outlined in Sec. II, which can account for the observed value of ∆a µ . As studied in great detail in Sec. III, the model is significantly constrained by experimental measurements of the T -parameter, by flavor changing neutral decays B s → µ + µ − and b → sγ, by precision measurement of the Zττ coupling and τ decays, by the measured production and decay channels for the SM-like Higgs boson, and by LHC searches for exotic Higgs decays. We find in Sec. IV that it is possible to satisfy all our requirements in certain parameter regions with characteristic features of the new Higgs bosons. Our results can be summarized as follows: • There are three distinctive mass regions that are viable as given in Eq. (29). These are distinguished by the mass values of H and A and by their dominant contributions to (g − 2) µ , as illustrated in Fig. 1.
• The Yukawa couplings of the new scalars to top quarks can play a leading role in generating the necessary enhancement. This in turn allows for a wider range of masses compared to the Type X model, including heavy pseudo-scalars, as well as very light scalars, with m H 10 GeV or m A ∼ 10 GeV.
• The new Higgs states all have large couplings to τ 's due to the necessarily large value of A l . This can have observable effects in precision measurements such as Z → τ τ , and universality in (semi-)leptonic decays, and it dominates the phenomenology of collider searches (Sec. V).
• Light new states (h 1 = A, H) decay to tau pairs. Heavier states will typically decay to either a pair h 1 h 1 or to h 1 Z, although a significant fraction may decay directly to a pair of taus or to tops. Thus, further searches for exotic Higgs decaying to taus or to taus plus Z in the final state can potentially discover the new particles of this model (Sec. VI).
• The leading production at the LHC may still be associated with the heavy quarks, such as gg → H/A, bb → H/A, gb → tH ± . New particles can also be pair produced through Z * → AH and W * → H + (A/H) mainly via gauge interactions. They can be competitive in the search when A/H is light (Sec. VI).
• The obtained solutions above may lead to significant modifications to the SM measurements, such as the deviation of the SM-like Higgs boson couplings, and a new decay channel to 4 taus (Sec. V A).
• Current LHC searches have already begun to impact the allowed parameter space and the model is sensitive to future searches.
Our work demonstrates the viability of a much larger class of theory related to the 2HDM, which can explain (g − 2) µ , with the Type X model as a special case within this framework.
The model presents interesting coupling patterns, new production channels for both the SM-like Higgs boson and the new states at hadron colliders, and a richer phenomenology at low energies.