Precision studies for string derived Z (cid:48) dynamics at the LHC

We consider Z (cid:48) s in heterotic string derived models and study Z (cid:48) resonant production at the TeV scale at the Large Hadron Collider (LHC). We use various kinematic diﬀer-ential distributions for the Drell-Yan process at NNLO in QCD to explore the parameter space of such models and investigate Z (cid:48) couplings. In particular, we study the impact of Z - Z (cid:48) kinetic-mixing interactions on forward-backward asymmetry ( A FB ) and other distributions at the LHC.


Introduction
With the observation of a scalar particle compatible with the Standard Model electroweak Higgs doublet, the last piece in the Standard Model has been confirmed experimentally, and in the coming years its basic properties will be probed and elucidated.Furthermore, the discovery of the Higgs particle at 125 GeV suggests that the electroweak symmetry breaking mechanism is perturbative, rather than non-perturbative.This supports the old expectation that the Standard Model may provide viable parameterisation of particle physics data up to the Planck scale, where gravitational effects break the effective field theory description.Alas, the question remains how the Higgs mass and the ensuing electroweak symmetry breaking scales are protected from radiative corrections from the higher energy scale.
While not yet seen at low energy experiments, supersymmetry remains one of the most appealing extensions of the Standard Model to address this puzzle.Furthermore, many supersymmetric models in fact predict that the Higgs particle should exist below 135 GeV.The observed Higgs particle is therefore compatible with supersymmetry.However, while supersymmetry protects the Higgs mass at the multi-loop level in the perturbative expansion, it introduces a problem at tree level, which is known as the µ-problem.Supersymmetry mandates the existence of a pair of Higgs doublets, which can receive a bilinear mass term in the superpotential.In the Minimal Supersymmetric Standard Model, the µ-parameter is therefore set by hand to be of the order of the electroweak scale, but nothing prevents it from being of the order of the Planck scale.Furthermore, absence of global symmetries in quantum gravity and the experience from quasirealistic string derived models, do in fact suggest that generically the µ-parameter will be of that order.A possible solution to the problem may be obtained if the Higgs states are chiral under an additional U (1) symmetry, which remains unbroken down to TeV scale.In this case, the bi-linear mass term for the supersymmetric Higgs pair can only be generated by the vacuum Expectation Value (VEV) that breaks the extra U (1) symmetry.
String theory provides the most developed framework to study the incorporation of the Standard Model in a perturbatively consistent theory of quantum gravity.Indeed, over the past few decades detailed phenomenological models have been built that produce the spectrum of the Minimal Supersymmetric Standard Model (MSSM) in the effective low energy field theory limit below the Planck scale (for review and references see e.g.[1]).The heterotic-string models constructed in the free fermionic formulation [2][3][4][5], which are Z 2 ×Z 2 orbifolds of six dimensional toroidal spaces [6], are among the most realistic models constructed to date, and give rise to an abundance of three generation models with different unbroken SO (10) subgroups below the string scale.However, the construction of string models with additional U (1) gauge symmetries, that may stay unbroken down to low scales, has proven to be a difficult task [7,8].The reason is that the extra U (1) gauge symmetries that are often discussed in the context of Grand Unified Theory (GUT) extensions of the Standard Model are anomalous in explicit string models and therefore cannot remain unbroken down to low scales.In fact to date, there exist a single string derived model in which the extra U (1) is anomaly free and can remain unbroken down to the TeV scale [9].As a bonus the string derived model predicts the existence of electroweak doublet pairs at the Z breaking scale, that are chiral under the extra U (1) symmetry.The bilinear mass term for these electroweak Higgs doublets can therefore be generated only by the VEV that breaks the extra U (1) symmetry.
In this paper, we therefore study some of the properties of the extra U (1) symmetry predicted by the string model and explore the parameter space of the Z by using various kinematic distributions of resonant lepton-pair production in the Drell-Yan (DY) process at the LHC.Other related literature can be found in refs.[10][11][12].In addition, we study the effect of gauge kinetic-mixing interactions χF µν F µν in the lagrangian [13].The impact of such a renormalizable operator is interesting in this context because it may in general be produced (at both field theory or string theory level) by physics at high scales above the Z breaking scale, with no suppression by the large-mass scale [14,15] and can impact the couplings of Z s at the TeV scale [16].The interplay between the gauge kinetic-mixing parameter and coupling of the Z , as well as other Z properties, can be explored by using forward-backward asymmetry (A F B ) distributions in DY, where this observable can be advantageous in extra-resonance searches and as a model discrimination tool [17][18][19][20][21][22][23][24][25][26][27].
This paper is organized as follows.In Sec. 2 we describe the structure of the string derived Z model.In Sec. 3 we discuss the phenomenological aspects of the model and analyze kinetic mixing.In Sec. 4 we illustrate and discuss our results.Finally, in Sec. 5 we present our conclusions.

The String Model
We elaborate in this section on the structure of the string derived Z model [9].The model was constructed in the free fermionic formulation [28][29][30].Only the important points relevant for the properties of the extra Z gauge boson at low scales are discussed here and further details can be found in the original literature.Particularly relevant for the low energy phenomenology is the precise combination of worldsheet U (1) currents predicted by the string model.In this respect we note that while similar aspects can be discussed in field theory models, in the string model the U (1) combination is forced due to the available scalar multiplets to break the GUT symmetry, whereas in field theory models the scalar multiplets are a matter of choice.
In the free fermionic formulation all the worldsheet degrees of freedom required to cancel the conformal anomaly are represented in terms of free fermions propogating on the string worldsheet.In the standard notation the 64 worldsheet fermions in the lightcone gauge are denoted as: Left-Movers: where the {y, ω|ȳ, ω} 1,••• ,6 correspond to the six dimensions of the internally compactified manifold; ψ1,••• ,5 generate the SO(10) GUT symmetry; φ1,••• ,8 generate the hidden sector gauge group; and η1,2,3 generate three U (1) gauge symmetries.Models in the free fermionic formulation are defined in terms of boundary condition basis vectors, which specify the transformation properties of the worldsheet fermions around the noncontractible loops of the vacuum to vacuum amplitude, and the Generalised GSO projection coefficients of the one loop partition function [28][29][30].The free fermion models are Z 2 × Z 2 orbifold of six dimensional toroidal manifolds with discrete Wilson lines [6].
The Standard Model particle charges and data motivates the embedding of the Standard Model states in representations of Grand Unified Theories, like SO( 10) and E 6 .The rank of these groups exceeds that of the Standard Model, and therefore they predict the existence of additional gauge symmetries, beyond the Standard Model.Compactifications of the heteroticstring to four dimensions do in fact give rise to E 6 GUT like models [31], which gives rise to string inspired Z models.The additional U (1) symmetries in these string inspired models have an E 6 embedding and produced numerous papers since the mid-eighties (for reviews and references see e.g.[10,19,[32][33][34]).The construction of string derived models that allow for an unbroken extra U (1) symmetry to remain unbroken down to low scales proves, however, to be very difficult.The symmetry breaking pattern E 6 → SO(10)×U (1) A in the string derived models entails that U (1) A is anomalous and cannot be part of an unbroken U (1) Z at low energy scales [35].String derived models with anomaly free U (1) Z / ∈ E 6 were analysed in [7,[36][37][38], but agreement with the measured values of sin 2 (θ) W (M Z ) and α s (M Z ) favours Z models with E 6 embedding [8].We remark that the anomaly free U (1) combination of U (1) B−L and U (1) T 3 R ∈ SO (10) may in principle stay unbroken down to low energy scales [39].However, adequate suppression of the left-handed neutrino masses is facilitated if this U (1) symmetry is broken at a high energy scale [40].Constructing string models that enable an extra U (1) ∈ E 6 symmetry to stay unbroken down to low energy scales necessitates the construction of string models in which U (1) A is anomaly free.One method of achieving this outcome is to enhance U (1) A to a non-Abelian gauge symmetry á la ref. [41].An alternative is the string derived model of ref. [9] that uses the Spinor-Vector Duality (SVD) that was observed in Z 2 × Z 2 orbifolds [42][43][44].The SVD is under the exchange of the total number of (16 + 16) representations of SO (10) with the total number of 10 representations, and is readily understood if we consider the enhancement of SO(10) × U (1) A to E 6 .The chiral and anti-chiral multiplets of E 6 decompose under SO(10) × U (1) as 27 = 16 + 10 + 1 and 27 = 16 + 10 + 1.In this case the # 1 of (16 + 16) and # 2 of 10 multiplets are equal.The E 6 symmetry point in the moduli space corresponds to a self-dual point under the exchange of the total number of SO (10) spinorial plus antispinorial, with the total number of vectorial, multiplets.Breaking the E 6 gauge symmetry to SO(10) × U (1) A results in the projection of some of the spinorial and vectorial multiplets, which results in U (1) A being anomalous.However, there may exist vacua with equal numbers of (16 + 16) spinorial, and 10 vectorial, multiplets, and traceless U (1) A , without enhancement of the SO(10) × U (1) A symmetry to E 6 [9].In such models the chiral spectrum still forms complete E 6 representations, but the gauge symmetry is not enhanced to E 6 .In this cases U (1) A is anomaly free and may remain unbroken down to low scales.
A fishing algorithm to extract models with specified physical properties was developed by using the free fermionic model building rules [5,[45][46][47][48][49][50][51].In ref. [9], using the free fermionic fishing algorithm, such a spinor-vector self dual model was obtained with subsequent breaking at the string scale of the SO(10) symmetry to the SO(6) × SO(4) subgroup, and preserves the spinor-vector self-duality.This model is a string derived model in which an extra U (1) with E 6 embedding may remain unbroken down to low scales.The full massless spectrum of the string derived model is given in ref. [9].The observable and hidden gauge groups at the string scale are produced by untwisted sector states and are given by: observable : The massless string spectrum contains the fields required to break the GUT symmetry to the Standard Model.The string model contains two anomalous U (1) symmetries TrU (1) 1 = 36 and TrU (1 The E 6 combination, given by, is anomaly free and can be part of an unbroken U (1) symmetry below the string scale.The SO(6) × SO(4) observable gauge symmetry is broken by the VEVs of the heavy Higgs fields H and H .The charges of these fields under the Standard Model gauge group factors is given by: The VEVs along the N and N directions preserves N = 1 spacetime supersymmetry along a flat direction and leave the U (1) combination unbroken below the string scale.This U (1) combination is anomaly free provided that U (1) ζ is anomaly free, as is the case in the string derived model [9].We emphasise that this symmetry breaking pattern, and the U (1) Z combination in eq. ( 3) is enforced in the string derived model due the available scalar states in the string spectrum to break the non-Abelian SO(6) × SO(4) gauge symmetry, and due to a doublet-triplet missing partner mechanism [52] that gives heavy mass to coloured scalar states that arise in the untwisted sector of the string model [9].Thus, the combination given in eq. ( 3) is the extra U (1) combination that can arise in the string derived model with E 6 embedding of the charges.Contrary to the situation in string inspired model, where any combination of U (1) C − U (1) L and U (1) ζ is possible, the U (1) combination U (1) Z is the uniquely predicted combination in the string derived model.
Anomaly cancellation of the U (1) Z gauge symmetry down to low scales, requires the vectorlike leptons {H i vl , Hi vl }, and quarks {D i , D i }, that are obtained from the vectorial 10 multiplets of SO (10), and the SO(10) singlets S i in the 27 representation of E 6 .The supermultiplet1 states below the SO(6) × SO(4) breaking scale are displayed schematically in Table 1.The three right-handed neutrino N i L states obtain heavy mass at the SU (2) R breaking scale, which produces the seesaw mechanism [53].The spectrum below the SU (2) R breaking scale is assumed to be supersymmetric.We include in the spectrum an additional pair of vector-like electroweak Higgs doublets, that facilitate gauge coupling unification at the GUT scale.This is justified in the string inspired models due to the string doublet-triplet splitting mechanism [54,55].We note from table 1 that this extra Higgs pair is not chiral with respect to U (1) Z , contrary to the three pairs, {H i vl , Hi vl }, that are chiral with respect to U (1) Z .Therefore, mass terms for these three pairs can only be generated by the breaking of U (1) Z , whereas the bilinear mass term for the extra vector-like Higgs pair can, in principle, be generated at a high scale.The states φ and φ are exotic Wilsonian states [9,56].Additionally, the existence of light states ζ i , that are neutral under the SU (3) C × SU (2) L × U (1) Y × U (1) Z low scale gauge group, is allowed.The U (1) Z gauge symmetry can be broken at low scale by the VEV of the SO(10) singlet fields S i and/or φ 1,2 .
Z quantum numbers, with i = 1, 2, 3 for the three light generations.The charges are displayed in the normalisation used in free fermionic heterotic-string models.

Phenomenological aspects of the model
In this section we illustrate details of the phenomenological analysis for the model described in the previous section whose charge assignment is given in Table 1.The analysis is based on previous work of two of the authors, published in refs.[38,[57][58][59] which is used as a reference for the notation.Some definitions are also imported for consistency.

Neutral gauge bosons sector
We start by analyzing the neutral gauge bosons sector.The covariant derivative relative to the low-scale gauge group Z is defined as where g 2 is the SU (2) coupling and W a µ and τ a are the SU (2) gauge fields and generators respectively; g Y is the U (1) Y coupling, while A Y µ is the gauge field and Y is the hypercharge.g Z is the U (1) Z coupling and B µ and z are the gauge field and charge of the U (1) Z respectively.g s is the strong coupling constant, and A a and T a are the SU (3) C gauge fields and generators respectively.
The lagrangian of the Higgs sector is given where the scalar components of the supermultiplets which give the two Higgs doublets are represented by the H 1 and H 2 fields, that have hypercharge Y H1 = −1, Y H2 = 1 and U (1) Z charges z H1 and z H2 .The singlets are the S i fields, and as in [57] we consider only one S field for simplicity.The hypercharge of S is Y S = 0 and its charge under U (1) Z is z s .The Higgs fields are parametrized as and the vacuum expectation values (VEVs) are defined as Expanding the Higgs lagrangian and collecting quadratic terms, we obtain the mass matrix in the (W 3 µ , A Y µ , B µ ) basis as where for convenience we have defined the following quantities The zero eigenvalue of the mass matrix corresponds to the photon.The other two eigenvalues correspond to the square mass of the Z and Z physical gauge bosons where From Eq. 10 we see that the presence of the Z induces corrections on the mass of the SM Z boson.In Fig. 1 we show the impact of these corrections on the SM Z boson mass within the current uncertainties reported in the PDG [60] as a function of g Z for different values of M Z .We note that for g Z 0.4 deviations are well within the current uncertainties.In the limit where v s is very large, the Z mass decouples and one obtains simplified expressions for the gauge boson masses, From the mass matrix diagonalization we obtain the normalized orthogonal eigenvector matrix which is used to rotate the gauge field basis to the physical field basis.The exact form of the elements of this matrix is given in Appendix A. When v s is large, the rotation matrix can be simplified as 91.17 where θ W is the Weinberg angle, sin θ W = g 2 /g and cos θ W = g Y /g as in the Standard Model, and the small mixing between the Z and Z is parametrized by the angle δ, defined by To simplify the notation, we define a perturbative parameter ε = sin δ.The rotation matrix in Eq. 12 to first order in ε can now be written as According to the parameter space explored in this work, ε is typically 10 −4 ε 10 −3 .After fields are rotated to the physical basis, the covariant derivative for the neutral current (NC) sector of Eq. 4 in terms of the mass eigenstate fields is given by The neutral current contributions in the Lagrangian, that are of the form ψγ µ c Z ψZ µ and ψγ µ c Z ψZ µ , can be rewritten using the relations where the charges τ 3 , Q, z correspond to the fermion and ε is the (small) mixing parameter.Using this notation, the interactions between fermions and electroweak neutral gauge bosons are expressed in terms of the vector and axial vector couplings, i.e., ψγ µ g 4 g Z V − g Z A γ 5 ψZ µ and ψγ µ g 4 g Z V − g Z A γ 5 ψZ µ .These couplings are explicitly given as where τ 3 L , z L/R refer to left/right-handed fermions.

Kinetic mixing
In this section we generalize the results obtained in the previous section to the case of kineticmixing interactions.The impact of kinetic mixing is assessed in Sec. 4, where applications to hadron collider phenomenology are discussed.
Since the low-energy gauge group includes two abelian groups, the Lagrangian for the kinetic terms can include a contribution which directly couples the A Y and B fields and does not break either U (1) Y or U (1) Z gauge invariance.This term corresponds to gauge kinetic mixing and is of the form where A Y µν and B µν are the field-strength tensors for the U (1) Y and U (1) Z symmetry groups respectively, and χ is a parameter such that |χ| < 1.This renormalizable operator plays an interesting role in the context of our model because it can be produced at scales much higher than the Z breaking scale and can affect the couplings of Z s at the TeV scale.It is therefore interesting to study the interplay between the gauge kinetic-mixing parameter and coupling of the Z .The expression for the electroweak NC covariant derivative is generalized as where Q T = (Y /2, z/2) is a vector containing the charges for the two U (1) groups, G is the mixing matrix which is defined as and contains the gauge couplings, and A = (A Y µ , B µ ) T contains the U (1) gauge fields.The off-diagonal elements g AB and g BA contain the kinetic mixing between the U (1) fields.
It is convenient to perform the rotation A → A = RA, G → G = GR −1 where R is an orthogonal 2 × 2 matrix, and choose g BA = 0 in order to eliminate the z s g BA A Y µ term in D N C µ .This way, the Higgs singlet does not affect the SU (2) × U (1) Y sector in the mass matrix.In the new basis, the mixing matrix reads where The electroweak NC covariant derivative becomes where A Y µ and B µ are the rotated U (1) gauge fields.Therefore, the inclusion of gauge kinetic mixing can be studied by using a single parameter g KM .The mass matrix in the (W where x z and N z have been replaced by their primed expressions given by The calculation for the gauge boson masses and the rotation matrix to the physical gauge boson basis proceeds identically to the previous section, with the replacements x z → x z and N z → N z .With these modifications, the general form of the fermionic couplings to the Z and Z in presence of kinetic mixing is where the perturbative parameter from the previous section is modified as The vector and axial vector couplings of the gauge bosons to the fermions are now

Hadron Collider Phenomenology Applications
In this section we explore the parameter space of the model and perform a detailed analysis for proton-proton collisions at the LHC using Drell-Yan (DY) kinematic distributions calculated at the next-to-next-to-leading order (NNLO) in the strong coupling constant α s of Quantum Chromodynamics (QCD).The contribution of electroweak corrections is not considered here and their impact will be analyzed in a forthcoming work.We explore the sensitivity of the forwardbackward asymmetry (A F B ) distribution to the Z and its parameters.In particular, we use A F B distributions to study the interplay between the gauge kinetic-mixing parameter and coupling of the Z .It has been pointed out in several works (e.g., see [17][18][19][20][21][22][23][24][25][26][27] and references therein) that A F B distributions are very sensitive to SM deviations in the electroweak sector, and in particular to the presence of extra gauge vector bosons.Moreover, the ATLAS, LHCb, and CMS collaborations at the LHC have performed high-precision measurements of A F B distributions at 7, 8, and 13 TeV collision energy respectively [61][62][63].The CMS collaboration [63] has recently set lower limits on the mass of additional gauge bosons from sequential SM extensions at around 4.4 TeV.Our theory predictions for the DY cross section are complemented by the calculation of uncertainties induced by parton distribution functions (PDFs) in the proton, as well as scale uncertainties for some of the distributions.PDFs represent one of the major sources of uncertainty in cross section calculations and complicate model validation and discrimination.The expression for the DY cross section in QCD factorization can be written as where f H 1 i and f H 2 j are the proton PDFs, which depend on the longitudinal momentum fraction x 1 and x 2 of parton i and j respectively, and on the factorization scale µ F .dσ is the hard scattering cross section, which is perturbatively calculable in QCD, µ R is the renormalization scale, and Λ 2 QCD /Q 2 represents power-suppressed contributions where Λ QCD is the QCD scale.The calculation of the differential distributions at NNLO in QCD which we present in this work has been performed by using an amended version of the MCFM-v9.0 computer program [64][65][66][67][68], which has been modified to incorporate the string-derived model with charge assignments in Table 1, the Z contribution, as well as the interference terms.We validated this implementation against other computer programs such as DY-Turbo [69] and FEWZ [70,71] and found agreement within 1%.Our results are presented at √ S = 13 TeV of center-of-mass energy and we used the CT18NNLO [72] PDFs with conservative uncertainties evaluated at the 90% confidence level (CL).As a case study, we have chosen a string-derived Z with M Z = 5 TeV.The impact of other recent PDF determinations [73][74][75] on A F B distributions in the context of extra resonance searches is studied in refs.[26,27].Scale uncertainties are obtained by using the 7-point variation, that is, varying µ R and µ F up and down independently by a factor of 2, and then taking the envelope.
For the scope of this analysis, it is sufficient to include in the total decay rate of the Z entering the cross section calculation, only the major decay channels which we report in the expression below where f runs over the quarks and leptons, W ± are the charged gauge vector bosons, and H ± are the Higgs bosons of the charged sector.The expression for the Γ Z →H + H − channel in presence of gauge kinetic mixing is given in Appendix B. The calculation of the partial rates relative to the other channels proceeds similarly to that in refs.[57,59]: in the expression for the Z decay rate in the fermion channel, the vector and axial vector couplings are replaced by those in Eq. 28, while in the expression for the rate in the W + W − channel ε is replaced by ε from Eq. 27.

Kinematic distribution results
In Fig. 2 we show the string-derived Z results for the dilepton invariant mass (m ll ) DY spectrum at the LHC 13 TeV in the electron channel, for a Z of mass M Z = 5 TeV and different values of the coupling g Z and kinetic-mixing parameter g KM .These theory predictions are compared to the SM and the cross sections are calculated in the full phase space.The error bands represent the CT18NNLO induced PDF uncertainty evaluated at the 90% CL.Central predictions are represented by lines with different dashing.We observe how the interplay between g KM and g Z modifies the shape of the resonance and in particular the width for g Z = 0.2 and g Z = g Y , when the strength of g KM is varied.The inclusion of kinetic mixing has a significant effect on the width of the resonance in the invariant mass distribution.
It is interesting to explore the impact of varying these parameters on more differential observables such as the Z transverse momentum p T spectrum and the Z rapidity distribution y Z , which we illustrate in Fig. 3. Here, the invariant mass of the final-state dilepton pair is chosen to be 4.6 < m l l < 5.4 TeV.In the left-column insets, where g KM = 0, PDF uncertainties are represented by red-hatched bands while scale uncertainties are represented by blue bands.In the right-column insets we illustrate for clarity the same distributions with no uncertainties but with g KM > 0. We observe that kinetic mixing has negligible impact on these two distributions for this choice of the parameters in the kinematic region 4.6 < m l l < 5.4 TeV.Again, we observe that PDF uncertainty dominates.δ PDF , SM δ PDF , g Z′ = 0.2 δ PDF , g Z′ = g Y δ PDF , g Z′ = 0.5 δ PDF , SM δ PDF , g KM = 0 δ PDF , g KM = 0.1 δ PDF , g KM = 0.2    Finally, in Fig. 4 we show the impact of varying the g Z and g KM parameters on the cos θ distribution, and on the pseudorapidity η e spectrum of the final-state electron.The angle θ is defined in the Collins-Soper frame [76].While there is almost no impact on the η e spectrum, a distortion of the central value of the cos θ angular distribution is observed when g KM is progressively increased from 0 to 0.2.As we shall see in the next section, this is reflected by the A F B distribution.However, all these effects in the cos θ distribution are buried by the almost complete overlap of the PDF uncertainty bands for g KM = 0, 0.1, 0.2.

Forward-Backward Asymmetry A F B distribution results
As anticipated in the previous sections, a more suitable variable to explore the Z parameters in DY is the forward-backward asymmetry A F B distribution.This observable is a function of the chiral quark and lepton couplings to the mediating gauge bosons and is very sensitive to Z properties.A F B is expressed in terms of the difference between the forward and backward angular contributions to the cross section, normalized to the total cross section.The forward and backward directions are represented by the angle θ between the negatively charged finalstate lepton and the initial-state quark in the dilepton center-of-mass frame.A F B is defined as where the forward (dσ F ) and backward (dσ B ) contributions are constructed by integrating the differential cross section over the forward and backward halves of the angular phase space θ l is the angle between the final-state lepton and initial-state quark in the Collins-Soper frame.
In our analysis, we study the A F B distribution as a function of m ll .In Fig. 5 we show the A F B over the same invariant mass range adopted before.The presence of a Z would lead to a distortion in the A F B which becomes wider as g Z increases.We calculated both PDF and scale uncertainties for the A F B distribution.Scale uncertainties are essentially invisible in Fig. 5, while PDF uncertainties are attenuated around the peak as compared to the distributions studied in the previous section, but are still large at high m ll .Fig. 5 shows the sensitivity of A F B to simultaneous variations of the g Z and g KM parameters.The distortions induced on the asymmetry tend to become shallower and wider as the kinetic mixing increases.Moving away from the resonance, when g KM = 0, the discriminatory power of A F B appears to be attenuated by PDF uncertainties.

Conclusions
In this work, we studied properties of heterotic string derived Z s in the TeV scale and explored their dynamics at the LHC.As a case study, we selected a Z of mass M Z = 5 TeV.We analyzed the impact of a gauge kinetic-mixing term in the Lagrangian and its interplay with the Z coupling.We explored the Z parameter space the model in presence of kinetic mixing and performed a thorough phenomenological analysis in which we used precision calculations at NNLO in QCD to predict kinematic distributions for the Z in the Drell-Yan process at the LHC.
In particular, we exploited the sensitivity of forward-backward asymmetry A F B distributions in Drell-Yan to further explore the parameter space and the interplay between the Z coupling and the kinetic-mixing parameter.Moreover, we estimated theory uncertainties on the distributions from perturbative (scale dependence) and nonperturbative (PDFs in the proton) sources in QCD.Proton PDFs still remain one of the major sources of uncertainty in Z searches.This is ascribed to the fact that the kinematic domain of heavy gauge boson production is sensitive, and is impacted by, PDFs at large x where uncertainties are still large due to lack of robust constraints from experimental data [77,78].The motivation to consider the particular U (1) Z combination in Eq. 3 stems from its extraction from a string derived heterotic-string model [9], in which it is the unique combination that may remain unbroken down to low scales.Its preservation as an unbroken gauge symmetry down to the TeV scale emanates from the role that it can play in suppressing some dangerous operators, like proton decay mediating operators, and furthermore by the role it can play in generating the electroweak symmetry breaking itself.Further analysis in that direction mandates the development of the interpolation tools between the string and electroweak scales and the extraction of further data from the string models.We note that while in this paper we treated the kinetic-mixing term as a free parameter, its computation, as that of many of the related parameters, can be obtained directly from the string models [15], hence increasing their predictive power.In this paper, we studied the case of Z vector bosons that have an E 6 embedding.The advantage of this choice is that it facilitates gauge unification at the unification scale [8].Such extra string inspired Z models have been discussed extensively in the literature since the mid-eighties [19,[32][33][34], but the combination given in Eq. 3 is obtained in the string derived model of ref. [9].Naturally, the string constructions can give rise to U (1) symmetries with different characteristics and these δ PDF , SM δ PDF , g KM = 0 δ PDF , g KM = 0.1 δ PDF , g KM = 0.2 have been of some interest in the literature [7,[36][37][38].The range of possibilities is more model dependent.The extra Z s in these cases can arise from flavour dependent U (1) symmetries, as well as U (1) symmetries that arise from the hidden sector of the heterotic-string.Naturally, the limits that we discussed in this paper will not apply to those cases (see for instance the discussion relative to Z boson searches in the PDG [60] and references therein), as they do not couple universally to all the three families.On the other hand, there may exist other constraints in these cases arising from flavour non-universality and the need to break flavour dependent symmetries to produce viable fermion masses.Nevertheless, these cases do represent interesting alternatives to the family universal U (1) and we will return to them in future work.

Figure 1 :
Figure 1: Z -induced corrections on the SM Z boson mass as a function of g Z for fixed values of M Z .The hatched band represents the current measured M Z uncertainty as reported in the PDG [60].

Figure 2 :
Figure 2: Z invariant mass distributions for varying g Z and g KM .The error bands represent the CT18NNLO induced PDF uncertainty on the cross section at 90% CL.Central predictions are represented by lines with different dashing.

Figure 3 :δδ
Figure 3: Transverse momentum p T (upper) and rapidity y Z of the Z (lower).The CT18NNLO PDF uncertainty is shown as a hatched red band and scale uncertainty is shown as a solid blue band.Central predictions are represented by a solid black line.

Figure 4 :
Figure 4: Left: Differential cross-section as a function of cos θ, where the θ angle is defined in the Collins-Soper frame.Right: Pseudorapidity distribution η e of the outgoing electron.Error bands with different hatching represent the CT18NNLO PDF uncertainties induced on the cross section, evaluated at the 90% CL.Central predictions are represented by lines with different dashing.

Figure 5 :
Figure 5: Forward-backward asymmetry (A F B ) shown as a function of the dilepton invariant mass (m ll ) for different values of g Z and g KM .The black hatched band represents the Standard Model prediction.The solid black line at A F B = 0 is introduced for reference.The CT18NNLO PDF uncertainties at the 90% CL are represented by hatched bands.Central predictions are shown as lines with different dashing.