## Preliminaries

An additional U(1)_{D}gauge symmetry added to the SM is theoretically well-motivated and occurs in many top-down and bottom-up extensions of the SM. The U(1)

_{D}vector boson (the "dark photon" or the "dark-Z") is usually referred to as A′, Z′, γ

_{D}, or Z

_{D}in the literature and various possibilities exist to connect the additional U(1)

_{D}to the SM (see e.g. [1,2,3,4] for reviews). The U(1)

_{D}can couple to the SM sector via a small gauge kinetic mixing term [1/2] ϵF′

_{μν}B

^{μν}[5,6,7] between the dark photon and the hypercharge gauge boson. This renormalizable interaction can be generated at a high scale in a grand unified theory or in the context of string theory with a wide range of ϵ ∼ 10

^{−17}− 10

^{−2}. This term effectively gives SM matter a dark milli-charge and allows for dark photon decay to SM particles and possible experimental detection. To avoid the tight constraints on new long-range forces, a `dark Higgs' S with a non-zero vacuum expectation value can give a non-zero mass to the A′. An A′ with a sub-GeV mass can be probed at beam dumps and colliders, and with measurements of the muon anomalous magnetic moment, supernova cooling, and rare meson decays, see Fig. 1 and e.g. [4] for a recent review. A broken U(1)

_{D}can also lead to exotic Higgs decays, especially if there is mixing between the two Higgs sectors. In this context we refer to the corresponding vector field as Z

_{D}. The possibility of h → Z

_{D}Z

_{D}through Higgs-to-dark-Higgs mixing or h → Z Z

_{D}through Z-Z

_{D}mass mixing (which is also induced by the above-mentioned kinetic mixing) was discussed in [8] and [9,10], respectively, with both occurring, for example, in hidden valley models [11,12]. To examine the range of possible exotic Higgs phenomena due to a U(1)

_{D}sector we examine the model of [8], but with m

_{h}set to 125 GeV and allowing for the full range of dark Higgs and dark-Z masses relevant to exotic Higgs decay phenomenology.

^{1}This includes Higgs-to-dark-Higgs mixing and kinetic mixing between the B boson and the dark vector Z

_{D}, but no explicit mass mixing between the Z and Z

_{D}.

^{2}We will assume prompt Z

_{D}decays, which requires m

_{ZD}>~10 MeV given the current constraints shown in Fig. 1.

Figure 1: Constraints on ϵ, m

_{ZD}for pure kinetic mixing (no additional source of Z-Z_{D}mass mixing) for m_{ZD}∼ MeV-10 GeV. The black dashed line separates prompt (cτ < 1 μm) from non-prompt decays. The three blue lines are contours of Br(h → Z Z_{D}) of 10^{−4}, 10^{−5}, 10^{−6}respectively. Shaded regions are existing experimental constraints, see e.g. [4] for a recent review.## Model Details

The model is defined by a U(1)_{D} gauge sector and a SM singlet S that has unit charge under the U(1)_{D}. The kinetic terms of the hypercharge and U(1)_{D} gauge bosons (adopting mostly the notation of [9]) are

(1)

with ∧B_{μν}= ∂

_{μ}∧B

_{ν}− ∂

_{ν}∧B

_{μ}, ∧Z

_{Dμν}= ∂

_{μ}∧Z

_{Dν}− ∂

_{ν}∧Z

_{Dμ}, and cosθ

_{W}= g/√{g

^{2}+ g′

^{2}} is the usual Weinberg mixing angle. The hatted quantities are fields before diagonalizing the kinetic term. The Higgs potential is

(2)

The dark Higgs S acquires a vacuum expectation value and gives Z_{D}, which `eats' the pseudoscalar component of S, some mass m

_{ZD}. There are two connections between the dark and the SM sectors: the gauge kinetic mixing ϵ and the Higgs mixing ζ. The phenomenology depends on which one dominates. The SM gauge boson Z

_{μ}= − sinθ

_{W}B

_{μ}+ cosθ

_{W}W

^{3}

_{μ}is no longer a mass eigenstate and to leading order in ϵ the mass eigenstates with masses m

_{Z}, m

_{ZD}+

*O*(ϵ

^{2}) are

(3)

(Henceforth, we omit the tildes and will refer to the mass eigenstates unless otherwise noted.) Therefore, there are interaction terms of the form 2 ϵ_{Z}[(m

_{ZD}

^{2})/v] h Z

_{μ}Z

_{D}

^{μ}and ϵ

_{Z}

^{2}[(m

_{ZD}

^{4})/( m

_{Z}

^{2}v)] h Z

_{Dμ}Z

_{D}

^{μ}which lead to h → Z

_{D}Z and h → Z

_{D}Z

_{D}decays (though the latter is strongly suppressed), see Figure 3. If Z

_{D}is the lightest state in the dark sector it will decay to SM particles. For m

_{ZD}>~ GeV the Z

_{D}, branching ratios are easily computed to lowest order and without QCD corrections, and are shown in Figure 2 (a). For m

_{ZD}<~ GeV, non-perturbative QCD effects are important. They can be computed from the QCD contribution to the imaginary part of the electromagnetic two-point function, which in turn is determined from cross-section measurements of e

^{+}e

^{−}→ hadrons [14]. The resulting branching ratios are shown in Figure 2 (b). The most important qualitative difference to the scalar decays considered in the SM+S and 2HDM+S scenarios is that branching ratios are ordered by gauge coupling instead of Yukawa coupling, meaning decays to e

^{+}e

^{−}and μ

^{+}μ

^{−}remain large above the τ thresholds. Prompt Z

_{D}decay requires ϵ >~10

^{−5}− 10

^{−3}, as indicated in Figure 1, which summarizes the constraints on Z

_{D}kinetic mixing for our regime of interest.

Figure 2:

(a) Branching ratios for Z

The Higgs potential is minimized by vacuum expectation values of H_{D}decay, to lowest order and without QCD corrections, assuming decays to the dark sector are kinematically forbidden.*l*stands for either e or μ. Hadronization effects likely invalidate our simple calculation in the shaded region. (b) Branching ratios for Z_{D}decay for m_{ZD}<~3 GeV, including non-perturbative QCD effects.^{0}and S

(4)

The mass eigenstates(5)

have masses(6)

(Again we drop the tildes from now on and always refer to the mass eigenstates.) The effective Lagrangian contains terms of the form κh s s where κ = ζ(m_{h}

^{3}+ 2 m

_{h}m

_{s}

^{2})/(√{16 λ} (m

_{h}

^{2}− m

_{s}

^{2})), and 2 ϵ

_{h}[(m

_{ZD}

^{2})/w] h Z

_{D μ}Z

_{D}

^{μ}, which lead to exotic Higgs decays h → s s and h → Z

_{D}Z

_{D}, see Figure 3. The vertex h s Z

_{D}is present but is suppressed by both mixings.

Figure 3: The dominant exotic Higgs decays in the SM+V model. The h → Z Z

We can now discuss the relevant limits of this theory for exotic Higgs phenomenology:
_{D}matrix element is proportional to the gauge kinetic mixing ϵ, while h → Z_{D}Z_{D}and h → s s are controlled by the Higgs mixing parameter ζ. The vertex h s Z_{D}is present but suppressed by both mixings.- Gauge mixing dominates:For ϵ >> ζ the dominant exotic Higgs decay is h → Z Z
_{D}. To leading order in m_{ZD}^{2}/m_{Z}^{2}the partial width is

(7)

This agrees with the full analytical expression to ∼ 10% for m_{h}− m

_{Z}− m

_{ZD}> 1 GeV. Figure 1 shows contours of Br(h → Z Z

_{D}) = 10

^{−4}, 10

^{−5}, 10

^{−6}. The largest Br allowed by indirect electroweak precision constraints is ∼ 10

^{−4}. In this limit, the SM+V theory leads to ff + Z exotic Higgs signatures , and exotic branching ratio sensitivities better than ∼ 10

^{−3}are possible at LHC8. Bounds from LHC14 could become competitive with indirect limits. For very light Z

_{D}above the electron threshold this would also lead to lepton-jets + Z signatures, see here and here [15].

- Higgs mixing dominates:When ζ >> ϵ and Higgs mixing dominates then h → Z
_{D}Z_{D}, s s are both possible, depending on the spectrum of the dark sector. (We still assume that ϵ is large enough for Z_{D}to decay promptly.) The partial decay widths to leading order in ζ are

(8)

Different regions of of the (m_{ZD}, m

_{s}) mass plane are shown in Figure 4, along with the size of the Higgs mixing ζ ∼ 10

^{−3}− 10

^{−2}required for Br(h → Z

_{D}Z

_{D}, s s) = 10% and the relative rates of h → ss vs h → Z

_{D}Z

_{D}decays when both are allowed. In Region A (m

_{s}> m

_{h}/2, m

_{ZD}< m

_{h}/2) the only relevant exotic Higgs decay is h → Z

_{D}Z

_{D}. This allows for spectacular h → 2

*l*2

*l*′ decays (

*l*,

*l*′ = e or μ) with a reconstructed Z

_{D}resonance above the τ- or b-thresholds. Region B allows exotic Higgs decays both to Z

_{D}Z

_{D}and ss. The presence of two resonances below half the Higgs mass gives a rich exotic decay phenomenology. h → s s → 4 Z

_{D}occurs with roughly equal probability as h → Z

_{D}Z

_{D}and can result in spectacular final states with as many as 8 leptons. Note that, in this simplified model, there is no corresponding Z

_{D}→ ss decay in the lower right hand corner of that mass plane. However, a (pseudo)scalar pair could be produced from dark vector decay in e.g. a 2HDM+V framework, resulting in final states with as many as 8 b-quarks. Already with current data, limits of Br(h→ Z

_{D}Z

_{D}) <~10

^{−4}can be achieved. Each of the above cases may, for suitable masses, also lead to interesting `lepton jet' signatures, see here.

- Intermediate Regime:Here the decays induced by kinetic and Higgs mixing are comparable. For example, Figure 1 shows that ϵ ∼ 10
^{−2}is not excluded for some values of m_{ZD}, allowing Br(h→ Z Z_{D}) ∼ 10^{−4}. The branching ratios for h→ Z_{D}Z_{D}, s s will be similar if ζ ∼ 10^{−4}.

Figure 4:

**Left:**Regions of the (m_{ZD}, m_{s}) mass plane in the SM+V model with different exotic Higgs decays for ζ >> ϵ (i.e. when the mixing between the Higgs and dark-Higgs dominates over the kinetic mixing). The black contours are the values of ζ×10^{3}required for Br(h → Z_{D}Z_{D}, s s) = 10%. Region A is the case examined by [8] (the dotted red line indicates m_{h}= m_{s}). Region C has no exotic Higgs decays. Region D reproduces the SM+S model. Region B has both h → s s and h → Z_{D}Z_{D}decays, with the h → ss fraction of exotic decays shown on the**right**. In the upper left shaded region, s → Z_{D}Z_{D}is the dominant decay mode of the dark scalar. This allows the Higgs to decay to up to 8 SM fermions.## Summary

In summary, the SM+V setup allows for many different kinds of exotic Higgs decays, including h → Z Z_{D}, h→ Z

_{D}Z

_{D}, and h→ ss, with Z

_{D}→ f―f , and s → f―f or s → Z

_{D}Z

_{D}→ (f―f)(f―f). This leads to final states of Z + (f―f), (f―f)(f―f), and ((f―f)(f―f))((f―f)(f―f)), where parentheses around a set of particles denotes a resonance (all final-state particles combined will form the Higgs resonance). Since the Z

_{D}(although not the s) couples to the fermions' gauge charges, final states with several light leptons have sizable branching fractions over the entire kinematically permitted mass range. Certain spectra can produce interesting lepton-jet signatures.

## MadGraph Model

## MadGraph Model

We have constructed a validated MadGraph 5 model for the SM + vector + dark higgs scenario outlined here. This can be used for signal generation in experimental and phenomenological studies.

Link to MadGraph5 model for SM + S + V

## Branching Ratio Tables

In "Illuminating Dark Photons with High-Energy Colliders" (Curtin, Essig, Gori, Shelton, arXiv:1412.0018) we compute the dark photon total width and branching ratio to leptons Br(ZD->ll) including QCD corrections, using PDG experimental data for mZD < 12 GeV, and by repurposing the calculations of hep-ph/0005139 for mZD > 12 GeV. We also compute the exotic higgs decay branching ratios Br(h->ZdZ*->4l) and Br(h->ZdZd->4l). All of these computed widths and branching ratios can be downloaded as a text table here. (See 1412.0018 for more details.)## References

[1]P. Langacker,

File translated from T

*The Physics of Heavy Z′ Gauge Bosons*,*Rev.Mod.Phys.*81 (2009) 1199-1228, [arXiv:0801.1345]. [2]J. Jaeckel and A. Ringwald,*The Low-Energy Frontier of Particle Physics*,*Ann.Rev.Nucl.Part.Sci.*60 (2010) 405-437, [`arXiv:1002.0329`]. [3]J. Hewett, H. Weerts, R. Brock, J. Butler, B. Casey,*et. al.*,*Fundamental Physics at the Intensity Frontier*, [`arXiv:1205.2671`]. [4]R. Essig, J. A. Jaros, W. Wester,*et. al.*,*Dark Sectors and New, Light, Weakly-Coupled Particles*, [`arXiv:1311.0029`]. [5]B. Holdom,*Two U(1)'s and Epsilon Charge Shifts*,*Phys.Lett.*B166 (1986) 196. [6]P. Galison and A. Manohar,*TWO Z's OR NOT TWO Z's?*,*Phys.Lett.*B136 (1984) 279. [7]K. R. Dienes, C. F. Kolda, and J. March-Russell,*Kinetic mixing and the supersymmetric gauge hierarchy*,*Nucl.Phys.*B492 (1997) 104-118, [`hep-ph/9610479`]. [8]S. Gopalakrishna, S. Jung, and J. D. Wells,*Higgs boson decays to four fermions through an abelian hidden sector*,*Phys.Rev.*D78 (2008) 055002, [`arXiv:0801.3456`]. [9]H. Davoudiasl, H.-S. Lee, and W. J. Marciano,*'Dark' Z implications for Parity Violation, Rare Meson Decays, and Higgs Physics*,*Phys.Rev.*D85 (2012) 115019, [`arXiv:1203.2947`]. [10]H. Davoudiasl, H.-S. Lee, I. Lewis, and W. J. Marciano,*Higgs Decays as a Window into the Dark Sector*, [`arXiv:1304.4935`]. [11]M. J. Strassler,*Why Unparticle Models with Mass Gaps are Examples of Hidden Valleys*, [`arXiv:0801.0629`]. [12]M. J. Strassler and K. M. Zurek,*Echoes of a hidden valley at hadron colliders*,*Phys.Lett.*B651 (2007) 374-379, [`hep-ph/0604261`]. [13]C.-F. Chang, E. Ma, and T.-C. Yuan,*Multilepton Higgs Decays through the Dark Portal*, [`arXiv:1308.6071`]. [14]Particle Data Group, J. Beringer*et. al.*,*Review of Particle Physics (RPP)*,*Phys.Rev.*D86 (2012) 010001. [15]N. Arkani-Hamed and N. Weiner,*LHC Signals for a SuperUnified Theory of Dark Matter*,*JHEP*0812 (2008) 104, [`arXiv:0810.0714`].### Footnotes:

^{1}Ref. [13] appeared while this work was being completed, performing a similar analysis with a different focus on constraining the couplings of the extended Higgs potential for relatively low m_{ZD}< 5 GeV.^{2}The constraints shown in Fig. 1 are altered in the presence of such pure mass mixing, which requires additional Higgs doublets that also carry dark charge. The resulting Z_{D}→ SM decays would be more Z-like and lead to additional constraints from rare meson decays as well as new parity-violating interactions [9]. However, we stress that the exotic Higgs phenomenology would not be qualitatively different.File translated from T

_{E}X by T_{T}H, version 4.03. On 17 Dec 2013, 16:29.