quantum chromodynamics , the theory describing the interactions of quarks and gluons , has been shown to provide excellent descriptions of a wide variety of phenomena ranging from hadron spectroscopy via lattice field theory to hard scattering phenomena via perturbation theory .
yet qcd is a theory whose degrees of freedom are the quarks and gluons  quanta that can not be observed directly due to their confinement in hadrons . on the other hand , experiments deal with hadrons and leptons . in order to calculate observables measured in high energy scattering
processes a technique is needed that allows conversion from a description in terms of hadrons to one in terms of quarks and gluons ( collectively referred to as partons ) . for hadrons in the initial state
this conversion is provided by parton distribution functions ( pdfs ) , while for final state hadrons a similar role is played by parton fragmentation functions .
the pdfs , denoted by @xmath0 , for each type of quark and antiquark and for the gluon ( @xmath1 ) , allow one to essentially describe a beam of hadrons as an effective beam of quarks and gluons .
in the infinite momentum frame , the pdfs can be interpreted as probability densities describing how the parent hadron s momentum is shared amongst the different types of partons , as a function of the hadron s momentum fraction @xmath2 . with knowledge of the pdfs and the feynman rules for qcd
, one can calculate hard scattering cross sections and differential distributions that can then be compared to data .
pdfs are traditionally determined by simultaneously fitting a wide variety of data for large momentum transfer processes .
the parameters of the fits describe the pdfs at some initial momentum transfer scale , while evolution equations are then used to calculate the pdfs at all other scales needed for the calculations .
the data sets used for such fits often include deep inelastic scattering ( dis ) of charged leptons on proton and deuterium targets , or neutrinos on heavy nuclear targets , lepton pair production on proton and deuterium targets , and the production of photons , vector bosons , or jets at large values of transverse momentum .
while the pdfs provide us with a detailed description of the partonic substructure of hadrons , they contain only partial information .
partons ( both quarks and gluons ) have nonzero spin , so that the fundamental distributions in nature are the pdfs for a specific helicity ( spin projection along the direction of motion ) , @xmath3 and @xmath4 , corresponding to parton spins aligned and antialigned with that of the hadron , respectively .
experiments with unpolarized beams and targets are therefore sensitive only to the sums of the helicity pdfs , @xmath5 , while information on the differences , @xmath6 , can be obtained from measurements involving polarized beams and/or targets .
it is the purpose of this topical review to give a current `` snapshot '' of the status of our knowledge of pdfs , both unpolarized and polarized . to this end
, brief reviews of both the theory and the types of data used in the fits will be presented , and theoretical issues that need to be addressed in the perturbative calculations will be outlined . throughout this review
we shall attempt to answer the question `` what do we know and how do we know it ? ''
recently there have been several excellent reviews of pdfs @xcite , updating the early treatises of pdfs in refs .
@xcite , which have focused on various aspects of unpolarized pdfs , such as the impact on the phenomenology of the large hadron collider ( lhc ) .
the spin structure of the nucleon , including spin - dependent pdfs , was also reviewed in refs .
the present pedagogical review should be viewed as complementary to these efforts , emphasizing the different aspects of hadron structure that can be revealed through the study of unpolarized and polarized pdfs , at both small and large parton momentum fractions @xmath2 .
there is an old saying in experimental particle physics that `` today s discovery is tomorrow s calibration . ''
a similar effect occurs for theoretical calculations .
what was once touted as a test of qcd , now often plays the role of providing a description of backgrounds for searches for new types of particles . of course , one wants these calculations to be as accurate as possible and this , in turn , requires precise knowledge of the pdfs that enter into the calculation .
the error estimates on pdfs will be discussed at some length for this reason . on the other hand , there are many aspects of hadron structure that are currently not understood well , or at all , either theoretically or experimentally , and clearly identifying and delineating the limits of our knowledge in the context of global pdf fits serves a valuable purpose .
the plan of the review is as follows .
sec .  2 contains an overview of the necessary qcd theory required to understand the results of global fits .
sec .  3 contains the review of unpolarized pdfs , while sec .
4 summarizes the status of polarized pdfs .
finally , sec .
5 provides some speculations as to how our knowledge of pdfs may be improved through anticipated results of ongoing and future experiments .
as noted in the introduction , pdfs are necessary ingredients for obtaining predictions for hard scattering hadron
hadron and lepton
hadron processes . the cross section for a typical hadron  hadron process involving collisions of hadrons @xmath7 and @xmath8 , with momenta @xmath9 and @xmath10 , respectively , producing a state @xmath11 in addition to other hadrons ( collectively denoted by @xmath12 ) , can be written in the form @xmath13 where @xmath14 denotes an @xmath15-th order parton
parton `` cross section '' that produces the desired final state , and the functions @xmath16 and @xmath17 are the pdfs of flavor @xmath18 in hadron @xmath7 and flavor @xmath19 in hadron @xmath8 , respectively .
the state @xmath11 might denote a high-@xmath20 jet , a lepton pair , a photon or weak vector boson , for example , and @xmath21 is a scale that characterizes the hard scattering . for high-@xmath20 jet or photon production the scale @xmath22 .
the form ( [ eq : xsec ] ) can also be applied to the case when the hadrons @xmath7 and @xmath8 are polarized ; in the remainder of this section , however , we shall for illustration focus primarily on the unpolarized case
.    the parton - level cross sections will generally possess ultraviolet singularities that must be renormalized .
doing so leads to the introduction of the running coupling @xmath23 which depends on a renormalization scale @xmath24 .
there can also be infrared singularities associated with loop graphs ; these will cancel corresponding singularities from the emission of real soft gluons , provided that the observable is suitably defined , _
i.e. , _ that it is `` infrared safe . ''
finally , there will be collinear singularities associated with , in this case , the initial partons emitting additional partons at zero angle .
these collinear configurations correspond to internal propagators going on - shell and , as such , correspond to long - distance physics .
such singular terms can be _
factorized _ and absorbed into the pdfs .
this process introduces a factorization scale , @xmath25 , that separates the long - distance and short - distance hard scattering physics .    as shown in eq .
( [ eq : xsec ] ) , the partonic hard scattering cross section can be expanded in powers of the running coupling .
the dependence on @xmath24 that enters via the running coupling cancels that which appears in the partonic cross sections .
similarly , the dependence on @xmath25 in the pdfs cancels against the @xmath25 dependence in the partonic cross sections .
these cancellations represent the fact that the physical cross section , if calculated to all orders in perturbation theory , should not depend on the two scales introduced to control the ultraviolet and collinear singularities .
however , at any fixed order in perturbation theory these cancellations will only be approximate .
indeed , at order @xmath26 one has the relations @xmath27 these results suggest that the dependence on the renormalization and factorization scale choices will decrease as one carries out the perturbative calculations to higher order .
this is indeed the case for the processes typically used in global fits for pdfs .
( if a new kinematic configuration opens up at a higher order , then the corrections at that order can be large .
this happens when going from lo to nlo for heavy quark production , for example , in which case the scale dependence can actually increase at that order .
the reduction then occurs at the next order . )
these processes are all known to at least nlo and some to nnlo , as will be discussed below . in principle , the two scales @xmath24 and @xmath25 can be chosen independently . however , in processes where there is one large scale characterizing the hard scattering it is often the case that the choice @xmath28 is used .
this simplifying choice will be utilized in the results discussed in the sections to follow .
the dependence on the renormalization scale of the running coupling is governed by the qcd @xmath29 function via the renormalization group equation , @xmath30 the qcd @xmath29 function is calculable in perturbation theory and can be written as @xmath31 where the coefficients for @xmath32 flavors are @xmath33 the strong coupling @xmath34 calculated using the expression for @xmath29 in eq .
( [ eq : alphas ] ) with terms up to @xmath35 is referred to as the @xmath36loop running coupling .
lo calculations retain only @xmath37 , while nlo calculations retain @xmath38 , and nnlo calculations retain also @xmath39 . note that the @xmath29 function is negative , resulting in a decrease of @xmath34 with increasing scale , a phenomenon known as asymptotic freedom .
this behavior is essential for the applicability of perturbation theory for hard scattering processes .      in eq .
( [ eq : xsec ] ) both the pdfs and the partonic cross sections depend on the factorization scale @xmath25 .
as seen in eq .
( [ eq : scale ] ) , this scale dependence cancels up to the order of perturbation theory used .
this fact allows one to calculate the @xmath41 dependence of the pdfs .
indeed , the singularities associated with collinear parton emission are universal and factorize from the hard scattering subprocesses .
the factorization scale dependence of the parton distributions is then also universal , reflecting the independence of physical quantities on the scale @xmath25 .
this leads to the @xmath40 evolution , or renormalization group equations ( rge ) , @xmath42 where the bjorken @xmath2 variable is ( related to ) the longitudinal momentum fraction of the partons .
the evolution kernels @xmath43 or `` splitting functions '' stem from the collinear divergences absorbed in the distributions and represent the ( collinear ) resolution of a parton @xmath44 in a parton @xmath45 .
they are calculated perturbatively as a series expansion in @xmath46 , @xmath47 at present these have been computed up to 3 loops ( @xmath48 ) , which is the necessary accuracy for a nnlo analysis .
the evolution of spin - dependent pdfs @xmath49 follows in a similar manner , with a corresponding set of spin - dependent splitting functions @xmath50 . in this case
the splitting functions have been computed to 2 loops , enabling analyses to be performed at nlo accuracy .
traditionally , information on the pdfs of the nucleon has come from the process of deep - inelastic scattering ( dis ) of leptons from protons or nuclei , beginning with the pioneering experiments at slac in the late 1960s . in the one - boson exchange approximation
, the differential dis cross section can be written as a product of lepton and hadron tensors , @xmath51 where @xmath52 is the fine structure constant , @xmath53 is the laboratory solid angle of the scattered lepton , and @xmath54 is the incoming ( outgoing ) lepton energy . for neutral - currents ,
the summation is over @xmath55 and the interference @xmath56 , while for charged - currents only @xmath57 exchange contributes .
the leptonic tensor @xmath58 depends on the charge @xmath59 and helicity @xmath60 of the lepton , @xmath61 where @xmath62 and @xmath63 are the initial and final electron momenta , and @xmath64 and @xmath65 are the electron vector and axial - vector charges , respectively . the factors @xmath66 in eq .
( [ eq : sigtens ] ) denote the ratios of the corresponding propagators and couplings to the photon propagator and coupling squared @xcite @xmath67 where @xmath68 is the fermi weak interaction coupling constant , and @xmath69 is the @xmath70 boson mass .
the hadronic tensor @xmath71 contains all of the information about the structure of the hadron target .
using constraints from lorentz and gauge invariance , together with parity conservation , it can be decomposed into spin - independent and spin - dependent contributions , @xmath72 , \label{eq : wmunu}\end{aligned}\ ] ] where @xmath73 and @xmath74 are the nucleon and exchanged boson four - momenta , @xmath75 , and @xmath76 .
the nucleon polarization four - vector @xmath77 satisfies @xmath78 and @xmath79
.    for spin - averaged scattering , the nucleon structure is parametrized in terms of the vector @xmath80 and @xmath81 structure functions , and the vector - axial vector interference @xmath82 structure function , which requires weak currents .
these are generally functions of two variables ( such as @xmath2 and @xmath40 ) , but become functions of @xmath2 only in the bjorken limit , in which both @xmath40 and @xmath83 , but @xmath2 is fixed . in this limit
the @xmath80 and @xmath81 structure functions become proportional , according to the callan - gross relation , @xmath84 , and in the parton model are given in terms of quark @xmath85 and antiquark @xmath86 distribution functions , @xmath87,\ ] ] where @xmath88 is interpreted as the probability to find a quark of flavor @xmath85 and charge @xmath89 in the nucleon with light - cone momentum fraction @xmath2 . at finite energies , the logarithmic @xmath40 dependence from the evolution equations described in sec .
[ ssec : q2evol ] , as well as residual @xmath40 dependence associated with power corrections ( see sec .  [
ssec : power ] below ) , give corrections to the simple parton model expectations .
the spin - dependent structure functions @xmath90 and @xmath91 can be extracted from measurements where longitudinally polarized leptons are scattered from a target that is polarized either longitudinally or transversely relative to the electron beam . for longitudinal beam and target polarization , the difference between the cross sections for spins aligned and antialigned is dominated at high energy by the @xmath90 structure function .
the @xmath91 structure function can be determined with additional measurement of cross sections for a nucleon polarized in a direction transverse to the beam polarization . in practice one
often measures the polarization asymmetry @xmath92 , which is given as a ratio of spin - dependent and spin - averaged structure functions , @xmath93.\end{aligned}\ ] ] at small values of @xmath94 , one then has @xmath95 .
if the @xmath40 dependence of the polarized and unpolarized structure functions is similar , the polarization asymmetry @xmath92 will be weakly dependent on @xmath40 .    in analogy with the unpolarized @xmath80 structure function , in the parton modelthe structure function @xmath90
can be expressed at lo in terms of differences between quark distributions with spins aligned ( @xmath96 ) and antialigned ( @xmath97 ) relative to that of the nucleon , @xmath98 , @xmath99 .
\label{eq : g1quark}\ ] ] the @xmath91 structure function , on the other hand , does not have a simple parton model interpretation .
however , its measurement provides important information on the subleading , higher twist contributions which parametrize long - range nonperturbative parton - parton correlations in the nucleon ( see sec .
[ ssec : power ] ) .      as discussed in sec .
[ sec : qcd ] , the pdfs are universal , the collinear singularities associated with initial parton collinear emission having been absorbed into the pdfs .
therefore , the parton distributions appearing , for example , in the expressions for dis structure functions , are the same that describe the structure of the incoming hadrons in hadronic production . repeating the general form given in eq .
( [ eq : xsec ] ) ( with @xmath100 ) , one has , for spin - unpolarized scattering , @xmath101 and an analogous expression for the spin - dependent cross section , involving the difference of cross sections with hadrons @xmath7 and @xmath8 polarized in the same and opposite directions , @xmath102 here @xmath103 and @xmath104 are the spin - dependent pdfs for flavor @xmath18 in a hadron @xmath7 and flavor @xmath19 in a hadron @xmath8 , respectively , and @xmath105 is the corresponding spin - dependent partonic cross section . as indicated in eqs .
( [ eq : xsec2 ] ) and ( [ eq : xsec_pol ] ) , the partonic cross sections are calculable in fixed
order perturbation theory as series expansions in @xmath106 , which starts with different powers depending on the process .
it will be useful for subsequent sections to describe the parton kinematics for several examples of hadron
hadron processes .
although all processes used in global fits are known at least to nlo , the lo parton kinematics nevertheless serves as a useful guide to what region of parton momentum fraction and which combinations of parton flavors will be constrained by a given set of data .    consider high-@xmath20 dijet production in lowest order where two jets are produced with approximately balancing transverse momenta @xmath20 and rapidities denoted by @xmath107 .
two - body phase constrains the momentum fractions to be @xmath108 where the dimensionless ratio @xmath109 . for centrally produced dijets
one sees that the parton momentum fractions are both approximately @xmath110 .
however , if one or both of the jets is produced at far forward or backward rapidity values , then the @xmath2 range can be considerably expanded towards large or small values thereby allowing a `` tuning '' of the probed @xmath2 range .    for the case of inclusive production of a jet with transverse momentum @xmath20 and rapidity @xmath111
there is , at lo , an integration over one of the parton momentum fractions , say @xmath112 , given by @xmath113 again , one sees that for @xmath111 near 0 the ranges of the momentum fractions are centered near @xmath110 .
another important class of experiments involves the production of a system of mass @xmath114 , rapidity @xmath111 , and transverse momentum @xmath20 .
examples of such systems include @xmath115 , @xmath116 , and charged lepton pairs @xmath117 .
the four  momentum of the produced system is conveniently given as @xmath118 where @xmath119 , and @xmath120 is the azimuthal angle .
the rapidity @xmath111 and , equivalently , the feynman@xmath2 variable @xmath121 are defined by @xmath122 which lead to the relations @xcite @xmath123 where @xmath124 .    depending on the detected final state it is possible to define different cross sections for hadron
hadron collisions . in the drell - yan mechanism
@xcite a quark from one hadron and an antiquark from the other one annihilate into an intermediate vector boson ( @xmath125 , @xmath116 or @xmath57 ) which subsequently decays into a lepton pair . for lepton pairs of invariant mass @xmath126 ,
the process is dominated by virtual photon exchange .
the experiments usually consist of proton(beam)nucleon ( proton or some other nucleus ) collisions and the cross sections are extracted from the detection of muon pairs ( dimuon production ) produced in the decay of the virtual photons .
the double differential cross section for lepton pair production at lo is given by @xmath127 , \label{eq : dy}\end{aligned}\ ] ] where the two parton momentum fractions are given by @xmath128 from eq .
( [ eq : dy ] ) one can see that for @xmath129 collisions this process is sensitive separately to sea and valence distributions , in contrast to the neutral current dis cross sections , where they enter as @xmath130 .
inclusion of the drell  yan data in global qcd analyses of pdfs is in fact instrumental in fixing ratio @xmath131 for @xmath132 at lo ( see eq .
( [ eq : sig_pd_pp ] ) below ) .
another alternative is the use of charged current neutrino dis structure functions , although this process entails the use of model dependent corrections for the heavy nuclear targets typically employed in such measurements .
nuclear corrections can be avoided if one instead uses the time reversed processes such as @xmath133 .
a process closely related to drell  yan dilepton production is the hadronic production of electroweak bosons @xcite . for @xmath134 scattering , the inclusive production cross section for @xmath57 bosons , for example ,
is given at lo by @xmath135 where @xmath136 are the ckm matrix elements , and the sum runs over all light quark and antiquark flavors in both hadrons .
the production of @xmath137 bosons is then sensitive primarily to the products @xmath138 , while @xmath139 production is sensitive to @xmath140 . at large @xmath70 boson rapidity , or
equivalently @xmath132 , @xmath57 production is sensitive to the @xmath141 and @xmath142 quark pdfs in the proton , respectively . provided the antiquark distributions at @xmath143 are known to sufficiently accuracy , the @xmath144 cross section ratio would then be a clean probe of the @xmath145 pdf ratio , which is presently not well determined at large @xmath112 .
conversely , if the quark pdfs are known in particular regions of @xmath146 , then the @xmath147 ratios can also be used to probe the antiquark distributions .
similar considerations apply when one or both of the protons are polarized , as for the case of the spin program at rhic . for collisions of longitudinally polarized protons from unpolarized protons , @xmath148 , the generalization of eq .
( [ eq : wprod ] ) simply involves replacing one of the quark or antiquark distributions by the corresponding polarized quark ( antiquark ) pdf , @xmath149 .
single - spin asymmetries can then constrain ratios of polarized to unpolarized pdfs , @xmath150 , for various flavors , over specific regions of @xmath146 ( see sec .
[ ssec : pol_sea ] ) .
double - spin asymmetries , formed by combinations of inclusive hadron or jet production cross sections in the scattering of polarized protons from polarized protons , @xmath151 , have also been studied in order to constrain the polarized gluon distribution @xmath152 ( see sec .  [
ssec : pol_gluon ] ) .
the treatment of heavy quarks in high - energy processes generally requires more care than for light quarks .
the production mechanism of heavy quarks in dis is predominantly the photon
gluon fusion ( pgf ) process , which makes it particularly sensitive to the gluon distribution in the nucleon .
although in principle there could exist an `` intrinsic '' ( initial state ) heavy quark content of the nucleon , existing measurements are well described through `` extrinsic '' ( generated in the hard scattering ) heavy quark production only ( _ cf .
_ @xcite and references therein ) , so that any intrinsic heavy quark content of the nucleon is small and restricted to large values of bjorken @xmath2 ( see sec .  [
ssec : hqresults ] below ) .    the natural framework for describing the pgf mechanism
is known as the `` fixed flavor number scheme '' ( ffns ) . in this approach the light - quark flavors ( @xmath153 ) and
the gluon are considered as massless partons within the nucleon .
this scheme allows one to compute the contributions of heavy quarks ( @xmath154 ) perturbatively as final - state quantum fluctuations , taking into account the full dependence of the production cross section on the mass @xmath155 of the heavy quark .
as is common in the @xmath156 renormalization scheme , the evaluation of the strong running coupling @xmath106 is nevertheless based on the usual variable flavor number scheme for the @xmath29function governing its scale dependence .
this procedure automatically resums the contributions to the strong coupling from heavy quarks , and consequently improves the stability of the perturbative expansion @xcite .
complete calculations of ( heavy  quark ) dis structure functions are presently known at lo @xcite and at nlo @xcite . beyond this
, approximations based on threshold resummation have been made @xcite ( nnlo * in sec .
[ ssec : hqresults ] ) , but thus far the nnlo coefficient functions are only partially known @xcite , which constitutes a major drawback of any precision qcd analysis at nnlo accuracy .
another aspect of these calculations is the mass definition used ; while the original calculations of the amplitudes employ the pole mass definition , it has been reported @xcite that the use of the running mass definition improves the stability of the perturbative series .    in many situations calculations within the ( fully massive ) ffns scheme become unduly complicated , and in many cases
even impossible due to the unknown massive ( @xmath157 ) matrix elements at nnlo or even nlo .
for this reason it is common to generate parton distributions in the so - called `` variable flavor number scheme '' ( vfns ) , where the heavy quarks are also considered to be massless constituents of the nucleon . here , the required nlo and nnlo cross sections have been computed for a variety of important production processes , in particular at hadron colliders . in the vfns scheme the effective heavy quark distributions
are generated perturbatively from the nonperturbative distributions of light quarks and gluons using the boundary conditions of ref .
@xcite , typically at the unphysical `` thresholds '' @xmath158 . from there
one proceeds to their renormalization group evolution with an increased number of flavors @xmath159 and ( eventually ) 6 .    for situations where the ( threshold ) invariant mass of the produced system far exceeds the mass of the interacting heavy flavor in the ffns , the vfns predictions
have been shown @xcite to deviate from the ffns typically by about 10% , although the exact amount depends on the particular process and energy scale .
this is usually within the margins of the renormalization and factorization scale uncertainties and other theoretical ambiguities related to pdfs .
eventually , one nevertheless has to _ assume _ that these massless `` heavy '' quark distributions are relevant asymptotically , in that they can correctly describe the relevant cross sections at scales much larger than the heavy quark masses involved .
unfortunately this is not the case for dis structure functions , for which the strictly massless approach , known in this context as the zero  mass vfns , or zm - vfns , is well known to be experimentally inadequate , particularly near the heavy quark production thresholds . to remedy this the heavy - quark mass effects are re  inserted in what are known as general - mass vfnss ( gm - vfns ) , for which several prescriptions are available in the literature @xcite
( see also ref .
these schemes interpolate ( in a model - dependent way ) between the ffns results near production threshold and the asymptotic results of the zm - vfns .
however , since the interpolating schemes are based on the same massive matrix elements as the ffns , they do not have complete information on the @xmath160 heavy flavor wilson coefficients ( even though they are sometimes referred to as `` nnlo gm - vfnss '' ) .
the importance for global pdf analysis of finite heavy quark mass effects in the calculation of dis structure functions had been previously emphasized in refs .
@xcite , but was only universally recognized after the analysis of @xcite .
currently the abm @xcite and jr @xcite collaborations use the ffns for the calculation of dis structure functions in their global analyses , while the mstw @xcite , ct @xcite , herapdf @xcite and nnpdf @xcite groups use variants of gm - vfns .
the calculations of hadron collider cross sections are carried out in the vfns in all cases .
the elegant machinery that has been developed within the framework of perturbative qcd to analyze leading twist pdfs is , strictly speaking , valid only at high values of @xmath40 and @xmath70 where all hadron mass scales are suppressed , @xmath161 . in real experiments performed at a finite beam energy @xmath162 , however , the maximum values of @xmath40 and @xmath70 are limited , which inevitably restricts the available coverage in bjorken @xmath2 .
this is especially relevant at large @xmath2 , where in dis the invariant mass squared of the produced hadronic system is given by @xmath163 with @xmath164 the mass of the target nucleon . for fixed @xmath40 , as @xmath165 the final state hadron mass @xmath70 decreases as one descends into the region dominated by nucleon resonances at @xmath166  gev .
the resonance region may be treated using the concept of quark - hadron duality @xcite , although this goes beyond the scope of the usual pqcd analysis .    in the region of low @xmath40 ,
power corrections to the bjorken limit results that scale as powers of @xmath167 become increasingly important . in the operator product expansion ,
these are associated with higher twist corrections , which arise from multi - parton correlations and characterize the long - range nonperturbative interactions between quarks and gluons . of tremendous interest in their own right as providing glimpses into the dynamics of quark confinement
, the power corrections are viewed as troublesome backgrounds to efforts aimed solely at extracting leading twist pdfs .    to avoid the complications from the higher twist corrections ,
the usual strategy in global pdf analyses is to apply cuts specifying minimum values of @xmath40 and @xmath168 . in many analyses of unpolarized scattering data ,
the cuts are of the order @xmath169  gev@xmath170 and @xmath171  gev@xmath170 @xcite , which in practice restrict the range of @xmath2 that can be accessed to @xmath172 . for spin - dependent pdfs , the scarcity of high - energy data forces most global analyses to use less restrictive cuts , of the order
gev@xmath170 and @xmath174  gev@xmath170 .
on the other hand , there are a number of important reasons for needing to know the large-@xmath2 behavior of pdfs , and several recent unpolarized pdf analyses @xcite have been performed with relaxed cuts of @xmath175 and @xmath176  gev@xmath170 which have allowed an expanded reach into large @xmath2 .
the improvement in the large-@xmath2 kinematic coverage amounts to about 1300 more data points for proton and deuteron targets , representing some @xmath177 increase in the total number of dis data points compared with the more restrictive cuts .
most importantly , the resulting leading twist pdfs fits have proved to be very stable even when such low cuts have been applied @xcite .    with the inclusion of the kinematic regions in which subleading @xmath178 effects play a non - negligible role , it is obviously crucial to account for the power corrections that might otherwise obfuscate the leading twist pdfs . among the different categories of @xmath178 effects , the simplest are the target mass corrections ( tmcs ) , which are formally associated with matrix elements of leading twist operators @xcite and are of kinematical origin .
others include genuine higher twist corrections , which arise from dynamical , multi - parton correlations , as well as higher order perturbative qcd corrections , which can also resemble power suppressed contributions at low @xmath40 .
the standard method to compute tmcs in dis is based on the operator product expansion , and was first formulated by georgi and politzer @xcite , and expressions for all unpolarized and polarized structure functions now exist in both @xmath2 and mellin space @xcite . an alternative method based on collinear factorization ( cf ) in momentum space
was developed by ellis , furmanski and petronzio @xcite , and extended by various authors @xcite ( see also the recent reviews of tmcs in refs .
the advantage of the ope method is that tmcs can be calculated to all orders in @xmath178 in dis , whereas tmcs in the cf approach have only been computed to @xmath179 . on the other hand ,
the ope is limited to inclusive dis , while the cf framework can be applied to computing tmcs also in other processes @xcite . since
typically the non - dis data are taken at very high @xmath40 where tmc effects are very small , the ope method is usually adopted .
all of the tmc methods also suffer to some extent from the threshold problem , whereby the target mass corrected structure function remains nonzero as @xmath165 @xcite .
for the purposes of global fits , however , the region where the threshold effects become problematic is @xmath180  gev @xcite , which is mostly outside of where even the most liberal cuts in @xmath70 and @xmath40 are made @xcite . from the seminal work of de  rjula _ et al . _
@xcite , the appearance of the threshold problem in the analysis of tmcs is attributed to the neglect of dynamical higher twist corrections , both of which scale as powers in @xmath178 .
( in fact , higher order perturbative qcd corrections can also resemble power suppressed contributions at low @xmath40 . )
regardless of their origin , the various power suppressed corrections that are not included in a leading twist calculation can be absorbed into phenomenological functions ; for example , for an unpolarized structure function @xmath181 , @xmath182 where @xmath183 denotes the leading twist contribution including tmcs .
the higher twist corrections are sometimes assumed to be multiplicative , with the functions @xmath184 proportional to the leading twist contribution , @xmath185 .
possible additional @xmath40 dependence of the higher twist contributions , from radiative @xmath186 corrections , is usually neglected .
the leading twist pdfs were found in ref .
@xcite to be essentially independent of the tmc prescription adopted , with the ht parameters able to compensate for the variations due to the different tmc formulations .
the isospin dependence of the ht corrections was also studied in refs .
existing data do not allow for an accurate determination of the @xmath187 corrections @xmath188 , and attempts to include them in global fits produce anomalously small values of @xmath189 and @xmath190 in conjuction with large compensating twist4 and twist6 contributions @xcite , possibly due to overfitting and blurring of the scaling violations .    in the polarized case
, the greater scarcity of data means that typically one can not afford the luxury of @xmath40 and @xmath70 cuts as stringent as those applied in some of the unpolarized pdf analyses .
most global analyses of spin - dependent pdfs therefore include structure function measurements down to @xmath191  gev@xmath170 , where higher twist corrections are believed to be important .
the higher twist contributions to @xmath90 and @xmath91 may be treated in an analogous way to the unpolarized @xmath181 structure functions , with the important difference that the twist3 contributions to @xmath91 are not @xmath40-supressed .
phenomenological pdf analyses exist which include ht contributions to @xmath90 @xcite , @xmath91 @xcite , and to both functions simultaneously @xcite .
while the focus in the present work is on the leading twist pdfs , the higher twist contributions to @xmath90 and @xmath91 are also of intrinsic interest in themselves , containing information , for example , about the correlations of color electric and magnetic fields in the nucleon with the nucleon s spin @xcite .
since nucleons bound in a nucleus are not free , the parton distributions @xmath192 in a nucleus @xmath7 deviate from a simple sum of pdfs in the free proton and neutron , @xmath193 , where @xmath194 is the number of protons .
this phenomenon is especially relevant at small values of @xmath2 , where nuclear shadowing ( or screening ) effects suppress the nuclear to free isoscalar nucleon ( @xmath195 ) ratio , @xmath196 , and at large @xmath2 , where the effects of fermi motion , nuclear binding , and nucleon off - shellness give rise to the `` nuclear emc effect '' @xcite .
in addition , for spin - dependent pdfs , the different polarizations of the bound nucleons and nuclei need to be taken into account .    in the nuclear impulse approximation , where scattering is assumed to take place incoherently from partons inside individual nucleons , the pdf in a nucleus can be expressed as a convolution of the pdf in a bound nucleon and a momentum distribution function @xmath197 of nucleons in the nucleus @xcite .
coherent rescattering effects involving partons in two or more nucleons give rise to nuclear shadowing corrections to the impulse approximation . in general , at large @xmath40 the pdfs in a nucleus and in a nucleon are related by @xmath198 where the additive term @xmath199 represents nucleon off - shell or relativistic corrections , and @xmath200 parametrizes the shadowing corrections .
a similar expression can be written for spin - dependent pdfs .
the momentum distribution , or `` smearing function '' , @xmath197 , can be computed from nuclear wave functions , incorporating nuclear binding and fermi motion effects . at @xmath201 the smearing function has a simple probabilistic interpretation in terms of the light - cone momentum fraction of the nucleus carried by the struck nucleon , where @xmath202 and @xmath203 are the four - momenta of the nucleon and nucleus , respectively , and @xmath204 is the nuclear mass . in this case
the smearing function is normalized to unity , @xmath205 . at finite @xmath40
, however , the smearing function depends in addition on the parameter @xmath206 , where @xmath207 and @xmath208 are the energy and three - momentum transfer , respectively , which characterizes the deviation from the bjorken limit @xcite .
typically , the function @xmath197 is steeply peaked around @xmath209 , becoming broader with increasing mass number @xmath7 as the effects of binding and fermi motion become more important . in the limit of zero binding , @xmath210 , and
one recovers the free - nucleon case .
this is the usual assumption made in most global pdf analyses .
recently , several analyses @xcite have accounted for the nuclear effects by explicitly calculating the corrections from microscopic nuclear models , or attempted to constrain them phenomenologically @xcite .
the most straightforward calculation is for the simplest nucleus  the deuteron , for which both nonrelativistic and relativistic wave functions are available @xcite , constrained by high - precision nucleon - nucleon scattering data .
experiments with deuterium targets play a vital role , in fact , in determining the flavor decomposition of the proton pdfs .
traditionally , the standard method for disentangling the @xmath141 and @xmath142 pdfs has been through charged lepton dis , which for a proton target is sensitive at large @xmath2 to the combination @xmath211 .
dis from a neutron would constrain @xmath212 ; however , the absence of free neutron targets has necessitated the use of deuterium as effective neutron targets . while the nuclear corrections in the deuteron are typically a few percent , in some regions of kinematics , most notably at large @xmath2 , they can give rise to large uncertainties in the extracted @xmath142 quark pdf in particular .    an alternative ,
purely phenomenological approach @xcite has attempted to constrain the nuclear corrections in deuterium directly from the data .
however , without an independent , high - precision measurement of the @xmath142 quark pdf from processes other than inclusive dis , it is difficult to unambiguously separate the effects of the nuclear corrections from uncertainties in the @xmath142 quark pdf , especially at high values @xmath2 . until future experiments ( see ref .
@xcite for a discussion ) are able to determine the @xmath142 quark pdf independent of uncertainties in nuclear models , a more practical approach adopted by the cj collaboration @xcite has been to produce a set of global pdfs for a range of nuclear models , corresponding to mild ( cj12min ) , medium ( cj12mid ) , and strong ( cj12max ) nuclear corrections in the deuteron .
the current uncertainties in pdfs from nuclear corrections can have important consequences far beyond the fixed - target experiments where they are encountered most directly .
in fact , the effects of the nuclear smearing corrections are not suppressed at large @xmath40 , and must be considered at all scales wherever data at @xmath213 are used @xcite . through @xmath40 evolution ,
pdfs at large @xmath2 and small @xmath40 evolve to lower @xmath2 and higher @xmath40 , so that large-@xmath2 uncertainties in fixed - target experiments can have significant consequences for collider measurements @xcite , examples of which were discussed for selected observables at the tevatron and lhc in ref .
@xcite .    for neutrino scattering , in order to increase the relatively low rates and obtain sufficient statistics ,
experiments have often resorted to using heavier nuclear targets , such as iron or lead .
this is particularly relevant for determinations of the strange quark pdf , which is typically extracted from opposite - sign dimuon events in @xmath207 and @xmath214 charm production , @xmath215 or @xmath216 @xcite .
such extractions are complicated by the presence of nuclear corrections in neutrino structure functions @xcite , as well as effects of the nuclear medium on the charm quark propagation in the final state @xcite .
uncertainties in the strange quark pdf can have significant impact on @xmath70 and @xmath194 boson measurements at the lhc , for example , so that understanding of the nuclear effects will have impact far beyond lepton
nucleus dis @xcite .
in order to systematically study the nuclear dependence of pdfs , without assuming specific relations between pdfs in nuclei and nucleons , several groups have parametrized nuclear pdfs directly @xcite .
one approach has been to parametrize the @xmath7 dependence of the initial proton pdf parameters , and then perform a global fit of the available hard scattering data on nuclear targets .
this implicitly assumes a smooth @xmath7 dependence for the pdfs , which is reasonable for large @xmath7 , but may break down for light nuclei such as deuterium .
an alternative method is to fit to one type of data set and examine how the results differ from those on proton targets .
this can be useful for testing whether the nuclear corrections are consistent with those obtained from various data sets .
such comparisons have revealed , for instance , a controversy about whether or not the @xmath7 dependence of the neutrino - nucleus data from the nutev collaboration is compatible with that observed in charged lepton dis @xcite .
differences between neutrino and electromagnetic nuclear interactions can arise , for instance , from the presence of the parity - odd @xmath82 structure function in @xmath217 scattering , which does not contribute to charged - lepton scattering @xcite .    for spin - dependent scattering , the scarcity of data and their larger uncertainties at very small @xmath2 and at high @xmath2 , where nuclear corrections are most prominent , has meant that almost all global analyses have thus far relied exclusively on the effective polarization _ ansatz _ , in which the polarized pdf in the nucleus @xmath218 is related to the polarized pdfs in the proton and neutron as @xmath219 , with @xmath220 the average polarizations of the proton ( neutron ) in the nucleus . in practice ,
only polarized deuterium and @xmath221he nuclei have been used in dis experiments , in addition to protons .
the new global analysis of helicity pdfs by the jam collaboration @xcite is the first systematic attempt to incorporate the effects of nuclear smearing in dis structure functions of the deuteron and @xmath221he .
this will be important as future spin - dependent scattering experiments probe the nucleon spin structure at increasingly large values of  @xmath2  @xcite .
the concept of global fitting for pdfs relies on the formalism outlined in the previous sections .
partonic cross sections for various hard scattering processes are convoluted with scale - dependent pdfs to generate results for physical observables which can then be compared to data .
the pdfs are parametrized at some convenient input scale @xmath222 and then evolved using the appropriate evolution equations to the scales needed for each calculation .
the values of the input parameters are estimated using a @xmath189 minimization technique .
all of the global pdf fitting groups use some variation of this technique . a typical parametrization at the input scale @xmath222 for a generic polarized or unpolarized parton ( quark , antiquark or gluon ) distribution function @xmath223 is @xmath224 where @xmath225 represents some smoothly varying function usually chosen as a polynomial in @xmath2 or @xmath226 , although exponential functions are used as well .
in addition , there are parametrizations inspired by a statistical model of the nucleon @xcite , as well as those based on neural networks @xcite and self - organizing maps @xcite
.    some of the parameters in the input distributions can be determined from physical constraints .
for example , in the unpolarized case valence quark number is conserved , @xmath227 while the momentum sum rule requires @xmath228 = 1 , \label{eq : mom_sum_rule}\ ] ] where we use the notation @xmath229 , and the number of flavors at the input scale @xmath230 is usually taken to be @xmath231 ( see sec .
[ ssec : hq ] ) . in the polarized case
the first moments of the charge - conjugation even ( or @xmath11-even ) distributions can be related to octet baryon weak decay constants . for the isovector combination , corresponding to the bjorken sum rule , @xmath232 where @xmath233 is the nucleon axial charge , while for the su(3 ) octet one has @xmath234 where the octet axial charge @xmath235 is extracted from hyperon @xmath29-decays assuming su(3 ) flavor symmetry @xcite .
note that the sum rules ( [ eq : sum_rules])([eq : hyperon ] ) are preserved under @xmath40 evolution .
while most of the pdf groups use similar procedures and data , one can nevertheless obtain pdf results that can be rather different .
some of the reasons for this include the following :    * differences in the specific parametrizations and input scale @xmath230 . * differences in data selection , choices of data sets , kinematic cuts , and the specific treatment of the correlated errors of the data .
this typically limits the amount of data available at large values of @xmath2 .
* differences in the theoretical framework used , including the particular form of the solutions of the rge employed , the treatment of heavy quarks , inclusion of higher twist contributions , target mass effects , and nuclear corrections .    these differences should be born in mind when comparing the results from different pdf groups .
pdfs are an essential ingredient for producing predictions for processes in high energy experiments .
the uncertainties in these predictions depend , in part , on how well determined the pdfs are themselves .
it is important to bear in mind that the predictions that are compared to data are convolutions of pdfs with partonic hard scattering cross sections .
there are thus three main sources of pdf uncertainties : the fitted data , the partonic cross sections , and the parametrizations used to describe the pdfs .
the following list describes the main sources of uncertainty in the determination of pdfs .    *
the experimental errors on the fitted data can be directly propagated to the fitted pdfs .
standard techniques include the hessian @xcite , lagrange @xcite , and monte carlo @xcite methods . *
uncertainties due to the use of perturbation theory : these can be estimated to some extent by doing lo , nlo , and nnlo fits , although not all processes are known to nnlo accuracy .
* scale dependence : the perturbative predictions depend to some extent on the choices made for the renormalization and factorization scales for each process .
these choices will change the results for each process and the fitted pdfs must compensate these changes ( see sec .
iii of ref .
@xcite for an example of the effects of choosing different scales ) .
* choice of the value of the running coupling @xmath236 : some pdf determinations fit @xmath236 while others use the global average value .
this is mostly a question of philosophy . on the one hand , since the strong coupling is a parameter of qcd , then there is no freedom to choose a different value for each process .
thus , information from data types not included in the global fits ( such as data from @xmath237 processes ) provides valid constraints on the value of the coupling and are included in determining the global average . on the other hand
, fitting the value of the strong coupling and comparing it to the global average provides an interesting consistency check on the description of the data provided by qcd .
* choice of data sets and kinematic cuts : these choices can affect the fitted pdfs and the user should be aware of these differences . *
treatment of heavy quarks : the various schemes that are currently used differ at higher orders and such differences can affect the fitted pdfs . *
parametrization dependence : fitted pdfs can differ simply through the choice of the initial parametrizations .
extensive efforts are made to choose flexible parametrizations which are well - constrained by data , but there is no control over the behavior outside the kinematic region covered by the data .
this can lead to different extrapolations at very large or very small values of @xmath2 .
a method for estimating this remaining uncertainty has been suggested in ref .
@xcite .    of all these sources of error ,
the easiest to treat is the first  the propagation of the experimental errors .
the hessian method is described in detail in ref .
the elements of the hessian matrix are given by @xmath238 where @xmath239 denotes the @xmath240 pdf parameter .
the hessian matrix is generated during the actual minimization procedure and its inverse is the error matrix .
the eigenvectors of the error matrix can then be used to define eigenvector parameter sets which , in turn , can be used to calculate error bands for the pdfs or for specific processes .
one particular subtlety is that the error bands generally depend on a @xmath189 tolerance .
mathematically , one expects the @xmath241 parameter errors to correspond to an increase of @xmath189 by one unit from the minimum value .
however , it has been suggested @xcite that inconsistencies between different data sets may require a larger value to be used .
this `` @xmath189 tolerance '' varies between groups and allowance must be made for this when comparing the resulting error bands .
the lagrange multiplier @xcite method is useful when one wants to determine the pdf error on a specific observable @xmath12 such as the cross section for @xmath70 or @xmath194 production .
let @xmath242 denote the @xmath189 for the global data set .
then one minimizes the function @xmath243 where @xmath18 denotes the set of pdf parameters and minimization is done for a range of values of @xmath244 .
the end result is a relation between @xmath242 and the value @xmath12 of the chosen observable .
once one specifies the @xmath189 tolerance , _
i.e. _ , the range of allowable @xmath189 values , this maps out a range of values @xmath245 where @xmath246 is the value of @xmath12 at the global minimum .
an alternative to the usual linear propagation of errors ( hessian method ) which can be useful for minima that are not well behaved or defined , is the so  called monte carlo method . in order to propagate the experimental errors a number of ` replica '
data sets are generated by using random numbers and the original errors to generate new data points @xcite .
these replica data sets are then fitted and the resulting replica pdf sets are treated using standard statistics ; the central values are given by the average over replicas , and the uncertainties by the envelope of predictions .
the main motivation for pursuing global fitting as a technique for determining pdfs is that the use of a wide range of data types with different kinematic coverage places many constraints simultaneously on the pdfs .
each type of observable depends on a particular linear combination of pdfs or products of pdfs .
obtaining the best fit to all of the observables simultaneously has proven to be an efficient method for extracting pdfs .
nevertheless , it remains true that specific observables may be sensitive to a specific pdf or combination of pdfs . in this section some examples using simplified lo kinematics will be presented since these can be useful in understanding how pdf errors can be reduced by future measurements .      *
deep - inelastic scattering * experiments provide direct information on pdfs since both the structure functions and the cross sections depend linearly on the pdfs .
the following list summarizes the main dependences .    1 .
charged lepton neutral current measurements on a proton target constrain the combination @xmath247 where we have assumed that only single photon exchange contributes ( @xmath194 boson exchange involves a different linear combination of pdfs ) . at large values of @xmath2 the antiquark pdfs become negligible leaving the combination @xmath248 .
if one had a neutron target , then the linear combination would be @xmath249 which becomes @xmath250 at large values of @xmath2 .
of course , neutron targets are not available , so deuterium is often used for this purpose requiring that nuclear corrections be made in order to extract the neutron target information as discussed in sec .
[ ssec : nuclear ] .
2 .   charged current neutrino interactions constrains the combinations @xmath251 where the plus ( minus ) sign corresponds to @xmath252 ( @xmath253 ) .
the actual linear combinations depend on the type of target used ; high statistics neutrino experiments employ different types of heavy targets for which model - dependent nuclear corrections must be made .
note that one can , in principle , isolate specific combinations of @xmath85 or @xmath86 pdfs , depending on the type of target used .
a special case is provided by charm production in neutrino and antineutrino scattering as these are proportional to the @xmath254 and @xmath255 pdfs , respectively .
the experimental signal is the production of opposite sign muon pairs requiring one to take into account the charm fragmentation and decay as well as nuclear corrections , depending on the type of target employed .
the gluon pdf enters in dis at @xmath256 .
thus , it is mainly constrained through the @xmath40 dependence of the structure functions and rosenbluth ( or longitudinal / transverse ) separated @xmath257 data ( see sec .
[ ssec : gluon ] ) .
* vector boson production * ( lepton pairs , @xmath115 , and @xmath116 ) in lo proceed through @xmath258 annihilation .
a few key examples are :    1 .
lepton pair production in proton  proton and proton  neutron collisions depend on the combinations @xmath259 where @xmath260 , and the  @xmath261indicate contributions from the @xmath254 , @xmath262 and @xmath19 quarks .
of course , the @xmath263 cross section is obtained from data from a deuterium target .
these cross sections can be used to constrain the @xmath264 ratio , for example .
2 .   @xmath115 production constrains products of the form @xmath265 with specific weights given by the appropriate ckm matrix elements , whereas @xmath116 production constrains @xmath258 . for @xmath266 collisions at the tevatron at large values of rapidity ,
one has approximately , @xmath267 where @xmath268 denotes the proton ( antiproton ) , and the @xmath261 represents contributions from heavier quarks .
note that these results are written in terms of proton pdfs . if the rapidity is large and positive , then @xmath269 and one can neglect the antiquark terms so that these cross sections directly constrain the @xmath270 pdfs .
due to the missing neutrino from the @xmath70 decays , one can not directly reconstruct the rapidity distributions .
what is usually presented is the charged lepton rapidity asymmetry for @xmath115 production . in this case
the decay process means that the constraints on the pdfs are less direct , but such measurements still provide useful constraints on the @xmath145 ratio at moderate values of @xmath2 .
of particular interest is the model dependent determination of the @xmath70 rapidity asymmetry by the cdf collaboration @xcite , which directly constrains the @xmath145 ratio at large values of @xmath2 since the asymmetry is for the @xmath70 and not the charged lepton from the decay .
recently the atlas collaboration has presented an analysis @xcite using @xmath70 decay lepton and @xmath116 rapidity distribution data which shows potential for constraining the strange quark pdf .
this analysis makes use of the fact that the @xmath116 cross section receives a significant contribution from @xmath271 annihilation while the @xmath70 cross sections help constrain the @xmath272 and @xmath273 pdfs .
the other major class of observables includes inclusive * jet or photon , dijet * , and * photon + jet production*. each of these constrain products of pdfs summed over all flavors .
nevertheless , there are certain pdfs that can be constrained by these observables .
this follows since the @xmath141 and @xmath142 pdfs are well constrained by the dis data and the vector boson and lepton pair production data provide strong constraints on the @xmath272 and @xmath273 pdfs .
hence , these data have the greatest impact on the gluon distribution .    1 .
direct photon production in lo has two subprocesses : @xmath274 and @xmath275 .
hence , this process was though to be a good candidate for constraining the gluon , especially at fixed target energies where the available @xmath110 range extended out to about 0.6 .
however , photons can also be created by bremsstrahlung from the charged quarks
 often referred to as the fragmentation process .
indeed , one can define photon fragmentation functions to help describe this production component .
it turns out that at fixed target energies this component receives very large soft gluon corrections , requiring threshold resummation techniques @xcite .
this , combined with some apparent disagreements between experimental data sets @xcite meant that photon production has not fulfilled the original promise of constraining the gluon pdf . recently , however , an analysis @xcite suggests that the use of _ isolated _ photon collider data may well help constrain the gluon pdf .
the use of isolation cuts reduces the fragmentation contribution and the use of higher energy collider reduces the need for threshold resummation .
high-@xmath20 jet production has played a significant role in constraining the gluon pdf .
even though the quark and antiquark pdfs are well constrained by other types of data , there are still significant contributions from @xmath276 and @xmath277 subprocesses @xcite .
dijet production triple differential cross sections yield more information than single jet cross sections because the rapidity of the second jet is also constrained , thereby helping to constrain the momentum fractions of the pdfs . by tuning the rapidity ranges for the two jets one can explore different regions for these momentum fractions .
an example is given in ref .
photon + jet production offers similar constraints , but now the subprocesses are weighted by the squared charge of the parton to which the photon couples .      as with unpolarized measurements ,
historically most constraints on spin - dependent pdfs have come from * polarized charged - lepton dis * experiments .    1 .   in the one - photon exchange approximation , the difference of inclusive dis cross sections for leptons polarized longitudinally with spin parallel and antiparallel to the target hadron polarization measures @xmath11-even quark combinations @xmath278 . for proton targets ,
one has the combination @xmath279 while for the neutron the combination would be @xmath280 with contributions from heavy quark polarization expected to be negligible . in practice ,
polarized @xmath221he targets are usually used as effective sources of polarized neutron , since the neutron carries almost 90% of the spin of @xmath221he , while polarized deuterons , to which the proton and neutron spins contribute equally , are used to provide the isoscalar combination @xmath281 because of the greater sensitivity of the proton @xmath90 structure function measurements to the polarized @xmath141 quark pdf , and uncertainties associated with nuclear corrections when extracting the neutron from @xmath221he or deuterium data , the uncertainty on the @xmath282 pdf is significantly larger at high @xmath2 than on @xmath283 .
2 .   at nlo
, the polarized gluon distribution @xmath152 also enters in the @xmath90 structure function .
the @xmath40 evolution of the flavor singlet contribution to @xmath90 can then be used to constrain @xmath152 , as for the unpolarized gluon pdf discussed above , although the constraints are weaker in the polarized because @xmath284 at small @xmath2 .
* semi - inclusive dis * provides additional independent combinations of spin - dependent pdfs that can be used to reconstruct individual quark and antiquark flavors .    1 .
semi - inclusive production of hadrons @xmath285 in the current fragmentation region , primarily pions or kaons , is proportional to products of spin - dependent pdfs and quark @xmath286 hadron fragmentation functions , @xmath287 where @xmath288 is the fraction of the quark s energy carried by the hadron @xmath285 .
the fragmentation functions @xmath289 can be determined from other reactions , such as inclusive hadron production in @xmath290 annihilation .
one can weight particular quark or antiquark flavors by selecting favored ( such as @xmath291 or @xmath292 ) or unfavored ( @xmath293 or @xmath294 ) fragmentation functions for specific hadrons ( in this case @xmath295 ) .
information on the polarized strange quark pdf @xmath296 in particular can be obtained from data on @xmath297 production .
2 .   the polarized gluon distribution @xmath152 can be constrained by semi - inclusive dis data on charmed or high-@xmath20 hadron production through the photon
gluon fusion process .
recently the analysis of this process has been performed at nlo , extending earlier lo extractions of @xmath152 .
inclusive particle production in * polarized proton
proton collisions * provides an additional method of determining spin - dependent sea quark and gluon distributions .    1 .
@xmath57 production cross sections for scattering longitudinally polarized protons from unpolarized protons , @xmath148 , depend on products of spin - dependent and spin - averaged pdfs , for example , @xmath298 at large positive or negative rapidities , @xmath132 or @xmath299 , the cross sections ( or asymmetries ) are dominated by a single flavor , while at mid - rapidities both @xmath141 and @xmath142 flavors contribute .
inclusive jet or @xmath300 production in double - polarized proton
proton scattering , @xmath301 , is sensitive to the polarized gluon pdf . the first evidence for a small , but nonzero @xmath152 was recently observed by the star collaboration at rhic in jet data at @xmath302  gev .
direct photon production in polarized @xmath134 scattering has also been suggested @xcite as a means of probing @xmath152 .
using the technology outlined in sec .
[ sec : qcd ] , a number of global qcd analyses of the world s high - energy scattering data have produced sets of proton pdfs , up to next - to - leading order ( nlo ) or next - to - next - to - leading order ( nnlo ) accuracy .
the efforts have indeed been global , with groups in europe and the us in the forefront of the data analyses .
the pdfs sets include parametrizations from the mstw @xcite and abm @xcite groups , both of which use standard global fitting methodology ; the nnpdf @xcite collaboration , which uses a newer approach based on neural networks ; the jr @xcite pdfs , which are dynamically generated through @xmath40 evolution from a low @xmath40 input scale ; and the herapdf @xcite group , which includes only data from the h1 and zeus experiments at hera .
the us - based efforts have been centered around the cteq collaboration @xcite , which at present involves two derivative analyses of nucleon pdfs , by the ct @xcite and cj ( cteq - jefferson lab ) @xcite groups , as well as the ncteq @xcite analysis of nuclear pdfs .
the cj collaboration in particular has focused on developing the methodologies needed for describing data over a broad energy range including the low-@xmath40 and @xmath70 domain @xcite .    in this section
we summarize the results for unpolarized pdfs from the various pdf fitting groups , discussing their similarities and contrasting their differences .
we will focus mostly on the physics issues , rather than on technical aspects of pdf fitting .
of course with new data arriving or soon anticipated at the lhc , jefferson lab , and other facilities , the global fitting efforts are constantly evolving , so that the information presented here can only be viewed as a snapshot of the field at the present time .
( red ) , @xmath303 ( blue ) , @xmath304 ( green ) , @xmath305 ( violet ) , @xmath306 ( orange ) and @xmath307 ( black ) pdfs for the cj12mid fit @xcite at @xmath308  gev@xmath170 , shown on logarithmic _ ( left ) _ and linear _ ( right ) _ scales in @xmath2 .
note that in the left panel the gluon is scaled by a factor 1/10.,title="fig:",width=302 ]   ( red ) , @xmath303 ( blue ) , @xmath304 ( green ) , @xmath305 ( violet ) , @xmath306 ( orange ) and @xmath307 ( black ) pdfs for the cj12mid fit @xcite at @xmath308  gev@xmath170 , shown on logarithmic _ ( left ) _ and linear _ ( right ) _ scales in @xmath2 .
note that in the left panel the gluon is scaled by a factor 1/10.,title="fig:",width=272 ]    as a typical example of modern pdfs and their uncertainties , fig .
[ fig : allpdfs ] shows the @xmath309 , @xmath303 , @xmath304 , @xmath305 , @xmath306 and @xmath307 pdfs ( with the gluon scaled by a factor 1/10 ) for the cj12 fit @xcite ( for the case of moderate nuclear corrections , cj12mid ) at a scale of @xmath308  gev@xmath170 .
the general behavior of the pdfs is similar for all the parametrizations @xcite , particularly where sufficient data exist to constrain the distributions in specific regions of @xmath2 and @xmath40 . in regions where data are scarce , such as at very low @xmath2 ( @xmath310 ) or high @xmath2 ( @xmath311 ) for certain observables , or where data on specific pdfs ( such as the strange quark ) are subject to large experimental or theoretical uncertainties ,
the details of the pdfs can depend strongly on the assumptions adopted in the extrapolations .     and @xmath312 , and the @xmath305 and @xmath313 sea distributions , for the available 3-flavor nnlo pdf sets from jr09 @xcite ( red solid line ) , mstw08 @xcite ( blue dashed ) and abm11 @xcite ( green dotted ) , at @xmath308  gev@xmath170.,width=453 ]    , but for the singlet - sector distributions @xmath314 , @xmath315 and @xmath316 and the gluon @xmath307 pdf.,width=453 ]    some of these differences are evident in figs .  [ fig : nspdf ] and [ fig : spdf ] , which illustrate the flavor nonsinglet combinations @xmath317 , and the singlet - sector distributions @xmath318 , respectively , for the nnlo mstw08 @xcite , abm11 @xcite , nnpdf @xcite and jr09 @xcite pdf sets .
while the valence @xmath319 and @xmath320 pdfs are fairly well constrained in the intermediate-@xmath2 region , at large @xmath2 there are significant uncertainties on the @xmath320 distribution in particular .
the su(2 ) nonsinglet distribution @xmath321 is well determined by data from the drell - yan reaction over the range @xmath322 , but is not constrained at higher @xmath2 .
the strange nonsinglet distribution is the most difficult to determine , requiring a combination of neutrino and antineutrino scattering data , which are subject to large experimental and nuclear uncertainties .    for the @xmath11-even distributions in fig .
[ fig : spdf ] , the @xmath323 and @xmath324 are the dominant quark distributions at large @xmath2 , but the strange quark pdf becomes increasing more important at small @xmath2 .
numerically , the gluon distribution dominates all other pdfs at small @xmath2 , but has sizable uncertainties at large @xmath2 values . in the following sections we discuss features of each of the pdfs in more detail .
valence quarks give the global properties of the nucleon , such as its charge and baryon number .
knowledge of their momentum distributions is important for many reasons , especially at high values of @xmath2 where a single quark carries most of the nucleon s momentum .
the large-@xmath2 region is in fact a unique laboratory for studying the nonperturbative flavor and spin dynamics of quarks @xcite , as well as testing predictions from perturbative qcd for the behavior of pdfs in the limit @xmath165 @xcite .
reliable determination of pdfs at large @xmath2 is also important for searches for new physics beyond the standard model in collider experiments at the lhc . through perturbative qcd evolution , uncertainties in pdfs from fixed - target experiments at high @xmath2 and low @xmath40
can propagate to larger @xmath40 to affect cross sections at smaller @xmath2 values @xcite .
this is especially true for the forward production of particles of mass @xmath155 at large rapidities @xmath111 , whose cross sections are given by products of pdfs with one evaluated at small @xmath325 and the other at large @xmath326 .
the production of heavy @xmath327 and @xmath328 bosons , for example , is sensitive to @xmath142 quark pdf uncertainties at high rapidities , exceeding 100% in the @xmath329 channel , which places limits on the accuracy of cross section measurements for masses near the kinematic thresholds @xcite .
furthermore , understanding pdfs at large @xmath2 is vital for the analysis of neutrino oscillation experiments @xcite , where one of the most significant uncertainties comes from neutrino  nucleus cross sections at the interface of the dis and resonance regions .
the growing need to better understand large-@xmath2 pdfs and their uncertainties has been reflected in the greater attention being paid recently to the physics of the large-@xmath2 region , with dedicated experiments planned at jefferson lab following its 12  gev upgrade @xcite , as well as at proposed new facilities such as the lhec @xcite and the electron - ion collider @xcite .
it has also been a catalyst for the recent concerted theoretical efforts by the cj collaboration @xcite in their global qcd analyses of pdfs extending into the lower @xmath40 and @xmath168 regions , with the aim of providing better constraints on pdfs at large values of @xmath2 @xcite .     and @xmath142 quark distributions from the cj12mid pdfs @xcite , with pdf errors only ( solid , red shaded ) with combined pdf and nuclear uncertainties ( dashed , blue shaded ) , and with a fit excluding all deuterium data ( dot - dashed , green shaded ) , relative to the reference cj12mid set at an arbitrary scale @xmath330  gev@xmath170 .
note the different vertical scales for the @xmath141 and @xmath142 quark pdfs.,title="fig:",width=283 ]   and @xmath142 quark distributions from the cj12mid pdfs @xcite , with pdf errors only ( solid , red shaded ) with combined pdf and nuclear uncertainties ( dashed , blue shaded ) , and with a fit excluding all deuterium data ( dot - dashed , green shaded ) , relative to the reference cj12mid set at an arbitrary scale @xmath330  gev@xmath170 .
note the different vertical scales for the @xmath141 and @xmath142 quark pdfs.,title="fig:",width=283 ]    the current constraints on the @xmath141 and @xmath142 quark distributions are illustrated in fig .
[ fig : dunod ] for the cj12mid pdf set @xcite , where the uncertainties arising from nuclear corrections and pdf ( experimental ) errors are indicated separately .
since the @xmath141 quark pdf is relatively well constrained by the proton dis data at large @xmath2 , the effect of the nuclear uncertainties is minor , increasing the total uncertainty at intermediate @xmath2 by a few percent . for the @xmath142 quark distribution
, on the other hand , whose determination at present requires both proton and deuterium dis data , the nuclear correction uncertainties in the deuteron significantly increase the overall error for @xmath331 .
however , even though the use of deuterium data introduces the need to confront the problem of nuclear effects and their uncertainties , without these data the error bands on the @xmath142 quark pdf would be even larger over most of the @xmath2 range , as fig .
[ fig : dunod ] demonstrates .
in addition to accounting for the finite-@xmath40 ( tmc and higher twist ) and nuclear effects which become increasingly important as @xmath2 tends to 1 at fixed @xmath40 , the behavior of pdfs in the @xmath165 limit also depends critically on the functional form used for the @xmath2 dependence .
the conventional parametrizations utilized in most global pdf analyses assume a @xmath332 dependence for both the @xmath141 and @xmath142 quark pdfs , as in eq .
( [ eq : param ] ) above , so that in the limit as @xmath165 , the @xmath145 ratio tends either to zero or infinity .
@xcite allowed for a more flexible parametrization of the valence @xmath142 quark pdf , in which the @xmath320 distribution receives a small admixture from the @xmath319 pdf , @xmath333 with @xmath19 and @xmath262 as free parameters . provided the @xmath142 quark has a softer momentum dependence ( larger exponent @xmath334 ) than the @xmath141 quark , which agrees with phenomenology
, the ratio of the distributions will approach a constant value , @xmath335 , as @xmath165 .
a  nonzero , finite value in this limit is in fact expected in several different nonperturbative models of nucleon structure @xcite .
ratio for the _
cj12min ( dot - dashed ) , cj12mid ( solid ) and cj12max ( dashed ) pdfs @xcite , and _
( right ) _
cj12mid ( solid ) compared with the ct10 @xcite ( dot - dashed ) , mstw08 @xcite ( dashed ) and abkm09 @xcite ( dot - dot - dashed ) parametrizations at @xmath330  gev@xmath170.,title="fig:",width=283 ]   ratio for the _ ( left ) _ cj12min ( dot - dashed ) , cj12mid ( solid ) and cj12max ( dashed ) pdfs @xcite , and _ ( right ) _
cj12mid ( solid ) compared with the ct10 @xcite ( dot - dashed ) , mstw08 @xcite ( dashed ) and abkm09 @xcite ( dot - dot - dashed ) parametrizations at @xmath330  gev@xmath170.,title="fig:",width=283 ]    the ratio of the @xmath142 to @xmath141 quark distributions is illustrated in fig .
[ fig : dupdf ] for several modern nlo pdf sets , including the cj12 @xcite , ct10 @xcite , mstw08 @xcite ( dashed ) and abkm09 @xcite distributions .
as expected , the total uncertainty on the @xmath145 ratio at large @xmath2 increases when different nuclear correction models are considered , with the central values extrapolated to @xmath336 increasing from @xmath337 for the cj12min pdfs , with the minimum nuclear correction model , to @xmath338 for the cj12max set with the maximum nuclear correction model from ref .
@xcite . combining all uncertainties ,
the cj12 analysis found @xmath339 } \pm 0.10 \,\rm{[nucl]}$ ] as @xmath165 , where the first error is from the pdf fits and the second is from the nuclear correction models .
in contrast , using the more restrictive , conventional parametrizations , the @xmath145 ratios for the mstw08 and abkm09 distributions tend to zero in the @xmath165 limit , while that for ct10 tends to infinity .
this behavior distorts somewhat the error bands , making the mstw08 ratio appear to have an anomalously small uncertainty , while the ct10 result has a much larger one .
the relative errors , however , are similar for the mstw08 and ct10 pdfs , while the cj12 and abkm09 uncertainties are reduced at larger @xmath2 because of the additional high-@xmath2 data used in those analyses .
unfortunately , the extrapolated values of @xmath145 at @xmath336 span most of the range of available predictions @xcite , so that discrimination between the physical mechanisms that lead to the different @xmath165 behaviors will only be possible with constraints from new experiments that do not involve deuterium data and the nuclear model uncertainties @xcite .
inclusive electromagnetic dis probes @xmath11-even combinations of pdfs , @xmath340 , weighted by their squared charges , or in terms of valence and sea quark ( or antiquark ) distributions , @xmath341 . the shapes of the quark valence and sea distributions differ markedly as a function of @xmath2 .
as we have seen in the preceding section , the valence quark pdfs dominate at intermediate and large values of @xmath2 , but vanish as @xmath342 for @xmath343 .
the sea , in contrast , dominates at small @xmath2 , but is strongly suppressed as @xmath165 . for the light antiquarks @xmath272 and @xmath273 , structure function data with proton targets are sensitive at small @xmath2 to the combination @xmath344 , while deuterium targets ( nuclear corrections aside ) probe the isosinglet combination @xmath345 . under @xmath40 evolution ,
the singlet quark distribution mixes with the gluon , which leads to the rapid growth in the antiquark pdfs with decreasing @xmath2 .
with increasing @xmath40 , the uncertainties on @xmath345 also decrease .
inclusive dis measurements can therefore provide important constraints on the behavior of both the @xmath346 and @xmath86 distributions in the regions where their contributions are dominant .
however , in the intermediate-@xmath2 region , @xmath347 , where the magnitude of valence and sea quark pdfs is comparable , additional information is required to uniquely determine the distributions individually . here data from neutrino scattering can allow access to the parity - violating @xmath82 structure function , which measures @xmath11-odd combinations of pdfs , @xmath348
. combined with the @xmath11-even pdfs from the parity - conserving @xmath80 and @xmath81 structure functions , one can then uniquely determine the @xmath85 and @xmath86 distributions individually .
unfortunately , neutrino scattering cross sections are usually statistics limited compared with electromagnetic data , and with few exceptions , have necessitated the use of nuclear targets such as iron or lead as a means of increasing the rates .
this subsequently brings with it the inherent complications associated with accounting for nuclear corrections in relating the nuclear structure functions with those in a free nucleon @xcite .
more direct constraints on @xmath86 distributions can be obtained from lepton pair production in inclusive hadron
hadron scattering , which involves products of pdfs evaluated at different values of beam ( @xmath112 ) and target ( @xmath143 ) momentum fractions . as discussed in sec .
[ sssec : h - h ] , the drell - yan process or @xmath57-boson production in @xmath134 scattering at nonzero rapidity , for example , allows the antiquark distribution at small @xmath112 to be extracted with knowledge of the corresponding quark distribution at @xmath143 .
in particular , the drell - yan reaction @xcite has been used by the na51 collaboration at cern @xcite and the e866/nusea collaboration at fermilab @xcite to determine the ratio of @xmath273 to @xmath272 distributions in the proton .
naive expectations from perturbative qcd based on @xmath349 splitting suggest that since the masses of the @xmath141 and @xmath142 quarks are very similar , the production of @xmath350 and @xmath351 pairs should be nearly identical , even with the inclusion of higher - order radiative corrections @xcite .
the first indication of a significant asymmetry in the proton sea was obtained by the new muon collaboration ( nmc ) @xcite at cern , who performed an accurate measurement of the ratio of dis cross sections for hydrogen and deuterium , from which the gottfried sum @xmath352 was determined  @xcite , @xmath353 where charge symmetry is assumed ( see sec .  [ ssec : csv ] below ) . for a flavor symmetric su(2 )
proton sea , the prediction for the gottfried sum would be @xmath354 . in contrast , the nmc measurement found @xmath355 @xcite , indicating a strong violation of flavor symmetry in the light antiquark sea , @xmath356 following the suggestion by ellis and stirling @xcite that the light antiquark sea could be more directly probed in the drell - yan process , the na51 experiment used the ratio of dimuon cross sections for @xmath357 and @xmath134 scattering to extract the ratio @xmath358 at an average @xmath359 at a feynman-@xmath2 of @xmath360 .
subsequently , the e866/nusea experiment at fermilab measured the @xmath361 ratio at large @xmath121 , where at leading order @xmath362,\ \ \ \ \
x_a \gg x_b .
\label{eq : sig_pd_pp}\end{aligned}\ ] ] the results , illustrated in fig .
[ fig : dbar_ubar ] , confirmed a significant deviation of the cross section ratio from the naive perturbative qcd expectation , indicating a large asymmetry between @xmath273 and @xmath272 .     and @xmath134 scattering as a function of the fractional momentum @xmath143 of the target parton .
the data points from the fermilab e866/nusea collaboration @xcite are compared with the cross section ratio calculated using the jr09 @xcite ( red solid line ) , mstw08 @xcite ( blue dashed ) , and abm11 @xcite ( green dotted ) , and nnpdf23 @xcite ( yellow dot - dashed ) pdf parametrizations.,width=453 ]    the earliest efforts to explain an enhancement of @xmath273 quarks over @xmath272 in the proton invoked the effects of pauli blocking , which arise from the presence of different numbers of valence @xmath141 and @xmath142 quarks in the proton @xcite .
estimates of this effect are very difficult to quantify reliably @xcite , with some calculations @xcite even suggesting that it produces a small excess of @xmath272 over @xmath273 .
alternative explanations have focused on the role of chiral symmetry breaking and the pion cloud of the nucleon @xcite . because of the heisenberg uncertainty principle , a component of the proton s wave function involves its dissociation into a nucleon and a virtual pion state .
the preferential coupling of a proton to a @xmath363 , with the @xmath364 containing a valence @xmath365 pair , implies a natural mechanism for generating an excess of @xmath273 over @xmath272 . while many models have been proposed to compute this effect quantitatively ( see refs .
@xcite for reviews ) , the underlying physics principles are model independent and simply the chiral symmetry properties of qcd @xcite .
while most of the modern pdf parametrizations are able to accommodate an enhanced @xmath273 sea over the range @xmath366 , there is some uncertainty in the trend of the @xmath264 ratio at larger @xmath2 .
the new drell - yan experiment e906/seaquest at fermilab @xcite will extend the @xmath2 coverage up to @xmath367 , and a proposal at the j - parc facility in japan @xcite would extend the range even further .
most pdf analyses have been performed under the assumption of charge symmetry , or independence of interactions under rotations in isospin space .
charge symmetry implies that the @xmath141 quark distribution in the proton is identical to the @xmath142 quark distribution in the neutron , and _
vice versa _ , and is believed to be accurate in nature at the @xmath368 level .
this symmetry is broken explicitly by light quark mass differences , @xmath369 , as well as by electromagnetic corrections .
such effects lead to nonzero values of the `` majority '' and `` minority '' charge symmetry violating asymmetries in the nucleon , defined as @xmath370 respectively , and similarly for the antiquarks @xmath371 .
although there is currently no direct evidence for charge symmetry violation ( csv ) in pdfs , it has been predicted in several nonperturbative models for both valence and sea quark pdfs @xcite .
in addition , electromagnetic radiative effects can be included in the @xmath40 evolution equations ( sec .
[ ssec : q2evol ] ) , in analogy with usual gluon radiation , by adding a term that accounts for the emission of photons by quarks @xcite , @xmath372 here @xmath373 are the new qed splitting functions ( expanded in terms of the electromagnetic coupling @xmath52 ) , and the sum over partons includes also the photon pdf @xmath374 . even with charge symmetric pdfs at the input scale , the different electric couplings of the photon to @xmath141 and @xmath142 quarks in the modified evolution equations will induce nonzero csv distributions ( [ eq : csv_u ] ) and ( [ eq : csv_d ] ) at a higher scale .
the valence quark asymmetries resulting from the qed - modified @xmath40 evolution ( [ eq : qeddglap ] ) are illustrated in fig .
[ fig : csv_ssbar ] at a scale @xmath375 gev@xmath170 .
the perturbatively generated csv distributions from glck _ et al .
_ @xcite are calculated with an average isoscalar quark mass @xmath376  mev , with the qed radiative corrections computed after fixing the qcd evolution effects .
the majority valence @xmath377 distribution is negative at intermediate @xmath2 values , and is of opposite sign to the minority valence distribution @xmath378 , which is smaller in magnitude .
qualitatively , these features are similar to the nonperturbative csv valence asymmetries computed in the bag model @xcite , which arise from differences between the spectator @xmath379 and @xmath380 diquark masses , as well as the proton  neutron mass difference .
the mrstqed @xcite analysis of csv was similar to that in ref .
@xcite , but involved unequal photon pdfs in the proton and neutron resulting from qed radiation with quark masses @xmath381  mev . since the photon now also contributes to the momentum sum rule , eq .
( [ eq : mom_sum_rule ] ) , the csv photon distributions induce charge symmetry violation in the valence quark pdfs , assuming the sea and gluon are charge symmetric . despite the different assumptions , the resulting distributions turn out to be similar to those in fig .
[ fig : csv_ssbar ] .     and `` minority '' @xmath378 valence quark distributions for the radiative qed fit of ref .
@xcite at @xmath375  gev@xmath170 , compared with the bag model estimates have been taken from ref .
_ ( right ) _ strange  antistrange asymmetry at @xmath382  gev@xmath170 ( appropriate for the nutev experiment ) for the asymmetric gjr08 @xcite ( red solid ) and the mstw2008 @xcite ( blue dashed ) pdfs.,title="fig:",width=253 ]    and `` minority '' @xmath378 valence quark distributions for the radiative qed fit of ref .
@xcite at @xmath375  gev@xmath170 , compared with the bag model estimates have been taken from ref .
_ ( right ) _ strange  antistrange asymmetry at @xmath382  gev@xmath170 ( appropriate for the nutev experiment ) for the asymmetric gjr08 @xcite ( red solid ) and the mstw2008 @xcite ( blue dashed ) pdfs.,title="fig:",width=364 ]    while the overall magnitude of the possible csv effects is rather small , they can nevertheless play an important role in understanding the discrepancy between the nutev measurement of @xmath383 @xcite and the standard model value @xcite .
the nutev collaboration extracted the paschos - wolfenstein ( pw ) ratio @xcite of total neutral to charged current cross sections , involving the difference between @xmath207 and @xmath384 cross sections .
in addition to assuming small uncertainties on the nuclear corrections necessary for translating the measured iron target data to an isoscalar nucleon target , the analysis also neglected any csv effects .
including csv corrections , the experimental pw ratio would receive a contribution proportional to @xmath385 @xcite . depending on the sign of the csv distributions ,
the asymmetry can therefore either remove some of the discrepancy with the nutev result , or exacerbate it .
the signs and magnitude of the predicted perturbative qed and nonperturbative quark model asymmetries in fig .
[ fig : csv_ssbar ] would in fact remove approximately 2/3 of the discrepancy .
an additional correction to the pw ratio can arise from differences between strange and antistrange distributions in the nucleon , @xmath386 , which we discuss in the next section .
traditionally the strange quark distribution has been more difficult to determine experimentally than the non - strange sea . in practice one
has often assumed that the @xmath254 and @xmath255 pdfs ( which are often assumed to be identical ) are proportional to the @xmath345 distributions , @xmath387 with @xmath388 ranging between @xmath389 and 0.5 .
the @xmath254 and @xmath255 distributions can be directly determined , however , through the measurement of dimuon pairs in neutrino dis , which are produced through the charged current by scattering from strange quarks with subsequent generation of charm .
in particular , while neutrino scattering produces charmed hadrons ( @xmath390 ) , antineutrino scattering results in anticharmed hadrons ( @xmath391 ) .
the semileptonic decay of the charmed hadrons ( such as @xmath392 and @xmath393 mesons ) then yields @xmath394 pairs along with the associated hadronic debris . because of
the low rates involved in neutrino scattering experiments , in order to increase the statistics these have typically made use of heavy nuclear targets , such as iron . since the nuclear effects for neutrino and charged lepton scattering are in general different @xcite , the nuclear corrections introduce an additional source of uncertainty in extractions of the strange quark pdf .
a number of nonperturbative models have in fact predicted an asymmetry between the @xmath254 and @xmath255 distributions in the nucleon , which can be probed by taking the difference between the @xmath207 and @xmath384 cross sections .
for example , in analogy with the pion cloud models invoked to explain the @xmath273 excess over @xmath272 in the proton in sec .
[ ssec : lightsea ] , the dissociation of a nucleon to a hyperon and a virtual kaon naturally produces unequal @xmath254 and @xmath255 distributions , although the magnitude and even the sign is difficult to determine unambiguously @xcite .
in addition , an asymmetry can also be generated through perturbative qcd evolution at nnlo because of unequal valence @xmath141 and @xmath142 distributions @xcite .
the ccfr @xcite and nutev @xcite collaborations at fermilab have collected the highest statistics data on inclusive charm production in @xmath207 and @xmath384 scattering on an iron target , enabling the most direct constraints to be placed on the strange quark pdfs as a function of @xmath2 .
the latter in particular was the first analysis using a complete nlo description of charm production .
the cross section for ( anti)neutrino charm production was computed to nlo in refs .
@xcite within the ffns , in which , apart from the gluon , only the light quark flavors ( @xmath395 ) are included as massless partons ( see secs .  [
ssec : hq ] and [ ssec : hqresults ] ) . the fully differential nlo calculation ,
including acceptance corrections , was performed in ref .
( note , however , that the use of the expressions in @xcite together with vfns distributions , as implemented in the nutev analysis @xcite , is strictly not consistent . )
comparison of the @xmath207 and @xmath384 events found a preferred fit for the @xmath396 distribution that peaked at @xmath397 , with the compensating negative contribution required by the sum rule @xmath398 to be restricted to the unmeasured region at @xmath399 @xcite . for the first moment , the value obtained at @xmath382  gev@xmath170 was @xmath400 , where the errors are statistical , systematic , and external .
the latter includes a dominant contribution from the charm semileptonic branching ratio uncertainty , as well as from nuclear corrections @xcite .    the strangeness asymmetry in the nucleon
is shown in fig .
[ fig : csv_ssbar ] at @xmath382  gev@xmath170 , as relevant for the nutev experiment , for several different pdf parametrizations .
the mstw08 @xcite strange asymmetry is fitted to the ccfr and nutev data assuming the form @xmath401 where the factor @xmath402 ensures zero net strangeness in the nucleon .
the strange asymmetry in ref .
@xcite is generated radiatively from the boundary condition @xmath403 at @xmath404  gev@xmath170 , but allowing @xmath405 to have the form in eq .
( [ eq : sminus ] ) .
the resulting @xmath405 asymmetry at the experimental scale is similar to the mstw08 fit , and to the phenomenological extraction in ref .
@xcite , within the currently sizable uncertainties .
the values for the momentum weighted asymmetry are the order @xmath406 for the mswt08 fit @xcite and @xmath407 for the dynamically generated pdfs @xcite at scales @xmath408  gev@xmath170 .
the cteq collaboration found a strange asymmetry in the range @xmath409 @xcite .
these values have the correct sign and similar magnitude to that required to reconcile the nutev @xmath383 measurement with the standard model value @xcite , once the charge symmetry violating effects in sec .
[ ssec : csv ] are also included @xcite .    finally , a recent measurement by the atlas collaboration @xcite of inclusive @xmath57 and @xmath194 boson production in @xmath134 collisions at the lhc found a significantly larger value of the strange to non - strange sea ratio in eq .
( [ eq : kappa ] ) , @xmath410 at @xmath411 for @xmath412  gev@xmath170 , than in previous analyses
. a reanalysis of the atlas and hera data by the nnpdf group @xcite found , however , that the uncertainty on the extracted @xmath388 was significantly underestimated , and that a more complete global analysis gives @xmath388 values consistent with the traditional values .
heavy quark production contributes considerably to inclusive dis ; for example , the charm ( bottom ) contributions to the dominant @xmath81 structure function constitute up to 30% ( 3% ) of the total in the small@xmath2 region . besides these contributions , heavy quark electroproduction is directly accessed experimentally in semi  inclusive dis , typically by detecting charm ( bottom ) mesons in the final state . because the main production mechanism is photon
gluon fusion , data on this process provide valuable constraints on the gluon pdf in the nucleon .
furthermore , they provide an appropriate context for testing the different approaches to these calculations that have been adopted in various pdf analyses ( see sec .  [
ssec : hq ] ) .         since it is by far the most relevant , we focus the discussion on the charm flavor , although similar considerations often also apply to bottom ( top contributions to this process , on the other hand , are negligible at accessible kinematics ) .
the most relevant data on heavy quark electroproduction are the ( h1 + zeus ) combination of hera measurements @xcite . a comparison of these data with ffns calculations at lo , nlo and nnlo * ( see sec .  [
ssec : hq ] ) in fig .
[ fig : charm ]  ( left ) shows excellent agreement between the ffns calculations .
in fact , the experimental analysis of the hera data is based primarily on the ffns fully differential ( exclusive ) calculations of ref .  @xcite . despite this fact ,
the constructs in the gm - vfnss are often motivated by arguing that pseudo - mass - divergences ( non - collinear logarithms of the form @xmath413 ) which appear in the wilson coefficients of the massive calculations need to be resummed . however , the ffns results are notably stable even for very large values of @xmath414 @xcite , in particular at the largest @xmath40 values accessible at hera , suggesting little need to resum these supposedly `` large logarithms '' in the context of global pdf analyses .    finally , in addition to the perturbative generation of charm , it is also possible for charm to be produced through nonperturbative mechanisms , such as those discussed in refs .
@xcite ( and references therein ) . first postulated to account for some early charm production data in hadronic reactions , these `` intrinsic '' charm contributions share several characteristic features , even though their details depend somewhat on the models .
in particular , they are typically valence - like , with significantly harder @xmath2 distributions than those generated perturbatively , and can produce large asymmetries between the @xmath262 and @xmath415 distributions .
measurements of the charm structure function @xmath416 by emc @xcite , especially the two data points at the highest @xmath40 in fig .
[ fig : charm ]  ( right ) , have been interpreted as indicating a nonzero intrinsic charm contribution at large @xmath2 .
the evidence is not conclusive , however , in view of the disagreement evident in fig .
[ fig : charm ]  ( right ) between the emc data and theoretical predictions at smaller @xmath2 , where the predictions are in agreement with the more recent hera results in fig .
[ fig : charm ]  ( left ) .    the current limits on the normalization of the nonperturbative charm quark distribution range from @xmath417 to 2% , depending on which data sets
are fitted and which models of intrinsic charm considered @xcite . although current global pdf analyses do not require an intrinsic charm component , the existence of nonperturbative charm remains an intriguing issue which can be addressed by high - precision measurements of charm production at future facilities such as j - parc @xcite , gsi - fair @xcite , jefferson lab at 12  gev @xcite , after@cern @xcite , or an electron - ion collider @xcite .
in addition to the heavy quark production discussed above , there are three primary sources of information on the gluon distribution : the @xmath40 dependence of the dis structure function @xmath81 at low values of @xmath2 ; the longitudinal structure function @xmath257 at all @xmath2 values ; and jet ( both at hera and at hadron colliders ) and photon production cross sections at moderate to high values of @xmath2 .
furthermore , the momentum sum rule ( [ eq : mom_sum_rule ] ) provides an additional constraint on the gluon pdf . since the normalizations of the valence pdfs are fixed by the flavor sum rules in eq .
( [ eq : sum_rules ] ) and the contributions of the light antiquarks are relatively well determined by dis and lepton pair production data , the momentum sum rule effectively constrains the gluon contribution .
the effect is an anticorrelation between the low- and high-@xmath2 behavior of the gluon pdf : a decrease ( increase ) of the gluon distribution at low @xmath2 leads to an increase ( decrease ) at high values of @xmath2 .
interestingly , this can result in more constrained gluon distributions in the medium-@xmath2 region ( see fig .
[ fig : gluon ] ) .
gev@xmath170.,width=377 ]    to understand the relationship between the scale dependence of @xmath81 and the gluon pdf , consider eq .
( [ eq : rge ] ) for the scale dependence of the quark pdf .
multiply both sides by @xmath2 and by the squared parton charge and sum over the active flavors . in this way one derives ( using the lo expression for @xmath81 ) , @xmath418 \label{eq : f2q2}\ ] ] where @xmath419 and @xmath420 are the quark  quark and quark
gluon splitting functions . at small values of @xmath2
the gluon term dominates the integrand so that the scale dependence of @xmath81 places a constraint on the gluon pdf .
a second constraint on the gluon from dis comes from the longitudinal structure function @xmath257 , which is defined as @xmath421 in lowest order , where only quarks contribute to @xmath80 and @xmath81 , one has @xmath422 so that @xmath423 .
the first nonzero contribution to @xmath257 starts at @xmath256 where , at small @xmath2 , there is a significant contribution from the subprocess @xmath424 , which again imposes strong constraints on the gluon pdf at small @xmath2 . at larger @xmath2 values
the valence quark contributions increase , the gluon pdf decreases , and the subprocess @xmath425 dominates .
however , since the same combinations of quark distributions are well determined from @xmath81 , where they enter with one power of @xmath34 less and the gluon does not contribute , precision measurements of the longitudinal structure function provide valuable information on the gluon pdf at large @xmath2 @xcite .
note also that the gluon pdf is accompanied by the same powers of @xmath34 in all dis structure functions , starting with @xmath426 .
this results in a correlation between the value of @xmath34 obtained in global pdf analysis and the size and shape of the gluon distribution , with larger @xmath34 leading to a smaller gluon pdf in the small-@xmath2 region and ( via the momentum sum rule correlation ) a larger gluon pdf at large @xmath2 .
this can be seen in fig .
[ fig : gluon ] , where the jr09 @xcite and abm11 @xcite analyses obtain @xmath34 values and small-@xmath2 gluon pdfs of similar size , while the mstw08 fit @xcite finds a larger @xmath34 value and a smaller gluon pdf in this region ( which even turns negative ) .
both possibilities describe well the inclusive cross sections measurements at hera .
this ambiguity can be reduced by using processes with a different leading @xmath34 power , typically jet production at the tevatron , which , as noted in sec .
[ sssec : unpol ] , can provide constraints on the gluon pdf . to understand this , recall that the calculation of the jet cross section involves a sum over @xmath427 , @xmath276 and @xmath277 induced subprocesses where `` @xmath85 '' stands for any flavor of quark or antiquark . at low values of jet @xmath110 ,
the @xmath277 subprocesses dominate , at intermediate values all three contribute , and eventually as @xmath110 approaches one only the @xmath427 terms survive  see fig .  3 in ref .
@xcite , for example .
jet data at the tevatron extend to @xmath110 values of approximately 0.5 , so these data provide some constraints on the gluon pdf , at least in the mid-@xmath2 range where the @xmath276 subprocesses make substantial contributions .    as discussed in ref .
@xcite , isolated photon production at collider energies has the potential to tighten the error bands in the @xmath2 range below about 0.2 . it remains a challenge to find processes that will constrain the gluon beyond @xmath428 .
the fundamental problem is that the valence pdfs here are larger than the gluon pdf .
most hard processes that receive contributions from @xmath276 subprocesses also receive contributions from @xmath427 processes , which reduce the sensitivity to the gluon pdf .
one counterexample is provided by direct photon production in @xmath134 collisions .
the two dominant subprocesses are initiated by @xmath258 and @xmath276 initial states . at large values of @xmath110 the fragmentation process
is suppressed ; it can be further reduced by photon isolation cuts .
the single photon cross section involves an integration over the awayside jet rapidity which tends to smear out the @xmath2 values of the contributing pdfs .
there are regions where one @xmath2 is small and one is large  one might think that this could constrain the large-@xmath2 gluon
. however , the favored configuration will be where the gluon or @xmath86 is at small values of @xmath2 while the valence quark is at the larger values .
instead , consider the photon + jet cross section with equal and opposite photon and jet rapidities , in which case both momentum fractions will be equal to @xmath429 using lo kinematics .
going to large values of @xmath110 and @xmath430 would then yield information on the gluon at large @xmath2 , since the @xmath86 contribution is much smaller there than that of the gluon .
this would also be the case for lepton pair @xcite and vector boson production in @xmath134 collisions , provided that one considers @xmath431 in order for the parton momentum fraction to be large . the analysis in ref .
@xcite suggests that photon + jet data with increased statistics from the lhc may help constrain the gluon and light quark pdfs .
another source of information on the gluon pdf is top quark pair production , for which the lowest order subprocess is @xmath432 . in ref .
@xcite the constraints on the gluon pdf provided by @xmath433 data from the tevatron and lhc were examined , with the results indicating that these data place strong constraints on the large-@xmath2 gluon pdf . also of note
is the analysis of ratios and double ratios of cross sections measured at different lhc energies @xcite .
such ratios have reduced systematic errors and have the potential to provide strong constraints on large-@xmath2 pdfs .
in addition to the issues discussed above , an obstacle for the determination of the gluon distribution , and the highly correlated @xmath34 values from pdf analysis , is that the calculations for the relevant processes are not available beyond nlo , with the important exceptions of the ( light quark contributions to the ) inclusive dis structure functions and drell - yan cross sections .
for example , the relative uncertainties of the jr09 gluon pdf in fig .  [ fig : gluon ] are below about 10% for @xmath434 , and increase as one proceeds to higher values of @xmath2 , where the constraints of the @xmath81 structure function from fixed target experiments used in that fit are weaker .
the uncertainty at large @xmath2 could in principle be reduced by including jet data . however
, this would also affect the central values obtained for the gluon pdf , changing them in the relevant region by an amount proportional to the missing nnlo corrections , which in the case of jet production are expected to be large @xcite .
global nnlo fits which include jet data , such as the mstw08 @xcite parametrization in fig .
[ fig : gluon ] , generally find larger gluon pdfs in the high-@xmath2 region and correspondingly also larger values of @xmath34 .
yet another factor which can alter the pdfs and the value of @xmath34 obtained in a particular analysis are the kinematic cuts applied ; of particular relevance are the cuts on the dis structure function data and the related treatment of contributions beyond leading twist @xcite .
although data selection ( data sets and cuts ) bring about certain arbitrariness in global pdf analysis , it can not account for all the existing differences between the various pdf sets , which are also affected by theoretical issues such as heavy quark schemes ( sec .
[ ssec : hq ] ) and the particular solutions of the rge equations used . even within a given framework
there are uncertainties due to the inability of the estimation procedure to find the optimal solution , for example , shortcomings of the parametrization to reproduce the optimal shape of the distributions , and the statistical estimation procedure .
these effects , referred in general as ` procedural bias ' , induce a dependence of the results on the choice of the scale at which the input distributions are parametrized , which in principle should not depend on these choices . for particularly sensitive quantities such as the gluon distribution and @xmath34
, this uncertainty can be comparable to the experimental uncertainty @xcite , as can be seen in fig .
[ fig : alphas ] .
values obtained for different parametrizations in ref .
the band indicates the uncertainties due to propagation of the experimental errors for a particular parametrization ( nnlo20).,width=377 ]    another example of this is the differences between `` standard '' distributions , generated from input pdfs at values for the input scale @xmath435  gev@xmath170 , and the so  called dynamical distributions @xcite , for which input scale is chosen to be @xmath436  gev@xmath170 . at these low scales
the input distributions naturally tend to valence  like ( positive definite ) shapes @xcite that vanish in the small@xmath2 limit , not only for valence but also for the sea quark and gluon input densities .
it was shown in refs .
@xcite , where also `` standard '' distributions were generated , that this implies that the dynamical distributions at small @xmath2 are more restricted and have smaller uncertainties than their standard counterparts .
furthermore , this procedure typically results in smaller @xmath34 values and steeper gluons in the small@xmath2 region , although both approaches provide comparable descriptions of data , in particular , of the inclusive dis cross sections from hera .
an alternative approach to fitting @xmath34 as part of the global fitting is to use the world average value for @xmath236 , thereby avoiding the correlations during the fitting process .
since @xmath34 is a parameter of the qcd lagrangian , it should be the same for all processes , not just those in the global fit .
of course , interesting information could be obtained by determining the pdfs both with a fixed value of @xmath236 and with a fitted value .
comparing the results could indicate which processes in the global fit are chiefly responsible for any differences .
the moments of leading twist pdfs can be evaluated from first principles using lattice qcd . within the operator product expansion ,
the moments of pdfs are related to matrix elements of local operators between hadron states , which can be computed numerically on the lattice .
the lattice pdf moments can be related to the quark and antiquark distributions in the nucleon as @xmath437 , \label{eq : mom_unp }                                  \\ \langle x^n \rangle_{\delta q } & = & \int_0 ^ 1 dx\ , x^n\ [ \delta q(x)+(-1)^n
\delta \bar   q(x ) ] , \label{eq : mom:_pol}\end{aligned}\ ] ] for the unpolarized and helicity distributions ( see sec .
[ sec : spin ] ) , respectively .
the lattice moments alternate between the total ( or @xmath11-even ) @xmath438 and valence ( @xmath11-odd ) @xmath348 distributions , depending on whether @xmath15 is even or odd .
while reconstructing the @xmath2 dependence of the pdfs from a finite number of calculated moments is challenging @xcite , comparison of the lattice moments with the phenomenological pdf moments can provide important tests of lattice techniques and provide constraints on specific pdfs which may be difficult to measure experimentally .
because of the discretization of space - time on the lattice , the rotational symmetry of the continuum is broken to the hyper - cubic subgroup , which introduces mixing with lower dimensional operators under renormalization . as a result , only the lowest few moments , corresponding to @xmath439 , can at present be reliably computed . as well as
the usual approximations of working at a finite lattice spacing @xmath18 , in a finite lattice volume @xmath440 , and at unphysically large values of the quark mass , simulations usually include contributions only from `` connected '' diagrams , representing operator insertions on quark lines connected to the source .
meaningful comparisons of lattice data can therefore only be made with nonsinglet combinations , such as @xmath441 , where the `` disconnected '' contributions cancel .
significant progress has been made in recent years in the computation of matrix elements of twist - two operators on the lattice , using dynamical ( unquenched ) configurations and quark ( or pion , by the gell - mann  oaks  renner relation ) masses approaching the physical limit .
unfortunately , the results for the lowest nontrivial ( @xmath442 ) nonsinglet moment of the unpolarized distributions , @xmath443 , lie systematically in the range @xmath444 at a scale @xmath445  gev@xmath170 , which is some 40% above the phenomenological pdf moments of @xmath446 @xcite .
a smaller , but persistent 10% overestimate of the lowest ( @xmath447 ) moment of the nonsinglet spin - dependent distribution , @xmath448 , has posed a serious challenge to lattice simulations @xcite .
it has been suggested @xcite that the discrepancy between the calculated and experimental pdf moments could be resolved by the nontrivial chiral behavior of the moments , when extrapolating the results from the lowest @xmath449 values at which lattice data exist to the physical mass .
verification of this behavior will require lattice simulations at very small quark masses , as well as at large enough volumes to fully include the physics of the pion cloud @xcite ( see also sec .
[ ssec : lightsea ] ) . while it may take some time before the lattice parameters can become sufficiently close to those in the physical realm
, such comparisons will continue to reflect the important synergy between lattice qcd and pdf phenomenology .
considerable progress has been made in understanding the spin structure of the nucleon since the first precision polarized dis experiments at cern in the late 1980s @xcite indicated an anomalously small fraction of the proton spin carried by quarks .
a rich program of spin - dependent inclusive and semi - inclusive dis , as well as polarized proton
proton scattering experiments has followed , vastly improving our knowledge of spin - dependent pdfs of the nucleon over the last two decades ( for a recent review , see ref .
@xcite and references therein ) .    while the spin - dependent data have not been as abundant as those available for constraining spin - averaged pdfs , the existing empirical information from polarized lepton and hadron facilities at cern @xcite , slac @xcite , desy @xcite , jefferson lab @xcite and rhic @xcite has enabled several dedicated global qcd analyses of spin - dependent parton distributions to be performed .
these include the nlo analyses from the more established dssv @xcite group , the european - based lss @xcite and bb @xcite groups , as well as the japanese asymmetry analysis collaboration ( aac ) @xcite , using the standard global fitting methodology outlined in sec .
[ sec : qcd ] .
more recent efforts have been made by the nnpdf collaboration @xcite , extending the neural network approach to the polarized sector , and the jefferson lab angular momentum ( jam ) collaboration @xcite , which has focused on utilizing data over a large kinematic range that includes the large-@xmath2 and low-@xmath40 regions , with a simultaneous fit of unpolarized pdfs .     and
@xmath450 , the antiquark @xmath451 , @xmath452 and @xmath453 , and polarized gluon @xmath454 distributions at @xmath191  gev@xmath170 , for the dssv09 @xcite , lss10 @xcite , bb10 @xcite , aac09 @xcite , and jam13 @xcite pdf sets.,width=548 ]    representative examples of the spin - dependent pdfs are shown in fig .
[ fig : spinpdfs ] , illustrating the total ( or @xmath11-even ) @xmath455 and @xmath456 pdfs , as well as the sea quark @xmath457 , @xmath458 and @xmath459 and polarized gluon @xmath152 distributions , for the dssv09 @xcite , lss10 @xcite , bb10 @xcite , aac09 @xcite , and jam13 @xcite pdf sets at @xmath191  gev@xmath170 .
the bb10 @xcite and jam13 @xcite analyses are based on inclusive dis data only , while the lss10 @xcite fit includes also semi - inclusive dis data .
the dssv09 pdfs are constrained in addition by polarized @xmath134 scattering data .
the latest aac analysis @xcite is based on inclusive dis data and a @xmath297-factor approximation for the nlo corrections for the @xmath134 data .    as for the unpolarized case in fig .
[ fig : nspdf ] , the @xmath455 distribution is the best constrained polarized pdf , mostly by measurements of the proton @xmath90 structure function over a relatively broad range of @xmath2 . the corresponding @xmath456 distribution , which has the opposite sign , is smaller in magnitude compared with the @xmath455 with somewhat larger uncertainties , especially at larger @xmath2 values .
the uncertainty band for the dssv09 parametrization of @xmath455 and @xmath456 is smaller than the variation between the different pdf sets , which reflects the fact that the systematic uncertainties associated with the choice of data sets and parametrization assumptions , as well as other theoretical inputs , are currently larger than the experimental errors .
the polarization of the sea is considerably smaller , and more strongly dependent upon assumptions about the flavor dependence in the analysis of semi - inclusive dis data . at present
there is no conclusive experimental evidence of a nonzero light quark sea , @xmath457 and @xmath458 , although there is a slight trend towards a more negatively polarized @xmath273 distribution than @xmath272 , with a positive @xmath460 .
the polarization of the strange sea is also very small , in contrast to suggestions from the early spin - dependent dis data analyses of a large negative @xmath459 .
the polarized gluon distribution is essentially unconstrained by existing inclusive and semi - inclusive dis data , and almost all information on @xmath152 comes from measurements of @xmath461 production in semi - inclusive dis @xcite , and inclusive pion and jet production in polarized @xmath134 scattering @xcite .
the data are consistent with a small value of @xmath462 , consistent with zero , although new measurements from rhic have the promise of resolving a small nonzero distribution .      at large values of @xmath2 ,
the spin - dependent pdfs are even more sensitive to the quark - gluon dynamics responsible for spin - flavor symmetry breaking than the spin - averaged valence quark pdfs discussed in sec .
[ ssec : valence ] .
however , while considerable progress has been made in determining spin - dependent pdfs over the last two decades , especially in the small-@xmath2 region , relatively little attention has been paid to structure function measurements at large @xmath2 .
the dearth of data in the valence region is especially striking for the neutron , where , with the exception of several data points from jefferson lab hall  a @xcite for the polarization asymmetry @xmath463 extending to @xmath464 , there are essentially no constraints for @xmath465 . on the other hand , there are a number of predictions for the behavior of polarized pdfs and ratios of polarized to unpolarized pdfs in the limit as @xmath165 that differ even in sign .
at large @xmath2 there are dramatically different predictions for the valence @xmath283 and @xmath282 distributions from nonperturbative and perturbative calculations . in the simple su(6 ) spin - flavor symmetric quark model , the equal probabilities of the spin-0 and spin-1 spectator diquark configurations lead to the simple relation @xmath466 , and corresponding polarized to unpolarized ratios @xmath467 and @xmath468 . in this limit the polarization asymmetry of the proton , which at lo can be written @xmath469 is predicted to be @xmath470 , while that of the neutron , @xmath471 would vanish , @xmath472 .    while spin - flavor symmetry is of course broken in nature ( in some cases badly ) , with the unpolarized @xmath142 quark distribution significantly softer than the @xmath141 , for example , the details of this breaking can affect the spin - dependent pdfs in strikingly different ways .
the mechanism of spin-1 diquark suppression , which leads to the vanishing of the @xmath145 ratio as @xmath165 , also predicts that @xmath473 in this limit , while @xmath474 remains unchanged from the su(6 ) expectation .
arguments based on helicity conservation in perturbative qcd , on the other hand , predict that for valence quarks in the ground state of the nucleon the scattering at large @xmath2 is predominantly from configurations in which the quark is coupled to a diquark of zero helicity @xcite . in this case
the pdfs associated with a quark polarized in the same direction as the nucleon are favored over the helicity - antialigned configurations , @xmath475 as @xmath165 .
this model leads to the expectation that @xmath473 and @xmath476 in the @xmath165 limit , and consequently a rapid approach to unity of both the proton and neutron @xmath92 asymmetries .
an additional feature of the helicity conservation model is that it correlates the @xmath165 behavior of the pdfs with the large-@xmath40 behavior of elastic form factors @xcite , which can be measured in elastic lepton
nucleon scattering experiments at jefferson lab and elsewhere @xcite .
this is also closely related to the notion of quark - hadron duality , which connects averages over structure functions in the resonance region at low @xmath70 with those extrapolated from the deep - inelastic continuum at higher @xmath70 ( lower @xmath2 ) @xcite .
since the existing dis data are for the most part consistent with a vanishingly small neutron asymmetry at @xmath477 , a very dramatic change in @xmath463 would be expected at larger @xmath2 .
an even more rapid transition at high @xmath2 may occur if quark orbital angular momentum plays an important role in nucleon structure @xcite .     and
@xmath478 distributions and their uncertainties for the jam13 pdf set @xcite at @xmath191  gev@xmath170 , illustrating the effects of the nuclear smearing ( green dotted ) , target mass ( blue dashed ) and higher twist ( red solid ) corrections , relative to the reference fit ( yellow dot - dashed ) .
note that the yellow and green bands overlap for @xmath479 and the green and blue bands overlap for @xmath478 .
_ ( right ) _ corresponding ratios of polarized to unpolarized @xmath480 and @xmath481 distributions.,title="fig:",width=302 ]   and @xmath478 distributions and their uncertainties for the jam13 pdf set @xcite at @xmath191  gev@xmath170 , illustrating the effects of the nuclear smearing ( green dotted ) , target mass ( blue dashed ) and higher twist ( red solid ) corrections , relative to the reference fit ( yellow dot - dashed ) .
note that the yellow and green bands overlap for @xmath479 and the green and blue bands overlap for @xmath478 .
_ ( right ) _ corresponding ratios of polarized to unpolarized @xmath480 and @xmath481 distributions.,title="fig:",width=302 ]    experimentally , the extraction of leading twist spin - dependent pdfs at high @xmath2 faces similar challenges to those discussed in sec .
[ ssec : valence ] for the unpolarized distributions .
in particular , information on the structure of the neutron is typically obtained from experiments involving polarized deuterium or @xmath221he nuclei , requiring the nuclear corrections to be understood .
analysis of high-@xmath2 data also necessitates careful treatment of finite-@xmath40 corrections such as target mass and higher twist contributions .
the effects of these corrections are illustrated in fig .
[ fig : pol_val ] , which shows the total @xmath455 and @xmath456 pdfs from the jam13 analysis @xcite for various fits with or without the finite-@xmath40 and nuclear corrections . at small values of @xmath2
their effects are negligible ; however , at @xmath482 they can give up @xmath483 corrections for @xmath455 and more @xmath177 corrections for @xmath456 ( the relative correction can be even larger at @xmath484 , although the pdfs are not constrained in this region ) .
the same effects are more clearly illustrated through the polarization ratios @xmath480 and @xmath481 in fig .
[ fig : pol_val ] , where the unpolarized distributions are fitted simultaneously within the same analysis as the polarized @xcite . this in principle eliminates any bias arising from the use of spin - averaged pdfs taken from analyses performed under different sets of assumptions . in the intermediate-@xmath2 region
the ratios for both @xmath141 and @xmath142 quarks are generally consistent with the symmetric quark model expectations , with the @xmath480 increasing towards unity at larger @xmath2 .
the @xmath481 ratio , on the other hand , remains negative for all @xmath2 where it is constrained , and shows no indication of the upturn predicted by the helicity conservation models .
the nuclear and finite-@xmath40 effects can significantly impact the limiting behavior as @xmath165 , and clearly additional data are needed in order to constrain their @xmath2 dependence at high @xmath2 .
a dedicated program of large-@xmath2 structure function measurements is planned at the 12  gev energy upgraded jefferson lab facility @xcite .
more complete information on spin - dependent pdfs will also be available once the @xmath91 structure function of the nucleon is measured more accurately @xcite .
this will provide additional information on possible higher twist corrections that enter in the extraction of leading twist pdfs from dis polarization asymmetries .
as for the unpolarized inclusive charged lepton dis measurements , inclusive @xmath90 structure function experiments measure @xmath11-even combinations of pdfs , @xmath278 . the individual quark and antiquark distributions can be separated with the help of semi - inclusive charged lepton scattering data , with coincident measurement of fast pions or kaons in the final state . typically such experiments measure the semi - inclusive polarization asymmetry , which at lo can be written @xmath485 where @xmath486 is the fragmentation function for the scattered quark or antiquark to produce a hadron @xmath285 ( @xmath487 ) in the current fragmentation region with a fraction @xmath488 . for large @xmath288 , the produced hadron has a high probability of containing the scattered parton , hence providing a tag on the initial state distribution of quarks and antiquarks .
the fragmentation functions @xmath289 are extracted from independent fits to single hadron production cross sections in @xmath290 , @xmath134 and other reactions .
, from the compass semi - inclusive dis measurement at @xmath489  gev@xmath170 @xcite , compared with several parametrizations and model calculations , and the unpolarized @xmath305 asymmetry ( solid curve ) .
( figure from ref .  @xcite . )
_ ( right ) _ single polarization asymmetry @xmath490 for inclusive @xmath57 production in @xmath491 scattering from the star experiment at rhic , at central electron pseudorapidity @xmath492 .
( figure from ref .
copyright ( 2011 ) by the american physical society.),title="fig:",width=340 ] , from the compass semi - inclusive dis measurement at @xmath489  gev@xmath170 @xcite , compared with several parametrizations and model calculations , and the unpolarized @xmath305 asymmetry ( solid curve ) .
( figure from ref .
_ ( right ) _ single polarization asymmetry @xmath490 for inclusive @xmath57 production in @xmath491 scattering from the star experiment at rhic , at central electron pseudorapidity @xmath492 . ( figure from ref .
copyright ( 2011 ) by the american physical society.),title="fig:",width=226 ]    the program of semi - inclusive dis measurements from proton and deuteron targets at hermes @xcite , compass @xcite , and earlier smc @xcite , have produced intriguing glimpses into the flavor content of the polarized sea .
the extracted sea quark and antiquark asymmetries are generally rather small , and consistent with zero within the current uncertainties . for the light quark polarized sea
, there is a slight trend towards a positive @xmath457 ( with an @xmath2-integrated value of @xmath493 from the latest compass data @xcite ) and a negative @xmath458 ( with an @xmath2-integrated value of @xmath494 ) , so that the difference @xmath460 is slightly positive , at the @xmath495 level @xcite , as fig .
[ fig : pol_sea ] indicates .
this is in contrast to the large negative values of the unpolarized @xmath496 asymmetry ( see sec .  [
ssec : lightsea ] ) , which would disfavor models such as the chiral quark soliton model in the large-@xmath497 limit @xcite that predict large @xmath498 .
very small antiquark spin asymmetries would be expected in models based on the pion cloud of the nucleon , such as those discussed in sec .
[ ssec : lightsea ] .
one should caution , however , that the experimental distributions in fig .
[ fig : pol_sea ] were extracted assuming lo expressions for the semi - inclusive asymmetries , which , strictly speaking , can not be directly compared with nlo pdfs @xcite .
the comparison should therefore be viewed as a qualitative one , although the general indication of a small polarized non - strange sea is unlikely to be qualitatively modified in an nlo treatment .
the semi - inclusive dis measurements are also consistent with a very small strange quark polarization , @xmath296 .
the hermes data @xcite give a slightly positive value when integrated between @xmath499 and 0.6 , while the compass data @xcite are consistent with zero .
this is in contrast with the analysis of inclusive dis data on protons , deuteron and @xmath221he , which suggests a small negative value .
the extractions of @xmath296 from the inclusive and semi - inlcusive data are somewhat dependent on the assumptions of su(3 ) flavor symmetry for the axial charges of hyperons , and on the specific fragmentation functions used in the semi - inclusive analysis .
future measurements of kaon fragmentation functions @xmath500 should improve the accuracy of the flavor decomposition of the proton s helicity distributions .    a complementary method , discussed by bourrely and soffer @xcite , of constraining the @xmath457 and @xmath458 distributions , which does not depend on knowledge of fragmentation functions , is through longitudinal polarization asymmetries in inclusive @xmath57 boson production from @xmath134 scattering with one proton polarized and the other unpolarized , @xmath501 .
the asymmetry for @xmath139 production , for example , can then be written at lo as @xcite @xmath502 where @xmath503 is the scattering angle of the electron in the partonic center - of - mass frame , while the corresponding @xmath137 asymmetry at lo is @xmath504 at large negative rapidity , where @xmath299 , the @xmath139 asymmetry then directly probes the @xmath457 pdf , @xmath505 , while at large positive rapidity one has @xmath132 , and the asymmetry is sensitive to the @xmath282 pdf at large @xmath2 , @xmath506 . for @xmath137 production
, the term proportional to @xmath507 will be similarly enhanced at small @xmath299 ; however , the angular factor will produce a suppression of this contribution at backward angles relative to the @xmath508 term , so that both terms will give competing contributions .
the first data from the star experiment ( run  9 ) at rhic @xcite are shown in fig .
[ fig : pol_sea ] for zero electron rapidity , and are generally in agreement with nlo calculations using existing pdf parametrizations . data from run  12 at rhic have subsequently been taken at both forward and backward rapidities @xcite , and will provide important constraints on the antiquark polarization .
information on the polarized gluon distribution @xmath152 can be obtained in several complementary ways .
firstly , from the @xmath40 evolution equations , the evolution of the quark singlet contribution to the @xmath90 structure function mixes with the gluon contribution , so that in principle the study of scaling violations in @xmath90 can constrain @xmath152 . in practice
, however , the @xmath40 range and precision of the available inclusive dis data are not sufficient to provide meaningful constraints , and dis - only fits of spin - dependent pdfs effectively leave @xmath152 undetermined .
values extracted at lo from charm and high-@xmath20 hadron production data from compass , smc and hermes .
( figure from ref .
copyright ( 2013 ) by the american physical society . )
_ ( right ) _ longitudinal double - spin asymmetry @xmath509 for inclusive jet production , as a function of the jet transverse momentum , in polarized @xmath510 scattering from the star experiment at rhic @xcite , compared with calculations based on several pdf parametrizations .
( figure from ref .
@xcite.),title="fig:",width=321 ]   values extracted at lo from charm and high-@xmath20 hadron production data from compass , smc and hermes .
( figure from ref .
copyright ( 2013 ) by the american physical society . )
_ ( right ) _ longitudinal double - spin asymmetry @xmath509 for inclusive jet production , as a function of the jet transverse momentum , in polarized @xmath510 scattering from the star experiment at rhic @xcite , compared with calculations based on several pdf parametrizations .
( figure from ref .
@xcite.),title="fig:",width=321 ]    more direct information on gluon polarization has come from the production of charm mesons and high-@xmath20 hadrons in semi - inclusive dis , through the process of photon
gluon fusion ( see sec .
[ ssec : hqresults ] ) .
data on hadron production from hermes , smc and compass , and open charm data from compass have been analyzed over a range of @xmath2 values , averaging from @xmath511 to @xmath512 , and @xmath40 between @xmath513  gev@xmath170 and 13  gev@xmath170 .
most of the analyses have been performed at lo , finding values of @xmath462 consistent with zero , albeit with large uncertainties , as illustrated in fig .
[ fig : pol_g ] .
the recent compass experiment @xcite was the first to be analyzed at nlo , yielding @xmath514 over the range @xmath515 with an average @xmath516 .
an alternative reaction that can be used to constrain @xmath152 is the inclusive production of jets or neutral pions in polarized proton
proton scattering .
the longitudinal double - spin asymmetry for this process , @xmath517 where @xmath518 is the differential cross section for longitudinally polarized protons with equal ( opposite ) helicities , is then sensitive to the polarized quark
gluon scattering cross section at the parton level .
following an early measurement of @xmath509 for multi-@xmath519 pairs from two - jet type events by the fermilab e581/704 collaboration @xcite , which found a small @xmath462 at @xmath520 , most of the experimental effort at measuring @xmath152 in @xmath134 scattering has been by the phenix and star collaborations at rhic .
the star experiment measures inclusive jet production in @xmath521 at an invariant center of mass energy @xmath302  gev .
the @xmath509 data from the 2006 and 2009 runs @xcite shown in fig .  [ fig : pol_g ] provide the clearest evidence to date for a small but nonzero @xmath152 .
a nonzero asymmetry is also observed by star in preliminary data on dijet events at mid - rapidity @xcite , which may provide even better constraints on the shape of @xmath152 .
the phenix data on inclusive neutral pion production , @xmath522 , from 2005 to 2009 suggest a very small asymmetry up to @xmath523  gev @xcite . taken together , the inclusive rhic data from several different channels suggest a small but non - vanishing gluon polarization in the proton .
though not as large as had been proposed in some early explanations of the small contribution of quarks to the proton spin , the precise role played by @xmath152 in the proton spin budget will be intriguing to understand from future measurements .
in this topical review we have presented a summary of the current knowledge of spin - averaged and ( longitudinal ) spin - dependent pdfs .
while great progress has been made in recent years , there are still regions of bjorken-@xmath2 where the pdf uncertainty bands are uncomfortably large .
this is primarily due to either lack of data in kinematic regions that would otherwise provide constraints on the pdfs , or to the particular pdf being smaller in that region than those of other flavors , or both .
one might naively conclude that if existing data do not constrain particular pdfs in certain domains , then it may not be too important to seek ways to constrain them there .
however , this viewpoint overlooks two important points : ( i ) there is intrinsic interest in the behavior of the different pdfs and relations to the properties of the parent hadrons , and ( ii ) it is possible that new phenomena may occur in extended kinematic regions which would require precise knowledge of the pdfs there .    as an example , suppose one were searching for a massive state produced at large rapidity in hadron
hadron collisions .
this would typically require one parton to have a large value of @xmath2 and the other a small value .
the errors at large @xmath2 may well dominate the uncertainty in the resulting prediction .
the key here is to measure standard model processes in expanded kinematic regions and include these processes in the global fits .
the resulting reduction in the errors can then be used to improve searches for new phenomena in these regions . moreover ,
if one reduces the errors at large values of @xmath2 , this will necessarily lead to reductions in the errors at intermediate @xmath2 .
since deuterium targets in dis are a prime source of information on the @xmath142 quark pdf at large values of @xmath2 , nuclear corrections are a significant source of uncertainty in its determination .
several planned experiments at the 12  gev energy upgraded jefferson lab @xcite will use techniques that reduce the effects of nuclear corrections on the extraction of @xmath145 at large @xmath2 .
constraints on the @xmath145 ratio in the proton can also be obtained , free of nuclear effects , through the measurement of @xmath57 charge asymmetries in @xmath134 collisions at very high rapidity @xcite , although this would require reconstruction of the @xmath57 distributions themselves from the decay lepton distributions and , unlike in @xmath524 collisions at the tevatron , would require precise knowledge of the @xmath272 and @xmath273 distributions @xcite .
another important area of investigation is the flavor separation of the pdfs .
for example , currently our knowledge of the @xmath254 and @xmath255 pdfs comes primarily from one neutrino experiment which used an iron target @xcite , thus requiring the use of model - dependent nuclear corrections .
other sources , for example weak vector boson production at the lhc @xcite , can provide additional constraints , especially as the statistics of the measurements improve . in a complementary effort ,
the hermes collaboration has measured @xcite @xmath525 production in semi - inclusive dis and used these data to place constraints on the strange quark pdf .
while an interesting analysis , it does depend on knowledge of the @xmath525 fragmentation functions , and a recent , more complete analysis @xcite of the hermes data suggests a somewhat smaller @xmath254 pdf than that extracted in ref .
better knowledge of the fragmentation functions would clearly help in the interpretation of the semi - inclusive dis measurements .
ideally , a high statistics neutrino experiment on _ hydrogen _ could further provide knowledge of both the @xmath254 and @xmath255 pdfs through the measurement of charm production ( using the @xmath526 and @xmath527 subprocesses ) without having to make model - dependent nuclear corrections .
note that such an experiment could also yield information on the @xmath145 ratio at high values of @xmath2 .
information on the su(2 ) flavor asymmetry of the proton sea will be provided by experiment e906/seaquest @xcite at fermilab , which is currently taking data on muon pair production using a 120  gev proton beam on a variety of targets .
it is designed to improve our knowledge of @xmath528 over the @xmath2 range out to @xmath529 , whereas the currently available data run out of statistics at @xmath530 .
the experiment will further enable a study of the effects of nuclear binding on the sea quark pdfs .
there is also the possibility of follow - on measurements of lepton pair production at the nuclear and particle physics facility at j - parc @xcite in japan .
the question of charge symmetry violation in pdfs can be addressed in future through study of drell - yan lepton pair production in @xmath531 scattering from the deuteron , which could be measured at the fermilab main injector or with the compass experiment at cern @xcite , or via measurement of @xmath532 electroproduction ratios in semi - inclusive dis from isoscalar targets @xcite .
another avenue to explore will be parity - violating dis on a deuterium target @xcite , for which the parity - violating asymmetry arising from @xmath533 interference is expected to be independent of hadronic structure at large @xmath2 in the absence of csv effects @xcite . beyond searches for csv effects ,
a program of parity - violating dis studies at the 12  gev jefferson lab @xcite could open up a new window on the vector @xmath534 interference structure functions , as well as provide glimpses into the parity - odd @xmath535 structure function .    for heavier quarks , the issue of intrinsic charm will require measuring the charm structure function at large values of @xmath2 , which may be possible at a future electron - ion collider facility @xcite .
the semi - inclusive production of prompt photons and charmed jets in @xmath134 collisions at the lhc , @xmath536 , may produce an enhancement at large transverse momenta from intrinsic charm @xcite .
measurements of jet , dijet , isolated @xmath519 , and @xmath519 + jet production at the lhc will all play a role in helping to further constrain the gluon pdf .
of particular interest would be @xmath519 + jet production where both the @xmath519 and the jet have equal and opposite large rapidities .
the dominant @xmath274 subprocess would thus require both a high-@xmath2 gluon and a high-@xmath2 quark .
furthermore , continued progress on resummation , especially for the fragmentation contribution at fixed target energies may provide additional information on the gluon at high values of @xmath2 .    in the polarized sector
, several experiments planned at jefferson lab will measure the longitudinal polarization asymmetries for hydrogen , deuterium and @xmath221he to unprecedented large values of @xmath2 @xcite .
these will constrain the spin - dependent valence quark pdfs , especially the @xmath282 distribution , to values of @xmath537 .
this should help answer long - standing questions about the behavior of the @xmath538 ratios as @xmath165 , and reveal insights into the role of quark orbital angular momentum in the nucleon .
the flavor asymmetry of the polarized light quark sea , @xmath460 , will be further probed via @xmath70 boson production in polarized @xmath134 scattering at rhic at @xmath539  gev . both the star and phenix experiments anticipate significantly improved statistics and kinematic coverage in their 2012 and 2013 runs @xcite . the strange quark polarization , despite appearing to be smaller than initially surmised from early polarized dis experiments ,
is neverthess important to pin down , as it allows a portal to nonperturbative qcd effects in the nucleon .
improved constraints on kaon fragmentation functions should lead to more reliable determinations of @xmath296 in semi - inclusive dis .
finally , the polarized gluon distribution is being actively investigated in semi - inclusive deep - inelastic hadron production at compass , and through double - spin asymmetries for pion and jet production in polarized @xmath134 collisions at rhic .
preliminary results from star on jet and dijet cross sections hint at nonzero values for @xmath152 , which will need to be confirmed by further measurements at small @xmath2 .
access to @xmath152 below @xmath540 will be possible with future running at @xmath541  gev and at large forward rapidity @xcite , although reaching values down to @xmath542 will only be possible with an electron - ion collider @xcite .
on the theory front , continued progress on nnlo calculations will eventually result in global fits where all the processes are treated consistently at nnlo , for both unpolarized and polarized observables .
the increased precision and decreased scale dependence will help to reduce the uncertainty bands on the pdfs .
the parametrization dependence of global fits will continue to be explored , allowing for greater flexibility and reducing bias introduced with the use of specific forms , such as for the @xmath145 ratio at large @xmath2 .
improvements in electroweak radiative corrections are being sought , incorporating , for example , pdfs of a photon with @xmath543 @xmath40 evolution .
progress on computing moments of pdfs in lattice qcd is also anticipated , with more reliable simulations performed on larger lattice volumes and at quark masses near the physical limit .
a new approach @xcite to calculating the @xmath2 dependence of the pdfs directly in the infinite momentum frame may provide complementary constraints on pdfs that are difficult to access experimentally .
we thank a.  accardi and k.  griffioen for helpful comments and a careful reading of the manuscript .
the work of p.j .- d . and
was supported by doe contract no .
de - ac05 - 06or23177 , under which jefferson science associates , llc operates jefferson lab . the work of j.f.o .
was supported by doe contract no .
de - fg02 - 13er41942 .
10        a.  accardi , _ large - x connections of nuclear and high - energy physics _ , to appear in mod . phys .
a ; pos ( confinement x ) 227 , + http://pos.sissa.it/archive/conferences/171/227/confinement%20x_227.pdf .
jefferson lab experiment e12 - 06 - 109 , _ the longitudinal spin structure of the nucleon _ ,
d.  crabb _ et al .
_ , spokespersons .
jefferson lab experiment e12 - 06 - 110 , _ measurement of the neutron spin asymmetry @xmath463 in the valence quark region using an 11  gev beam in hall  c _ , j .- p .
chen _ et al .
_ , spokespersons .
jefferson lab experiment e12 - 06 - 122 , _ measurement of the neutron asymmetry @xmath463 in the valence quark region using 8.8 and 6.6  gev beam energies and the bigbite spectrometer in hall  a _ , b.  wojtsekhowski _ et al .
_ , spokespersons .
a.  accardi _ et al .
_ , _ electron ion collider : the next qcd frontier - understanding the glue that binds us all _ , bnl-98815 - 2012-ja , jlab - phy-12 - 1652 , arxiv:1212.1701 [ nucl - ex ]
.                                                              jefferson lab experiment e12 - 10 - 102 [ bonus12 ] , _ the structure of the free neutron at large @xmath2-bjorken _ , s.  bltmann , m.  e.  christy , h.  fenker , k.  griffioen , c.  e.  keppel , s.  kuhn and w.  melnitchouk , spokespersons .
jefferson lab experiment e12 - 10 - 103 [ marathon ] , _ measurement of the @xmath544 , @xmath145 ratios and @xmath545 emc effect in dis off the tritium and helium mirror nuclei _ ,
g.  g.  petratos , j.  gomez , r.  j.  holt and r.  d.  ransome , spokespersons .
fermilab e906 experiment , _ drell - yan measurements of nucleon and nuclear structure with the fermilab main injector _ , d.  f.  geesaman and p.  e.  reimer , spokespersons ; http://www.phy.anl.gov/mep/seaquest/index.html .
a.  airapetian _ et al . _
[ hermes collaboration ] , in print ; h.  jackson , talk given at the _ xxi international workshop on deep
inelastic scattering and related subjects ( dis2013 ) _ , marseilles , france , april 22 - 26 , 2013 .
jefferson lab experiment e12 - 09 - 002 , _ precise measurement of @xmath532 ratios semi - inclusive dis : charge symmetry violating quark distributions _
, k.  hafidi _ et al .
_ , spokespersons .
c.  y.  prescott _
et al . _ ,
b * 77 * , 347 ( 1978 ) .