Unitarized Diffractive Scattering in QCD and Application to Virtual Photon Total Cross Sections

The problem of restoring Froissart bound to the BFKL-Pomeron is studied in an extended leading-log approximation of QCD. We consider parton-parton scattering amplitude and show that the sum of all Feynman-diagram contributions can be written in an eikonal form. In this form dynamics is determined by the phase shift, and subleading-logs of all orders needed to restore the Froissart bound are automatically provided. The main technical difficulty is to find a way to extract these subleading contributions without having to compute each Feynman diagram beyond the leading order. We solve that problem by using nonabelian cut diagrams introduced elsewhere. They can be considered as colour filters used to isolate the multi-Reggeon contributions that supply these subleading-log terms. Illustration of the formalism is given for amplitudes and phase shifts up to three loops. For diffractive scattering, only phase shifts governed by one and two Reggeon exchanges are needed. They can be computed from the leading-log-Reggeon and the BFKL-Pomeron amplitudes. In applications, we argue that the dependence of the energy-growth exponent on virtuality $Q^2$ for $\gamma^*P$ total cross section observed at HERA can be interpreted as the first sign of a slowdown of energy growth towards satisfying the Froissart bound. An attempt to understand these exponents with the present formalism is discussed.Comment: 41 pages in revtex preprint format, with 10 figure

Saturation and Wilson Line Distributions

We introduce a Wilson line distribution function bar{W}_tau(v) to study gluon saturation at small Feynman x_F, or large tau=ln(1/x_F). This new distribution can be obtained from the distribution W_tau(alpha) of the Color Glass Condensate model and the JIMWLK renormalization group equation. bar{W}_tau(v) is physically more relevant, and mathematically simpler to deal with because of unitarity of the Wilson line v. A JIMWLK equation is derived for bar{W}_tau(v); its properties are studied. These properties are used to complete Mueller's derivation of the JIMWLK equation, though for bar{W}_tau(v) and not W_tau(alpha). They are used to derive a generalized Balitsky-Kovchegov equation for higher multipole amplitudes. They are also used to compute the unintegrated gluon distribution at x_F=0, yielding a completely flat spectrum in transverse momentum squared k^2, with a known height. This is similar but not identical to the mean field result at small k^2.Comment: One reference and two short comments added. To appear in Physical Revies

Nonlinear QCD Evolution: Saturation without Unitarization

We consider the perturbative description of saturation based on the nonlinear QCD evolution equation of Balitsky and Kovchegov (BK). Although the nonlinear corrections lead to saturation of the scattering amplitude locally in impact parameter space, we show that they do not unitarize the total cross section. The total cross section for the scattering of a strongly interacting probe on a hadronic target is found to grow exponentially with rapidity. The origin of this violation of unitarity is the presence of long range Coulomb fields away from the saturation region. The growth of these fields with rapidity is not tempered by the nonlinearity of the BK equation.Comment: 4 pages, RevTe

Production of gluons in the classical field model for heavy ion collisions

The initial stages of relativistic heavy ion collisions are studied numerically in the framework of a 2+1 dimensional classical Yang-Mills theory. We calculate the energy and number densities and momentum spectra of the produced gluons. The model is also applied to non central collisions. The numerical results are discussed in the light of RHIC measurements of energy and multiplicity and other theoretical calculations. Some problems of the present approach are pointed out.Comment: 9 pages, 11 figures, RevTeX; error in eq. (11) corrected, figures clarified, published in Phys. Rev.

Spatially explicit estimates of N2O emissions from croplands suggest climate mitigation opportunities from improved fertilizer management

With increasing nitrogen (N) application to croplands required to support growing food demand, mitigating N2O emissions from agricultural soils is a global challenge. National greenhouse gas emissions accounting typically estimates N2O emissions at the country scale by aggregating all crops, under the assumption that N2O emissions are linearly related to N application. However, field studies and meta-analyses indicate a nonlinear relationship, in which N2O emissions are relatively greater at higher N application rates. Here we apply a super-linear emissions response model to crop-specific, spatially-explicit synthetic N fertilizer and manure N inputs to provide subnational accounting of global N2O emissions from croplands. We estimate 0.66 Tg of N2O-N direct global emissions circa 2000, with 50% of emissions concentrated in 13% of harvested area. Compared to estimates from the IPCC Tier 1 linear model, our updated N2O emissions range from 20-40% lower throughout Sub-Saharan Africa and Eastern Europe, to >120% greater in some Western European countries. At low N application rates, the weak non-linear response of N2O emissions suggests that relatively large increases in N fertilizer application would generate relatively small increases in N2O emissions. Since aggregated fertilizer data generate underestimation bias in nonlinear models, high-resolution N application data are critical to support accurate N2O emissions estimates

Higher twists in the pion structure function

We calculate the QCD moments of the pion structure function using Drell-Yan data on the quark distributions in the pion and a phenomenological model for the resonance region. The extracted higher twist corrections are found to be larger than those for the nucleon, contributing around 50% of the lowest moment at Q^2=1 GeV^2.Comment: 8 pages, 3 figures, to appear in Phys. Rev.

Proton strangeness form factors in (4,1) clustering configurations

We reexamine a recent result within a nonrelativistic constituent quark model (NRCQM) which maintains that the uuds\bar s component in the proton has its uuds subsystem in P state, with its \bar s in S state (configuration I). When the result are corrected, contrary to the previous result, we find that all the empirical signs of the form factors data can be described by the lowest-lying uuds\bar s configuration with \bar s in P state that has its uuds subsystem in $S$ state (configuration II). Further, it is also found that the removal of the center-of-mass (CM) motion of the clusters will enhance the contributions of the transition current considerably. We also show that a reasonable description of the existing form factors data can be obtained with a very small probability P_{s\bar s}=0.025% for the uuds\bar s component. We further see that the agreement of our prediction with the data for G_A^s at low-q^2 region can be markedly improved by a small admixture of configuration I. It is also found that by not removing CM motion, P_{s\bar s} would be overestimated by about a factor of four in the case when transition dominates over direct currents. Then, we also study the consequence of a recent estimate reached from analyzing the existing data on quark distributions that P_{s\bar s} lies between 2.4-2.9% which would lead to a large size for the five-quark (5q) system, as well as a small bump in both G^s_E+\eta G^s_M and G^s_E in the region of q^2 =< 0.1 GeV^2.Comment: Prepared for The Fifth Asia-Pacific Conference on Few-Body Problems in Physics 2011 in Seoul, South Korea, 22-26 August 201

Nuclear effects in the Drell-Yan process at very high energies

We study Drell-Yan (DY) dilepton production in proton(deuterium)-nucleus and in nucleus-nucleus collisions within the light-cone color dipole formalism. This approach is especially suitable for predicting nuclear effects in the DY cross section for heavy ion collisions, as it provides the impact parameter dependence of nuclear shadowing and transverse momentum broadening, quantities that are not available from the standard parton model. For p(D)+A collisions we calculate nuclear shadowing and investigate nuclear modification of the DY transverse momentum distribution at RHIC and LHC for kinematics corresponding to coherence length much longer than the nuclear size. Calculations are performed separately for transversely and longitudinally polarized DY photons, and predictions are presented for the dilepton angular distribution. Furthermore, we calculate nuclear broadening of the mean transverse momentum squared of DY dileptons as function of the nuclear mass number and energy. We also predict nuclear effects for the cross section of the DY process in heavy ion collisions. We found a substantial nuclear shadowing for valence quarks, stronger than for the sea.Comment: 46 pages, 18 figures, title changed and some discussion added, accepted for publication in PR

Virtual photon fragmentation functions

We introduce operator definitions for virtual photon fragmentation functions, which are needed for reliable calculations of Drell-Yan transverse momentum ($Q_T$) distributions when $Q_T$ is much larger than the invariant mass $Q$. We derive the evolution equations for these fragmentation functions. We calculate the leading order evolution kernels for partons to fragment into a unpolarized as well as a polarized virtual photon. We find that fragmentation functions to a longitudinally polarized virtual photon are most important at small $z$, and the fragmentation functions to a transversely polarized virtual photon dominate the large $z$ region. We discuss the implications of this finding to the J/$\psi$ mesons' polarization at large transverse momentum.Comment: Latex, 19 pages including 6 figures. An error in the first version has been corrected, and references update

Neutral weak currents in pion electroproduction on the nucleon

Parity violating asymmetry in inclusive scattering of longitudinally polarized electrons by unpolarized protons with $\pi^0$ or $\pi^+$ meson production, is calculated as a function of the momentum transfer squared $Q^2$ and the total energy $W$ of the $\pi N$-system. This asymmetry, which is induced by the interference of the one-photon exchange amplitude with the parity-odd part of the $Z^0$-exchange amplitude, is calculated for the $\gamma^*(Z^*)+p\to N+\pi$ processes ($\gamma^*$ is a virtual photon and $Z^*$ a virtual Z-boson) considering the $\Delta$-contribution in the $s-$channel, the standard Born contributions and vector meson ($\rho$ and $\omega$) exchanges in the $t-$channel. Taking into account the known isotopic properties of the hadron electromagnetic and neutral currents, we show that the P-odd term is the sum of two contributions. The main term is model independent and it can be calculated exactly in terms of fundamental constants. It is found to be linear in $Q^2$. The second term is a relatively small correction which is determined by the isoscalar component of the electromagnetic current. Near threshold and in the $\Delta$-region, this isoscalar part is much smaller (in absolute value) than the isovector one: its contribution to the asymmetry depend on the polarization state (longitudinal or transverse) of the virtual photon.Comment: 30 pages 9 figure