We update our estimates of charged and neutral current neutrino total cross sections on isoscalar nucleons at ultrahigh energies using a global (x, Q^2) fit, motivated by the Froissart bound, to the F_2 (electron-proton) structure function utilizing
the most recent analysis of the complete ZEUS and H1 data sets from HERA I. Using the large Q^2, small Bjorken-x limits of the wee parton model, we connect the ultrahigh energy neutrino cross sections directly to the large Q^2, small-x extrapolation of our new fit, which we assume saturates the Froissart bound. We compare both to our previous work, which utilized only the smaller ZEUS data set, as well as to recent results of a calculation using the ZEUS-S based global perturbative QCD parton distributions using the combined HERA I results as input. Our new results substantiate our previous conclusions, again predicting significantly smaller cross sections than those predicted by extrapolating pQCD calculations to neutrino energies above 10^9 GeV.
Using repeated Laplace transform techniques, along with newly-developed accurate numerical inverse Laplace transform algorithms, we transform the coupled, integral-differential NLO singlet DGLAP equations first into coupled differential equations, th
en into coupled algebraic equations, which we can solve iteratively. After Laplace inverting the algebraic solution analytically, we numerically invert the solutions of the decoupled differential equations. Finally, we arrive at the decoupled NLO evolved solutions F_s(x,Q^2)=calF_s(F_{s0}(x),G_0(x)) and G(x,Q^2)=calG(F_{s0}(x),G_0(x)), where calF_s and calG are known functions - determined using the DGLAP splitting functions up to NLO in the strong coupling constant alpha_s(Q^2). The functions F_{s0}(x)=F_s(x,Q_0^2) and G_0(x)=G(x,Q_0^2) are the starting functions for the evolution at Q_0^2. This approach furnishes us with a new tool for readily obtaining, independently, the effects of the starting functions on either the evolved gluon or singlet structure functions, as a function of both Q^2 and Q_0^2. It is not necessary to evolve coupled integral-differential equations numerically on a two-dimensional grid, as is currently done. The same approach can be used for NLO non-singlet distributions where it is simpler, only requiring one Laplace transform. We make successful NLO numerical comparisons to two non-singlet distributions, using NLO quark distributions published by the MSTW collaboration, over a large range of x and Q^2. Our method is readily generalized to higher orders in the strong coupling constant alpha_s(Q^2).
We recently derived an explicit expression for the gluon distribution function G(x, Q^2) = xg(x, Q^2) in terms of the proton structure function F_2^{gamma p} (x, Q^2) in leading-order (LO) QCD by solving the the LO DGLAP equation for the Q^2 evolutio
n of F_2^{gamma p} (x, Q^2) analytically, using a differential-equation method. We re-derive and extend the results here using a Laplace-transform technique, and show that the singlet quark structure function F_S(x,Q^2) can be determined directly in terms of G from the DGLAP gluon evolution equation. To illustrate the method and check the consistency of existing LO quark and gluon distributions, we used the published values of the LO quark distributions from the CTEQ5L and MRST2001LO analyses to form F_2^{gamma p} (x, Q^2), and then solved analytically for G(x,Q^2). We find that the analytic and fitted gluon distributions from MRST2001LO agree well with each other for all x and Q^2, while those from CTEQ5L differ significantly from each other for large x values, x>~0.03 - 0.05 at all Q^2. We conclude that the published CTEQ5L distributions are incompatible in this region. Using a non-singlet evolution equation, we obtain a sensitive test of quark distributions which holds in both LO and NLO perturbative QCD. We find in either case that the CTEQ5 quark distributions satisfy the tests numerically for small x, but fail the tests for x>~0.03 - 0.05 - their use could potentially lead to significant shifts in predictions of quantities sensitive to large x. We encountered no problems with the MRST2001LO distributions or later CTEQ distributions. We suggest caution in the use of the CTEQ5 distributions.
We derive a second-order linear differential equation for the leading order gluon distribution function G(x,Q^2) = xg(x,Q^2) which determines G(x,Q^2) directly from the proton structure function F_2^p(x,Q^2). This equation is derived from the leading
order DGLAP evolution equation for F_2^p(x,Q^2), and does not require knowledge of either the individual quark distributions or the gluon evolution equation. Given an analytic expression that successfully reproduces the known experimental data for F_2^p(x,Q^2) in a domain x_min<=x<=x_max, Q_min^2<=Q^2<=Q_max^2 of the Bjorken variable x and the virtuality Q^2 in deep inelastic scattering, G(x,Q^2) is uniquely determined in the same domain. We give the general solution and illustrate the method using the recently proposed Froissart bound type parametrization of F_2^p(x,Q^2) of E. L. Berger, M. M. Block and C-I. Tan, PRL 98, 242001, (2007). Existing leading-order gluon distributions based on power-law description of individual parton distributions agree roughly with the new distributions for x>~10^-3 as they should, but are much larger for x<~10^-3.
Estimates are made of the ultra-high energy neutrino cross sections based on an extrapolation to very small Bjorken x of the logarithmic Froissart dependence in x shown previously to provide an excellent fit to the measured proton structure function
F_2^p(x,Q^2) over a broad range of the virtuality Q^2. Expressions are obtained for both the neutral current and the charged current cross sections. Comparison with an extrapolation based on perturbative QCD shows good agreement for energies where both fit data, but our rates are as much as a factor of 10 smaller for neutrino energies above 10^9 GeV, with important implications for experiments searching for extra-galactic neutrinos.
