#### Abstract

New Abelian vector bosons can kinetically mix with the hypercharge gauge boson of the Standard Model. This letter computes the model-independent limits on vector bosons with masses from 1 GeV to 1 TeV. The limits arise from the numerous experiments that have been performed in this energy range and bound the kinetic mixing by for most of the mass range studied, regardless of any additional interactions that the new vector boson may have.

#### 1. Introduction

The Standard Model (SM) successfully describes all known interactions of SM fermions and gauge bosons; however, there are several phenomena that motivate physics beyond the SM. Chief among these open questions is the identity of dark matter and its interactions with the SM. Recent anomalies in cosmic ray and direct detection experiments have motivated the exploration of new gauge interactions in a putative dark sector [1, 2]. New Abelian vector bosons provide one of the most robust portals for dark matter—SM interactions. The new vector boson can interact with the SM, even if no SM fermions are directly charged under the additional gauge symmetry. This interaction occurs via mixed kinetic terms between the SM's hypercharge field strength and the new Abelian field strength [3].

The Lagrangian for a kinetically mixed theory is where is the field strength for the new vector boson, is the field strength for the SM hypercharge, and encapsulates the interactions of the with fields in the dark sector. The mixed kinetic term is interesting for several reasons. First, it is a dimension 4 operator, meaning that it can be generated at high energies without decoupling. Second, this coupling allows communication with a secluded sector that otherwise has no interactions with SM fields. These considerations have motivated a dedicated program to search for kinetically mixed vector bosons [4–7].

Most of the current program for discovering a kinetically mixed vector boson involves producing the state and searching for its subsequent decays. This method is promising, but has the drawback that it assumes that the searches can recognize the decay products of the . If the dark sector has states lighter than the , then the will preferentially decay to the dark sector over SM states because the kinetic mixing parameter almost always satisfies . Searching for the by looking for the dark sector final states requires a wide-ranging search program because the dark sector may decay back to the SM in a variety of different ways, for example, lepton jets [8–12]. Model-independent searches are possible using completely inclusive searches, for example, , but these are challenging and few of these searches have actually been performed.

At low masses, the best model-independent bounds arise from the measurements of the electron and muon [13]; however, the power of begins to weaken for . At masses far above collider energies, the can be integrated out and its effects can be encapsulated in higher-dimension operators, most importantly and [14–16]. colliders have probed up to GeV and therefore, the effects of the cannot be parameterized as local higher-dimension operators for masses less than this energy scale.

This letter computes the model-independent constraints on the kinetic mixing parameter, , for masses between 1 GeV and 1 TeV by looking for the effects of virtual s on precision SM observables. This approach has the benefit of not requiring any knowledge of the decay modes of the and sets an upper limit on regardless of the behavior of the decay modes to the dark sector.

#### 2. Kinetic Mixing

Kinetic mixing changes the mass eigenstates and interactions of the vector bosons. What follows is a brief synopsis of the results in [17], see also [18–21]. After diagonalizing the kinetic terms and going to the mass eigenstate basis, the SM neutral current interactions are modified. Absorbing the gauge coupling constants into the definition of the currents, the neutral current interactions are where the notation for the gauge and mass eigenstates is and the currents are with the diagonalization matrix stand for cosine, sine, and tangent, respectively, and and are the cosine and sine of the weak mixing angle. The photon's interactions are unaltered due to its residual gauge invariance. The angle is defined as with . After changing to the mass eigenstate basis, the physical mass of the , , is and the physical mass of the new vector boson is

These corrections to the SM neutral currents and to the mass of the place model-independent bounds on . The next section describes the SM measurements that are sensitive to these modified neutral current interactions.

#### 3. Precision SM Measurements

Virtual exchange modifies measured observables such as Bhabha scattering, forward-backward asymmetry measurements, , and the total hadronic cross sections. The mass of the is the most powerful single measurement but the constraint is augmented by other measurements at and above the pole. Additionally, if the has a sizeable branching ratio back to the Standard Model, resonant production of the bounds the parameter space at specific energies.

The strongest constraint on comes from the shift of the mass [22]. Notice that in (2.6) changes sign as goes through , meaning that the corrections to the mass vanish at this point. Defining , the correction to the mass is given by so that as , there is no bound on resulting from the measurement of . There is a reduction in the limits on for where other measurements must take over for the mass measurement.

colliders measure the SM neutral current interactions and when , the couples dominantly to the electromagnetic current, which causes a kink to appear in the running of the fine structure constant. Differential Bhabha scattering measures and there is a wealth of data from experiments such as OPAL [23], DELPHI [24], SLD [25], TASSO [26], CELLO [27], and TRISTAN [28]. As a result, the new vector boson changes the predictions for differential Bhabha scattering. All the experiments above have a large forward bin of , so only a small range of is probed at each experiment. The forward bin normalizes the luminosity and cannot be used as a constraint, thereby limiting the power of these measurements. Differential Bhabha scattering for is not useful but provides additional constraints at and above the pole, where corrections to are less powerful.

In addition to differential Bhabha scattering, the forward-backward asymmetries for the bottom, charm, muon, and tau are measured at , effectively fixing and [29]. The modification to the SM neutral currents alters and and leads to a conflict with other SM predictions, most notably and , that is, .

Resonant and on-shell production of the can be relevant even if there is a small width directly back to the SM. The has a decay width into the SM and dark sectors given by The width of the into the dark sector is unknown; however, given that bounds from the mass set and is bounded by , there can be a detectable width for the back into the SM. As a way to parameterize these effects, two different dark sector widths are used in setting limits On-shell production of the is calculated using MadGraph 4.4.32 [30]. Only the interference between the and the SM is explicitly computed by zeroing out the and squared matrix elements. This results in a deviation from the SM that scales as and the calculations can be compared with measurements using the methods described in the next section.

The total hadronic cross sections, , are measured at LEP2 with GeV [23, 24] as well as at many other experiments with 22 GeV ≤ GeV [31]. These measurements provide additional bounds because the results from differential Bhabha scattering are not reported at every energy. While the error bars are large compared to the differential Bhabha scattering, resonant production enhances sensitivity if . Radiative return processes involving the could in principle constrain the theory for away from ; however, these never provide competitive measurements.

The can have exotic decays into the hidden sector, and assuming that there are no mass thresholds in the hidden sector between and , then The line shape measurement constrains in a model-independent manner giving a bound on [22]. The bound on the width of the is shown in Figure 1.

In addition to precision measurements, direct searches at the Tevatron can produce on-shell s. This letter finds that the Tevatron's sensitivity is just beneath precision electroweak results even assuming [17, 32].

#### 4. Results

The regions in the parameter space consistent with precision SM measurements are found by performing a global fit to the SM parameters. This letter uses a test that is a function of , and the Standard Model parameters, , and . is fit by decay and does not vary in practice. Most of the data is statistics dominated so that in the case where the signal predicts an excess of events, the simplifies into Due to the high statistics, these Poisson summations can be approximated as gaussian integrals.

The advantage of using the method is that it is not diluted by superfluous measurements that have no * a priori* possibility of constraining a theory with a given , only experiments that have a significant impact.

While superfluous measurements are ignored, measurements that only slightly affect the data can have a significant influence if there are enough of them. This can be illustrated by considering experiments which all give the same result, only slightly different from the SM. A typical analysis will never exclude the SM because the /d.o.f. is small. The method will eventually exclude the SM because will be small for large .

For GeV, the effects of the new vector boson can be encapsulated in terms of local operators and coincide with the precision electroweak analyses, for example, the , parameters [14–16] or more recently [33–35]. For , the bounds are close to those from [17–21], which only use the constraint from .

Figure 2 shows the 95% confidence level (CL) excluded regions in the plane obtained in this study. Wide s are best constrained by the mass of the for most of the parameter space. The exception occurs near the mass. The forward-backward asymmetries, , and augment the limits when the corrections to the mass vanish and also for GeV where LEP2 forward-backward measurements are more constraining than the mass. Limits on narrow s are enhanced for for the numerous experiments. The forward-backward asymmetries, hadronic cross section, and differential Bhabha scattering measurements provide the additional constraints. The peaks appearing in the exclusion region can be traced to experimental energies at which various experiments were conducted. The constraint on near the is illustrated in the inset of Figure 2. For comparison, the bounds from , and model-dependent BaBar searches from [4, 36] are shown.

The main improvement of this work over previous papers is the use of lower energy experiments SLD, TASSO, CELLO, and TRISTAN to measure the running of and how it changes due to the addition of a new boson. These low energy experiments augment the sensitivity to a new boson, although the single strongest constraint is still the change in mass of the boson. These low energy experiments are sensitive to channel production of the and hence the width plays an important role. The width of the depends on interactions with a hidden sector and cannot be a determined model, independently. Two different scenarios of the width of the are considered. If the is narrow, then there are many new small dips which are excluded (see Figure 2). If the has a larger width, then the effects are washed out and the results are in agreement with [21].

The model-independent limits on kinetic mixing were computed in this letter and found to be for most of the mass range studied, 1 GeV < GeV. The possible use of radiated return to place tighter constraints on was investigated at both LEP1 and LEP2 energies; however, this channel did not help place tighter bounds on kinetic mixing. Even with the constraints found in this letter, there still is a vast parameter space available for a kinetically mixed vector boson to mediate interactions between a dark sector and the SM. The current program of searching for model-dependent decay modes at low energy experiments will augment these model-independent limits for GeV. For higher energies, only the LHC will provide additional information for 200 GeV 3 TeV [17]. The relatively weak limits for GeV motivate new high intensity experiments to potentially discover new interactions of this form.

#### Acknowledgments

The authors would like to thank D. E. Kaplan, M. Lisanti, and M. Peskin for numerous useful conversations. A. Hook, E. Izaguirre, and J. G. Wacker are supported by the US DOE under contract number DE-AC02-76SF00515. A. Hook, E. Izaguirre, and J. G. Wacker receive partial support from the Stanford Institute for Theoretical Physics. J. G. Wacker is partially supported by the US DOE's Outstanding Junior Investigator Award.