Scrutinizing a hidden SM-like gauge model with corrections to oblique parameters

In view of the recent high precision measurement of the Standard Model W boson mass at the CDF II detector, we compute the contributions to the oblique parameters S, T and U coming from the two additional Higgs doublets (one inert and one hidden) as well as the hidden neutral dark gauge bosons and extra heavy fermions in the gauged two-Higgs-doublet model (G2HDM). While the effects from the hidden Higgs doublet and new heavy fermions are found to be minuscule, the hidden gauge sector SU(2)H×U(1)X\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$SU(2)_H \times U(1)_X$$\end{document} with gauge coupling strength ≳10-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gtrsim 10^{-2}$$\end{document} and gauge boson mass ≳100\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gtrsim 100$$\end{document} GeV can readily explain the W boson mass anomaly but is nevertheless excluded by the dilepton high-mass resonance searches at the Large Hadron Collider. On the other hand, the new global fits to the oblique parameters due to the new W boson mass measurement can give discernible impacts on the mass splitting and mixing angle for the inert Higgs doublet in G2HDM. We also study the impact on the signal strength of diphoton mode of the 125 GeV Higgs boson h→γγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h \rightarrow \gamma \gamma $$\end{document} and the detectability of the yet to be observed process h→Zγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h \rightarrow Z \gamma $$\end{document} at the High Luminosity Large Hadron Collider. Current constraints for the dark matter candidate W′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W^\prime $$\end{document} including the dark matter relic density, dark matter direct detections and invisible Higgs decays are also taken into account in this study.


Introduction
Based on the data with an integrated luminosity of 8.
which is ∼ 7σ away from the SM prediction from electroweak (EW) global fits [2] m W (SM − Global Fits) = 80, 359.1 ± 5.2 MeV/c 2 . (2) The CDF result (1) also represents a ∼ 3σ deviation from other direct measurements from the more recent ATLAS [3] and LHCb [4] experiments. This immediately stirs a great deal of excitement in the field and triggers many subsequent studies. This recent measurement has significant implications for the electroweak precision global fit [5][6][7]. Specifically, the fit can accommodate the recent CDF m W measurement by incorporating a large and nonzero value of the U parameter or a large positive value of T parameter in the case of U = 0. These two parameters T and U , together with a third parameter S, are well-known in the community of electroweak precision test as oblique parameters [8]. We will briefly review them later in Sect. 3. We note that the deviation of the CDF W boson mass measurement with the global fit only shifts slightly from 7σ to 6σ if a theoretical calculation tool used by CDF is updated by a more recent version [9]. More important and optimistic view is that the global fit could entail new physics (NP) beyond the SM (BSM). While combined result of the measurements from LEP, Tevatron and ATLAS are still lacking, pending on evaluation of uncertainty correlations [1,10], and the new CDF result of m W is needed to be independently confirmed, BSM enthusiast has already offering various NP interpretations of the new CDF result. See for example, Refs. [11][12][13] for extra U (1) implications, [14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29] for extended or composite Higgs sectors, [30,31] for SMEFT, [32,33] for lepto-quark, [34,35] for dark matter (DM) models, [36][37][38][39][40] for low energy supersymmetry, and [41][42][43][44] for grand unification, etc. Another interesting point to support the BSM physics as pointed out in [45] is that the hadronic uncertainties in the fine structure constant and hadronic vacuum polarization that affects the SM W boson mass and muon anomalous magnetic dipole moment respectively are anti-correlated with each other.
In this work, we study the impact of the new CDF result to the extra particle mass spectra in the gauged two-Higgs-doublet model (G2HDM) first proposed in [46] and explored further in [47][48][49][50][51][52][53][54][55][56]. G2HDM is a gauged DM model based on the extended electroweak gauge group G G2HDM = SU (2) The extra gauge group SU (2) H × U (1) X represents the dark gauge sector interacting feebly with the visible SM sector, in the sense that the new gauge couplings g H and g X are much smaller than the SM electroweak gauge couplings g and g , as suggested by our recent detailed studies [54,55]. Additional Higgses and heavy fermions must be included in G2HDM to make it phenomenologically viable and free from gauge and gravitational anomalies. The crucial idea of G2HDM is that the usual two Higgs doublets (H 1 and H 2 ) in general two-Higgs-doublet model (2HDM) is lumped together in a 2 dimensional spinor representation H = (H 1 H 2 ) T of a hidden local SU (2) H gauge group. Hence we christened the model as gauged 2HDM. Nevertheless, a hidden doublet H of SU (2) H is necessarily introduced so as to provide realistic Yukawa couplings and mass spectra for the extra vector bosons and fermions in the model. Effectively we have a tailor-made 3HDM in the scalar sector.
One distinctive feature of the model is that there is no need to impose an ad hoc discrete Z 2 symmetry to stabilize the DM. There is a hidden h-parity in the model [52], admitted readily once one writes down all possible gauge invariant and renormalizable interactions, that will guarantee the lightest h-parity odd particle to be the DM candidate, provided that it is not broken spontaneously. 1 Another interesting feature is that the new gauge bosons W ( p,m) ≡ (W 1 ∓ i W 2 )/ √ 2 are electrically neutral and don't mix with the SM W ± bosons. They are also h-parity odd and can be the DM candidate. A scenario of sub-GeV low mass W ( p,m) as DM was studied in [54,55]. In this work, we will turn our attention to the scenario of W ( p,m) as DM candidate with a wider mass range.
The G2HDM model is renormalizable and anomaly-free. Furthermore, the model has a range of intriguing phenomenological implications, including the emergence of new signals from extra gauge bosons at the LHC [50], the enhancement of di-Higgs production [51], and the possibility of charged lepton flavor violation processes [56]. Recent studies also suggest that the G2HDM model has great potential to be probed in both dark matter and dark photon search experiments [54,55]. In addition, the DM candidate W ( p,m) has non-abelian and complex nature which can make it to be both asymmetric and self-interacting DM. These features of the W ( p,m) in G2HDM as DM candidate may be interesting to explore in the future so as to resolve some well-known issues such as the too-big-to-fail and core-cusp problems encountered in the conventional WIMP cold dark matter scenario [57] in astrophysics and cosmology.
In the next Sect. 2, we will briefly review the particle content in the minimal G2HDM [54,55], their masses and interactions that are relevant to our study. In Sect. 3, we remind ourselves by reviewing the Peskin-Takeuchi oblique parameters [8] that not only entered in the art of global fit analysis which correlates all electroweak observables but also provided important constraints on NP models. In Sect. 4, we compute the new contributions to the oblique parameters from the mixings of the neutral gauge bosons (Sect. 4.1), extended Higgs sector (Sects. 4.2 and 4.3) and extra heavy fermions (Sect. 4.4) in G2HDM. In Sect. 5, we present and discuss our numerical results by including the new contributions to the oblique parameters to compare with the updated global fits [6]. In addition, we will take the opportunity in this section to explore the detectability of the process h → Z γ in the High Luminosity Large Hadron Collider (HL-LHC). Conclusions are given in Sect. 6. We reserve an appendix for the analytical expressions of the one-loop amplitudes of the two processes h i → γ γ (i = 1, 2) and h i → Z j γ (i = 1, 2; j = 1, 2, 3) in G2HDM, where h 1 and Z 1 are identified as the SM 125 GeV Higgs scalar (h) and 91 GeV Z vector boson respectively.

Model setup
In this section, we will briefly review the minimal G2HDM. The original model was introduced in Ref. [46], and various refinements [48,49,52] and collider implications [47,50,51] were pursued subsequently with the same particle content as the original model. As advocated recently in [54,55], we will drop the triplet field H of the extra SU (2) H . This is mainly due to the fact that realistic mass spectra for all particles in G2HDM can be achieved without including the triplet. This important fact has been overlooked in the original work [46] because a couple of terms were missing in the scalar potential. Additionally, removing the triplet field results in a reduction of 6 arbitrary parameters in the scalar potential.
We will refer this as the minimal G2HDM in what follows. The quantum numbers of the matter particles in G2HDM Spin 1/2 Fermions: Quarks In G2HDM, the SM W boson does not mix with W and its mass is the same as in SM: m W = gv/2. However in general the SM Z SM will mix further with the gauge field W 3 associated with the third generator of SU (2) H and the U (1) X gauge field X via the following mass matrix: where The real and symmetric mass matrix M 2 Z in (11) can be diagonalized by a 3 by 3 orthogonal matrix O G , i.e.
where m Z i is the mass of the physical fields Z i for i = 1, 2, 3. We will identify Z 1 ≡ Z to be the neutral gauge boson resonance with a mass of 91.1876 GeV observed at LEP [58]. The lighter/heavier of the other two states is the dark photon (γ )/dark Z (Z ). These neutral gauge bosons are h-parity even in the model, despite the adjective 'dark' are used for the other two states. We note that these neutral gauge bosons can decay into SM particles and thus they can be constrained by experimental data, including the electroweak precision measurement at the Z pole physics from LEP, searches for dark Z and dark photon at colliders, beam-dump experiments, and astrophysical observations. The DM candidate considered in this work is W ( p,m) , which is electrically neutral but carries one unit of dark charge and chosen to be the lightest h-parity odd particle in the parameter space. In The mixing angle θ 1 is given by The masses of h 1 and h 2 are given by The mixing angle θ 2 satisfies and the mass of D is In the Feynman-'t Hooft gauge the Goldstone fieldG * (G) has the same mass as the W p (W m ) which is given by (13). Finally the charged Higgs H ± is also h-parity odd and has a mass One can do the inversion to express the fundamental parameters in the scalar potential in terms of the particle masses [54,55]: From (13), we also have Thus one can use m h 2 , m W , m D , m H ± , θ 1 and θ 2 as input in our numerical scan. The Feynman rules can be straightforwardly derived by rewriting the two Lagrangians (6) and (7) using the physical fields mentioned above. The h-parity odd particles in G2HDM are W ( p,m) , D ( * ) ,G ( * ) , H ± , and all new heavy fermions f H . Among them, W ( p,m) , D ( * ) , and ν H are electrically neutral and hence any one of them can be a DM candidate. The connector sector linking the SM particles and the DM W ( p,m) consists of the h-parity even or odd particles. Specifically, we have γ , Z i (i = 1, 2, 3) and h n (n = 1, 2) in the s-channel, and new heavy fermions f H s or even the DM W ( p,m) itself in the t-channel and/or u-channel. Phenomenology of a complex scalar D ( * ) as DM was studied in detail in [52,53] and for low mass W ( p,m) as DM, see [54,55]. For further details of G2HDM, we refer our readers to the earlier works [48,49]. Phenomenology of a new heavy neutrino ν H as DM in the model, which is necessarily implying both DM and neutrino physics, has yet to be explored.

Oblique parameters and W boson mass shift
As is well known, the oblique parameters S, T , and U [8] represent the most important electroweak radiative corrections since they are defined by the transverse pieces of the vacuum polarization tensors of the SM vector gauge bosons. They are process independent whereas the other vertex and box corrections are necessarily attached to the particles in the initial and final states in the elementary processes in high precision experiments.
The vacuum polarization tensor i μν I J (q) involving the SM gauge bosons I and J has the following decomposition The form factor I J (q 2 ) needs no concern to us since at high energy experiments like LEP I and II where electroweak precision measurements were carried out, q μ will dot into the helicity spinors of light leptons and will give vanishing results. The vacuum polarization amplitude I J (q 2 ) has the following expansion The oblique parameters S, T and U are defined with an overall factor ofα =ê 2 /4π extracted out in front as [8] The W boson mass shift can be related to the oblique parameters according to [8]

New contributions to the oblique parameters in G2HDM
In this section, we compute the contributions to the oblique parameters from all the new particles introduced in G2HDM. First, we will handle the tree level contributions to the oblique parameters coming from the mass mixing matrix for the three massive neutral gauge bosons given in (11).

Contributions from the tree level mixings of Z SM , W 3 and X
The kinetic and mass mixings of an extra U (1) boson with the SM B and W 3 gauge fields violate the custodial symmetry in the SM at the tree level and can give rise to the effective shift of the oblique parameters [60] (see also [13]). The neutral gauge boson mass mixings (11) of Z SM , W 3 and X in G2HDM also violate the custodial symmetry in the SM and can give rise to a non-vanishing vacuum polarization amplitude Z Z (q) at the tree level. The other 3 vacuum polarization amplitudes vanish at tree level in the model. For g H , g X g, g , from the Feynman diagram in Fig. 1, one Tree level diagrams (to second orders of g H and g X ) that contribute to the oblique parameters from the mass mixings (11) between the Z SM , W 3 and X in G2HDM The W mass shift from the tree level mixing of the neutral gauge bosons spanned on the (m W , g H ) plane. Here we set g X = 0. The solid green, solid blue and solid red lines represent the W mass shift of 1 MeV, 20 MeV and the one measured at the CDF II, respectively. The light gray shaded region is the excluded region ( m Z > 2.1 MeV) from the measurement of Z boson mass at LEP [58]. The orange region is the excluded region from the di-lepton high mass resonance search from ATLAS [62] obtains where m 2 W is given by (13) and m 2 X is the 33 entry of (11), namely The Z boson propagator is modified as This modified propagator implies a mass correction δm 2 Z and a wave function renormalization constant Z for the Z boson, viz.
The mass shift for the physical Z field is then given by 3 One can then deduce the effective oblique parameters by using the EFT approach with the following effective Lagrangian [61] where L SM is the SM Lagrangian and L NP is the Lagrangian from NP effects which is defined as HereF μν andẐ μν represent the usual Abelian field strengths, is the W (Z ) mass in the SM and the coefficients A, B, C, G, w and z represent either the wave function renormalizations or mass corrections from NP. Since neutral gauge boson mass mixings are tree-level effects, this effective Lagrangian approach seems more appropriate to deduce the oblique parameters. We therefore identify 4 All the other effective parameters A, B, G and w in Eq. (43) are zeros at tree level in G2HDM. Using Eq. (2) in [61], we then obtain the tree level oblique parameters as Eq. (12) and Z Z (m 2 Z ) by Eq. (37). 4 One can also extract C and z directly from the Z mass and neutral current interaction terms in the Lagrangian. In particular, where the rotation matrix O G and physical Z boson mass m Z1 are obtained by diagonalizing the mass matrix in Eq. (11). 2 (lower right panel) from the Higgs signal strength measurement at the LHC. We note this upper bound for sin 2 θ 1 (or cos 2 θ 1 ) changes if m h2 is close to m h1 [65] The corresponding W boson mass shift at the tree level can be computed using (36).
In Fig. 2, we show the contours of W boson mass shift from the contribution of the tree level mixings between the neutral gauge bosons projected on the plane of m W and g H . Here we have fixed g X = 0 for simplicity. Including contributions from both W and X do not change our conclusions in a significant way. In m W m Z regions and for relative large values of g H , one can reach the W boson mass shift ( m W ≈ 75 MeV) measured at the CDF II. However these large values of g H also yield a large Z boson mass shift which is in conflict with the great precision measurement of the Z mass at LEP [58] m Z (LEP) = 91.1876 ± 0.0021 GeV. (46) In the heavy mass region of W boson, the gauge coupling g H is also constrained to be g H 10 −2 from the di-lepton high mass resonance searches at ATLAS [62]. In the allowed region, the W boson mass shift is relatively small from these extra massive neutral gauge bosons in G2HDM. In other words, the mass mixing effects from (11) to the W boson mass shift is not significant. We note that for lower mass regions of W and X bosons (m W ,X < m Z ), the oblique parameters in Eqs. (44,45) calculated from the EFT approach may not be justified. However in these low mass regions, the gauge couplings are constrained to be smaller (g H,X 10 −3 ) [54,55], one expects the contribution to the W boson mass shift from these mixing effects is also negligible.
In the next three subsections, we will turn to one-loop contributions to the oblique parameters from all other new particles in G2HDM. We will treat Z 1 = Z ≈ Z SM and work in the Feynman-'t Hooft gauge.

Contributions from the SU (2) H doublet H
Since H is a SU (2) H doublet but a SU (2) L singlet, its contribution can only arise from the mixing effects between h SM and φ 2 and therefore similar to the singlet extension of the SM. The relevant Feynman diagrams are depicted in Fig. 3.
For the T parameter, we find The above result agrees with [63,64]. Compact expressions for the S( H ) and U ( H ) parameters can also be obtained using their definitions given in (33) and (35). and As one expects, all three S( H ), T ( H ) and U ( H ) vanish as m h 2 → m h 1 or θ 1 → 0. In this limit, from Eq. (16), we must require λ H = 0, λ H = λ and μ 2 H = μ 2 . The scalar potential with λ H = 0, λ H = λ and μ 2 H = μ 2 for the three doublets H 1 , H 2 and H can then be simplified to manifest a custodial symmetry of SO (12), preventing the deviation of the oblique parameters from nil for the degenerated h 1 and h 2 . The other two couplings λ H and λ H in the scalar potential break this custodial symmetry but they are not relevant for the masses of h 1 and h 2 .
We note that the above expressions (47), (48), and (49) are written down assuming the lighter h 1 is identified as the observed Higgs boson. If instead the heavier h 2 is identified to be the observed Higgs boson, analogous expressions can be obtained by setting sin 2 θ 1 → cos 2 θ 1 and m h 1 ↔ m h 2 in (47), (48), and (49).
We will first study the case where h 1 is the observed Higgs boson. In the upper left panel of Fig. 4, we show the oblique parameters divided by sin 2 θ 1 calculated from (47), (48), and (49) as a function of m h 2 . One can see that S is positive value while T is opposite, and | U | is much smaller than | S| and | T |. In the upper right panel of Fig. 4, we show the contours for the W mass shift projected on the (m h 2 , sin 2 θ 1 ) plane. The W boson mass shift is seen to be negative for Similar plots in the lower left and right panels of Fig. 4 are made for the case where the heavier h 2 is the observed Higgs boson. In this case, the oblique parameters flip their signs as compared with the previous case as shown in the lower left panel in Fig. 4. As a result, a positive W boson mass shift is obtained as shown in the lower right panel in Fig. 4. The mass shift can be large in the large cos 2 θ 1 and small m h 1 regions.
Due to the Higgs data at the LHC which require sin 2 θ 1 (or cos 2 θ 1 ) 0.2 (from the Higgs boson coupling modifier κ Z = 0.99 ± 0.06 [66]) for the case where h 1 (or h 2 ) is the observed Higgs boson, the W boson mass shift from this contribution is constrained to be relatively small in both cases and thus makes it difficult for the hidden SU (2) H doublet H to explain the CDF W boson mass anomaly.

Contributions from the inert Higgs doublet H 2 in G2HDM
The second Higgs doublet H 2 in G2HDM plays a similar role as the inert Higgs doublet in I2HDM [67][68][69][70]. However the neutral component H 0 2 of H 2 can be treated as a complex scalar field in G2HDM instead of decomposing into (S + i A)/ √ 2 in I2HDM, where S and A are the scalar and pseudoscalar fields. In addition, as already given in (17) For the T parameter, we got where m W is the mass the Goldstone bosonG which is absorbed by the new gauge boson W ( p,m) of SU (2) H , and the function F(m 1 , m 2 ) is defined as Again, analytical formulas for the S(H 2 ) and U (H 2 ) parameters can be obtained using their definitions given in (33) and (35). and where Note that in the limit of m 2 → m 1 , G(m 1 , m 1 ) = 12.
In the pure inert limit of θ 2 → 0, D → H 0 2 and we simply have lim lim The above expressions of (55) and (56) are consistent with the inert Higgs results [68].

Contributions from the new heavy fermions in G2HDM
Since all the new heavy fermions f H in G2HDM are SU (2) L singlets, they don't interact with the charge W ± bosons. Under the assumption that g H , g X g, g , the heavy fermions interacts with both the SM γ and Z are vector-like described by the following Lagrangian where we have dropped terms that are proportional to g H or g X . Thus we have where Q Q is the oblique loop amplitude γ γ with both the color and electric charge factors trimming off [71], i.e.
with E ≡ 2 −γ E +log 4π . Using the above expressions (59)-(63), we can demonstrate easily that all the oblique parameters from the heavy fermions f H s in G2HDM vanish: The non-trivial leading contributions from the heavy fermions in G2HDM to the oblique parameters start at order g 2 H /16π 2 and g 2 X /16π 2 , which we are neglecting in this study. The rationale behind the vector-like results in Eqs. (59)- (62) and (64)  To summarize, we have computed all possible nontrivial sources of new physics effects to the oblique parameters S, T and U in G2HDM under the approximation of g H , g X g, g . While the tree level mixings of the neutral gauge bosons and the extra heavy fermions in G2HDM contribute to the oblique parameters are of order g 2 H,X and g 2 H,X /16π 2 respectively and therefore not significant, there are new contributions from the inert doublet H 2 and the H 2 -H mixings, as well as from the hidden doublet H through the H 1 -H mixings. These new contributions from the extended scalar sector in G2HDM are of order g 2 /16π 2 and g 2 /16π 2 which are the same order as the SM one-loop contributions. In the next section, we will focus on the detailed numerical analysis of these new contributions.

Numerical results
In this section we present the numerical results in the light of new W boson mass measurement at the CDF II. Among the h-parity odd particles, we require W to be the lightest particle so that it can be a DM candidate in this model. We propose two setups of scan based on the mass of the DM, one is the light DM mass scenario and another is the heavy DM mass scenario. The parameter space setup for these two scenarios is given in Table 1. We sample the parameter space in the model using MCMC scans emcee [72]. For the light DM mass scenario, m W , M X and g X are scanned in the log scale, while the rest are in the linear scale in both scenarios except M X and g X are fixed in the heavy DM mass scenario. We also assume that the lighter h 1 is the observed Higgs boson and all heavy fermion masses are degenerated and fixed to be 3 TeV.
We closely follow the analysis for the current constraints in the model from [54,55]. In particular, we take into account  [75] for a recent review). To take into account the new W boson mass measurement at the CDF II, the recent global fit values for the oblique parameters due to NP contributions is adapted from [6], which are given as S = 0.005 ± 0.096, We also take into account the constraints from DM searches including the DM relic density h 2 = 0.120 ± 0.001 measured from Planck collaboration [76], DM direct detections from CRESST III [77], DarkSide-50 [78], XENON1T [79,80], PandaX-4T [81] and LZ [82], and the Higgs invisible decays constraint from the LHC [83]. We note that, due to the kinematical forbiddance, the Higgs invisible decays constraint is not applied for the heavy DM mass scenario. The branching ratio of invisible Higgs decay is given in [54,55]. We use micrOMEGAs package [84] to calculate the DM relic density and the DM-proton scattering cross section. We note that the production cross section of the monojet signals pp → W p W m j in the model is small due to the smallness of the gauge coupling g H [54,55] and hence evading the current constraint from the LHC [85,86]. Thus we do not include the mono-jet constraint in our analysis.
Hereafter, we denote the scan without the DM constraints as CDF-2022 and with the DM constraints as CDF-2022+DM. To see the impacts due to the CDF II W boson mass measurement, we employ other scans with the old global fit values for the oblique parameters taken from the Particle Data Group (PDG) [58] which are given as S = −0.01 ± 0.1,

Heavy DM mass scenario
In Fig. 6, we show first the favored regions from CDF-2022 (orange regions) and PDG-2021 (gray regions) spanned on the parameter space without imposing the DM constraints.
On the left panel of Fig. 6, we project the favored regions on the plane of the mass splitting ( m ≡ m H ± − m D ) and the mixing angle θ 2 . The mixing angle θ 2 is allowed to be either a nearly maximal mixing region (θ 2 ∼ −π/2) or a tiny mixing region (θ 2 5 × 10 −3 ). Since the relation between θ 2 and g H is given as the upper bound on θ 2 is due to the upper bound on the gauge coupling g H which is from the constraints of the Z mass shift and the di-lepton high mass resonance search at the LHC as shown in Fig. 2. For the tiny mixing θ 2 region, the impact from the new W mass measurement at CDF II is significant. In particular, within 2σ favored region, CDF-2022 prefers a large mass splitting between the charged Higgs and dark Higgs, while the PDG-2021 prefers a smaller mass splitting and even allows the degenerated case. On the other hand, for the nearly maximal mixing region, both CDF-2022 and PDG-2021 allow the degenerated mass between the charged Higgs and dark Higgs. However CDF-2022 still allows for a larger mass splitting as compared with the region favored by PDG-2021. On the right panel of Fig. 6, we project the favored regions on the plane of the charged Higgs and W masses. The favored regions from CDF-2022 and PDG-2021 are almost the same. For the region of m W 1.4 TeV, the favored region has a thin cigar shape indicating the charged Higgs and W Similar to the well known WIMP DM scenario, the DM candidate W is kept in the chemical equilibrium with the SM thermal bath via its 2 ↔ 2 annihilations before starting to freeze-out due to the expansion of the universe. The DM relic abundance can be determined by solving the Boltzmann equation for the evolution of the DM number density which heavily influenced by the 2 ↔ 2 annihilations between the DM and SM particles.
In G2HDM, beside the standard annihilation of the W ( p,m) pair to pairs of SM particles, the coannihilationmutual annihilation of multiple h-parity odd species-to pairs of SM particles can be also occurred. The Feynman diagrams for the DM annihilation and coannihilation processes are depicted in Fig. 7. In the coannihilation limit, the χ i ↔ χ j with χ i = {W , H ± , D} and i = j, is rapid until after annihilations decouple. The effective annihilation rate in this limit can be given by [87] where σ i j is the annihilation cross section, x = m W /T with T is the temperature, i = (m i − m W )/m W , g i counts the number of internal degrees of freedom and For the heavy DM mass scenario, we found out that the main contribution that yields the DM relic density observed at Planck Collaboration is the coannihilation channels. The annihilation processes are subdominant because the cross sections are suppressed due to the smallness of the new gauge couplings g H and g X . Figure 8 shows the mass difference between the DM W and h-parity odd particles H ± and D when the DM constraints are included. The mass differences (m H ± − m W )/m W and (m D − m W )/m W are required to be O(10 −2 − 10 −1 ) within 2σ region. The CDF-2022+DM prefers a larger region of the mass differences while the PDG-2021+DM can extend to a lower region. Figure 9 shows the favored regions on ( m, θ 2 ) (left panel) and (m W , m H ± ) (right panel) planes when the DM constraints are taken into account. As compared with the results without the DM constraints (shown in Fig. 6), the nearly maximal mixing regions of the angle θ 2 are not favored anymore. The degeneracy of the charged Higgs and dark Higgs (i.e. m = 0) favored by the PDG-2021 is no longer favored when the DM constraint is included. Within   Fig. 10 The 1σ (dark blue) and 2σ (light blue) favored regions of the data from CDF-2022+DM projected on the plane of the DM mass and spin independent DM-proton scattering cross section for the heavy DM mass scenario. The dark and light gray regions are the exclusion from PandaX 4T [81] and LZ [82] experiments. The dashed green, red and purple lines represent future sensitivities from the DM direct detections at XENONnT [88], DarkSide-20k [89] and DARWIN [90], respectively. Green region is the neutrino floor background 2σ region, it requires the mass splitting to be in the range of 33 GeV m 78 GeV for PDG-2021+DM and 72 GeV m 100 GeV for CDF-2022+DM. The mixing angle is also required to be smaller θ 2 2 × 10 −3 . The cigar shape of the favored region on the right panel of Fig. 9 is due to the happenstance of DM coannihilation i.e. (m H ± − m W )/m W ∼ O(10 −2 − 10 −1 ) as suggested already in the right panel of Fig. 6 for m W 1.4 TeV even before imposing DM constraints. The DM mass is predicted to be m W > 700 GeV while the charged Higgs mass is m H ± > 800 GeV within 2σ region. The CDF-2022+DM prefers a bit higher in the charged Higgs mass region as compared with the result from PDG-2021+DM.
In Fig. 10, we show the CDF-2022+DM favored region on the (m W , σ SI W p ) plane for the DM direct detection. Due to the constraint from the di-lepton high mass resonance search at the LHC [62], the gauge coupling is required to be g H 10 −2 in the favored DM mass region. It results in a small DM-proton spin-independent scattering cross section σ SI W p and thus the favored region lies far below the current limits from PandaX 4T [81] and LZ [82] (gray regions). Most of the CDF-2022+DM favored region overlaps with the neutrino floor background region. However a portion of the favored region predicted by the model can be probed by future DM direct detections at DarkSide-20k [89] and DARWIN [90].

Light DM mass scenario
Again we show first the favored regions from CDF-2022 and PDG-2021 for the light DM mass scenario in Fig. 11 without the DM constraints imposed. Similar to the heavy DM mass scenario, the significant difference between CDF-2022 and PDG-2021 is the favored region projected on the mass splitting between the charged Higgs and dark Higgs which is shown in the left panel of Fig. 11. In particular, within 2σ region, CDF-2022 prefers a mass splitting 60 GeV m 100 GeV, while PDG-2021 favors a smaller region 10 GeV m 75 GeV. However, unlike the heavy DM mass scenario, both CDF-2022 and PDG-2021 in this scenario disfavor the degeneracy of the charged Higgs and dark Higgs masses. The favored regions projected on (m W − θ 2 ) plane is shown in the right panel of Fig. 11. The results between CDF-2022 and PDG-2021 projected on this plane is slightly different. CDF-2022 prefers a region with the DM mass (m W 20 GeV within 2σ region) and the mixing angle (θ 2 0.15 rad within 2σ region) while PDG-2021 favors a bit lighter DM mass and larger mixing angle region. There is a lower bound on the DM mass (m W 0.02 GeV within 2σ region) which is due to the constraints from the dark photon searches [54,55]. In particular, because the relation between m W and m Z shown in (11), the lower bound on m Z due to NA64 [91], E141 [92] ν-CAL I [93] experiments can be translated to a lower bound on the DM mass. The constraints from the dark photon searches also put a strong upper limit on the gauge couplings g H and g X [54,55]. Since the relation between g H and θ 2 given in (67), the upper limit on g H induces the upper limit on θ 2 as shown in the right panel of Fig. 11.
Next, we take into account the DM constraints for this light DM mass scenario. Due to the scanned mass ranges of D, H ± and W , the mass splitting between DM and other h−parity odd particles H ± and D are relatively large and hence forbids the DM coannihilation channels. The annihilations of W p W m to the SM fermion pairs via s-channel with the mediation of Z , γ are dominant processes. However the cross sections of such processes are suppressed due to the smallness of the gauge couplings g H and g X except near the resonant regions where the mass of mediators are about twice DM mass. Because the mass of the lighter γ can not satisfy the resonant condition i.e. m γ ≈ 2 × m W , only the annihilation process with a heavier mediator Z near the resonance can provide a DM relic density observed at the Planck Collaboration.
We show the results with the DM constraints included in Fig. 12. In this scenario, since the DM physics is not significantly impacted by the charged Higgs and dark Higgs, the favored regions projected on ( m, m H ± ) plane as shown in Fig. 12 Same as Fig. 11 for the light DM mass scenario but with the DM constraints included. The dark (light) blue region represents the 1σ (2σ ) favored by the CDF-2022+DM, while the dark (light) gray region indicates the 1σ (2σ ) favored by the PDG-2021+DM Fig. 13 The 1σ (dark blue) and 2σ (light blue) favored regions of the data from CDF-2022+DM projected on the plane of the DM mass and spin independent DM-proton scattering cross section for the light DM mass scenario. The gray regions are the exclusion from CRESST-III [77], DarkSide-50 [78] and XENON1T [80] experiments. The dashed green, red, purple and orange lines represent the future sensitivities from DM direct detection experiments at NEWS-G [94], SuperCDMS [95], CDEX [96] and HeRALD [97], respectively. Green region is the neutrino floor background the left panel of Fig. 12 are almost unchanged as compared with the one without DM constraints as shown in the left panel of Fig. 11. On the other hand, due to the DM direct detection constraints from CRESST-III [77], DarkSide-50 [78] and XENON1T [80] experiments, the DM mass is required to be m W 2 GeV within 2σ region as shown in the right panel of Fig. 12. In the same panel, we also see that it is required to have a lower bound on θ 2 1.8 × 10 −2 rad. This is because of the DM relic density observed at the Planck Collaboration. A smaller θ 2 results in a smaller g H v/m W value as shown in (67), which can cause a smaller DM anni-hilation cross section and eventually lead to an overabundant DM in the universe.
In Fig. 13, we project the favored regions from the CDF-2022+DM data on the plane of the DM mass and spin independent DM-proton scattering cross section. The current constraints from CRESST-III [77], DarkSide [78] and XENON1T [80] experiments are shown as the gray shaded regions. These constraints put an upper limit on the DM mass in this scenario. We find out that a portion of the favored region at m W ∼ 1 GeV and σ SI W p ∼ [10 −44 − 10 −39 ] cm 2 can be probed by future DM direct detection experiments at NEWS-G [94], SuperCDMS [95], CDEX [96] and HeRALD [97]. Note also a small fraction of the 2σ favored region lies below the neutrino floor background (light green). s

5.3
The processes h → γ γ and h → Z γ at the LHC Given the favored parameter space obtained in previous two subsections for the heavy and light DM scenarios, one would like to study their impacts to the diphoton and Z γ channels from the 125 GeV Higgs decays and the future searches for these processes at the hadron collider. A precise measurement of the signal strength of these processes at the colliders could reveal the existence of new particles coupled to the Higgs (see e.g. [98,99]). For the diphoton channel, measurements with the current 13 TeV LHC data yield a signal strength via the gluon-gluon fusion production mechanism: μ ggh γ γ = 0.96 ± 0.14 [73]. For 14 TeV LHC with 3 ab −1 luminosity, the signal strength is expected to be measured with a ±0.04 uncertainty [100]. On the other hand, the Z γ channel has not yet been observed with the current data at the LHC [101,102]. With an expected uncertainty of 24% for the signal strength measurement at the 14 TeV LHC with 3 ab −1 luminosity, Fig. 14 The one-loop form factors for the individual particle running inside the loop of h → γ γ (left panel) and h → Z γ (right panel) processes as a function of | sin θ 1 |. The data points are taken in χ 2 < 5.99 region from CDF-2022+DM for the heavy DM mass scenario. The orange, green and red points represent the form factors for the SM fermions, charged Higgs, new heavy fermions respectively. The blue points represent the negative value of form factor for the W boson. Note that F γ γ (X ) ≡ F 1 (X ) defined in Eqs. (93,94,95,96) and F Z γ (X ) ≡ F 11 (X ) defined in Eqs. (73,80,81,87) in the Appendix. this Z γ channel can be observed with 4.9 σ significance at ATLAS experiment [100].
The analytical expressions for the one-loop amplitudes of h → γ γ and h → Z γ in the G2HDM are given in the Appendix. As compared with the SM predictions, the production rates of these two channels can be modified due to new contributions from the charged Higgs and new charged heavy fermions running inside the loops as well as effects from the mixing angle θ 1 between the SM Higgs boson h SM and the hidden scalar φ H in the model. We note again that in this analysis, we assume all new heavy fermions to be degenerated and fixed their masses to be 3 TeV. We select the data points from CDF-2022+DM with χ 2 = χ 2 − χ 2 min < 5.99 where χ 2 is the total χ 2 calculated from the Higgs data at the LHC, Z -boson mass measurement from LEP II, oblique parameters from Ref. [6], DM relic density from Planck collaboration and the one-side limits from Z searches, DM direct detection and Higgs invisible decay. Figure 14 shows the one-loop form factor for the individual particle including W ± , H ± , charged SM fermions (mainly top quark) and new charged heavy fermions running inside the loop of h → γ γ (left panel) and h → Z γ (right panel) processes as a function of | sin θ 1 | in the heavy DM scenario. The current data from ATLAS [66] puts an upper limit on the mixing angle | sin θ 1 | 0.35. In both h → γ γ and h → Z γ processes, the W ± form factor gets negative values but its magnitude is dominant. From the right panel of Fig. 9 and the left panel of Fig. 12, we learnt that our scan results prefer a heavy charged Higgs mass in both heavy and light DM mass scenarios. One expects the charged Higgs form factor is rather small. From Fig. 14, we have F γ γ,Z γ (H ± ) ∼ O(10 −7 − 10 −5 )/ GeV. The charged Higgs form factor is smaller than the SM top quark form factor in the h → γ γ process but it can be larger in the h → Z γ process. We also see that the new charged heavy fermions are suppressed due to the large VEV v in this data. We note that the SM fermion and new heavy fermion form factors for the h → γ γ process are about one order of magnitude larger than those for the h → Z γ process, while the W boson (as well as the charged Higgs) form factors stay almost the same for these two processes. In the limit of no mixing between the SM Higgs boson h SM and the hidden scalar φ H , i.e. sin θ 1 = 0, the new heavy fermion form factor would be vanished.  W ( p,m) . Furthermore, a lighter DM mass also results in a smaller v value, this would enhance the contribution of new heavy fermions in the loop. As shown in the left panel of Fig. 15 for the γ γ case, the contribution from the new heavy fermions can be comparable with the charged Higgs when | sin θ 1 | 0.1.
In Fig. 16, we project the favored data points from CDF-2022+DM for the heavy DM mass scenario on the (μ ggh γ γ , |δ h Zγ |) plane, where is the normalized deviation of the h → Z γ decay width from its SM value. The signal strength of diphoton production from the gluon-gluon fusion is given by where the decay width of h into two gluons in the model can be found in [46]. The deviation |δ h Zγ | can be large if the mixing angle θ 1 is large as shown in the left panel of Fig. 16. Similarly, μ ggh γ γ can be also impacted by this mixing effects, a larger mixing angle θ 1 results in a wider region of μ ggh γ γ . Moreover, due to the destructive interference between the charged Higgs and W boson contributions, the diphoton signal strength becomes smaller when the charged Higgs has a lighter mass as shown in the right panel of Fig. 16. Future measurements for the diphoton and Z γ production from the HL-LHC can probe a large portion of the CDF-2022+DM favored region in this heavy DM scenario. The low mass region of charged Higgs can be probed by the expected diphoton signal strength measurement with a 4% uncertainty at the HL-LHC [100].
Similar plots as Fig. 16 but for the light DM mass scenario are shown in Fig. 17. Similar to the heavy DM mass scenario, both |δ h Zγ | and μ ggh γ γ can be impacted by the mixing effects and the charged Higgs mass. We note that, in this scenario, Higgs boson is kinematically allowed to decay into light neutral gauge bosons Z Z , γ γ , Z Z and Z γ , however the branching ratio of these processes are minuscule due to the smallness of the gauge couplings g H and g X . As also shown in Fig. 17, precision measurements on the diphoton and Z γ productions at the HL-LHC can further constraint the parameter space of the model especially on the mixing angle θ 1 and charged Higgs mass for the light DM mass scenario.

Conclusions
Motivated by the new high-precision measurement of the SM W boson mass by the CDF II [1], we scrutinize the minimal G2HDM by computing the contributions from the inert Higgs doublet, hidden Higgs doublet, hidden neutral dark gauge bosons and extra heavy fermions in the model to the oblique parameters S, T , and U . We found out that only inert Higgs doublet can give a significant contribution while contributions from other new particles in the model are minuscule.
The favored regions on the model parameter space are obtained using the updated electroweak global fits of these oblique parameters [2,6] and the current constraints in the model including the theoretical constraints on the scalar potential, the Higgs data from the LHC, Z , dark photon physics constraints as well as relic density, direct detection and Higgs invisible decays constraints for the dark matter candidate W . We initiated two scenarios based on the dark matter mass, one is heavy DM mass scenario and another is light DM mass scenario. By comparing with the old global fit values for the oblique parameters from PDG, we found out that the CDF W boson mass measurement significantly impacts the mass splitting m = (m ± H − m D ) in the inert h-parity odd Higgs doublet H 2 . In particular, the result using the updated global fit values [6] favors a larger mass splitting (72 GeV m 100 GeV for the heavy DM mass scenario and 60 GeV m 100 GeV for the light DM mass scenario) as compared with one using the old global fit values (33 GeV m 78 GeV for the heavy DM mass scenario and 10 GeV m 75 GeV for the light DM mass scenario) from the PDG. We also found that to satisfy the DM relic density observed at the Planck collaboration, it is required the W coannihilation processes to be occurred for the heavy DM mass scenario and the nearly resonant annihilation process with the dark photon as a mediator for the light DM mass scenario. The former requires the differences (m H ± − m W )/m W and (m D − m W )/m W to be in ∼ O(10 −2 , 10 −1 ), while the latter requires the mediator mass to be around twice the DM mass. For the heavy DM mass scenario, the favored dark matter mass is around 1 TeV and can be probed by future dark matter direct detection experiments at DarkSide-20k and DARWIN as shown in Fig. 10. Whereas, for the light DM mass scenario, the favored DM mass locates at sub-GeV ∼ O(0.1−1) GeV and also can be probed by future dark matter direct detections at CDEX, NEWS-G, SuperCDMS, and HeRALD as shown in Fig. 13.
Furthermore we studied the diphoton and Z γ productions from the Higgs decays and the detectability of these processes at the LH-LHC. We showed that the contribution from the charged Higgs can be larger than the SM fermions (mainly top quark) for the h → Z γ process but not for the diphoton channel. In the low charged Higgs mass region, the contribution from charged Higgs can be significant and thus reduces the decay widths of h → γ γ and h → Z γ due to the destructive interference with the W boson contribution. We also showed that the deviations of the diphoton and Z γ productions from the SM predictions are highly depending on the mixing angle θ 1 between the SM Higgs boson and the scalar from the hidden doublet. A large portion of the favored region, especially the large mixing angle θ 1 and low charged Higgs mass regions, can be probed by more precise measurement of the h → γ γ signal strength and the detection of h → Z γ process at the HL-LHC.
Other rare decay modes h → γ γ , γ Z (one-loop) and h → γ γ , Z Z , and γ Z (tree) in G2HDM may be kinematical allowed and thus have important impacts at collider physics. For instance, the dark photon (γ ) and dark Z (Z ) may be long-lived, traversing some macroscopic distances for small enough g X and g H before they decay into SM light fermion pairs. We reserve the exploration of these modes for future studies. For a recent analysis of h → γ γ at the LHC, see [103].
Before closing, we reiterate that the results of the present work are obtained under the approximation of the new gauge couplings of SU (2) H ×U (1) X are much smaller than the SM gauge couplings, as suggested from the analysis of Z mass shift and dark matter direct detection in our previous works [54,55]. Consequently we can ignore the loop contributions from the extra gauge bosons (W ( p,m) , Z , γ ) to the oblique parameters since the effects would be down by the loop factor 1/16π 2 compared with the tree level mass mixing effects studied in Sect. 4.1. Furthermore, in analogous to the SM case, there should be 8 more vacuum polarization amplitudes i and i μν Z Z that are needed to be worried about. The formalism of the oblique parameters by Peskin and Takeuchi [8] would have to extend properly so as to take into account all new hidden particle contributions to all possible oblique parameters in G2HDM. Such a task is very interesting but beyond the scope of this present work. We hope to return to this challenge in the future. where We note that the arguments of the functions f (x) and g (x) in (76) and (77) are defined to be the inverse of those in [46]. All the charged fermions in G2HDM, including both the SM fermions f SM and the new heavy fermions f H contribute to F 1/2 i j . Thus where and with N c f being the color factor and Q f the electric charge of f in unit of e > 0; C f V j is the vector coupling of Z j with fermion f given by 6 where 6 It is well known that the axial vector couplings C f A j = (−C f L j + C f R j )/2 do not contribute to the h i → Z j γ fermion loop amplitudes.  Table 3 Coupling coefficients C f L j and C f R j for leptons  Table 5 Coupling coefficients C f V j for leptons  Table 2 for quarks and Table 3 for leptons. The vector couplings C f V j of quarks and leptons are listed in Tables 4 and 5 respectively.
There is only one charged Higgs H ± in G2HDM. Thus where Q H + = +1, and g h i H + H − and g Z j H + H − are the h i H + H − and Z j H + H − couplings in the G2HDM respectively. Explicitly they are h i → γ γ The decay rate for h i → γ γ is with similar definitions of F s i like F s i j before. The factor of 1/64 instead of 1/32 is due to Bose statistics in the diphoton final state. The previous formulas of F i j (W ± ), F i j ( f SM ), F i j ( f H ) and F i j (H ± ) can be easily translated into the diphoton case by taking m Z j → 0 limit or equivalently λ jl = 4 m 2 l /m 2 Z j → ∞ and noting that The results are Note that both the charged Higgs and new heavy fermion contributions were not handled correctly in h i → γ γ and h i → Z j γ in [46].