arxiv:hep-ph/ v4 24 Nov 2005

Similar documents
Physics Oct Reading

Electric dipole moment of nuclei

x a x a Lecture 2 Series (See Chapter 1 in Boas)

The standard deviation of the mean

1 Adiabatic and diabatic representations

a b c d e f g h Supplementary Information

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

Analysis of Experimental Measurements

ANALYSIS OF EXPERIMENTAL ERRORS

Science & Technology Peer Review. Nucleon Spin Structure

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle

Linear Regression Demystified

Quantum Annealing for Heisenberg Spin Chains

R is a scalar defined as follows:

Kinetics of Complex Reactions

Statistical Analysis on Uncertainty for Autocorrelated Measurements and its Applications to Key Comparisons

Statistics 511 Additional Materials

Probability, Expectation Value and Uncertainty

Simona Malace University of South Carolina

Sequences. Notation. Convergence of a Sequence

The Born-Oppenheimer approximation

Coherent and incoherent DVCS and π 0 production on 4 He

GUIDELINES ON REPRESENTATIVE SAMPLING

Physics 232 Gauge invariance of the magnetic susceptibilty

The Riemann Zeta Function

Regression with an Evaporating Logarithmic Trend

Chapter 6 Sampling Distributions

The Random Walk For Dummies

Tests of Hypotheses Based on a Single Sample (Devore Chapter Eight)

1 6 = 1 6 = + Factorials and Euler s Gamma function

FIR Filter Design: Part II

Singular Continuous Measures by Michael Pejic 5/14/10

Beurling Integers: Part 2

Chapter 7: The z-transform. Chih-Wei Liu

CHAPTER 10 INFINITE SEQUENCES AND SERIES

Topic 9: Sampling Distributions of Estimators

ENGI 4421 Confidence Intervals (Two Samples) Page 12-01

Hydrogen (atoms, molecules) in external fields. Static electric and magnetic fields Oscyllating electromagnetic fields

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions

SECTION 2 Electrostatics

THE KALMAN FILTER RAUL ROJAS

Physics of n-a. Ulrich Mosel

True Nature of Potential Energy of a Hydrogen Atom

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Optimally Sparse SVMs

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion

Chapter 9 - CD companion 1. A Generic Implementation; The Common-Merge Amplifier. 1 τ is. ω ch. τ io

Statistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Infinite Sequences and Series

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT

6.3 Testing Series With Positive Terms

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Streamfunction-Vorticity Formulation

A Brief Introduction to the Physical Basis for Electron Spin Resonance

Properties and Hypothesis Testing

The Quark Puzzle A 3D printable model and/or paper printable puzzle that allows students to learn the laws of colour charge through inquiry.

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures

10. Comparative Tests among Spatial Regression Models. Here we revisit the example in Section 8.1 of estimating the mean of a normal random

Lecture 1 Probability and Statistics

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

L = n i, i=1. dp p n 1

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

11 Correlation and Regression

THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS

Lecture 19: Convergence

Discrete Mathematics for CS Spring 2008 David Wagner Note 22

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples:

PHY4905: Nearly-Free Electron Model (NFE)

ECON 3150/4150, Spring term Lecture 3

5.76 Lecture #33 5/08/91 Page 1 of 10 pages. Lecture #33: Vibronic Coupling

arxiv:hep-ph/ v1 13 Oct 2004

PSYCHOLOGICAL RESEARCH (PYC 304-C) Lecture 9

TEACHER CERTIFICATION STUDY GUIDE

7.1 Convergence of sequences of random variables

1 Inferential Methods for Correlation and Regression Analysis

A statistical method to determine sample size to estimate characteristic value of soil parameters

1 Review of Probability & Statistics

Name Solutions to Test 2 October 14, 2015

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =

Random Variables, Sampling and Estimation

Topic 9: Sampling Distributions of Estimators

Math 2784 (or 2794W) University of Connecticut

Measurement uncertainty of the sound absorption

Chapter 3. Strong convergence. 3.1 Definition of almost sure convergence

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1

Math 113 Exam 3 Practice

Topic 9: Sampling Distributions of Estimators

a for a 1 1 matrix. a b a b 2 2 matrix: We define det ad bc 3 3 matrix: We define a a a a a a a a a a a a a a a a a a

6 Sample Size Calculations

Recent Results from T2K

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model.

Geometry of LS. LECTURE 3 GEOMETRY OF LS, PROPERTIES OF σ 2, PARTITIONED REGRESSION, GOODNESS OF FIT

If, for instance, we were required to test whether the population mean μ could be equal to a certain value μ

Because it tests for differences between multiple pairs of means in one test, it is called an omnibus test.

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

Time-Domain Representations of LTI Systems

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

SNAP Centre Workshop. Basic Algebraic Manipulation

Transcription:

JLAB-THY-05-437 RM3-TH/05-10 Leadig twist momets of the eutro structure fuctio F 2 arxiv:hep-ph/0510189v4 24 Nov 2005 M. Osipeko 1,2, W. Melitchouk 3, S. Simula 4, S. Kulagi 5, G. Ricco 1,6 1 Istituto Nazioale di Fisica Nucleare, Sezioe di Geova, 16146 Geova, Italy 2 Skobeltsy Istitute of Nuclear Physics, 119992 Moscow, Russia 3 Jefferso Lab, Newport News, Virgiia 23606, USA 4 Istituto Nazioale di Fisica Nucleare, Sezioe Roma III, 00146 Roma, Italy 5 Istitute for Nuclear Research of Russia Academy of Sciece, 117312 Moscow, Russia 6 Dipartimeto di Fisica dell Uiversità di Geova, 16146 Geova, Italy Abstract We perform a global aalysis of eutro F2 structure fuctio data, obtaied by combiig proto ad deutero measuremets over a large rage of kiematics. From these data the lowest momets ( 10) of the leadig twist eutro F2 structure fuctio are extracted. Particular attetio is paid to uclear effects i the deutero, which become icreasigly importat for the higher momets. Our results for the osiglet, isovector p combiatio of the leadig twist momets are used to test recet lattice simulatios. We also determie the lowest few momets of the higher twist cotributios, ad fid these to be approximately isospi idepedet, suggestig the possible domiace of ud correlatios over uu ad dd i the ucleo.

1 Itroductio Iclusive lepto scatterig experimets have for some time bee the stadard tool with which to study the iteral quark-gluo structure of ucleos ad uclei. I upolarized scatterig, this structure is represeted through two structure fuctios, F 1 ad F 2, which parameterize the differetial cross sectio, d 2 σ dωde = σ Mott [ 2F1 (x,q 2 ) M ta 2 θ 2 + F 2(x,Q 2 ) ν ], (1) where σ Mott = (4α 2 E 2 /Q 4 )cos 2 θ 2 is the Mott cross sectio for the scatterig from poit particles, with Ω the scattered lepto solid agle (θ is the polar scatterig agle), ad E the recoil lepto eergy. Here ν is the eergy trasfer ad M is the ucleo mass. The structure fuctios F 1,2 are fuctios of the four-mometum trasfer squared q 2 Q 2, ad the Bjorke scalig variable x = Q 2 /2Mν. I the Bjorke limit (Q 2,ν, x fixed), the structure fuctios ca be related to parto distributio fuctios, which give the light-coe mometum distributios of quarks (ad gluos, or partos i geeral) i the hadro. At leadig order, oe has F 2 = 2xF 1 = x ( i e2 i qi (x,q 2 ) + q i (x,q 2 ) ). At fiite Q 2 radiative correctios ad higher twist effects give rise to a Q 2 depedece i the structure fuctios, leadig to scalig violatios. The coectio betwee the partoic structure of hadros ad QCD ca be most readily made through momets, or x-weighted itegrals, of structure fuctios, which through the operator product expasio (OPE) ca be related to hadroic matrix elemets of local operators. I the case of the F 2 structure fuctio, to which we restrict ourselves i the followig, the momet is defied by M (Q 2 ) = 1 0 dx x 2 F 2 (x,q 2 ). (2) While i the Bjorke limit the momets are give by real umbers, idepedet of Q 2, at fiite Q 2 the reormalizatio group equatios lead to a specific Q 2 depedece of the momets that ca be calculated from perturbative QCD. I fact, observatio of this Q 2 depedece was istrumetal i establishig QCD as the theory of the strog uclear iteractios. The momets ca also be calculated operturbatively usig lattice QCD techiques. A wealth of iformatio has bee accumulated over several decades o the structure of the proto, deutero, as well as heavier uclei. To completely describe the structure of the ucleo, however, requires kowledge of both the proto ad eutro structure fuctios. Typically, the eutro structure fuctio is obtaied by combiig data o the deutero ad proto. I most aalyses the deutero structure fuctio is assumed to be give by a simple sum of free proto ad eutro structure fuctios, F D 2 F p 2 +F 2, which is a good approximatio for some regios of kiematics (0.1 x 0.4). However, it is well kow that at small x (x 0.1) uclear shadowig reduces the deutero structure fuctio relative to the free ucleo [1, 2, 3, 4], while at large x (x 0.4) 2

uclear bidig ad Fermi motio correctios are expected to give rise to a small uclear EMC effect i the deutero. I terms of momets, ay uclear effect at large x will maifest itself i high momets (large ). It is crucial, therefore, that a momet aalysis of the eutro structure fuctio icludes a careful treatmet of uclear correctios. I this paper we for the first time perform such a aalysis, utilizig the recet high-precisio data from Jefferso Lab together with earlier data o F p 2 ad F 2 D (see Refs. [5, 6] ad refereces therei). The earlier world data o the F 2 structure fuctio cover a wide rage of kiematics, from Q 2 0 to 100 (GeV/c) 2. However, the large-x domai has bee almost completely uexplored, due to either experimetal limitatios or itepretatio difficulties. New, high precisio data from Jefferso Lab have allowed oe to fill this gap which has i the past preveted a complete momet based aalysis of data. Jefferso Lab data cover the etire resoace regio startig from the elastic peak (quasi-elastic for the deutero) ad extedig to ivariat masses of the produced hadroic system of about 3.5 (GeV/c). This regio is of particular importace for higher momets. I fact, already for the = 4 momet the cotributio of this kiematic regio is sigificat for Q 2 < 5 (GeV/c) 2. As a cosequece for Q 2 < 5 (GeV/c) 2 ad > 2 the ucertaities come mostly from Jefferso Lab data. The lowest momet ( = 2) is also sesitive to the small-x regio, where much data, i particular for Q 2 < 5 (GeV/c) 2, have bee accumulated i previous experimets. For the proto F 2 structure fuctio the lowest x were achieved with the HERA collider data which reached dow to x 10 5. For the deutero F 2 structure fuctio, the miimum x values are somewhat larger, x 10 4, ad were obtaied with muo beams i fixed target experimets at CERN ad Fermilab. The extrapolatio to x = 0 has bee carefully studied i Refs. [5, 6] ad the correspodig systematic ucertaities have bee estimated. At Q 2 values above those curretly accessible at Jefferso Lab (Q 2 > 5 (GeV/c) 2 ), the large-x extrapolatio was also studied by a compariso of various models ad the correspodig systematic ucertaities were icluded i the momet aalyses. Ulike previous studies, which have focussed o the x depedece, ad have igored or uderestimated the uclear correctios, we extract the lowest few momets of the eutro F2 takig ito accout the systematic ucertaities associated with our kowledge of uclear effects i the deutero. Furthermore, this aalysis directly uses the leadig twist parts of the momets extracted from the F p 2 ad F 2 D data, which allows a more straightforward implemetatio of the impulse approximatio. The geeral framework for aalyzig the data is preseted i Sec. 2, where we discuss various models of ucleo light-coe mometum distributios i the deutero. The extracted leadig twist eutro momets are preseted i Sec. 3, where we also discuss correctios to the simple oe-dimesioal covolutio formalism, i the form of ucleo off-shell effects, as well as uclear shadowig ad meso exchage currets. The results for the extracted momets are compared with those obtaied from several global parto distributio fuctios. Moreover we costruct the osiglet, isovector p combiatio of the leadig twist momets ad compare our results with available lattice oes. 3

As a iterestig applicatio of our extractio, i Sec. 4 we idetify the isospi depedece of the higher twist cotributios to the momets, ad speculate about possible iterpretatios of our fidigs withi a diquark picture. Fially, i Sec. 5 we summarize our results ad briefly outlie some future directios. 2 Neutro structure ad uclear effects i the deutero I the absece of free eutro targets, iformatio o the structure fuctio of the eutro is typically extracted from measuremets usig deuteros, together with kowledge of the correspodig proto structure fuctio. Sice the deutero is a ucleus, its ucleo costituets are boud, albeit weakly, ad i geeral the structure of the boud ucleos will differ from the free ucleo structure. I this sectio we describe the various uclear correctios which eter i the extractio of the free eutro F 2 structure fuctio, ad i particular its momets, from the structure fuctio of the deutero. We begi with the stadard uclear impulse approximatio, i which the virtual photo scatters icoheretly from idividual ucleos i the deutero, ad discuss several implemetatios of the covolutio formula. Followig this we describe correctios to covolutio, associated with ucleo off-shell effects, uclear shadowig, ad meso exchage currets. 2.1 Impulse approximatio For moderate to large values of x (x 0.1) iclusive scatterig off a deutero ca be described, usig the optical theorem, i terms of the amplitudes γ N γ N ad N d N d, where N refers to the fact that the ucleo is off its mass or eergy shell. I this approximatio iteractios betwee the itermediate states of the γ N ad N d subprocesses are eglected. Formally, the deutero F 2 structure fuctio ca the be writte as ] F2 D (x,q2 ) = d 4 p Tr [Ŵ(p,q) Â(p,p D ), (3) where p, p D ad q are the virtual ucleo, deutero ad virtual photo four-mometa, respectively, the Dirac matrix Ŵ(p,q) = I W 0 + p W 1 + q W 2 is the trucated (offshell) ucleo structure fuctio, ad Â(p,p D) = IA 0 +γ α A α 1 is the relativistic deutero spectral fuctio. The amplitudes A 0 ad A α 1 deped o the deutero wave fuctios. For a o-shell ucleo, the F 2 structure fuctio is proportioal to M W 0 +M 2 W 1 + p q W 2. I geeral the presece of several Dirac structures ad correspodig structure fuctios W i i the trace i the itegrad i Eq. (3) complicates the relatio betwee the deutero ad ucleo structure fuctios. There are various ways which the formal expressio i Eq. (3) ca be simplified, however, as we ow discuss. 2.1.1 Norelativistic expasio The formal expressio for the deutero structure fuctio F2 D i Eq. (3) is Loretzivariat ad ca be evaluated i ay referece frame. Without loss of geerality, how- 4

ever, oe ca choose to work i the deutero rest frame, ad take the z-axis such that q z = q. Performig a orelativistic reductio of the Dirac traces i Eq. (3), oe ca show that to order p 2 /M 2 i the Bjorke limit the product of the Dirac structures reduces to a sigle structure proportioal to a geeralized (off-shell) ucleo structure fuctio [4], F2 D (x,q 2 ) = d 3 p ( (2π) 3 1 + p ) z Ψ D (p) 2 F2 N M ( x z,q2,p 2). (4) Here F2 N = 1 2 (F p 2 +F 2 ) is the isoscalar ucleo structure fuctio, z = (p 0+p z )/M is the light-coe mometum fractio of the deutero carried by the iteractig ucleo, where p 0 = M + ε D p 2 /(2M) is the eergy of the struck ucleo, ad ε D = 2.22 MeV is the deutero bidig eergy. The deutero wave fuctio Ψ D cotais the usual orelativistic S- ad D-wave compoets, ψ 0 ad ψ 2, ormalized such that 1 2π 2 Ψ D(p) 2 = ψ 2 0( p ) + ψ 2 2( p ), (5) 1 2π 2 0 d p p 2 Ψ D (p) 2 = 1. (6) Note that i additio to the Bjorke variable ad Q 2, the ucleo structure fuctio i (4) also depeds o the ucleo off-shelless p 2 p 2 0 p2. If oe eglects the depedece of F2 N o p 2, the Eq. (4) leads to the simple oe-dimesioal covolutio formula i terms of the free ucleo structure fuctio [7, 8], F D(cov) MD 2 (x,q 2 /M ) = dz fr D (z) F 2 N x ( x z,q2), (7) where M D is the deutero mass, ad fr(z) D is the orelativistic ucleo (light-coe) mometum distributio i the deutero, d fr D 3 (z) = p ( (2π) 3 1 + p ) ( z Ψ D (p) 2 δ z p ) 0 + p z. (8) M M Physically the fuctio fr D (z) describes the Fermi motio ad bidig of ucleos i the deutero. Note that some implemetatios of the covolutio formula (7) do ot iclude the flux factor (1 + p z /M) [9]. However, sice this appears at O(p/M), cosistecy to order p 2 /M 2 requires its iclusio [4]. The wave fuctio ormalizatio (6) esures that MD /M 0 dz f D r(z) = 1, (9) which guaratees that (for the valece compoet) the covolutio term aloe preserves the baryo umber of the deutero. 5

2.1.2 Relativistic effects Because the deutero is a weakly boud ucleus, i most calculatios of deutero structure the orelativistic approximatio is adequate. I the large-x regio, however, oe is more sesitive to the large-mometum compoets of the deutero wave fuctio. I this case relativistic effects are expected to become more importat. Sice the distributio fuctio fr D i Eq. (8) was obtaied i the orelativistic approximatio, it applies i the regio z 1 1. Outside of this regio relativistic correctios may be importat. I order to test the sesitivity to these effects we also cosider a relativistic versio of the distributio i Eq. (8) which is calculated with the same orelativistic wave fuctio but usig relativistic kiematics, amely p 0 = M D E p, where E p = M 2 + p 2 is the eergy of the spectator ucleo. Futhermore, oe also replaces the ratio p z /M i the flux factor with p z /E p. While most deutero wave fuctios calculated to date have bee orelativistic, there have bee some models which have treated the two-ucleo system relativistically [10, 11, 12]. I order to cosistetly utilize the relativistic wave fuctios, oe eeds to work with the full expressio for the structure fuctio i Eq. (3) directly. This ecessarily ivolves modelig the off-shell ucleo structure fuctios, which itroduces additioal model depedece i the calculatio [13]. Oe ca show, however, that the exact expressio i Eq. (3) ca be reduced to a o-shell covolutio part, similar to that i Eq. (7), plus uclear off-shell ad relativistic correctios [14, 15], F D 2 (x,q2 ) = F D(cov) 2 (x,q 2 ) + δ (off) F D 2 (x,q2 ). (10) The covolutio compoet here ow ivolves a relativistic ucleo light-coe distributio fuctio [14], pmax frel D (z) = zm dp 2 E p D 16π 2 Ψ D (p) 2, (11) p 0 where E p = M 2 + p 2 = (M 2 p 2 + M 2 D )/2M D is the eergy of the (o-shell) ucleo spectator, ad the upper limit of the ucleo virtuality p 2 i the itegratio is give by p 2 max = zm2 D /2 zm2 /(2 z). I additio to the usual S- ad D-state wave fuctios u ad w, the wave fuctio Ψ D (p) also icludes the relativistic P-states [11, 12] 1 2π 2 Ψ D (p) 2 = ψ0 2 ( p ) + ψ2 2 ( p ) + ψ2 1t ( p ) + ψ2 1s ( p ), (12) where ψ 1s ad ψ 1t correspod to the siglet ad triplet P-states, respectively. Note that the fuctio frel D also icludes the relativistic flux factor, (E p /p 0 ) z [9], ad therefore satisfies the same ormalizatio coditio as i Eq. (9), dzfrel D = 1. Oe should poit out, however, that the separatio of F2 D ito the two terms i Eq. (10) is ot uique (the sum, of course, is uique). However, the choice (11) is the most atural oe, ad correspods to that which has bee followed i may pheomeological treatmets. 6

2.1.3 Light-coe formulatio A alterative to the covariat (istat-form) formulatios of scatterig from the deutero is the Hamiltoia light-coe (LC) approach. Here the itermediate state ucleos are o-mass-shell, but off-eergy-shell, so that oe avoids the explicit p 2 depedece i the boud ucleo structure fuctio i Eq. (4). Because the ucleos are o-mass-shell, the derivatio of the oe-dimesioal covolutio from Eq. (3) is trivial sice each of the three terms i the trace factor become proportioal [13]. The LC ucleo mometum distributio i the deutero is give by [16] flc(z) D = d p p M 2 + p 2 Ψ LC D (p) 2, (13) p mi where p mi = M z 1 / z(2 z), ad Ψ LC D is the LC deutero wave fuctio. Sice most of the two-ucleo pheomeology has traditioally bee formulated withi the istat-form approach, i practice for Ψ LC D oe has usually employed the S- ad D-state wave fuctios as calculated i orelativistic potetial models [see Eq. (5)]. Note that i the limit where the ucleos are free (i.e. ot boud), the deutero ca be described i terms of a NN Fock state compoet. Iclusio of bidig requires the additio of N N +meso Fock states. Techically, bidig or off-shell effects i the LC approach ca be icluded by costructig a wave fuctio which explicitly depeds o a additioal light-like four-vector related to the orietatio of the ull-plae [17]. This itroduces a umber of additioal wave fuctios, which would eed to be icluded i flc D. Recetly calculatios of deutero LC wave fuctios have bee performed withi this framework, ad some of the additioal wave fuctios foud to be importat at large mometum. While it is iterestig to explore the pheomeology usig these wave fuctios, i the preset study we shall compare the results of the LC distributio flc D usig the orelativistic S- ad D-state wave fuctios as used i the literature. I the ext sectio we compare the various light-coe mometum distributios discussed above, for differet model deutero wave fuctio models. 2.2 Nucleo mometum distributios Withi the covolutio approximatio there are two sources of uclear model depedece i the calculatio of the deutero structure fuctio. Firstly, as discussed above, differet approaches are used to formulate scatterig from the two-ucleo system orelativistic, relativistic, or light-coe. I additio, withi each approach there is depedece o the deutero wave fuctio or N N potetial. It is beyod the scope of this work to critically aalyze the merits of the various models; istead we survey a reasoably large sample of models ad estimate the systematic ucertaity from the spread i their predictios. I practice the deutero wave fuctios are determied by fittig N N scatterig data. Sice most of the data are at low mometum, the wave fuctios are relatively well costraied at small p, but the large-p tail of the wave fuctio is more model depedet. I terms of the light-coe variable z, the large-p regio correspods to 7

10-1 10-1 f D (y)y 7 10-2 f D (y)y 7 10-2 10-3 10-3 10-4 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 y 10-4 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 y Figure 1: Nucleo light-coe mometum distributio f D weighted by y 7, where y = z/2, for various models: (a) orelativistic distributio (8) (solid), distributio (8) with relativistic kiematics (dashed), relativistic model (11) (dotted), light-coe distributio (13), all usig the Paris wave fuctio [18]; (b) relativistic distributio (11) usig the Bo [19] (solid), Paris [18] (dashed), ad relativistic Gross wave fuctios from Refs. [11] (dotted) ad [12] (dot-dashed). large values of z. To emphasize the large-z regio of f D we plot i Fig. 1 the various ucleo mometum distributios, for differet wave fuctios, weighted by y 7, where y z/2. Note that the mometum itegratio i each of the curves has bee limited to p = 0.5 GeV (see below). As expected the model depedece is relatively weak for moderate y, ad oly becomes sigificat for y 0.6. Ucertaities due to the uclear models ca also be quatified i terms of momets. This is actually more relevat for the preset aalysis sice ultimately we extract the momets of the eutro F2 structure fuctio. If we defie the -th momet of the distributio f D (z) by F D = MD /M 0 dz z 1 f D (z), (14) the clearly the ucertaities i the large-z part of the distributio will be reflected i the ucertaities of the momets of f D at large-, which are illustrated i Fig. 2. As metioed above, the differeces betwee the mometum distributios at large z arise because the deutero wave fuctio is ot well costraied at large mometum. I fact, N N scatterig data which have bee used to costrai the deutero wave fuctio typically exted oly to p 0.4 GeV. At large mometum (p 0.3 0.4 GeV) the iter-ucleo distaces probed becomes smaller tha the ucleo radius, so that here a descriptio i terms of ucleo degrees of freedom aloe may be questioable. To estimate the ucertaity arisig from the large-p tail of the deutero wave fuctio, we 8

1.2 1.2 1.15 1.15 1.1 1.1 1.05 1.05 I D 1 I D 1 0.95 0.95 0.9 0.9 0.85 0.85 0.8 2 4 6 8 10 0.8 2 4 6 8 10 Figure 2: Momets F D ( = 2 10) of the ucleo light-coe mometum distributio fd (z) i Eq. (11) (a) for differet deutero wave fuctios: BoCD [20] (solid), Nijmege 93 [21] (dashed), Argoe v 18 [22] (dotted), Gross [11] (dot-dashed); (b) for the Gross wave fuctio [12] with the mometum itegratio restricted to p max = 0.5 GeV (solid) ad 1 GeV (dashed). examie the effect of trucatig the mometum itegratio at p = p max. The momets F D of the distributio f D (z) with p max = 0.5 ad 1 GeV are plotted i Fig. 2(b), where the distributios have bee reormalized to satisfy Eq. (9). For low momets ( = 2,4) the effect is egligible, while for the = 10 momet the effect of the p max = 0.5 (1) GeV cuts is 6% (1%). We estimate that the ucertaity itroduced by our poor kowledge of the large-p compoets of the deutero wave fuctio is of the same order of magitude as the wave fuctio model depedece. The maximum deviatio amog the momets obtaied with differet deutero wave fuctios will be take as a estimate of the systematic error due to the wave fuctio model depedece. A idepedet test of the ucleo distributio fuctios may be made by cosiderig quasi-elastic scatterig. Withi the impulse approximatio ad eglectig fiite-q 2 effects, the quasi-elastic cotributio to the deutero structure fuctio, F D(QE) 2, ca be writte i factorized form, i terms of the ucleo distributio fuctio f D ad the ucleo elastic form factors, F D(QE) 2 (x D,Q 2 ) 1 { G p 2 2 x Df D E (x D ) (Q2 ) + G E 2(Q2 ) + τ(g p 2 M (Q2 ) + G M 2(Q2 )) 1 + τ (15) where x D Q 2 /2M D ν x/2. Thus, oe ca defie pheomeological ucleo distributio momets F D (Q 2 ) i terms of the momets M D(QE) (Q 2 ) of the deutero quasi-elastic }, 9

structure fuctio F D(QE) 2 ad the ucleo form factors, F D (Q 2 ) = (1 + τ)m D(QE) (Q 2 ) G p 2 E (Q2 ) + G E 2(Q2 ) + τ(g p 2 M (Q2 ) + G M 2(Q2 )). (16) We expect that F D (Q 2 ) F D at large eough Q 2. Note that for the momet M D(QE) (Q 2 ) we use the Nachtma defiitio [23] to elimiate the fiite-q 2 effects due to the ozero target mass (i.e. the deutero mass), M D(QE) (Q 2 ) where r D = 1 0 ξd +1 dx D x 3 F D(QE) 2 (x D,Q 2 ) 3 + 3( + 1)r D + ( + 2)rD 2 D ( + 2)( + 3) 1 + 4M 2 D x2 D /Q2 ad ξ D = 2x D /(1 + r D )., (17) 1.2 1.2 1.15 1.15 1.1 1.1 I ~ 2 D 1.05 I ~ 6 D 1.05 1 1 0.95 0.95 0.9 1 10 10 2 0.9 1 10 10 2 Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 Figure 3: Momets F 2 D(Q2 ) (a) ad F 6 D(Q2 ) (b) extracted from the quasi-elastic peak momets defied i Eq. (16) (triagles). The hatched areas represet the correspodig momets F2 D ad F6 D calculated from Eq. (14), icludig the uclear model ucertaities show i Fig. 2. The deutero quasi-elastic cotributio to the structure fuctio F D(QE) 2 has bee calculated usig the approach developed i Ref. [24], which icludes fial state iteractios ad has bee tested agaist the recet deutero experimetal data of Ref. [6]. As a illustratio, the extracted F D 2 (Q2 ) ad F D 6 (Q2 ) momets are show i Fig. 3 as a fuctio of Q 2, ad compared with the correspodig momets F D 2 ad FD 6 calculated from the distributio fuctio f D (z) i Eq. (14). The results agree well withi the uclear model ucertaities. 10

2.3 Off-shell effects The expressios for the deutero structure fuctio thus far have bee i the framework of the oe-dimesioal covolutio formula, i which ucleo off-shell, relativistic, ad other effects beyod the impulse approximatio are eglected. While this costitutes a reasoable approximatio at moderate x, or for low momets, at large (ad very small) x values these effects will become more importat. Withi the orelativistic approach discussed i Sec. 2.1.1, the fact that the structure of boud ucleos ca differ from that of free ucleos is icorporated i the depedece of the boud ucleo structure fuctio i Eq. (4) o the ucleo virtuality, or offshelless, p 2. Sice the deutero is treated as a orelativistic system, the off-shelless parameter v (p 2 M 2 )/M 2 is o average a small umber. Expadig the leadig twist structure fuctio i v ad keepig oly the first order correctio i v, oe fids F N 2 (x,q2,p 2 ) = F N 2 (x,q2 )(1 + δf 2 (x)v), (18) δf 2 = l F N 2 / l p2, (19) where the first term o the right i Eq. (18) is the structure fuctio of the o-mass-shell ucleo, ad the derivative is evaluated at p 2 = M 2. The deutero structure fuctio F2 D ca the be writte as i Eq. (10), with d δ (off) F2 D (x,q 2 3 p ( ) = (2π) 3 1 + p ) z Ψ D (p) 2 δf 2 (x/z,q 2 ) v F2 N (x/z,q 2 ).(20) M The off-shell correctio δf 2 was studied pheomeologically i Ref. [25] by aalyzig data o the ratios of structure fuctios of differet uclei (EMC ratios). It was foud that to a good approximatio δf 2 was idepedet of Q 2, ad could be parameterized as δf 2 = C N (x x 1 )(x x 0 )(1 + x 0 x), (21) with the best fit parameters C N = 8.1 ± 0.3 stat ± 0.5 syst, x 1 = 0.05, ad x 0 = 0.448 ± 0.005 stat ± 0.007 syst. Futhermore, the off-shell correctio to the structure fuctio F 3 is give by the same fuctio δf 2, which is cosistet with the ormalizatio of uclear valece quark distributio (for more details see Ref. [25]). Expressig the off-shell correctio to F2 D i terms of momets, oe ca chage the order of the p ad x itegratios to obtai M D (Q2 ) = F D MN (Q2 ) + δc δm N (Q2 ), (22) where the first term o the right had side is the usual covolutio term ad the secod is the off-shell correctio, with d 3 p ( δc = (2π) 3 1 + p ) z z 1 v Ψ D (p) 2 (23) M ad δm N (Q2 ) = 1 0 dxx 2 δf 2 (x) F N 2 (x,q2 ). (24) 11

Writig the ucleo momet i terms of the proto ad eutro momets explicitly, M N = (Mp + M )/2, oe ca solve Eq. (22) to express the eutro momet i terms of the experimetal proto ad deutero momets, ad the theoretical F D ad off-shell correctios, ( ) M (Q2 ) = 2M D 1 (off) (Q2 ) F D M p (Q2 ), (25) where the off-shell ratio (off) ad δ (off) is defied as (off) (Q 2 ) = is give by 1 dx x 2 δ (off) (x) F2 D (x,q 2 ), (26) M D δ (off) (x,q 2 ) = δ(off) F D 2 (x,q2 ) F D 2 (x,q2 ). (27) From Eqs. (20), (23) ad (24) oe ca also write the momet correctio as (off) (Q 2 ) = δc δm N 2M D. (28) For the relativistic model of Sec. 2.1.2, the off-shell correctio δ (off) F2 D receives two cotributios: oe which depeds o the relativistic P-state wave fuctios ψ 1s ad ψ 1t, ad oe which is give by the off-shell part of the trucated ucleo structure fuctios W 0,1,2. The former was foud to be of the order 1% ad egative over the whole x rage, while the latter was somewhat smaller i magitude. The total off-shell ratio δ (off) was parameterized at Q 2 5 GeV 2 as [14] δ (off) (x,q 2 ) = a 0 (1 + a 1 x a 2 )(1 (a 3 x a 4 ) a 5 ), (29) with the parameters a 0 = 0.014, a 1 = 3.0, a 2 = 20.0, a 3 = 1.067, a 4 = 1.5, ad a 5 = 18.0. The Q 2 depedece of this ratio was foud to be weak. Note that the estimate of δ (off) i Ref. [14] was based o the set of wave fuctios i Ref. [11] derived from a NN potetial with a pseudoscalar πnn couplig, which are kow to give larger P-state cotributios tha more moder treatmets [10, 12]. The correctio δ (off) i Eq. (29) should therefore serve as a upper limit o the size of the off-shell correctio i this model. The x depedece of the correctios δ (off) i the models of Refs. [14, 25], evaluated usig the leadig twist parameterizatio [6] of F2 D at Q2 = 10 GeV 2, is displayed i Fig. 4. Both of the correctios are aroud 1 2% for x 0.6 0.7, but icrease sigificatly i magitude for x 0.8. The correctio at large x is egative i both cases, ad is are listed i Table 1 for = 2 10. The magitude of the correctio for the lowest momet is similar i the two models, but is of opposite sig. For higher, the relative importace of the especially dramatic for the model of Ref. [25]. The resultig momets (off) 12

0.02 0-0.02-0.04 δ (off) -0.06-0.08-0.1-0.12-0.14 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x Figure 4: Nucleo off-shell correctio δ (off) (x, Q 2 ) as a fuctio of x at Q 2 = 10 GeV 2, i the orelativistic [25] (dashed) ad relativistic [14] (solid) models. Table 1: Momets of the o-covolutio correctios to the deutero F2 D structure fuctio for = 2 10: the first two colums give off-shell correctios (off) i the relativistic [14] ad orelativistic [25] models, respectively; the third ad fourth colums represet the cotributios from meso exchage currets ad uclear shadowig [3], respectively. (off) [14] (off) [25] (MEC) [3] (shad) [3] 2-5.67 10 3 5.08 10 3 2.6 10 4-1.0 10 4 4-1.21 10 2-3.34 10 3 7.6 10 6 7.6 10 6 6-1.39 10 2-1.68 10 2 8.4 10 7 1.5 10 6 8-1.46 10 2-3.01 10 2 1.6 10 7 4.1 10 7 10-1.53 10 2-3.97 10 2 4.4 10 8 1.4 10 7 off-shell correctio icreases. Because of the steep decrease of δ (off) i the orelativistic model [25], the magitude of the correctio (off) at = 10 is about 2-3 times larger tha the oe i the model of Ref. [14]. Note that part of the off-shell correctios i the model of Ref. [14] are icluded i the relativistic distributio fuctio (11), which i the model of Ref. [25] are cotaied oly i the o-covolutio term. Therefore the defiitio of off-shell correctios ca differ betwee various models, ad oe should esure that off-shell effects i geeral are calculated self-cosistetly. The differece betwee the two models provides a estimate of the systematic error associated with the poor kowledge of the chage i the ucleo 13

structure whe the ucleo is boud i deutero (see also Ref. [26]). 2.4 Small-x effects: shadowig ad meso exchage currets For completeess, i this sectio we briefly discuss the uclear correctios i the deutero at small values of x. From the defiitio of the momets of the F 2 structure fuctio, it is clear that small-x effects will mostly be relevat for the lowest momet, ad will be egligible for 4. At small values of x, the coherece legth of a virtual q q fluctuatio of the photo probe i the laboratory frame is 1/Mx, which for x 0.1 becomes larger tha the iter-ucleo separatio. Cosequetly at small x the scatterig from the deutero ca take place via double scatterig from both ucleos, which gives rise to the pheomeo of uclear shadowig (see e.g. Refs. [4, 2]). At low Q 2 the q q fluctuatios are usually represeted i the form of the vector mesos, ρ, ω ad φ. Their cotributios scale as 1/Q 2, however, ad are formally of higher twist. Leadig twist shadowig cotributios arise from diffractive scatterig of the virtual photo from the boud ucleos, which is described through Pomero exchage, as well as cotributios i which the exchaged quata are mesos (most otably the pio). 0.002 0.001 0 δ (P,MEC) -0.001-0.002-0.003-0.004-0.005 0 0.1 0.2 0.3 0.4 0.5 x Figure 5: Shadowig (Pomero exchage) ad meso exchage curret correctios δ (P) (x, Q 2 ) (solid) ad δ (MEC) (x, Q 2 ) (dotted) as a fuctio of x at Q 2 = 10 GeV 2 [3]. Both the Pomero (P) exchage ad meso exchage curret (MEC) cotributios to the deutero F D 2 structure fuctio have bee estimated by several authors [1, 2, 3, 4, 25, 27], ad eter as additive correctios to F D 2, F D 2 (x,q2 ) F D 2 (x,q2 ) + δ (P,MEC) F D 2 (x,q2 ). (30) 14

The Pomero shadowig cotributio δ (P) F2 D is egative, while the meso (pio) exchage cotributio δ (MEC) F2 D is positive ad hece correspods to a atishadowig effect. To estimate their effect o the momet aalysis, we use the results from Ref. [3], which ca be coveietly parameterized at Q 2 5 GeV 2 as δ (P) F D 2 = b 0 x b 1 (1 x) b 2 (1 + b 3 x b 4 ), (31) for the P-exchage, where b 0 = 0.003, b 1 = 0.13, b 2 = 5.0, b 3 = 2.2 ad b 4 = 0.4, ad δ (MEC) F D 2 = c 0 x c 1 (1 x) c 2 (32) for meso exchage, with c 0 = 0.002, c 1 = 0.03 ad c 2 = 6.0. As a fractio of the empirical leadig twist F2 D, the ratios δ (P,MEC) (x,q 2 ) = δ(p,mec) F D 2 (x,q2 ) F D 2 (x,q2 ) (33) are plotted i Fig. 5. The correspodig momets, defied i aalogy with Eq. (26), (P,MEC) (Q 2 ) = 1 M D dx x 2 δ (P,MEC) (x,q 2 ) F2 D (x,q2 ), (34) are listed i Table 1, where it is clear that these correctios are egligible for 4. 3 Momets of the eutro structure fuctio I this sectio we preset the extractio of the lowest few momets of the eutro structure fuctio F2, takig ito accout the systematic ucertaities associated with our icomplete kowledge of uclear effects i the deutero. We shall focus o the leadig twist (LT) part of the structure fuctio, ad make use of the previous determiatio of the LT momets of the proto ad deutero structure fuctios obtaied by the twist aalyses of Refs. [5, 6]. Higher twist cotributios to the momets, i the form of subleadig 1/Q 2 power correctios, will be discussed i Sec. 4. Moreover, sice the twist extractio performed i Refs. [5, 6] is based o the Nachtma defiitio of the momets [23], our LT terms are determied without ay a priori kowledge of the x- shape of parto distributio fuctios (PDFs). The momets of the LT F 2 structure fuctio of the eutro are related to those of proto ad deutero by M(Q 2 ) = 2M D (Q 2 ) 1 (off) (Q 2 ) (P,MEC) (Q 2 ) F D M(Q p 2 ), (35) where F D /[1 (off) (Q 2 ) (P,MEC) (Q 2 )] represets the uclear correctio factor, or i other words the EMC effect of the deutero i momet space. From the values reported i Fig. 2 ad i Table 1, the EMC effect o the secod momet is expected to be below 1% [see also Fig. 3(a)]. 15

Usig the uclear igrediets discussed i the previous Sectio, the extracted LT momets M (Q2 ) are listed i Table 2 ad show i Fig. 6. The cetral values refer to the case of the relativistic distributio (11) with the deutero wave fuctio from Ref. [12]. For compariso we also plot the correspodig proto (stars) ad deutero (ope crosses) LT momets. It turs out that the eutro LT momets with 8 ca be determied with quite good precisio, amely 18% (statistical) ad 30% (systematic). Differet cotributios to the total systematic error are show i Fig. 7. It ca be see that all the ucertaities related to the uclear correctio factor steeply rise with, so that the total errors due to uclear effects become comparable with the experimetal oes aroud = 10. Thus, a improvemet i the measuremets of the proto ad deutero structure fuctios at large x would be beeficial ad also the momets at = 12 may be still reasoably accessible. I Fig. 8 our extracted eutro momets are compared with the correspodig momets calculated from various PDF sets available i the literature [28, 29, 30, 31]. For > 2 all PDF sets predict sigificatly lower momets at low Q 2. This differece is due to Sudakov effects at large-x; ideed, as was firstly ivestigated i Ref. [32], soft gluo resummatio plays a relevat role for > 2 already at Q 2 = 3 4 GeV 2. While resummatio effects are take ito accout i the aalysis of Refs. [5, 6], all the PDF sets of Refs. [28, 29, 30, 31] are based o ext-to-leadig (NLO) or NNLO treatmets. At higher Q 2 the PDF predictios are cosistet with our results for 10 withi the quoted statistical errors. For M 2 however all PDFs give somehow smaller values at large Q 2, though withi the statistical ucertaities. The large- behavior of the momets is also closely related to the asymptotic x 1 behavior of the structure fuctios. Ideed, ispectig the defiitio of the momets, Eq. (2), we coclude that at 1 the factor x 2 i the itegrad suppresses the cotributio from the regio of small x, leavig the momet effectively saturated by a regio close to x = 1. Furthermore, assumig that the structure fuctios F p 2 ad F 2 vaish as (1 x) β with the same expoet for the proto ad eutro (as oe would expect from perturbative QCD) [33], it is possible to derive the asymptotic relatio betwee the ratios of the momets ad the structure fuctios M lim M p F2 = lim (x). (36) x 1 (x) This relatio ca be useful i obtaiig some isights ito the large-x asymptotics of the d/u ratio from the large- limit of the eutro ad proto momets. Note that most of the PDF fits [28, 29, 30, 31] are characterized by the limitig value lim x 1 F 2 /F p 2 = 1/4, which is expected if soft, operturbative effects domiate the scatterig at large x (amely, if the eergy of the spectator uu system which accompaies the d quark is icreased relative to the ud system which is relevat for the u quark [34]). O the other had, argumets based o helicity coservatio i pqcd, whe combied with a uperturbed spi-flavor symmetric wave fuctio, suggest a asymptotic value of 3/7 [35]. The ratios of eutro to proto LT momets are show i Fig. 9 for several, where they are also compared to various PDF parameterizatios [28, 29, 30, 31]. The ratio F p 2 16

Table 2: Extracted leadig twist momets of the eutro F 2 structure fuctios, M (Q 2 ), for = 2, 4, 6, 8, show i Fig. 6, together with statistical ad systematic ucertaities. Q 2 [(GeV/c) 2 ] M2 (Q2 ) 10 1 M4 (Q2 ) 10 2 M6 (Q2 ) 10 3 M8 (Q2 ) 10 3 1.025 1.37 ± 0.06 ± 0.12 1.95 ± 0.20 ± 0.11 8.06 ± 0.45 ± 1.15 5.51 ± 0.98 ± 1.49 1.075 1.36 ± 0.06 ± 0.12 1.88 ± 0.19 ± 0.11 7.35 ± 0.41 ± 1.04 4.59 ± 0.81 ± 1.24 1.125 1.36 ± 0.06 ± 0.12 1.82 ± 0.18 ± 0.11 6.78 ± 0.38 ± 0.96 3.95 ± 0.70 ± 1.07 1.175 1.36 ± 0.06 ± 0.12 1.76 ± 0.18 ± 0.10 6.33 ± 0.35 ± 0.90 3.48 ± 0.62 ± 0.94 1.225 1.36 ± 0.05 ± 0.11 1.71 ± 0.17 ± 0.10 5.95 ± 0.33 ± 0.85 3.13 ± 0.56 ± 0.85 1.275 1.35 ± 0.05 ± 0.11 1.67 ± 0.17 ± 0.10 5.63 ± 0.31 ± 0.80 2.85 ± 0.51 ± 0.77 1.325 1.35 ± 0.05 ± 0.11 1.63 ± 0.16 ± 0.09 5.36 ± 0.30 ± 0.76 2.63 ± 0.47 ± 0.71 1.375 1.35 ± 0.05 ± 0.11 1.59 ± 0.16 ± 0.09 5.13 ± 0.28 ± 0.73 2.44 ± 0.43 ± 0.66 1.425 1.34 ± 0.05 ± 0.11 1.56 ± 0.16 ± 0.09 4.93 ± 0.27 ± 0.70 2.29 ± 0.41 ± 0.62 1.475 1.34 ± 0.05 ± 0.11 1.53 ± 0.15 ± 0.09 4.74 ± 0.26 ± 0.68 2.16 ± 0.38 ± 0.58 1.525 1.34 ± 0.05 ± 0.11 1.50 ± 0.15 ± 0.09 4.58 ± 0.25 ± 0.65 2.04 ± 0.36 ± 0.55 1.575 1.34 ± 0.05 ± 0.11 1.48 ± 0.15 ± 0.09 4.44 ± 0.25 ± 0.63 1.94 ± 0.35 ± 0.53 1.625 1.33 ± 0.05 ± 0.11 1.45 ± 0.15 ± 0.08 4.31 ± 0.24 ± 0.61 1.86 ± 0.33 ± 0.50 1.675 1.33 ± 0.05 ± 0.11 1.43 ± 0.14 ± 0.08 4.19 ± 0.23 ± 0.60 1.78 ± 0.32 ± 0.48 1.725 1.33 ± 0.05 ± 0.11 1.41 ± 0.14 ± 0.08 4.08 ± 0.23 ± 0.58 1.71 ± 0.30 ± 0.46 1.775 1.32 ± 0.05 ± 0.11 1.39 ± 0.14 ± 0.08 3.99 ± 0.22 ± 0.57 1.65 ± 0.29 ± 0.45 1.825 1.32 ± 0.05 ± 0.11 1.37 ± 0.14 ± 0.08 3.90 ± 0.22 ± 0.55 1.60 ± 0.28 ± 0.43 1.875 1.32 ± 0.05 ± 0.11 1.36 ± 0.14 ± 0.08 3.81 ± 0.21 ± 0.54 1.55 ± 0.27 ± 0.42 1.925 1.32 ± 0.05 ± 0.11 1.34 ± 0.13 ± 0.08 3.73 ± 0.21 ± 0.53 1.50 ± 0.27 ± 0.41 1.975 1.32 ± 0.05 ± 0.11 1.33 ± 0.13 ± 0.08 3.66 ± 0.20 ± 0.52 1.46 ± 0.26 ± 0.40 2.025 1.31 ± 0.05 ± 0.11 1.31 ± 0.13 ± 0.08 3.59 ± 0.20 ± 0.51 1.42 ± 0.25 ± 0.38 2.075 1.31 ± 0.05 ± 0.11 1.30 ± 0.13 ± 0.08 3.53 ± 0.20 ± 0.50 1.38 ± 0.25 ± 0.38 2.125 1.31 ± 0.05 ± 0.11 1.29 ± 0.13 ± 0.07 3.47 ± 0.19 ± 0.49 1.35 ± 0.24 ± 0.37 2.175 1.31 ± 0.05 ± 0.10 1.27 ± 0.13 ± 0.07 3.42 ± 0.19 ± 0.49 1.32 ± 0.23 ± 0.36 2.225 1.31 ± 0.05 ± 0.10 1.26 ± 0.13 ± 0.07 3.36 ± 0.19 ± 0.48 1.29 ± 0.23 ± 0.35 2.275 1.31 ± 0.05 ± 0.10 1.25 ± 0.13 ± 0.07 3.32 ± 0.18 ± 0.47 1.27 ± 0.22 ± 0.34 2.325 1.31 ± 0.05 ± 0.10 1.24 ± 0.12 ± 0.07 3.28 ± 0.18 ± 0.47 1.25 ± 0.22 ± 0.34 2.375 1.30 ± 0.05 ± 0.10 1.23 ± 0.12 ± 0.07 3.24 ± 0.18 ± 0.46 1.23 ± 0.22 ± 0.33 2.425 1.30 ± 0.05 ± 0.10 1.22 ± 0.12 ± 0.07 3.20 ± 0.18 ± 0.46 1.21 ± 0.21 ± 0.33 2.475 1.30 ± 0.05 ± 0.10 1.22 ± 0.12 ± 0.07 3.17 ± 0.18 ± 0.45 1.19 ± 0.21 ± 0.32 2.525 1.30 ± 0.05 ± 0.10 1.21 ± 0.12 ± 0.07 3.13 ± 0.17 ± 0.45 1.17 ± 0.21 ± 0.32 2.575 1.30 ± 0.05 ± 0.10 1.20 ± 0.12 ± 0.07 3.10 ± 0.17 ± 0.44 1.15 ± 0.21 ± 0.31 2.625 1.30 ± 0.05 ± 0.10 1.19 ± 0.12 ± 0.07 3.07 ± 0.17 ± 0.44 1.14 ± 0.20 ± 0.31 2.675 1.30 ± 0.05 ± 0.10 1.18 ± 0.12 ± 0.07 3.04 ± 0.17 ± 0.43 1.12 ± 0.20 ± 0.30 2.725 1.30 ± 0.05 ± 0.10 1.18 ± 0.12 ± 0.07 3.01 ± 0.17 ± 0.43 1.11 ± 0.20 ± 0.30 2.775 1.30 ± 0.05 ± 0.10 1.17 ± 0.12 ± 0.07 2.98 ± 0.17 ± 0.42 1.10 ± 0.19 ± 0.30 2.825 1.30 ± 0.05 ± 0.10 1.16 ± 0.12 ± 0.07 2.96 ± 0.16 ± 0.42 1.08 ± 0.19 ± 0.29 2.875 1.30 ± 0.05 ± 0.10 1.16 ± 0.12 ± 0.07 2.93 ± 0.16 ± 0.42 1.07 ± 0.19 ± 0.29 2.925 1.30 ± 0.05 ± 0.10 1.15 ± 0.12 ± 0.07 2.91 ± 0.16 ± 0.41 1.06 ± 0.19 ± 0.29 2.975 1.30 ± 0.05 ± 0.10 1.15 ± 0.12 ± 0.07 2.88 ± 0.16 ± 0.41 1.05 ± 0.19 ± 0.28 3.025 1.30 ± 0.05 ± 0.10 1.14 ± 0.11 ± 0.07 2.86 ± 0.16 ± 0.41 1.04 ± 0.18 ± 0.28 3.075 1.30 ± 0.05 ± 0.10 1.13 ± 0.11 ± 0.07 2.84 ± 0.16 ± 0.40 1.03 ± 0.18 ± 0.28 3.125 1.30 ± 0.05 ± 0.10 1.13 ± 0.11 ± 0.07 2.82 ± 0.16 ± 0.40 1.02 ± 0.18 ± 0.28 3.175 1.30 ± 0.05 ± 0.10 1.12 ± 0.11 ± 0.07 2.80 ± 0.16 ± 0.40 1.01 ± 0.18 ± 0.27 3.225 1.30 ± 0.05 ± 0.10 1.12 ± 0.11 ± 0.06 2.78 ± 0.15 ± 0.40 1.00 ± 0.18 ± 0.27 3.275 1.30 ± 0.05 ± 0.10 1.11 ± 0.11 ± 0.06 2.76 ± 0.15 ± 0.39 0.99 ± 0.18 ± 0.27 3.325 1.30 ± 0.05 ± 0.10 1.11 ± 0.11 ± 0.06 2.74 ± 0.15 ± 0.39 0.98 ± 0.17 ± 0.27 3.375 1.30 ± 0.05 ± 0.10 1.10 ± 0.11 ± 0.06 2.72 ± 0.15 ± 0.39 0.97 ± 0.17 ± 0.26 3.425 1.30 ± 0.05 ± 0.10 1.10 ± 0.11 ± 0.06 2.70 ± 0.15 ± 0.38 0.96 ± 0.17 ± 0.26 3.475 1.30 ± 0.05 ± 0.10 1.09 ± 0.11 ± 0.06 2.68 ± 0.15 ± 0.38 0.95 ± 0.17 ± 0.26 3.525 1.29 ± 0.05 ± 0.10 1.09 ± 0.11 ± 0.06 2.67 ± 0.15 ± 0.38 0.94 ± 0.17 ± 0.26 3.575 1.29 ± 0.05 ± 0.10 1.08 ± 0.11 ± 0.06 2.65 ± 0.15 ± 0.38 0.94 ± 0.17 ± 0.25 3.625 1.29 ± 0.05 ± 0.10 1.08 ± 0.11 ± 0.06 2.64 ± 0.15 ± 0.38 0.93 ± 0.17 ± 0.25 3.675 1.29 ± 0.05 ± 0.10 1.08 ± 0.11 ± 0.06 2.62 ± 0.15 ± 0.37 0.92 ± 0.16 ± 0.25 3.725 1.29 ± 0.05 ± 0.10 1.07 ± 0.11 ± 0.06 2.60 ± 0.14 ± 0.37 0.92 ± 0.16 ± 0.25 3.775 1.29 ± 0.05 ± 0.10 1.07 ± 0.11 ± 0.06 2.59 ± 0.14 ± 0.37 0.91 ± 0.16 ± 0.25 3.825 1.29 ± 0.05 ± 0.10 1.06 ± 0.11 ± 0.06 2.58 ± 0.14 ± 0.37 0.90 ± 0.16 ± 0.24 17

Q 2 [(GeV/c) 2 ] M2 (Q2 ) 10 1 M4 (Q2 ) 10 2 M6 (Q2 ) 10 3 M8 (Q2 ) 10 3 3.875 1.29 ± 0.05 ± 0.10 1.06 ± 0.11 ± 0.06 2.56 ± 0.14 ± 0.36 0.90 ± 0.16 ± 0.24 3.925 1.29 ± 0.05 ± 0.10 1.06 ± 0.11 ± 0.06 2.55 ± 0.14 ± 0.36 0.89 ± 0.16 ± 0.24 3.975 1.29 ± 0.05 ± 0.10 1.05 ± 0.11 ± 0.06 2.53 ± 0.14 ± 0.36 0.88 ± 0.16 ± 0.24 4.025 1.29 ± 0.05 ± 0.10 1.05 ± 0.11 ± 0.06 2.52 ± 0.14 ± 0.36 0.88 ± 0.16 ± 0.24 4.075 1.29 ± 0.05 ± 0.10 1.05 ± 0.11 ± 0.06 2.51 ± 0.14 ± 0.36 0.87 ± 0.15 ± 0.24 4.125 1.29 ± 0.05 ± 0.10 1.04 ± 0.10 ± 0.06 2.50 ± 0.14 ± 0.36 0.87 ± 0.15 ± 0.23 4.175 1.29 ± 0.05 ± 0.10 1.04 ± 0.10 ± 0.06 2.48 ± 0.14 ± 0.35 0.86 ± 0.15 ± 0.23 4.225 1.29 ± 0.05 ± 0.10 1.04 ± 0.10 ± 0.06 2.47 ± 0.14 ± 0.35 0.86 ± 0.15 ± 0.23 4.275 1.29 ± 0.05 ± 0.10 1.03 ± 0.10 ± 0.06 2.46 ± 0.14 ± 0.35 0.85 ± 0.15 ± 0.23 4.325 1.29 ± 0.05 ± 0.10 1.03 ± 0.10 ± 0.06 2.45 ± 0.14 ± 0.35 0.84 ± 0.15 ± 0.23 4.375 1.29 ± 0.05 ± 0.10 1.03 ± 0.10 ± 0.06 2.44 ± 0.14 ± 0.35 0.84 ± 0.15 ± 0.23 4.425 1.29 ± 0.05 ± 0.10 1.02 ± 0.10 ± 0.06 2.43 ± 0.13 ± 0.35 0.84 ± 0.15 ± 0.23 4.475 1.29 ± 0.05 ± 0.10 1.02 ± 0.10 ± 0.06 2.42 ± 0.13 ± 0.34 0.83 ± 0.15 ± 0.23 4.525 1.29 ± 0.05 ± 0.10 1.02 ± 0.10 ± 0.06 2.41 ± 0.13 ± 0.34 0.83 ± 0.15 ± 0.22 4.575 1.29 ± 0.05 ± 0.10 1.01 ± 0.10 ± 0.06 2.40 ± 0.13 ± 0.34 0.82 ± 0.15 ± 0.22 4.625 1.29 ± 0.05 ± 0.10 1.01 ± 0.10 ± 0.06 2.39 ± 0.13 ± 0.34 0.82 ± 0.15 ± 0.22 4.675 1.29 ± 0.05 ± 0.10 1.01 ± 0.10 ± 0.06 2.38 ± 0.13 ± 0.34 0.81 ± 0.14 ± 0.22 4.725 1.29 ± 0.05 ± 0.10 1.01 ± 0.10 ± 0.06 2.37 ± 0.13 ± 0.34 0.81 ± 0.14 ± 0.22 4.775 1.29 ± 0.05 ± 0.10 1.00 ± 0.10 ± 0.06 2.36 ± 0.13 ± 0.34 0.80 ± 0.14 ± 0.22 4.825 1.29 ± 0.05 ± 0.10 1.00 ± 0.10 ± 0.06 2.35 ± 0.13 ± 0.33 0.80 ± 0.14 ± 0.22 4.875 1.29 ± 0.05 ± 0.10 1.00 ± 0.10 ± 0.06 2.34 ± 0.13 ± 0.33 0.80 ± 0.14 ± 0.22 4.925 1.29 ± 0.05 ± 0.10 1.00 ± 0.10 ± 0.06 2.33 ± 0.13 ± 0.33 0.79 ± 0.14 ± 0.21 4.975 1.29 ± 0.05 ± 0.10 0.99 ± 0.10 ± 0.06 2.32 ± 0.13 ± 0.33 0.79 ± 0.14 ± 0.21 5.025 1.29 ± 0.05 ± 0.10 0.99 ± 0.10 ± 0.06 2.31 ± 0.13 ± 0.33 0.78 ± 0.14 ± 0.21 5.075 1.29 ± 0.05 ± 0.10 0.99 ± 0.10 ± 0.06 2.30 ± 0.13 ± 0.33 0.78 ± 0.14 ± 0.21 5.125 1.29 ± 0.05 ± 0.10 0.99 ± 0.10 ± 0.06 2.30 ± 0.13 ± 0.33 0.78 ± 0.14 ± 0.21 5.275 1.29 ± 0.05 ± 0.10 0.98 ± 0.10 ± 0.06 2.27 ± 0.13 ± 0.32 0.77 ± 0.14 ± 0.21 5.325 1.29 ± 0.05 ± 0.10 0.98 ± 0.10 ± 0.06 2.26 ± 0.13 ± 0.32 0.76 ± 0.14 ± 0.21 5.375 1.29 ± 0.05 ± 0.10 0.97 ± 0.10 ± 0.06 2.26 ± 0.13 ± 0.32 0.76 ± 0.13 ± 0.21 5.475 1.28 ± 0.05 ± 0.09 0.97 ± 0.10 ± 0.06 2.24 ± 0.12 ± 0.32 0.75 ± 0.13 ± 0.20 5.525 1.28 ± 0.05 ± 0.09 0.97 ± 0.10 ± 0.06 2.23 ± 0.12 ± 0.32 0.75 ± 0.13 ± 0.20 5.625 1.28 ± 0.05 ± 0.09 0.96 ± 0.10 ± 0.06 2.22 ± 0.12 ± 0.32 0.74 ± 0.13 ± 0.20 5.675 1.28 ± 0.05 ± 0.09 0.96 ± 0.10 ± 0.06 2.21 ± 0.12 ± 0.31 0.74 ± 0.13 ± 0.20 5.725 1.28 ± 0.05 ± 0.09 0.96 ± 0.10 ± 0.06 2.21 ± 0.12 ± 0.31 0.74 ± 0.13 ± 0.20 5.955 1.28 ± 0.05 ± 0.09 0.95 ± 0.10 ± 0.06 2.18 ± 0.12 ± 0.31 0.72 ± 0.13 ± 0.20 6.915 1.28 ± 0.04 ± 0.09 0.92 ± 0.09 ± 0.05 2.07 ± 0.12 ± 0.29 0.68 ± 0.12 ± 0.18 7.267 1.28 ± 0.04 ± 0.09 0.91 ± 0.09 ± 0.05 2.04 ± 0.11 ± 0.29 0.67 ± 0.12 ± 0.18 7.630 1.28 ± 0.04 ± 0.09 0.90 ± 0.09 ± 0.05 2.00 ± 0.11 ± 0.29 0.65 ± 0.12 ± 0.18 8.021 1.28 ± 0.04 ± 0.09 0.89 ± 0.09 ± 0.05 1.97 ± 0.11 ± 0.28 0.64 ± 0.11 ± 0.17 8.847 1.28 ± 0.04 ± 0.09 0.87 ± 0.09 ± 0.05 1.91 ± 0.11 ± 0.27 0.62 ± 0.11 ± 0.17 9.775 1.27 ± 0.04 ± 0.09 0.85 ± 0.09 ± 0.05 1.86 ± 0.10 ± 0.26 0.59 ± 0.11 ± 0.16 10.267 1.27 ± 0.04 ± 0.09 0.85 ± 0.09 ± 0.05 1.83 ± 0.10 ± 0.26 0.58 ± 0.10 ± 0.16 10.762 1.27 ± 0.04 ± 0.09 0.84 ± 0.08 ± 0.05 1.81 ± 0.10 ± 0.26 0.57 ± 0.10 ± 0.16 11.344 1.27 ± 0.04 ± 0.09 0.83 ± 0.08 ± 0.05 1.78 ± 0.10 ± 0.25 0.56 ± 0.10 ± 0.15 12.580 1.27 ± 0.04 ± 0.09 0.81 ± 0.08 ± 0.05 1.73 ± 0.10 ± 0.25 0.54 ± 0.10 ± 0.15 13.238 1.27 ± 0.04 ± 0.09 0.80 ± 0.08 ± 0.05 1.71 ± 0.10 ± 0.24 0.53 ± 0.09 ± 0.15 14.689 1.27 ± 0.04 ± 0.09 0.79 ± 0.08 ± 0.05 1.66 ± 0.09 ± 0.24 0.52 ± 0.09 ± 0.14 17.108 1.26 ± 0.04 ± 0.09 0.77 ± 0.08 ± 0.04 1.60 ± 0.09 ± 0.23 0.49 ± 0.09 ± 0.13 19.072 1.26 ± 0.04 ± 0.09 0.75 ± 0.08 ± 0.04 1.56 ± 0.09 ± 0.22 0.48 ± 0.08 ± 0.13 20.108 1.26 ± 0.04 ± 0.09 0.75 ± 0.08 ± 0.04 1.54 ± 0.09 ± 0.22 0.47 ± 0.08 ± 0.13 21.097 1.26 ± 0.04 ± 0.09 0.74 ± 0.07 ± 0.04 1.52 ± 0.08 ± 0.22 0.46 ± 0.08 ± 0.13 24.259 1.26 ± 0.04 ± 0.08 0.73 ± 0.07 ± 0.04 1.47 ± 0.08 ± 0.21 0.44 ± 0.08 ± 0.12 26.680 1.26 ± 0.04 ± 0.08 0.71 ± 0.07 ± 0.04 1.44 ± 0.08 ± 0.21 0.43 ± 0.08 ± 0.12 32.500 1.25 ± 0.04 ± 0.08 0.69 ± 0.07 ± 0.04 1.39 ± 0.08 ± 0.20 0.41 ± 0.07 ± 0.11 34.932 1.25 ± 0.04 ± 0.08 0.69 ± 0.07 ± 0.04 1.37 ± 0.08 ± 0.19 0.41 ± 0.07 ± 0.11 36.750 1.25 ± 0.04 ± 0.08 0.68 ± 0.07 ± 0.04 1.35 ± 0.08 ± 0.19 0.40 ± 0.07 ± 0.11 43.970 1.25 ± 0.04 ± 0.08 0.66 ± 0.07 ± 0.04 1.31 ± 0.07 ± 0.19 0.38 ± 0.07 ± 0.10 47.440 1.25 ± 0.04 ± 0.08 0.66 ± 0.07 ± 0.04 1.29 ± 0.07 ± 0.18 0.38 ± 0.07 ± 0.10 64.270 1.25 ± 0.04 ± 0.08 0.63 ± 0.06 ± 0.04 1.22 ± 0.07 ± 0.17 0.35 ± 0.06 ± 0.10 75.000 1.25 ± 0.04 ± 0.08 0.62 ± 0.06 ± 0.04 1.19 ± 0.07 ± 0.17 0.34 ± 0.06 ± 0.09 86.000 1.24 ± 0.04 ± 0.08 0.61 ± 0.06 ± 0.04 1.16 ± 0.06 ± 0.17 0.33 ± 0.06 ± 0.09 97.690 1.24 ± 0.04 ± 0.08 0.60 ± 0.06 ± 0.03 1.13 ± 0.06 ± 0.16 0.32 ± 0.06 ± 0.09 18

0.2 M 2 0.1 M 4 10-2 0.09 0.08 0.07 0.06 1 10 10 2 1 10 10 2 Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 10-2 10-2 M 6 M 8 10-3 10-3 1 10 10 2 1 10 10 2 Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 10-4 Figure 6: Leadig twist momets of the eutro structure fuctio for = 2, 4, 6 ad 8 as a fuctio of Q 2 (triagles), extracted usig the relativistic distributio i Eq. (11) with the deutero wave fuctio from Ref. [12]. For compariso we also show the proto (stars) ad deutero (ope crosses) momets. The errors show are statistical. depeds oly slightly o Q 2 for = 2, while it is almost completely Q 2 -idepedet for larger. The latter are i reasoable agreemet with our extracted ratios withi both statistical ad systematical errors. The /p ratio systematically decreases with ad reaches the value 0.34±0.12 stat ± 0.13 syst at = 10. To compare this with the theoretical predictios the mea value of x probed by our momets eeds to be estimated. This ca be doe i the followig way: 19

0.2 0.175 0.15 δm /M 0.125 0.1 0.075 0.05 0.025 0 2 4 6 8 10 Figure 7: Various cotributios to the systematic errors o the LT eutro momets for = 2 10: deutero LT momet error (solid); proto LT momet error (dashed); error due to the mometum cut (dotted); depedece o the relativistic treatmet (dot-dashed); offshelless ucertaity (left-hatched); wave fuctio model depedece (cross-hatched). for ay calculated -th eutro ad proto momet we estimate a correspodig average value of x defied as x (M +1 + Mp +1 )/(M + Mp ); the the ratio of eutro to proto structure fuctios is computed at the values x. The results of this procedure are show i Fig. 10 (dashed curve), where they are also compared with the ratio of the extracted LT momets (solid curve). The two ratios coicide with good accuracy, suggestig that the extracted ratio of momets is similar to the ratio of structure fuctios take at the average value x. Note that at = 10 oe has x 10 0.70, ad the F 2 /F p 2 ratio at x 0.70 is equal to 0.34 ± 0.12 stat ± 0.13 syst. As discussed above (see Fig. 7), momets for = 12 may still be accessible by improvig measuremets of the proto ad deutero structure fuctios at large x, ad with = 12 oe ca probe values of x as large as 0.75. However, higher x values are difficult to access by this method because of the rapid icrease of the ucertaities o the uclear correctios with the order. The sesitivity of the momet s ratio to the high-x tail of the PDFs is illustrated i Fig. 11, where our /p momet s ratio (triagles) is compared with the PDF parameterizatios from Ref. [30] with ad without modifyig the d-quark distributio as: d(x) d(x) + 0.1 x m (1 + x) u(x). (37) with m a parameter. Depedig o the value of m, the form i Eq. (37) allows the x 1 behavior of the d/u ratio to vary betwee 0 ad 1/5, which correspod respectively to the structure fuctio limits lim x 1 F2 (x)/f p 2 (x) of 1/4 ad 3/7. The dotted lie i 20

0.15 10-2 M 2 0.1 M 4 0.09 10-3 0.08 0.07 1 10 10 2 1 10 10 2 Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 M 6 M 8 10-3 10-3 1 10 10 2 1 10 10 2 Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 10-4 Figure 8: Extracted momets of the eutro structure fuctio for = 2, 4, 6 ad 8 as a fuctio of Q 2 i compariso to momets of differet PDFs: solid lie represets GRV [28], dashed lie represets Alekhi [31], dotted lie represets CTEQ [30], dot-dashed lie represets MRST [29]. The hatched area shows systematic errors. Fig. 11 correspods to the case m = 1 i Eq. (37), which was advocated i the NNLO aalysis of both electro ad eutrio DIS data of Ref. [36]. It appears that our results rule out such a ehacemet of the d quark distributio, sice the correctio is ot egligible eve at values of x as low as x 0.5. A larger value of the parameter m cofies the ehacemet oly to larger values of x, ad the results obtaied adoptig value m = 4, correspodig to the dashed lie i Fig. 11, are more cosistet with our 21

0.9 0.8 0.7 0.6 M /M p 0.5 0.4 0.3 0.2 0.1 0 2 4 6 8 10 Figure 9: Ratio of eutro to proto leadig twist momets for various (triagles), together with the statistical ad systematic ucertaities (hatched area). The horizotal lies deote the ratios obtaied with differet PDF sets: GRV [28] (solid), Alekhi [31] (dashed), CTEQ [30] (dotted), MRST [29] (dot-dashed). extracted /p ratio. Before closig this sectio, our results for the LT momets of the osiglet p part of the F 2 structure fuctio, M p (Q 2 ) = M(Q p 2 ) M (Q2 ), are listed i Table 3 together with their statistical ad systematic ucertaities. We poit out that a importat fractio of the total systematic error for the secod momets M p 2 (Q2 ) ad M2 (Q2 ) comes from the extrapolatio i the umeasured low-x regios. Such systematic ucertaity basically cacels out i the isovector p combiatio. This allows us to determie M p 2 (Q 2 ) withi 20% global accuracy. Low-order momets of the isovector u d quark distributios have bee recetly calculated i both uqueched [37] ad queched [38] lattice QCD. Our results for the osiglet LT momets M p ca be easily compared with those of the isovector momets ( x 1 u x 1 d ), takig ito accout that for eve values of the order oe has M p (Q 2 ) = 1 ( 3 CNS x 1 u x 1 ) d, (38) where C NS is the ( 1)-th momet of the quark coefficiet fuctio. At the scale Q 2 = 4 GeV 2 ad i the MS reormalizatio scheme at NLO, which are the scale ad the scheme adopted i Refs. [37, 38], oe gets C2 NS 1.01 ad C4 NS 1.15. The quark masses accessible i preset lattice calculatios are still relatively high ad thus the extrapolatio to the physical poit is oe of the ope issues i lattice applicatios to light hadro pheomeology. I Refs. [37, 38] the lattice data are extrapolated liearly 22