Standard Model Parton Distributions at Very High Energies

We compute the leading-order evolution of parton distribution functions for all the Standard Model fermions and bosons up to energy scales far above the electroweak scale, where electroweak symmetry is restored. Our results include the 52 PDFs of the unpolarized proton, evolving according to the SU(3), SU(2), U(1), mixed SU(2) x U(1) and Yukawa interactions. We illustrate the numerical effects on parton distributions at large energies, and show that this can lead to important corrections to parton luminosities at a future 100 TeV collider.


Introduction
Experiments at the Large Hadron Collider are now probing the structure of matter at scales comparable with, and even beyond, the characteristic scale of electroweak symmetry breaking. So far, no evidence has been found for a breakdown of the Standard Model (SM) in particle collisions. Indeed, there is a logical possibility that the SM remains a good description of hard scattering processes up to scales far beyond those of any conceivable particle colliders. It is therefore of interest to examine the features predicted by the SM for collider events well above the electroweak scale. For this purpose, Monte Carlo event generators including all the SM interactions on an equal footing are necessary. Such generators would be useful for investigating the limits of LHC searches, the potential of possible future colliders and cosmic processes at ultrahigh energies.
To construct a general-purpose SM event generator, 1 the three phases of a hard collision, namely initial-state parton showering, parton-parton collision and final-state showering, need to be simulated including all SM particles and interactions. For the initial-state showering, parton distribution functions (PDFs) for all the SM fermions and bosons need to be computed and tabulated beforehand, so that showering can be generated backwards from the hard process, guided by the scale dependence of the PDFs [2,3].
Recently, a final-state parton shower including emissions from all interactions in the Standard Model was developed [4], which illustrated the importance of electroweak splittings at high energies. For initial-state radiation the generalization of the DGLAP [5][6][7] evolution equations using all the Standard Model interactions has been worked out in [8], but so far no numerical implementation of these results has been published.
As already mentioned, understanding the DGLAP evolution of PDFs using all interactions of the SM is a required first step in developing a complete initial state parton shower. Moreover, it already allows us to study many new qualitative features of very high-energy processes, such as lepton-initiated processes in hadron collisions and the polarization induced by electroweak PDF evolution.
The PDF evolution equations for the full Standard Model have been presented in Ref. [8]. In the present paper we recast those equations in a form suitable for event generation and solve them numerically for a given set of input distributions at the electroweak scale. The resulting PDF set extends through the region of interest for future colliders and well beyond, so that we can study the onset of the regime where all the SM interactions start to become comparable.
Our solutions to the SM evolution equations are obtained in the approximation of exact SU(3)×SU(2)×U(1) symmetry. That is, we neglect fermion and Higgs masses and the Higgs vacuum expectation value, the effects of these being power-suppressed at high scales. We impose an infra-red cutoff m V on interactions that involve the emission of an electroweak vector boson, V = W i for SU (2) or B for U (1). Leading-order evolution kernels and one-loop running couplings are used. All the electroweak PDFs are generated dynamically from the QCD plus photon PDFs, starting from a matching scale q 0 ∼ m V . In practice we take q 0 = m V = 100 GeV. For the evolution of the photon, we decompose its PDF into W 3 , B and mixed B/W 3 components at the input scale, evolve these components, and reconstruct the photon PDF from them at higher scales using the running SU(2) and U(1) couplings. For the top quark, we set the PDF to zero below the top mass scale and then use the leading-order massless evolution kernels, as for other fermions. This treatment of the transition region around the electroweak scale is clearly over-simplified but it should give a reliable indication of the magnitude of electroweak effects at higher energies.
The accuracy of our resulting PDFs is leading logarithmic, with subleading logarithmic effects included where possible, but not in a complete way. Contributions to the evolution from the U(1), SU(3) and Yukawa interactions are therefore correct at the single logarithmic level. However, as mentioned above, the SU(2) interactions give rise to double logarithmic effects in the PDF evolution, such that single logarithmic effects in SU(2) non-singlet quantities are not fully under control.
The organization of the paper is as follows. In Sec. 2 we define the relevant parton distribution functions for unpolarized proton beams and the general form of their evolution equations, paying particular attention to the conservation of momentum in the presence of the cutoff m V for vector boson emission. After specifying all the necessary splitting functions and running couplings, we write the explicit evolution equations associated with the five interactions: SU(3), U(1), SU(2), Yukawa and mixed U(1)×SU (2), for all the SM partons in a flavor basis. As usual for DGLAP evolution, we do not include 4 point interactions which are suppressed at high energies.
For a numerical implementation, as described in Sec. 3, the flavor basis is not convenient, as too many coupled equations are involved. Instead we use the basis of conserved quantum numbers introduced in Ref. [8]. As shown there, the double-logarithmic evolution of SU(2) non-singlet PDFs can then be factored out, which stabilizes and accelerates the solution of the equations. In this way we are able to evolve all the SM PDFs to arbitrarily high scales with satisfactory speed and precision. In practice we evolve up to 10 8 GeV, where the approach to asymptotic behavior is well established.
In Sec. 4, we present a selection of results that illustrate the extent to which electroweak effects change the behavior of the various PDFs. In particular, we show changes in the PDFs of strongly interacting particles relative to pure QCD evolution, and show the size of the PDFs for electroweak gauge bosons relative to the gluon PDF. Finally, we present results of the associated changes in parton-parton luminosities at a 100 TeV pp collider. Our conclusions are presented in Sec. 5.

Definition of the parton distribution functions
The standard definition of an x-weighted parton distribution is given by the matrix element of a bi-local operator, separated along the lightcone. For fermions, one finds the standard definition, but without spin averaging as we are separating the fermions into left-and right-handed. Thus, each fermion has only one possible spin determined by its helicity and the sign of its momentum where µ is the renormalization scale. Since we have separate left-and right-handed PDFs, for each generation there are a total of 8 quark PDFs and 6 lepton PDFs to consider, giving a total of 42 fermion PDFs. Parton distributions functions of the vector bosons are given by Since SU(3) is unbroken, we consider a single PDF to describe the gluon field. For the SU(2) ⊗ U(1) symmetry, on the other hand, one needs to take the symmetry breaking into account. For the W + and W − boson we simply include separate PDFs for each of the two gauge bosons. For the B and W 3 , however, one needs to be more careful to take the mixed contributions of these two bosons into account. Such contributions arise from the fact that the left-handed fermions and Higgs carry both isospin and hypercharge. This implies that besides B and W 3 PDFs one needs to include a mixed PDF, which is given by 2 From these PDFs one can then construct the PDF for the photon, the transversely-polarized Z 0 and their mixed state as a transformation of the PDF for the B, the W 3 and their mixed state. Using For the electroweak input at scale µ = q 0 we have f γ = 0 and f Z = f γZ = 0, so the input conditions at that scale are After evolving these three unbroken PDFs to a higher scale q, the physical photon and Z 0 PDFs are reconstructed there using the corresponding running values of c W and s W . Finally, one needs to include PDFs for the scalar bosons. One writes and PDFs for each of the 4 Higgs fields H 0 ,H 0 , H + and H − are included. The relationship to the 4 Higgs fields in the unbroken basis to the physical Higgs and the longitudinal gauge bosons is as follows: The H ± PDFs correspond to those of the longitudinally polarized W ± . In the notation of Ref. [8], the neutral Higgs fields are where h and Z L represent the Higgs and the longitudinal Z 0 fields, respectively. The corresponding PDFs are and one can also define the mixed PDFs Both of these mixed PDF carry non-zero hypercharge, such that they are not produced by the DGLAP evolution in the unbroken gauge theory as considered in this paper 3 . Thus, one immediately finds and (2.14) In summary, there are a total of 52 parton distribution functions that need to be considered. Apart from the QCD quark and gluon distributions and the electroweak PDFs (2.6), all the other SM PDFs are set to zero at scale q 0 = m V and evolve according to the generalized DGLAP equations presented below.

General evolution equations
We consider the x-weighted PDFs of parton species i at momentum fraction x and scale q, f i (x, q). In general they satisfy evolution equations of the following forms: 3 They are only produced through insertions of the Higgs vacuum Here, the sum over I goes over the different interactions in the Standard Model and the notation [q ∂/∂qf i (x, q)] I implies that we only keep the terms proportional to the coupling α I when taking the derivative 4 . For the rest of the section, we will show the evolution of each f i (x, q). We choose I = 1, 2, 3 for the pure U(1), SU(2) and SU(3) gauge interactions, I = Y for Yukawa interactions, and I = M for the mixed interaction proportional to The first contribution, proportional to P V i,I , denotes the virtual contribution to the PDF evolution (the disappearance of a flavor i), while the second contribution is the real contribution (the appearance of flavor i due to the splitting of a flavor j). The maximum value of z in the integration of the real contribution depends on the type of splitting and interaction, and we choose that is, we apply an infrared cutoff m V , of the order of the electroweak scale, when a B or W boson is emitted. This regulates the divergence of the splitting function for those emissions as z → 1. Such a cutoff is mandatory for I = 2 because there are PDF contributions that are SU(2) non-singlets. The evolution equations for SU (3) are regular in the absence of a cutoff, as hadron PDFs are color singlets. Similarly for U(1), the unpolarized PDFs have zero hypercharge, 5 but we include the same cutoff for I = 1, since the B and W 3 are mixed in the physical Z and γ states.
Note that the precise choice of the cutoff is somewhat arbitrary, and as already mentioned, we choose m V = 100 GeV in this paper. Changing this value changes our results by subleading logarithmic effect, at the same level as other effects not included. However, given that the SU(2) evolution is double logarithmic, this implies that the ambiguity is single logarithmic for the SU(2) coupling. By matching our results to fixed order, one would account for these term at first order in α 2 . This is beyond the scope of this paper.
While the flavor basis chosen above is the most intuitive basis, the fact that all 52 PDFs are coupled to one another makes it quite difficult to solve the evolution equations. To decouple some of the equations, it helps to change the basis such that the ingredients have quantum numbers that are conserved in the Standard Model. Choosing the total isospin T and CP as the quantum numbers, the PDFs for each set of quantum numbers required are shown in Table 1.
Note that in general there can be additional mixed PDFs, which however are zero in our initial conditions and which are not generated in the evolution. In particular, there can be states mixing left-and right-handed fermions, but they are not present in the initial condition when only considering unpolarized beams because those states are not Lorentz scalar. Thus, we can drop these states from our evolution. 4 Note that [. . .] I is only introduced for notational convenience and should not be interpreted as setting all other couplings to zero. In particular, the PDFs appearing on the right-hand side of Eq. (2.15) still depend on the value of all coupling constants 5 Although there can be contributions with non-zero hypercharge for transversely polarized beams [8]. The sum of momenta of all non-mixed PDFs in the particle basis is conserved, since it is the momentum of the proton. Momentum conservation applies independently for each interaction This is equivalent to the sum over all T = 0, CP = + PDFs in the isospin and CP basis because only these states contribute to a sum over the PDFs in the particle basis. For the other values of T and CP, the PDFs correspond to differences of PDFs in the particle basis.
For example an isospin 1 PDF is added in PDF of an up-type fermion, but subtracted in the down-type PDF, thus it has no effect on the sum. Combining Eqs. (2.15) and (2.18) gives where we have defined the momentum averaged PDF Solving the equation for each of the f i (q) , since all the input particle PDFs can be set independently, we get Thus, momentum conservation determines the factor P V i,I for all non-mixed fields in the particle basis.
Note that the result from momentum conservation agrees up to power corrections with the more traditional definition of the virtual corrections as loop insertions on the fields of the PDF. Summing over possible loops, one has where and similarly for C V V,I and C HH,I . The sums in Eq. (2.22) extend over particles, and not their anti-particles. To see that Eqs. (2.21) and (2.22) agree with each other, we will work it out explicitly for the virtual contribution to a fermion. One uses for the fermions that and C f f,I = C V f,I to obtain the correct relation: where . . . denotes power corrections in 1 − z max . The argument is exactly the same for P V H,I (q), while for P V V,I (q) one simply uses that P R V V,I (z) and P R f V,I (z), and P R hV,I (z) and P R f H,I (z), are symmetric in z ↔ 1−z to write z dz = dz/2. In our implementation of the evolution equations, we use Eq. (2.21), to ensure exact momentum conservation without explicit power corrections.
Since the mixed PDF f BW is a pure T = 1 state, it does not contribute to the momentum sum. This implies that one cannot derive its associated virtual contribution from momentum conservation. However, using the traditional definition in terms of loops, one sees that in this case the U(1) and SU(2) virtual corrections each apply to only one of the two fields involved, and thereforẽ while the virtual contribution is zero for the other interactions. One can simplify the general evolution equations in Eq. (2.15) by defining a full Sudakov factor as well as a partial Sudakov factor for each interaction where q 0 is an arbitrary cutoff, which for convenience we set equal to m V . This allows us to write where again the notation [. . .] I implies that only terms from the interaction I are kept. This gives

Splitting functions
The splitting functions depend only on the type of particles, which for the Standard Model are the spin 1/2 fermions, denoted by f , spin 1 gauge bosons, denoted by V , as well as spin 0 Higgs bosons, denoted by H.
Denoting the three gauge interactions of the Standard Model collectively by I = G, the splitting functions involving gauge bosons are given by The factor of 1/2 in P f V has to be included since we are considering fermions with definite chirality. For the Yukawa interaction (Y ), one obtains

Running couplings
The one-loop running of the gauge couplings α I (I = 1, 2, 3) is given by where, for n g generations and n H Higgs doublets, At scale M Z = 91.2 GeV we take We set all Yukawa couplings to zero, except for the top Yukawa coupling α Y = y 2 t /4π. Its running receives significant Yukawa and QCD contributions: where β Y = 9/2 and β S = 8. The solution is We take m t (m t ) = 163 GeV, which implies α Y (m t ) = 0.0349, and α 3 (m t ) = 0.109.

I = 3: SU(3) interactions
We start by considering the well known case of SU(3) interactions. The relevant degrees of freedom are the gluon, as well as left and right-handed quarks. The coupling constants are (with C F = 4/3, C A = 3, T R = 1/2) This gives for the evolution of a quark or gluon 6 The Sudakov factor can be obtained from Eq. (2.21) using the coupling constants in Eq. (2.53). This gives where the hypercharges of the different fermions are given by where h here stands for any of the four Higgs boson PDFs. 6 From now on we omit the arguments of functions for brevity.
Plugging this into the general evolution equation gives The virtual splitting functions, required for the Sudakov factor are given by where we have used in the second line that for each generation there are 4 left-handed quarks (one needs to count particles and antiparticles separately), 2 right-handed up-type quarks, 2 right-handed down-type quarks, 4 left-handed leptons and 2 right-handed electrons, and that there are a total of 4 Higgs bosons.

I = 2: SU(2) interactions
The SU (2) interactions are more complicated, since the emission of W ± bosons changes the flavor of the emitting particle. This, combined with the SU(2) breaking in the input hadron PDFs, leads to double-logarithmic scale dependence in the DGLAP evolution, rather than only single-logarithmic dependence as in the evolution based on U(1) and SU(3). The relevant coupling constants are (where u L and d L denote any up-and down-type left-handed fermion) where as before the color factor N f = 3 for quarks, 1 for leptons. The couplings of the W 3 state to the Higgs are given by where again h stands for any of the 4 Higgs bosons, while those of the charged W states are given by The couplings for the charge-conjugate states are the same.
This gives for the evolution of the fermions For the W + and W 3 bosons we have where the sum in the last line is over all left-handed fermions and anti-fermions. The equation for the W − can be obtained from that of the W + by charge conjugation.
Finally, for the Higgs bosons we have The virtual splitting functions are from which the Sudakov factor can be constructed using Eq. (2.29). An important aspect of the SU(2) evolution equations is that, contrary to the other gauge groups, the dependence on the ratio m V /q does not cancel between the real and virtual splitting functions. As an example, consider the evolution equation for an up-type fermion, given on the first line of Eq. (2.74), with the virtual contribution given by the first line of Eq. (2.80). The sum of the contributions of real and virtual splitting functions is given by where . . . represents less singular terms. Thus, the SU(2) breaking in the proton, which renders f u (z) = f d (z), gives rise to a logarithmic dependence on m V /q, which leads to a double-logarithmic dependence upon integration over q. As we will see later, the effect of this dependence is to double-logarithmically suppress the SU(2) breaking effects at high energies.

I = Y : Yukawa interactions
The interaction of Higgs particles with fermions is described by the Yukawa interactions.
In this work we only keep the top Yukawa coupling, setting all others to zero. This gives the following couplings where q 3 L denotes either the left-handed top or bottom quark. We furthermore need This gives contributions to the top quark PDFs, as well as the left-handed bottom PDF: It also contributes to the evolution of the Higgs bosons: The Sudakov factors can be obtained using Eq. (2.29) with Finally, we need to consider the evolution involving the mixed BW boson PDF. The nonvanishing couplings are where, as for the fermions, the same relations hold for the charge-conjugate states.
Plugging these into the general evolution equation gives As already discussed, the mixed gauge field PDF f BW has Sudakov factors associated with the U(1) and SU(2) interactions, given by Eq. (2.29). Since there is no corresponding real emission term in the evolution equation for f BW , it evolves double-logarithmically and is suppressed at high scales relative to the unmixed PDFs.

Implementation details
Our treatment assumes that the SM PDFs at very high energies can be obtained by smoothly matching the broken and unbroken symmetry regimes at a matching scale q 0 ∼ m V , which in practice we take to be 100 GeV. Our input PDFs at 100 GeV are obtained as follows. We take the CT14qed PDF set [15] at 10 GeV and replace the photon PDF by that of the LUXqed set [16]. We do not use the CT14qed photon because the LUXqed photon, while being consistent with CT14qed, has much smaller uncertanties and a smoother x dependence. The LUXqed PDF set combines the PDF4LHC15 nnlo 100 parton set [31] with a determination of the photon PDF from structure function and elastic form factor fits in electron-proton scattering. However, we do not use the LUXqed partons, because being NNLO they are not positive-definite, which we require for our LO treatment and is satisfied by CT14qed. We evolve this hybrid CT14-LUX PDF set from 10 to 100 GeV using leading-order QCD plus QED evolution, which incidentally generates the charged leptons. The resulting parton, photon and lepton PDFs form our input to the unbroken SM evolution upwards from 100 GeV. The input left-and right-handed fermion PDFs are identical. The input W 3 , B and mixed B/W 3 PDFs are determined by the photon (and the absence of the Z 0 ) at the matching scale according to Eq. (2.6). The remaining vector boson, neutrino and Higgs PDFs are all generated dynamically starting from zero at the matching scale.
The equations given in Sections 2.5 to 2.9 completely define the evolution of all parton distribution functions in the unbroken symmetry regime. However, as already explained, one can rewrite the equations slightly to make them more amenable to a numerical implementation. First, switching to a basis of states with well-defined isospin decouples the set of 52 equations to some degree. In this new basis another transformation eliminates the double logarithmic sensitivity to the ratio m V /q. Second, by combining the virtual and real splitting functions into +-distributions, one can reduce numerical sensitivity to the cutoff of the z integrations. We will now discuss these simplifications in turn.

Switching to a basis of conserved quantum numbers
As we already explained in Section 2.2, the set of 52 evolution equations can be decoupled to some degree by switching to a basis of well-defined isospin T and CP. Writing a fermion PDF with T and CP as f TCP i , we write the left-handed fermions as where u L and d L refer to left-handed up-and down-type fermions. Right-handed fermions are given by The SU(3) and U(1) boson PDFs have T = 0, CP = + while the SU(2) boson PDFs can have T = 0, 1, 2 with respectively CP = +, −, + The mixed BW boson state is a combination of 0 − and 1 − and therefore its PDF has For the Higgs boson, one writes similarly to the fermions In terms of these states the longitudinal vector boson and Higgs PDFs are then, using Eq. (2.14),

Cancellation of double-logarithmic dependence in evolution equations
In the {T, CP} basis the singular contributions to the evolution equations (those that are proportional to the splitting functions P R f f,G (z), P R V V,G (z) and P R HH,G (z), which diverge in the limit z → 1) are diagonal, such that the PDF multiplying the divergent splitting function is the same as that appearing on the left-hand side. Here, as in f TCP i , the label i now refers to a parton species f, V, H rather than a particular parton. Recalling that the Sudakov factor takes the form where . . . represents less divergent terms, and where k and l are particular partons, we have where The +-prescription defined by Eq. (3.17) regulates the divergence in the integrand as z → 1 and therefore if we define the modifying factor then the evolution equation (3.13) becomes 19) with no logarithmic dependence on m V /q on the right-hand side. For all interactions except SU(2), one can show that D TCP i,I = C i,I , so that the modifying factor (3.18) is unity 7 . For SU(2) we have explicitly: 8 while so that For the mixed PDF f BW we have D 1+ BW,2 = 0 and therefore The equations finally used to evolve the PDFs in the conserved-quantum-number basis are given in Appendix A.

Results
We begin by showing how the PDFs of strongly interacting particles are changed by including the evolution of the full Standard Model. Figure 1 shows results on the evolution of left-and right-handed quark PDFs, shown solid and dashed respectively, normalized to their values assuming pure QCD evolution. In each plot we show the results at three different scales, namely q = 10 4 GeV, q = 10 6 GeV and q = 10 8 GeV. The values of 10 6 and 10 8 GeV are of course far away from energy scales one can reach at any collider in the near or distant future. However, showing the results at such unattainable values helps to illustrate their approach to asymptotic behavior.
All the light quarks (and antiquarks, not shown) evolve to lower values compared to pure QCD at small x, due to an overall loss of energy to the electroweak gauge bosons through the additional splittings q → qW and q → qB. At higher x values, the up and down quarks (top row) exhibit different behaviors, with the left-handed up PDF evolving more rapidly to lower values compared to pure QCD, while the down quark eventually evolves to higher values. This is because the left-handed up and down distributions evolve towards each other, their difference being double-logarithmically suppressed at high scales. The right-handed quark PDFs have no double-logarithmic component and evolve to slightly  lower values than pure QCD, due to energy loss through the additional splitting q R → q R B.
The asymmetry between left-handed charm and strange quarks also evolves doublelogarithmically towards zero, primarily through a more rapid decrease of the strange PDF.
At high x the behavior is more complicated because the input CT14qed charm PDF is larger than the strange above x ∼ 0.7. The right-handed quarks behave qualitatively the same as those of the first generation.
The left-handed top and bottom quarks also must evolve towards equal values, which in this case means that the top has higher values that in pure QCD, while the bottom evolution looks similar to strange, relative to pure QCD. The right-handed b-quark behaves qualitatively like the right-handed quarks of the first and second generation, while the righthanded top quark, being generated purely dynamically, behaves differently at large x. Since the right-handed top has vanishing initial condition, the splitting t R → t R B, which would decrease the PDF, is sub-dominant compared to the process B → t RtR . This means that at large x the right-handed top PDF is increased, rather than decreased.
The effect on the gluon PDF is shown in last row of Fig. 1. While the effects are quite small up to q ∼ 10 4 GeV, at larger scales the back-reaction from the changing quark PDFs is affecting the gluon PDF at an appreciable level. It is interesting to study how rapidly electroweak symmetry is restored. To illustrate this, we show the asymmetry compared to the result if only QCD evolution were turned on. This asymmetry ratio is shown in Figure 2 for the three generation of quarks as a function of q, for various values of x. For all generations the asymmetry decreases as q gets larger, driving the PDFs of the different isospin states towards each other. The onset of the deviation from pure QCD is in the range 1 − 10 TeV. The ratio between the full asymmetry and the result using only QCD evolution is given by where ∆ f,2 (q) is the fermion Sudakov factor, as given in Eq. (3.18), independent of the generation. Next, we study the size of the PDFs of particles not charged under the strong interaction. Since these PDFs are only generated by emissions due to the U(1), SU (2) or Yukawa interactions, they are vanishing at all scales if one is including only SU(3) evolution. The only exception is the photon, which has a non-vanishing initial condition at q = 100 GeV. Figure 3 shows results on the electroweak boson PDFs normalized to the gluon PDF, both evolved using the full Standard Model. One can see that the electroweak gauge boson PDFs become a significant fraction of the gluon PDF, especially at large values of x. The  photon PDF is the largest mainly because it has a non-zero input. The PDF for the W + boson is initially larger than the W − boson PDF at large x because the W + is mainly generated through emissions from the up-quark, whose PDF is larger than the down-quark which mainly generates the W − . Since the difference between W + and W − has isospin 1, the W + evolves more slowly and the W − more rapidly, so that they approach each other at high q. At low x they are more similar as are the up-quark and down-quark PDFs. The Z 0 PDF is similar to the W + but it is smaller at low x and larger at large x. The mixed γZ PDF is small and positive at small x and negative at large x. There is no constraint to be positive definite for a mixed PDF as it is the product of two amplitudes rather than the square modulus of one. Its absolute value becomes very large at large x and q.
We also show the PDFs for the longitudinally polarized gauge bosons, the Higgs boson, the mixed PDF between the Higgs and the Z L and the leptons. The Z L PDF is the same as the Higgs in our approximation, see Eq. (2.14), so we do not make a separate plot for it. The boson PDFs are shown in Fig. 4, and the leptons in Fig. 5, both normalized to the gluon. Both are expected to be much smaller than the transverse vector boson PDFs, because they are generated via a second order effect of emission from the vector bosons and via Yukawa emission from the top and bottom quarks, which are much smaller than the up and down quarks. The mixed PDF is even smaller because it is generated by the asymmetry between transverse W + and W − PDFs and the top and anti-top PDFs. The W + L and W − L PDFs are very similar, for the same reason.   As a final result, we study several parton luminosities, choosing a future 100 TeV pp collider as a reference. While the energy scales that can be reached at such a collider are not quite large enough to get O(1) effects, the effects of the full Standard Model evolution are still numerically relevant. In Figure 6 we show the q LqL luminosities for the six different quark flavors, normalized to their values if only QCD evolution is taken into account. One can see that all except the tt luminosity are reduced appreciably from their values if only QCD evolution were taken into account. This will affect searches for Z -like particles at a future 100 TeV collider. The dd luminosity is decreasing more slowly as the double-logarithmic evolution drives it larger than QCD at high x (see Fig. 1).
We also show selected luminosities of vector bosons combined with quarks, normalized to the average of the uū and dd luminosities. One can see that luminosities involving one transverse vector boson become of comparable magnitude to the qq luminosities. Luminosities involving the longitudinal gauge and Higgs bosons are much smaller.

Conclusions
The energy regime around and beyond the electroweak scale is currently being explored by the LHC experiments, and so far they have found no firm evidence for physics beyond the Standard Model. In the present paper, we have examined the consequences of assuming that the parton distributions of the proton continue to be described by the Standard Model up to very high energies, in the approximation that its symmetries are unbroken above the electroweak scale.
We have implemented numerically the full set of generalized DGLAP evolution equations for all the parton species and interactions of the unbroken SM in leading order. The input PDFs of 5 quark flavors, the gluon, photon and charged leptons at a starting scale q 0 = 100 GeV for the full SM evolution are obtained from parton and photon PDFs at 10 GeV by QCD plus QED evolution. The input left-and right-handed fermion PDFs are thus identical at scale q 0 but they evolve differently above that scale. The top quark PDFs (not present in the input) start to evolve from the top mass scale. The input photon is resolved into its U(1), SU(2) and mixed components, which are evolved independently from scale q 0 and reassembled into the photon and transversely polarized Z 0 at higher scales. The charged and longitudinal vector boson, Higgs and neutrino PDFs are generated dynamically starting from zero at scale q 0 . This simplified treatment misses some symmetry-breaking effects around the electroweak scale, but these are power-suppressed at higher scales and our results should provide a guide to the ways in which the PDFs deviate from pure QCD evolution.
Amongst the most interesting features of the SM is the distinction between left-and right-handed fermions. The evolution of the right-handed PDFs deviates little from pure QCD, owing to the weakness of the U(1) interaction. The left-handed PDFs generally deviate from pure QCD at the 5-10% level by 10 TeV.
Another important SM characteristic is the restoration of isospin symmetry at high scales. This is manifest in the decreasing asymmetry between the up-and down-type quark PDFs, which sets in at 1 − 10 TeV, the up-type being pulled down in the first generation and conversely in the third. The suppression of the asymmetry is a double-logarithmic effect that can be treated in fixed order at present energies but is resummed to all orders in the evolution.
The electroweak bosons are generated quite copiously, the W + in particular at high x due to splitting of valence up quarks into dW + . The photon and Z 0 PDFs also grow rapidly, eventually exceeding the gluon at high x. The PDFs of the longitudinal vector bosons, the Higgs boson and the leptons are generally much smaller as they arise from second-order splittings.
Finally, we have used the generated PDFs to present some parton-parton luminosities at a 100 TeV pp collider. These results are just an illustration of the size of the effects that can be expected at such a future collider, and a more detailed phenomenological analysis will be presented in a forthcoming publication.
In conclusion, we find a rich structure in the proton when probed beyond the electroweak scale. The associated PDFs are interesting and useful in their own right. They also represent a key component of event generators that aim to embody the full Standard Model in initial-state parton showering, a topic we plan to explore further.  • T = 1 and CP = −: