text:
# 1 Introduction

## 1 Introduction

In principle, quantum chromodynamics provides a fundamental description of hadron and nuclear physics in terms of quark and gluon degrees of freedom. QCD has been developed and successfully tested extensively, particularly in inclusive and exclusive processes involving collisions at large momentum transfer where factorization theorems and the smallness of the QCD effective coupling allow perturbative predictions. However, despite its empirical successes, many fundamental questions about QCD have not been resolved. These include a rigorous proof of color confinement, the behavior of the QCD coupling at small momenta, a rigorous description of the structure of hadrons in terms of their quark and gluon degrees of freedom, the problem of asymptotic $`n!`$ growth of the perturbation theory (renormalon phenomena), the nature of the pomeron and Reggeons, the nature of shadowing and anti-shadowing in nuclear collisions, the apparent conflict between QCD vacuum structure and the small size of the cosmological constant, and the problems of scale and scheme ambiguities in perturbative QCD expansion.

In these lectures I will focus on one of the central questions in QCD—the non-perturbative description of the proton and other hadrons as composites of confined, relativistic quark and gluon quanta. The goal is a frame-independent, quantum-mechanical representation of hadrons at the amplitude level capable of encoding multi-quark, hidden-color and gluon momentum, helicity, and flavor correlations in the form of universal process-independent hadron wavefunctions. For example, the measurement and interpretation of the basic parameters of the electroweak theory and $`CP`$ violation depends on an understanding of the dynamics and phase structure of $`B`$ decays at the amplitude level. As I will discuss in these lectures, light-front quantization allows a unifying wavefunction representation of non-perturbative hadron dynamics in QCD. Remarkably, it is now possible to measure the wavefunctions of a relativistic hadron by diffractively dissociating it into jets whose momentum distribution is correlated with the valence quarks’ momenta . It is also particularly important to understand the shape of the gauge- and process-independent meson and baryon valence-quark distribution amplitudes  $`\varphi _M(x,Q)`$, and $`\varphi _B(x_i,Q)`$. These quantities specify how a hadron shares its longitudinal momentum among its valence quarks; they control virtually all exclusive processes involving a hard scale $`Q`$, including form factors, Compton scattering and photoproduction at large momentum transfer, as well as the decay of a heavy hadron into specific final states .

What do we mean by a hadronic wavefunction, and what are its degrees of freedom? In the chiral effective Lagrangian approximation, a baryon can be represented as a classical soliton solution as in the Skyrme model. On the other hand, bag models, the observed baryon spectroscopy, and magnetic moment phenomenology suggests that baryons are composites of three “constituent” quarks or perhaps a quark bound to an effective spin-0 or spin-1 diquark. The cloudy bag model  takes into account the effect of meson-baryon fluctuations such as $`n\pi ^+`$ and $`\mathrm{\Lambda }K^+`$ in the proton. This in turn implies the existence of sea quarks with an $`\overline{u}(x)\overline{d}(x)`$ and $`s(x)\overline{s}(x)`$ asymmetries in the momentum and spin distributions of the quark sea.

The $`n`$-particle Schrödinger momentum space wavefunction $`\psi _N(\stackrel{}{p}_i)`$ of a composite system is the projection of the exact eigenstate of the equal-time Hamiltonian on the $`n`$-particle states of the non-interacting Hamiltonian, the Fock basis. It represents the amplitude for finding the constituents with three-momentum $`\stackrel{}{p}_i`$, orbital angular momentum, and spin, subject to three-momentum conservation and angular momentum sum rules. The constituents are on their mass shell, $`E_i=\sqrt{\stackrel{}{p}_i^2+m_i^2}`$ but do not conserve energy $`_{i=1}^nE_i>E=\sqrt{\stackrel{}{p}^2+M^2}`$. However, in a relativistic quantum theory, a bound-state cannot be represented as a state with a fixed number of constituents. For example, the existence of gluons which propagate between the valence quarks necessarily implies that the hadron wavefunction must describe states with an arbitrary number of gluons. Thus a hadronic wavefunction must describe fluctuations in particle number $`n`$, as well as momenta and spin. One has to take into account fluctuations in the wavefunction which allow for any number of sea quarks, as long as the total quantum numbers of the constituents are compatible with the overall quantum numbers of the baryon.

Gluon exchange between quarks of different nucleons in a nucleus will leave the nucleons in an exited octet state of $`SU(3)_C`$, thus requiring the existence of “hidden color” configurations in a nuclear wavefunction. What is worse is that the description of the state at a given time $`t_0`$ depends not only on the choice of gauge but also on the Lorentz frame of the observer since two observers will differ on their definition of simultaneity.

A number of non-perturbative formalisms have been developed which in principle could provide the exact spectra and wavefunctions of bound states in relativistic quantum field theory, including QCD. These include Euclidean lattice gauge theory, the covariant Bethe-Salpeter bound state equation, and Hamiltonian methods based on equal-time or light-front quantization. Lattice gauge theory, particularly in the heavy quark approximation, has been an enormously helpful guide to the hadronic and glueball spectrum of QCD; however, progress in evaluating hadron wavefunctions has been so far limited to the evaluation of certain moments of hadronic matrix elements .

The Bethe-Salpeter formalism in principle provides a rigorous and systematic, renormalizable, and explicitly covariant treatment of bound states of two interacting fields. There has been recent progress applying model forms of Dyson-Schwinger equations to the Bethe-Salpeter formalism in order to incorporate the effects of dynamical symmetry breaking in QCD . However, the actual interaction kernel in the bound state equation and the evaluation of current matrix elements requires a sum over an infinite number of irreducible graphs; the use of a single kernel (the ladder approximation) breaks gauge invariance, crossing symmetry, and has an improper heavy quark limit. It is also intrinsically difficult to apply the Bethe Salpeter to bound states of three or more constituents.

There has been progress applying equal-time Hamiltonian methods to QCD in Coulomb gauge ; however, the renormalization procedure in this formalism is not well-defined and the results are typically frame-dependent. The evaluation of current matrix elements in the equal-time formalism is especially problematic since one must include contributions to the current from vacuum excitations. For similar reasons, the expansion of an eigenstate on the equal-time Fock basis of free quanta is not well-defined.

In the light-front quantization, one takes the light-cone time variable $`t+z/c`$ as the evolution parameter instead of ordinary time $`t`$. (The $`\widehat{z}`$ direction is an arbitrary reference direction.) The method is often called “light-front” quantization rather than “light-cone” quantization since the equation $`x^+=\tau =0`$ defines a hyperplane corresponding to a light-front. The light-front fixes the initial boundary conditions of a composite system as its constituents are intercepted by a light-wave evaluated at a specific value of $`x^+=t+z/c`$. In contrast, determining an atomic wavefunction at a given instant $`t=t_0`$ requires measuring the simultaneous scattering of $`Z`$ photons on the $`Z`$ electrons.

An extensive review and guide to the light-front quantization literature can be found in Ref. . I will use here the notation $`A^\mu =(A^+,A,A_{}),`$ where $`A^\pm =A^0\pm A^z`$ and the metric $`AB=\frac{1}{2}(A^+B^{}+A^{}B^+)A_{}B_{}.`$

The origins of the light-front quantization method can be traced to Dirac  who noted that 7 out of the 10 Poincare’ generators, including a Lorentz boost $`K_3`$, are kinematical (interaction-independent) when one quantizes a theory at fixed light-cone time. This in turn leads to the remarkable property that the light-front wavefunctions of a hadron are independent of the hadron’s total momentum, whether it is at rest or moving. Thus once one has solved for the light-front wavefunctions, one can compute hadron matrix elements of currents between hadronic states of arbitrary momentum. In contrast, knowing the rest frame wavefunction at equal time, does not determine the moving hadron’s wavefunction.

Light-front wavefunctions are related to momentum-space Bethe-Salpeter wavefunctions by integrating over the relative momenta $`k^{}=k^0k^z`$ since this projects out $`x^+=0.`$ The light-front Fock space is the eigenstates of the free light-front Hamiltonian; i.e., it is a Hilbert space of non-interacting quarks and gluons, each of which satisfy $`k^2=m^2`$ and $`k^{}=\frac{m^2+k_{}^2}{k^+}0.`$ The light-front wavefunctions are the projections of the hadronic eigenstate on the light-front Fock basis, the complete set of color singlet states of the free Hamiltonian. An essentially equivalent approach, pioneered by Weinberg , is to evaluate the equal-time theory from the perspective of an observer moving in the negative $`\widehat{z}`$ direction with arbitrarily large momentum $`P_z\mathrm{}.`$ The light-cone fraction $`x=\frac{k^+}{p^+}`$ of a constituent can be identified with the longitudinal momentum $`x=\frac{k^z}{P^z}`$ in a hadron moving with large momentum $`P^z.`$

It is convenient to define the invariant light-front Hamiltonian: $`H_{LC}^{QCD}=P^+P^{}\stackrel{}{P}_{}^{}{}_{}{}^{2}`$ where $`P^\pm =P^0\pm P^z`$. The operator $`P^{}=i\frac{d}{d\tau }`$ generates light-cone time translations. The $`P^+`$ and $`\stackrel{}{P}_{}`$ momentum operators are independent of the interactions, and thus are conserved at all orders. The eigen-spectrum of $`H_{LC}^{QCD}`$ in principle gives the entire mass squared spectrum of color-singlet hadron states in QCD, together with their respective light-front wavefunctions. For example, the proton state satisfies: $`H_{LC}^{QCD}|\mathrm{\Psi }_p=M_p^2|\mathrm{\Psi }_p`$. The projection of the proton’s eigensolution $`|\mathrm{\Psi }_p`$ on the color-singlet $`B=1`$, $`Q=1`$ eigenstates $`\{|n\}`$ of the free Hamiltonian $`H_{LC}^{QCD}(g=0)`$ gives the light-front Fock expansion:

$`|\mathrm{\Psi }_p;P^+,\stackrel{}{P}_{},\lambda `$ $`=`$ $`{\displaystyle \underset{n3,\lambda _i}{}}{\displaystyle \mathrm{\Pi }_{i=1}^n\frac{d^2k_idx_i}{\sqrt{x_i}16\pi ^3}}`$
$`16\pi ^3\delta \left(1{\displaystyle \underset{j}{\overset{n}{}}}x_j\right)\delta ^{(2)}\left({\displaystyle \underset{\mathrm{}}{\overset{n}{}}}\stackrel{}{k}_{\mathrm{}}\right)`$
$`|n;x_iP^+,x_i\stackrel{}{P}_{}+\stackrel{}{k}_i,\lambda _i\psi _{n/p}(x_i,\stackrel{}{k}_i,\lambda _i).`$
It is especially convenient to develop the light-front formalism in the light-cone gauge $`A^+=A^0+A^z=0`$. In this gauge the $`A^{}`$ field becomes a dependent degree of freedom, and it can be eliminated from the gauge theory Hamiltonian, with the addition of a set of specific instantaneous light-cone time interactions. In fact in $`QCD(1+1)`$ theory, this instantaneous interaction provides the confining linear $`x^{}`$ interaction between quarks. In $`3+1`$ dimensions, the transverse field $`A^{}`$ propagates massless spin-one gluon quanta with polarization vectors  which satisfy both the gauge condition $`ϵ_\lambda ^+=0`$ and the Lorentz condition $`kϵ=0`$. Thus no extra condition on the Hilbert space is required.

The light-front Fock wavefunctions $`\psi _{n/H}(x_i,\stackrel{}{k}_i,\lambda _i)`$ thus interpolate between the hadron $`H`$ and its quark and gluon degrees of freedom. The light-cone momentum fractions of the constituents, $`x_i=k_i^+/P^+`$ with $`_{i=1}^nx_i=1,`$ and the transverse momenta $`\stackrel{}{k}_i`$ with $`_{i=1}^n\stackrel{}{k}_i=\stackrel{}{0}_{}`$ appear as the momentum coordinates of the light-front Fock wavefunctions. A crucial feature is the frame-independence of the light-front wavefunctions. The $`x_i`$ and $`\stackrel{}{k}_i`$ are relative coordinates independent of the hadron’s momentum $`P^\mu `$. The actual physical transverse momenta are $`\stackrel{}{p}_i=x_i\stackrel{}{P}_{}+\stackrel{}{k}_i.`$ The $`\lambda _i`$ label the light-front spin $`S^z`$ projections of the quarks and gluons along the $`z`$ direction. The physical gluon polarization vectors $`ϵ^\mu (k,\lambda =\pm 1)`$ are specified in light-cone gauge by the conditions $`kϵ=0,\eta ϵ=ϵ^+=0.`$ Each light-front Fock wavefunction satisfies conservation of the $`z`$ projection of angular momentum: $`J^z=_{i=1}^nS_i^z+_{j=1}^{n1}l_j^z.`$ The sum over $`S_i^z`$ represents the contribution of the intrinsic spins of the $`n`$ Fock state constituents. The sum over orbital angular momenta $`l_j^z=\mathrm{i}(k_j^1\frac{}{k_j^2}k_j^2\frac{}{k_j^1})`$ derives from the $`n1`$ relative momenta. This excludes the contribution to the orbital angular momentum due to the motion of the center of mass, which is not an intrinsic property of the hadron .

The interaction Hamiltonian of QCD in light-cone gauge can be derived by systematically applying the Dirac bracket method to identify the independent fields . It contains the usual Dirac interactions between the quarks and gluons, the three-point and four-point gluon non-Abelian interactions plus instantaneous light-front-time gluon exchange and quark exchange contributions

$`_{int}`$ $`=`$ $`g\overline{\psi }^i\gamma ^\mu A_{\mu }^{}{}_{}{}^{ij}\psi ^j`$ (3)
$`+{\displaystyle \frac{g}{2}}f^{abc}(_\mu A_{}^{a}{}_{\nu }{}^{}_\nu A_{}^{a}{}_{\mu }{}^{})A^{b\mu }A^{c\nu }`$
$`+{\displaystyle \frac{g^2}{4}}f^{abc}f^{ade}A_{b\mu }A^{d\mu }A_{c\nu }A^{e\nu }`$
$`{\displaystyle \frac{g^2}{2}}\overline{\psi }^i\gamma ^+(\gamma ^{^{}}A_{^{}})^{ij}{\displaystyle \frac{1}{i_{}}}(\gamma ^{}A_{})^{jk}\psi ^k`$
$`{\displaystyle \frac{g^2}{2}}j_{}^{+}{}_{a}{}^{}{\displaystyle \frac{1}{(_{})^2}}j_{}^{+}{}_{a}{}^{}`$
where

$$j_{}^{+}{}_{a}{}^{}=\overline{\psi }^i\gamma ^+(t_a)^{ij}\psi ^j+f_{abc}(_{}A_{b\mu })A^{c\mu }.$$
(4)
In light-cone time-ordered perturbation theory, a Green’s functions is expanded as a power series in the interactions with light-front energy denominators $`_{\mathrm{initial}}k_i^{}_{\mathrm{intermediate}}k_i^{}+iϵ`$ replacing the usual energy denominators. \[For a review see Ref. .\] In general each Feynman diagram with $`n`$ vertices corresponds to the sum of $`n!`$ time-ordered contributions. However, in light-cone-time-ordered perturbation theory, only those few graphs where all $`k_i^+0`$ survive. In addition the form of the light-front kinetic energies is rational: $`k^{}=\frac{k_{}^2+m^2}{k^+}`$, replacing the nonanalytic $`k^0=\sqrt{\stackrel{}{k}^2+m^2}`$ of equal-time theory. Thus light-cone-time-ordered perturbation theory provides a viable computational method where one can trace the physical evolution of a process. The integration measures are only three-dimensional $`d^2k_{}dx`$; in effect, the $`k^{}`$ integral of the covariant perturbation theory is performed automatically.

Alternatively, one derive Feynman rules for QCD in light-cone gauge, thus allowing the use of standard covariant computational tools and renormalization methods including dimensional regularization. Prem Srivastava and I  have recently presented a systematic study of light-front-quantized gauge theory in light-cone gauge using a Dyson-Wick S-matrix expansion based on light-cone-time-ordered products. The gluon propagator has the form

$$0|T(A_{}^{a}{}_{\mu }{}^{}(x)A_{}^{b}{}_{\nu }{}^{}(0))|0=\frac{i\delta ^{ab}}{(2\pi )^4}d^4ke^{ikx}\frac{D_{\mu \nu }(k)}{k^2+iϵ}$$
(5)
where we have defined

$$D_{\mu \nu }(k)=D_{\nu \mu }(k)=g_{\mu \nu }+\frac{n_\mu k_\nu +n_\nu k_\mu }{(nk)}\frac{k^2}{(nk)^2}n_\mu n_\nu .$$
(6)
Here $`n_\mu `$ is a null four-vector, gauge direction, whose components are chosen to be $`n_\mu =\delta _{\mu }^{}{}_{}{}^{+}`$, $`n^\mu =\delta _{}^{\mu }{}_{}{}^{}`$. Note also

$`D_{\mu \lambda }(k)D_{}^{\lambda }{}_{\nu }{}^{}(k)=D_\mu (k)D_{}^{}{}_{\nu }{}^{}(k)`$ $`=`$ $`D_{\mu \nu }(k),`$
$`k^\mu D_{\mu \nu }(k)=0,`$ $`n^\mu D_{\mu \nu }(k)D_\nu (k)=0,`$
$`D_{\lambda \mu }(q)D^{\mu \nu }(k)D_{\nu \rho }(q^{})`$ $`=`$ $`D_{\lambda \mu }(q)D^{\mu \rho }(q^{}).`$ (7)
The gauge field propagator $`iD_{\mu \nu }(k)/(k^2+iϵ)`$ is transverse not only to the gauge direction $`n_\mu `$ but also to $`k_\mu `$, i.e., it is doubly-transverse. This leads to appreciable simplifications in the computations in QCD. For example, the coupling of gluons to propagators carrying high momenta is automatic. The absence of collinear divergences in irreducible diagrams in the light-cone gauge greatly simplifies the leading-twist factorization of soft and hard gluonic corrections in high momentum transfer inclusive and exclusive reactions  since the numerators associated with the gluon coupling only have transverse components. The renormalization factors in the light-cone gauge are independent of the reference direction $`n^\mu `$. Since the gluon only has physical polarization, its renormalization factors satisfy $`Z_1=Z_3`$. Because of its explicit unitarity in each graph, the doubly-transverse gauge is well suited for calculations identifying the “pinch” effective charge .

The running coupling constant and QCD $`\beta `$ function have also been computed at one loop in the doubly-transverse light-cone gauge . It is also possible to effectively quantize QCD using light-front methods in covariant Feynman gauge .

A remarkable advantage of light-front quantization is that the vacuum state $`|\mathrm{\hspace{0.17em}0}`$ of the full QCD Hamiltonian evidently coincides with the free vacuum. The light-front vacuum is effectively trivial if the interaction Hamiltonian applied to the perturbative vacuum is zero. Note that all particles in the Hilbert space have positive energy $`k^0=\frac{1}{2}(k^++k^{})`$, and thus positive light-front $`k^\pm `$. Since the plus momenta $`k_i^+`$ is conserved by the interactions, the perturbative vacuum can only couple to states with particles in which all $`k_i^+`$ = 0; i.e., so called zero-mode states. In the case of QED, a massive electron cannot have $`k^+=0`$ unless it also has infinite energy. In a remarkable calculation, Bassetto and collaborators  have shown that the computation of the spectrum of $`QCD(1+1)`$ in equal time quantization requires constructing the full spectrum of non perturbative contributions (instantons). In contrast, in the light-front quantization of gauge theory, where the $`k^+=0`$ singularity of the instantaneous interaction is defined by a simple infrared regularization, one obtains the correct spectrum of $`QCD(1+1)`$ without any need for vacuum-related contributions.

In the case of QCD(3+1), the momentum-independent four-gluon non-Abelian interaction in principle can couple the perturbative vacuum to a state with four collinear gluons in which all of the gluons have all components $`k_i^\mu =0,`$ thus hinting at role for zero modes in theories with massless quanta. In fact, zero modes of auxiliary fields are necessary to distinguish the theta-vacua of massless QED(1+1) , or to represent a theory in the presence of static external boundary conditions or other constraints. Zero-modes provide the light-front representation of spontaneous symmetry breaking in scalar theories .

There are a number of other simplifications of the light-front formalism:

1. The light-front wavefunctions describe quanta which have positive energy, positive norm, and physical polarization. The formalism is thus physical, and unitary. No ghosts fields appear explicitly, even in non-Abelian theory. The wavefunctions are only functions of three rather than four physical momentum variables: the light-front momentum fractions $`x_i`$ and transverse momenta $`k_{}`$. The quarks and gluons each have two physical polarization states.

2. The set of light-front wavefunctions provide a frame-independent, quantum-mechanical description of hadrons at the amplitude level capable of encoding multi-quark and gluon momentum, helicity, and flavor correlations in the form of universal process-independent hadron wavefunctions. Matrix elements of spacelike currents such as the spacelike electromagnetic form factors have an exact representation in terms of simple overlaps of the light-front wavefunctions in momentum space with the same $`x_i`$ and unchanged parton number . In the case of timelike decays, such as those determined by semileptonic $`B`$ decay, one needs to include contributions in which the parton number $`\mathrm{\Delta }n=2.`$ . The leading-twist off-forward parton distributions measured in deeply virtual Compton scattering have a similar light-front wavefunction representation .

3. The high $`x1`$ and high $`k_{}`$ limits of the hadron wavefunctions control processes and reactions in which the hadron wavefunctions are highly stressed. Such configurations involve far-off-shell intermediate states and can be systematically treated in perturbation theory .

4. The leading-twist structure functions $`q_i(x,Q)`$ and $`g(x,Q)`$ measured in deep inelastic scattering can be computed from the absolute squares of the light-front wavefunctions, integrated over the transverse momentum up to the resolution scale $`Q`$. All helicity distributions are thus encoded in terms of the light-front wavefunctions. The DGLAP evolution of the structure functions can be derived from the high $`k_{}`$ properties of the light-front wavefunctions. Thus given the light-front wavefunctions, one can compute  all of the leading twist helicity and transversity distributions measured in polarized deep inelastic lepton scattering. For example, the helicity-specific quark distributions at resolution $`\mathrm{\Lambda }`$ correspond to

$`q_{\lambda _q/\mathrm{\Lambda }_p}(x,\mathrm{\Lambda })`$ $`=`$ $`{\displaystyle \underset{n,q_a}{}}{\displaystyle \underset{j=1}{\overset{n}{}}\frac{dx_jd^2k_j}{16\pi ^3}\underset{\lambda _i}{}|\psi _{n/H}^{(\mathrm{\Lambda })}(x_i,\stackrel{}{k}_i,\lambda _i)|^2}`$
$`\times 16\pi ^3\delta \left(1{\displaystyle \underset{i}{\overset{n}{}}}x_i\right)\delta ^{(2)}\left({\displaystyle \underset{i}{\overset{n}{}}}\stackrel{}{k}_i\right)\delta (xx_q)\delta _{\lambda ,\lambda _q}\mathrm{\Theta }(\mathrm{\Lambda }^2_n^2),`$
where the sum is over all quarks $`q_a`$ which match the quantum numbers, light-front momentum fraction $`x,`$ and helicity of the struck quark. Similarly, the transversity distributions and off-diagonal helicity convolutions are defined as a density matrix of the light-front wavefunctions. This defines the LC factorization scheme  where the invariant mass squared $`_n^2=_{i=1}^n(k_i^2+m_i^2)/x_i`$ of the $`n`$ partons of the light-front wavefunctions is limited to $`_n^2<\mathrm{\Lambda }^2`$

5. The distribution of spectator particles in the final state in the proton fragmentation region in deep inelastic scattering at an electron-proton collider are encoded in the light-front wavefunctions of the target proton. Conversely, the light-front wavefunctions can be used to describe the coalescence of comoving quarks into final state hadrons.

6. The light-front wavefunctions also specify the multi-quark and gluon correlations of the hadron. Despite the many sources of power-law corrections to the deep inelastic cross section, certain types of dynamical contributions will stand out at large $`x_{bj}`$ since they arise from compact, highly-correlated fluctuations of the proton wavefunction. In particular, there are particularly interesting dynamical $`𝒪(1/Q^2)`$ corrections which are due to the interference of quark currents; i.e., contributions which involve leptons scattering amplitudes from two different quarks of the target nucleon .

7. The higher Fock states of the light hadrons describe the sea quark structure of the deep inelastic structure functions, including “intrinsic” strangeness and charm fluctuations specific to the hadron’s structure rather than gluon substructure . Ladder relations connecting state of different particle number follow from the QCD equation of motion and lead to Regge behavior of the quark and gluon distributions at $`x0`$ .

8. The gauge- and process-independent meson and baryon valence-quark distribution amplitudes $`\varphi _M(x,Q)`$, and $`\varphi _B(x_i,Q)`$ which control exclusive processes involving a hard scale $`Q`$, including heavy quark decays, are given by the valence light-front Fock state wavefunctions integrated over the transverse momentum up to the resolution scale $`Q`$. The evolution equations for distribution amplitudes follow from the perturbative high transverse momentum behavior of the light-front wavefunctions .

9. The line-integrals needed to defining distribution amplitudes and structure functions as gauge invariant matrix elements of operator products vanish in light-front gauge.

10. Proofs of factorization theorems in hard exclusive and inclusive reactions are greatly simplified since the propagating gluons in light-cone gauge couple only to transverse currents; collinear divergences are thus automatically suppressed.

11. At high energies each light-front Fock state interacts distinctly; e.g., Fock states with small particle number and small impact separation have small color dipole moments and can traverse a nucleus with minimal interactions. This is the basis for the predictions for “color transparency” in hard quasi-exclusive  and diffractive reactions .

12. The Fock state wavefunctions of hadron can be resolved by a high energy diffractive interaction, producing forward jets with momenta which follow the light-front momenta of the wavefunction .

13. The deuteron form factor at high $`Q^2`$ is sensitive to wavefunction configurations where all six quarks overlap within an impact separation $`b_i<𝒪(1/Q).`$ The leading power-law fall off predicted by QCD is $`F_d(Q^2)=f(\alpha _s(Q^2))/(Q^2)^5`$, where, asymptotically, $`f(\alpha _s(Q^2))\alpha _s(Q^2)^{5+2\gamma }`$ . In general, the six-quark wavefunction of a deuteron is a mixture of five different color-singlet states. The dominant color configuration at large distances corresponds to the usual proton-neutron bound state. However at small impact space separation, all five Fock color-singlet components eventually evolve to a state with equal weight, i.e., the deuteron wavefunction evolves to 80% “hidden color” . The relatively large normalization of the deuteron form factor observed at large $`Q^2`$ hints at sizable hidden-color contributions . Hidden color components can also play a predominant role in the reaction $`\gamma dJ/\psi pn`$ at threshold if it is dominated by the multi-fusion process $`\gamma ggJ/\psi `$ . Hard exclusive nuclear processes can also be analyzed in terms of “reduced amplitudes” which remove the effects of nucleon substructure.

Light-front wavefunctions are thus the frame-independent interpolating functions between hadron and quark and gluon degrees of freedom. Hadron amplitudes are computed from the convolution of the light-front wavefunctions with irreducible quark-gluon amplitudes. More generally, all multi-quark and gluon correlations in the bound state are represented by the light-front wavefunctions. The light-front Fock representation is thus a representation of the underlying quantum field theory. I will discuss progress in computing light-front wavefunctions directly from QCD in Sections 9 and 10.

## 2 Other Theoretical Tools

In addition to the light-front Fock expansion, a number of other useful theoretical tools are available to eliminate theoretical ambiguities in QCD predictions:

(1) Conformal symmetry provides a template for QCD predictions , leading to relations between observables which are present even in a theory which is not scale invariant. For example, the natural representation of distribution amplitudes is in terms of an expansion of orthonormal conformal functions multiplied by anomalous dimensions determined by QCD evolution equations . Thus an important guide in QCD analyses is to identify the underlying conformal relations of QCD which are manifest if we drop quark masses and effects due to the running of the QCD couplings. In fact, if QCD has an infrared fixed point (vanishing of the Gell Mann-Low function at low momenta), the theory will closely resemble a scale-free conformally symmetric theory in many applications.

(2) Commensurate scale relations  are perturbative QCD predictions which relate observable to observable at fixed relative scale, such as the “generalized Crewther relation” , which connects the Bjorken and Gross-Llewellyn Smith deep inelastic scattering sum rules to measurements of the $`e^+e^{}`$ annihilation cross section. Such relations have no renormalization scale or scheme ambiguity. The coefficients in the perturbative series for commensurate scale relations are identical to those of conformal QCD; thus no infrared renormalons are present . One can identify the required conformal coefficients at any finite order by expanding the coefficients of the usual PQCD expansion around a formal infrared fixed point, as in the Banks-Zak method . All non-conformal effects are absorbed by fixing the ratio of the respective momentum transfer and energy scales. In the case of fixed-point theories, commensurate scale relations relate both the ratio of couplings and the ratio of scales as the fixed point is approached .

(3) $`\alpha _V`$ and Skeleton Schemes. A physically natural scheme for defining the QCD coupling in exclusive and other processes is the $`\alpha _V(Q^2)`$ scheme defined from the potential of static heavy quarks. Heavy-quark lattice gauge theory can provide highly precise values for the coupling. All vacuum polarization corrections due to fermion pairs are then automatically and analytically incorporated into the Gell Mann-Low function, thus avoiding the problem of explicitly computing and resumming quark mass corrections related to the running of the coupling . The use of a finite effective charge such as $`\alpha _V`$ as the expansion parameter also provides a basis for regulating the infrared nonperturbative domain of the QCD coupling. A similar coupling and scheme can be based on an assumed skeleton expansion of the theory .

(4) The Abelian Correspondence Principle. One can consider QCD predictions as analytic functions of the number of colors $`N_C`$ and flavors $`N_F`$. In particular, one can show at all orders of perturbation theory that PQCD predictions reduce to those of an Abelian theory at $`N_C0`$ with $`\widehat{\alpha }=C_F\alpha _s`$ and $`\widehat{N}_F=2N_F/C_F`$ held fixed . There is thus a deep connection between QCD processes and their corresponding QED analogs.

## 3 Applications of light-front Wavefunctions to Current Matrix Elements

The light-front Fock representation of current matrix elements has a number of simplifying properties. The space-like local operators for the coupling of photons, gravitons and the deep inelastic structure functions can all be expressed as overlaps of light-front wavefunctions with the same number of Fock constituents. This is possible since one can choose the special frame $`q^+=0`$  for space-like momentum transfer and take matrix elements of “plus” components of currents such as $`J^+`$ and $`T^{++}`$. No contributions to the current matrix elements from vacuum fluctuations occur. Similarly, given the local operators for the energy-momentum tensor $`T^{\mu \nu }(x)`$ and the angular momentum tensor $`M^{\mu \nu \lambda }(x)`$, one can directly compute momentum fractions, spin properties, and the form factors $`A(q^2)`$ and $`B(q^2)`$ appearing in the coupling of gravitons to composite systems .

In the case of a spin-$`\frac{1}{2}`$ composite system, the Dirac and Pauli form factors $`F_1(q^2)`$ and $`F_2(q^2)`$ are defined by

$$P^{}|J^\mu (0)|P=\overline{u}(P^{})\left[F_1(q^2)\gamma ^\mu +F_2(q^2)\frac{i}{2M}\sigma ^{\mu \alpha }q_\alpha \right]u(P),$$
(9)
where $`q^\mu =(P^{}P)^\mu `$ and $`u(P)`$ is the bound state spinor. In the light-front formalism it is convenient to identify the Dirac and Pauli form factors from the helicity-conserving and helicity-flip vector current matrix elements of the $`J^+`$ current :

$$P+q,\left|\frac{J^+(0)}{2P^+}\right|P,=F_1(q^2),$$
(10)
$$P+q,\left|\frac{J^+(0)}{2P^+}\right|P,=(q^1\mathrm{i}q^2)\frac{F_2(q^2)}{2M}.$$
(11)
The magnetic moment of a composite system is one of its most basic properties. The magnetic moment is defined at the $`q^20`$ limit,

$$\mu =\frac{e}{2M}\left[F_1(0)+F_2(0)\right],$$
(12)
where $`e`$ is the charge and $`M`$ is the mass of the composite system. We use the standard light-front frame ($`q^\pm =q^0\pm q^3`$):

$`q`$ $`=`$ $`(q^+,q^{},\stackrel{}{q}_{})=(0,{\displaystyle \frac{q^2}{P^+}},\stackrel{}{q}_{}),`$
$`P`$ $`=`$ $`(P^+,P^{},\stackrel{}{P}_{})=(P^+,{\displaystyle \frac{M^2}{P^+}},\stackrel{}{0}_{}),`$ (13)
where $`q^2=2Pq=\stackrel{}{q}_{}^2`$ is 4-momentum square transferred by the photon.

The Pauli form factor and the anomalous magnetic moment $`\kappa =\frac{e}{2M}F_2(0)`$ can then be calculated from the expression

$$(q^1\mathrm{i}q^2)\frac{F_2(q^2)}{2M}=\underset{a}{}\frac{\mathrm{d}^2\stackrel{}{k}_{}\mathrm{d}x}{16\pi ^3}\underset{j}{}e_j\psi _a^{}(x_i,\stackrel{}{k}_i^{},\lambda _i)\psi _a^{}(x_i,\stackrel{}{k}_i,\lambda _i),$$
(14)
where the summation is over all contributing Fock states $`a`$ and struck constituent charges $`e_j`$. The arguments of the final-state light-front wavefunction are

$$\stackrel{}{k}_i^{}=\stackrel{}{k}_i+(1x_i)\stackrel{}{q}_{}$$
(15)
for the struck constituent and

$$\stackrel{}{k}_i^{}=\stackrel{}{k}_ix_i\stackrel{}{q}_{}$$
(16)
for each spectator. Notice that the magnetic moment must be calculated from the spin-flip non-forward matrix element of the current. In the ultra-relativistic limit where the radius of the system is small compared to its Compton scale $`1/M`$, the anomalous magnetic moment must vanish . The light-front formalism is consistent with this theorem.

The form factors of the energy-momentum tensor for a spin-$`\frac{1}{2}`$ composite are defined by

$`P^{}|T^{\mu \nu }(0)|P`$ $`=`$ $`\overline{u}(P^{})[A(q^2)\gamma ^{(\mu }\overline{P}^{\nu )}+B(q^2){\displaystyle \frac{i}{2M}}\overline{P}^{(\mu }\sigma ^{\nu )\alpha }q_\alpha `$ (17)
$`+C(q^2){\displaystyle \frac{1}{M}}(q^\mu q^\nu g^{\mu \nu }q^2)]u(P),`$
where $`q^\mu =(P^{}P)^\mu `$, $`\overline{P}^\mu =\frac{1}{2}(P^{}+P)^\mu `$, $`a^{(\mu }b^{\nu )}=\frac{1}{2}(a^\mu b^\nu +a^\nu b^\mu )`$, and $`u(P)`$ is the spinor of the system. One can also readily obtain the light-front representation of the $`A(q^2)`$ and $`B(q^2)`$ form factors of the energy-tensor Eq. (17) . In the interaction picture, only the non-interacting parts of the energy momentum tensor $`T^{++}(0)`$ need to be computed:

$$P+q,\left|\frac{T^{++}(0)}{2(P^+)^2}\right|P,=A(q^2),$$
(18)
$$P+q,\left|\frac{T^{++}(0)}{2(P^+)^2}\right|P,=(q^1\mathrm{i}q^2)\frac{B(q^2)}{2M}.$$
(19)
The $`A(q^2)`$ and $`B(q^2)`$ form factors Eqs. (18) and (19) are similar to the $`F_1(q^2)`$ and $`F_2(q^2)`$ form factors Eqs. (10) and (11) with an additional factor of the light-front momentum fraction $`x=k^+/P^+`$ of the struck constituent in the integrand. The $`B(q^2)`$ form factor is obtained from the non-forward spin-flip amplitude. The value of $`B(0)`$ is obtained in the $`q^20`$ limit. The angular momentum projection of a state is given by

$$J^i=\frac{1}{2}ϵ^{ijk}d^3xT^{0k}x^jT^{0j}x^k=A(0)L^i+\left[A(0)+B(0)\right]\overline{u}(P)\frac{1}{2}\sigma ^iu(P).$$
(20)
This result is derived using a wave-packet description of the state. The $`L^i`$ term is the orbital angular momentum of the center of mass motion with respect to an arbitrary origin and can be dropped. The coefficient of the $`L^i`$ term must be 1; $`A(0)=1`$ also follows when we evaluate the four-momentum expectation value $`P^\mu `$. Thus the total intrinsic angular momentum $`J^z`$ of a nucleon can be identified with the values of the form factors $`A(q^2)`$ and $`B(q^2)`$ at $`q^2=0:`$

$$J^z=\frac{1}{2}\sigma ^z\left[A(0)+B(0)\right].$$
(21)
The anomalous moment coupling $`B(0)`$ to a graviton is shown to vanish for any composite system. This remarkable result, first derived by Okun and Kobzarev , is shown to follow directly from the Lorentz boost properties of the light-front Fock representation .

Dae Sung Hwang, Bo-Qiang Ma, Ivan Schmidt, and I  have recently shown that the light-front wavefunctions generated by the radiative corrections to the electron in QED provides a simple system for understanding the spin and angular momentum decomposition of relativistic systems. This perturbative model also illustrates the interconnections between Fock states of different number. The model is patterned after the quantum structure which occurs in the one-loop Schwinger $`\alpha /2\pi `$ correction to the electron magnetic moment . In effect, we can represent a spin-$`\frac{1}{2}`$ system as a composite of a spin-$`\frac{1}{2}`$ fermion and spin-one vector boson with arbitrary masses. A similar model has been used to illustrate the matrix elements and evolution of light-front helicity and orbital angular momentum operators . This representation of a composite system is particularly useful because it is based on two constituents but yet is totally relativistic. We can then explicitly compute the form factors $`F_1(q^2)`$ and $`F_2(q^2)`$ of the electromagnetic current, and the various contributions to the form factors $`A(q^2)`$ and $`B(q^2)`$ of the energy-momentum tensor.

Another remarkable advantage of the light-front formalism is that exclusive semileptonic $`B`$-decay amplitudes such as $`BA\mathrm{}\overline{\nu }`$ can also be evaluated exactly . The time-like decay matrix elements require the computation of the diagonal matrix element $`nn`$ where parton number is conserved, and the off-diagonal $`n+1n1`$ convolution where the current operator annihilates a $`q\overline{q^{}}`$ pair in the initial $`B`$ wavefunction. See Fig. 1. This term is a consequence of the fact that the time-like decay $`q^2=(p_{\mathrm{}}+p_{\overline{\nu }})^2>0`$ requires a positive light-front momentum fraction $`q^+>0`$. Conversely for space-like currents, one can choose $`q^+=0`$, as in the Drell-Yan-West representation of the space-like electromagnetic form factors. However, as can be seen from the explicit analysis of the form factor in a perturbative model, the off-diagonal convolution can yield a nonzero $`q^+/q^+`$ limiting form as $`q^+0`$. This extra term appears specifically in the case of “bad” currents such as $`J^{}`$ in which the coupling to $`q\overline{q}`$ fluctuations in the light-front wavefunctions are favored. In effect, the $`q^+0`$ limit generates $`\delta (x)`$ contributions as residues of the $`n+1n1`$ contributions. The necessity for such “zero mode” $`\delta (x)`$ terms has been noted by Chang, Root and Yan , Burkardt , and Ji and Choi .

The off-diagonal $`n+1n1`$ contributions give a new perspective for the physics of $`B`$-decays. A semileptonic decay involves not only matrix elements where a quark changes flavor, but also a contribution where the leptonic pair is created from the annihilation of a $`q\overline{q^{}}`$ pair within the Fock states of the initial $`B`$ wavefunction. The semileptonic decay thus can occur from the annihilation of a nonvalence quark-antiquark pair in the initial hadron. This feature will carry over to exclusive hadronic $`B`$-decays, such as $`B^0\pi ^{}D^+`$. In this case the pion can be produced from the coalescence of a $`d\overline{u}`$ pair emerging from the initial higher particle number Fock wavefunction of the $`B`$. The $`D`$ meson is then formed from the remaining quarks after the internal exchange of a $`W`$ boson.

In principle, a precise evaluation of the hadronic matrix elements needed for $`B`$-decays and other exclusive electroweak decay amplitudes requires knowledge of all of the light-front Fock wavefunctions of the initial and final state hadrons. In the case of model gauge theories such as QCD(1+1)  or collinear QCD  in one-space and one-time dimensions, the complete evaluation of the light-front wavefunction is possible for each baryon or meson bound-state using the DLCQ method. It would be interesting to use such solutions as a model for physical $`B`$-decays.

## 4 Light-front Representation of Deeply Virtual Compton Scattering

The virtual Compton scattering process $`\frac{d\sigma }{dt}(\gamma ^{}p\gamma p)`$ for large initial photon virtuality $`q^2=Q^2`$ has extraordinary sensitivity to fundamental features of the proton’s structure. Even though the final state photon is on-shell, the deeply virtual process probes the elementary quark structure of the proton near the light cone as an effective local current. In contrast to deep inelastic scattering, which measures only the absorptive part of the forward virtual Compton amplitude $`Im𝒯_{\gamma ^{}p\gamma ^{}p}`$, deeply virtual Compton scattering allows the measurement of the phase and spin structure of proton matrix elements for general momentum transfer squared $`t`$. In addition, the interference of the virtual Compton amplitude and Bethe-Heitler wide angle scattering Bremsstrahlung amplitude where the photon is emitted from the lepton line leads to an electron-positron asymmetry in the $`e^\pm pe^\pm \gamma p`$ cross section which is proportional to the real part of the Compton amplitude . The deeply virtual Compton amplitude $`\gamma ^{}p\gamma p`$ is related by crossing to another important process $`\gamma ^{}\gamma `$ hadron pairs at fixed invariant mass which can be measured in electron-photon collisions .

To leading order in $`1/Q`$, the deeply virtual Compton scattering amplitude $`\gamma ^{}(q)p(P)\gamma (q^{})p(P^{})`$ factorizes as the convolution in $`x`$ of the amplitude $`t^{\mu \nu }`$ for hard Compton scattering on a quark line with the generalized Compton form factors $`H(x,t,\zeta ),`$ $`E(x,t,\zeta )`$, $`\stackrel{~}{H}(x,t,\zeta ),`$ and $`\stackrel{~}{E}(x,t,\zeta )`$ of the target proton . Here $`x`$ is the light-front momentum fraction of the struck quark, and $`\zeta =Q^2/2Pq`$ plays the role of the Bjorken variable. The square of the four-momentum transfer from the proton is given by $`t=\mathrm{\Delta }^2=2P\mathrm{\Delta }=\frac{(\zeta ^2M^2+\stackrel{}{\mathrm{\Delta }}_{}^{}{}_{}{}^{2})}{(1\zeta )}`$ , where $`\mathrm{\Delta }`$ is the difference of initial and final momenta of the proton ($`P=P^{}+\mathrm{\Delta }`$). We will be interested in deeply virtual Compton scattering where $`q^2`$ is large compared to the masses and $`t`$. Then, to leading order in $`1/Q^2`$, $`\frac{q^2}{2P_Iq}=\zeta .`$ Thus $`\zeta `$ plays the role of the Bjorken variable in deeply virtual Compton scattering. For a fixed value of $`t`$, the allowed range of $`\zeta `$ is given by

$$0\zeta \frac{(t)}{2M^2}\left(\sqrt{1+\frac{4M^2}{(t)}}1\right).$$
(22)
The form factor $`H(x,t,\zeta )`$ describes the proton response when the helicity of the proton is unchanged, and $`E(x,t,\zeta )`$ is for the case when the proton helicity is flipped. Two additional functions $`\stackrel{~}{H}(x,t,\zeta ),`$ and $`\stackrel{~}{E}(x,t,\zeta )`$ appear, corresponding to the dependence of the Compton amplitude on quark helicity.

Recently, Markus Diehl, Dae Sung Hwang and I  have shown how the deeply virtual Compton amplitude can be evaluated explicitly in the Fock state representation using the matrix elements of the currents and the boost properties of the light-front wavefunctions. For the $`nn`$ diagonal term ($`\mathrm{\Delta }n=0`$), the arguments of the final-state hadron wavefunction are $`\frac{x_1\zeta }{1\zeta }`$, $`\stackrel{}{k}_1\frac{1x_1}{1\zeta }\stackrel{}{\mathrm{\Delta }}_{}`$ for the struck quark and $`\frac{x_i}{1\zeta }`$, $`\stackrel{}{k}_i+\frac{x_i}{1\zeta }\stackrel{}{\mathrm{\Delta }}_{}`$ for the $`n1`$ spectators. We thus obtain formulae for the diagonal (parton-number-conserving) contribution to the generalized form factors for deeply virtual Compton amplitude in the domain  $`\zeta x_11`$:

$`\sqrt{1\zeta }f_{1(nn)}(x_1,t,\zeta ){\displaystyle \frac{\zeta ^2}{4\sqrt{1\zeta }}}f_{2(nn)}(x_1,t,\zeta )`$ (23)
$`=`$ $`{\displaystyle \underset{n,\lambda }{}}{\displaystyle \underset{i=1}{\overset{n}{}}}{\displaystyle _0^1}𝑑x_{i(i1)}{\displaystyle \frac{d^2\stackrel{}{k}_i}{2(2\pi )^3}\delta \left(1\underset{j=1}{\overset{n}{}}x_j\right)\delta ^{(2)}\left(\underset{j=1}{\overset{n}{}}\stackrel{}{k}_j\right)}`$
$`\times \psi _{(n)}^{}(x_i^{},\stackrel{}{k}_i^{},\lambda _i)\psi _{(n)}^{}(x_i,\stackrel{}{k}_i,\lambda _i)(\sqrt{1\zeta })^{1n},`$
$`\sqrt{1\zeta }\left(\mathrm{\hspace{0.17em}1}+{\displaystyle \frac{\zeta }{2(1\zeta )}}\right){\displaystyle \frac{(\mathrm{\Delta }^1\mathrm{i}\mathrm{\Delta }^2)}{2M}}f_{2(nn)}(x_1,t,\zeta )`$ (24)
$`=`$ $`{\displaystyle \underset{n,\lambda }{}}{\displaystyle \underset{i=1}{\overset{n}{}}}{\displaystyle _0^1}𝑑x_{i(i1)}{\displaystyle \frac{d^2\stackrel{}{k}_i}{2(2\pi )^3}\delta \left(1\underset{j=1}{\overset{n}{}}x_j\right)\delta ^{(2)}\left(\underset{j=1}{\overset{n}{}}\stackrel{}{k}_j\right)}`$
$`\times \psi _{(n)}^{}(x_i^{},\stackrel{}{k}_i^{},\lambda _i)\psi _{(n)}^{}(x_i,\stackrel{}{k}_i,\lambda _i)(\sqrt{1\zeta })^{1n},`$
where

$$\{\begin{array}{ccc}x_1^{}=\frac{x_1\zeta }{1\zeta },\hfill & \stackrel{}{k}_1^{}=\stackrel{}{k}_1\frac{1x_1}{1\zeta }\stackrel{}{\mathrm{\Delta }}_{}\hfill & \text{for the struck quark,}\hfill \\ x_i^{}=\frac{x_i}{1\zeta },\hfill & \stackrel{}{k}_i^{}=\stackrel{}{k}_i+\frac{x_i}{1\zeta }\stackrel{}{\mathrm{\Delta }}_{}\hfill & \text{for the }(n1)\text{ spectators.}\hfill \end{array}$$
(25)
A sum over all possible helicities $`\lambda _i`$ is understood. If quark masses are neglected, the currents conserve helicity. We also can check that $`_{i=1}^nx_i^{}=1`$, $`_{i=1}^n\stackrel{}{k}_i^{}=\stackrel{}{0}_{}`$.

For the $`n+1n1`$ off-diagonal term ($`\mathrm{\Delta }n=2`$), consider the case where partons $`1`$ and $`n+1`$ of the initial wavefunction annihilate into the current leaving $`n1`$ spectators. Then $`x_{n+1}=\zeta x_1`$, $`\stackrel{}{k}_{n+1}=\stackrel{}{\mathrm{\Delta }}_{}\stackrel{}{k}_1`$. The remaining $`n1`$ partons have total momentum $`((1\zeta )P^+,\stackrel{}{\mathrm{\Delta }}_{})`$. The final wavefunction then has arguments $`x_i^{}=\frac{x_i}{1\zeta }`$ and $`\stackrel{}{k}_i^{}=\stackrel{}{k}_i+\frac{x_i}{1\zeta }\stackrel{}{\mathrm{\Delta }}_{}`$. We thus obtain the formulae for the off-diagonal matrix element of the Compton amplitude in the domain $`0x_1\zeta `$:

$`\sqrt{1\zeta }f_{1(n+1n1)}(x_1,t,\zeta ){\displaystyle \frac{\zeta ^2}{4\sqrt{1\zeta }}}f_{2(n+1n1)}(x_1,t,\zeta )`$ (26)
$`=`$ $`{\displaystyle \underset{n,\lambda }{}}{\displaystyle _0^1}𝑑x_{n+1}{\displaystyle \frac{d^2\stackrel{}{k}_1}{2(2\pi )^3}\frac{d^2\stackrel{}{k}_{n+1}}{2(2\pi )^3}\underset{i=2}{\overset{n}{}}_0^1𝑑x_i\frac{d^2\stackrel{}{k}_i}{2(2\pi )^3}}`$
$`\times \delta \left(1{\displaystyle \underset{j=1}{\overset{n+1}{}}}x_j\right)\delta ^{(2)}\left({\displaystyle \underset{j=1}{\overset{n+1}{}}}\stackrel{}{k}_j\right)[\sqrt{1\zeta }]^{1n}`$
$`\times \psi _{(n1)}^{}(x_i^{},\stackrel{}{k}_i^{},\lambda _i)\psi _{(n+1)}^{}(\{x_1,x_i,x_{n+1}=\zeta x_1\},`$
$`\{\stackrel{}{k}_1,\stackrel{}{k}_i,\stackrel{}{k}_{n+1}=\stackrel{}{\mathrm{\Delta }}_{}\stackrel{}{k}_1\},\{\lambda _1,\lambda _i,\lambda _{n+1}=\lambda _1\}),`$
$`\sqrt{1\zeta }\left(\mathrm{\hspace{0.17em}1}+{\displaystyle \frac{\zeta }{2(1\zeta )}}\right){\displaystyle \frac{(\mathrm{\Delta }^1\mathrm{i}\mathrm{\Delta }^2)}{2M}}f_{2(n+1n1)}(x_1,t,\zeta )`$ (27)
$`=`$ $`{\displaystyle \underset{n,\lambda }{}}{\displaystyle _0^1}𝑑x_{n+1}{\displaystyle \frac{d^2\stackrel{}{k}_1}{2(2\pi )^3}\frac{d^2\stackrel{}{k}_{n+1}}{2(2\pi )^3}\underset{i=2}{\overset{n}{}}_0^1𝑑x_i\frac{d^2\stackrel{}{k}_i}{2(2\pi )^3}}`$
$`\times \delta \left(1{\displaystyle \underset{j=1}{\overset{n+1}{}}}x_j\right)\delta ^{(2)}\left({\displaystyle \underset{j=1}{\overset{n+1}{}}}\stackrel{}{k}_j\right)[\sqrt{1\zeta }]^{1n}`$
$`\times \psi _{(n1)}^{}(x_i^{},\stackrel{}{k}_i^{},\lambda _i)\psi _{(n+1)}^{}(\{x_1,x_i,x_{n+1}=\zeta x_1\},`$
$`\{\stackrel{}{k}_1,\stackrel{}{k}_i,\stackrel{}{k}_{n+1}=\stackrel{}{\mathrm{\Delta }}_{}\stackrel{}{k}_1\},\{\lambda _1,\lambda _i,\lambda _{n+1}=\lambda _1\}),`$
where $`i=2,3,\mathrm{},n`$ label the $`n1`$ spectator partons which appear in the final-state hadron wavefunction with

$$x_i^{}=\frac{x_i}{1\zeta },\stackrel{}{k}_i^{}=\stackrel{}{k}_i+\frac{x_i}{1\zeta }\stackrel{}{\mathrm{\Delta }}_{}.$$
(28)
We can again check that the arguments of the final-state wavefunction satisfy $`_{i=2}^nx_i^{}=1`$, $`_{i=2}^n\stackrel{}{k}_i^{}=\stackrel{}{0}_{}`$.

The above representation is the general form for the generalized form factors of the deeply virtual Compton amplitude for any composite system. Thus given the light-front Fock state wavefunctions of the eigensolutions of the light-front Hamiltonian, we can compute the amplitude for virtual Compton scattering including all spin correlations. The formulae are accurate to leading order in $`1/Q^2`$. Radiative corrections to the quark Compton amplitude of order $`\alpha _s(Q^2)`$ from diagrams in which a hard gluon interacts between the two photons have also been neglected.

## 5 Applications of QCD Factorization to Hard QCD Processes

Factorization theorems for hard exclusive, semi-exclusive, and diffractive processes allow the separation of soft non-perturbative dynamics of the bound state hadrons from the hard dynamics of a perturbatively-calculable quark-gluon scattering amplitude. The factorization of inclusive reactions is reviewed in ref. For reviews and bibliography of exclusive process calculations in QCD (see Ref. ).

The light-front formalism provides a physical factorization scheme which conveniently separates and factorizes soft non-perturbative physics from hard perturbative dynamics in both exclusive and inclusive reactions .

In hard inclusive reactions all intermediate states are divided according to $`_n^2<\mathrm{\Lambda }^2`$ and $`_n^2>\mathrm{\Lambda }^2`$ domains. The lower mass regime is associated with the quark and gluon distributions defined from the absolute squares of the LC wavefunctions in the light cone factorization scheme. In the high invariant mass regime, intrinsic transverse momenta can be ignored, so that the structure of the process at leading power has the form of hard scattering on collinear quark and gluon constituents, as in the parton model. The attachment of gluons from the LC wavefunction to a propagator in a hard subprocess is power-law suppressed in LC gauge, so that the minimal quark-gluon particle-number subprocesses dominate. It is then straightforward to derive the DGLAP equations from the evolution of the distributions with $`\mathrm{log}\mathrm{\Lambda }^2`$. The anomaly contribution to singlet helicity structure function $`g_1(x,Q)`$ can be explicitly identified in the LC factorization scheme as due to the $`\gamma ^{}gq\overline{q}`$ fusion process. The anomaly contribution would be zero if the gluon is on shell. However, if the off-shellness of the state is larger than the quark pair mass, one obtains the usual anomaly contribution .

In exclusive amplitudes, the LC wavefunctions are the interpolating amplitudes connecting the quark and gluons to the hadronic states. In an exclusive amplitude involving a hard scale $`Q^2`$ all intermediate states can be divided according to $`_n^2<\mathrm{\Lambda }^2<Q^2`$ and $`_n^2<\mathrm{\Lambda }^2`$ invariant mass domains. The high invariant mass contributions to the amplitude has the structure of a hard scattering process $`T_H`$ in which the hadrons are replaced by their respective (collinear) quarks and gluons. In light-cone gauge only the minimal Fock states contribute to the leading power-law fall-off of the exclusive amplitude. The wavefunctions in the lower invariant mass domain can be integrated up to an arbitrary intermediate invariant mass cutoff $`\mathrm{\Lambda }`$. The invariant mass domain beyond this cutoff is included in the hard scattering amplitude $`T_H`$. The $`T_H`$ satisfy dimensional counting rules . Final-state and initial state corrections from gluon attachments to lines connected to the color-singlet distribution amplitudes cancel at leading twist. Explicit examples of perturbative QCD factorization will be discussed in more detail in the next section.

The key non-perturbative input for exclusive processes is thus the gauge and frame independent hadron distribution amplitude  defined as the integral of the valence (lowest particle number) Fock wavefunction; e.g. for the pion

$$\varphi _\pi (x_i,\mathrm{\Lambda })d^2k_{}\psi _{q\overline{q}/\pi }^{(\mathrm{\Lambda })}(x_i,\stackrel{}{k}_i,\lambda )$$
(29)
where the global cutoff $`\mathrm{\Lambda }`$ is identified with the resolution $`Q`$. The distribution amplitude controls leading-twist exclusive amplitudes at high momentum transfer, and it can be related to the gauge-invariant Bethe-Salpeter wavefunction at equal light-cone time. The logarithmic evolution of hadron distribution amplitudes $`\varphi _H(x_i,Q)`$ can be derived from the perturbatively-computable tail of the valence light-front wavefunction in the high transverse momentum regime . The conformal basis for the evolution of the three-quark distribution amplitudes for the baryons  has recently been obtained by V. Braun et al..

The existence of an exact formalism provides a basis for systematic approximations and a control over neglected terms. For example, one can analyze exclusive semi-leptonic $`B`$-decays which involve hard internal momentum transfer using a perturbative QCD formalism  patterned after the perturbative analysis of form factors at large momentum transfer. The hard-scattering analysis again proceeds by writing each hadronic wavefunction as a sum of soft and hard contributions

$$\psi _n=\psi _n^{\mathrm{soft}}(_n^2<\mathrm{\Lambda }^2)+\psi _n^{\mathrm{hard}}(_n^2>\mathrm{\Lambda }^2),$$
(30)
where $`_n^2`$ is the invariant mass of the partons in the $`n`$-particle Fock state and $`\mathrm{\Lambda }`$ is the separation scale. The high internal momentum contributions to the wavefunction $`\psi _n^{\mathrm{hard}}`$ can be calculated systematically from QCD perturbation theory by iterating the gluon exchange kernel. The contributions from high momentum transfer exchange to the $`B`$-decay amplitude can then be written as a convolution of a hard-scattering quark-gluon scattering amplitude $`T_H`$ with the distribution amplitudes $`\varphi (x_i,\mathrm{\Lambda })`$, the valence wavefunctions obtained by integrating the constituent momenta up to the separation scale $`_n<\mathrm{\Lambda }<Q`$. Furthermore in processes such as $`B\pi D`$ where the pion is effectively produced as a rapidly-moving small Fock state with a small color-dipole interactions, final state interactions are suppressed by color transparency. This is the basis for the perturbative hard-scattering analyses . In a systematic analysis, one can identify the hard PQCD contribution as well as the soft contribution from the convolution of the light-front wavefunctions. Furthermore, the hard-scattering contribution can be systematically improved.

Given the solution for the hadronic wavefunctions $`\psi _n^{(\mathrm{\Lambda })}`$ with $`_n^2<\mathrm{\Lambda }^2`$, one can construct the wavefunction in the hard regime with $`_n^2>\mathrm{\Lambda }^2`$ using projection operator techniques. The construction can be done perturbatively in QCD since only high invariant mass, far off-shell matrix elements are involved. One can use this method to derive the physical properties of the LC wavefunctions and their matrix elements at high invariant mass. Since $`_n^2=_{i=1}^n\left(\frac{k_{}^2+m^2}{x}\right)_i`$, this method also allows the derivation of the asymptotic behavior of light-front wavefunctions at large $`k_{}`$, which in turn leads to predictions for the fall-off of form factors and other exclusive matrix elements at large momentum transfer, such as the quark counting rules for predicting the nominal power-law fall-off of two-body scattering amplitudes at fixed $`\theta _{cm}`$  and helicity selection rules . The phenomenological successes of these rules can be understood within QCD if the coupling $`\alpha _V(Q)`$ freezes in a range of relatively small momentum transfer .

## 6 Two-Photon Processes

The simplest and perhaps the most elegant illustration of an exclusive reaction in QCD is the evaluation of the photon-to-pion transition form factor $`F_{\gamma \pi }(Q^2)`$  which is measurable in single-tagged two-photon $`eeee\pi ^0`$ reactions. The form factor is defined via the invariant amplitude $`\mathrm{\Gamma }^\mu =ie^2F_{\pi \gamma }(Q^2)ϵ^{\mu \nu \rho \sigma }p_\nu ^\pi ϵ_\rho q_\sigma .`$ As in inclusive reactions, one must specify a factorization scheme which divides the integration regions of the loop integrals into hard and soft momenta, compared to the resolution scale $`\stackrel{~}{Q}`$. At leading twist, the transition form factor then factorizes as a convolution of the $`\gamma ^{}\gamma q\overline{q}`$ amplitude (where the quarks are collinear with the final state pion) with the valence light-front wavefunction of the pion:

$$F_{\gamma M}(Q^2)=\frac{4}{\sqrt{3}}_0^1𝑑x\varphi _M(x,\stackrel{~}{Q})T_{\gamma M}^H(x,Q^2).$$
(31)
The hard scattering amplitude for $`\gamma \gamma ^{}q\overline{q}`$ is $`T_{\gamma M}^H(x,Q^2)=[(1x)Q^2]^1\times \left(1+𝒪(\alpha _s)\right).`$ The leading QCD corrections have been computed by Braaten . The evaluation of the next-to-leading corrections in the physical $`\alpha _V`$ scheme is given in Ref. . For the asymptotic distribution amplitude $`\varphi _\pi ^{\mathrm{asympt}}(x)=\sqrt{3}f_\pi x(1x)`$ one predicts $`Q^2F_{\gamma \pi }(Q^2)=2f_\pi \left(1\frac{5}{3}\frac{\alpha _V(Q^{})}{\pi }\right)`$ where $`Q^{}=e^{3/2}Q`$ is the BLM scale for the pion form factor. The PQCD predictions have been tested in measurements of $`e\gamma e\pi ^0`$ by the CLEO collaboration . See Fig. 4 (b). The observed flat scaling of the $`Q^2F_{\gamma \pi }(Q^2)`$ data from $`Q^2=2`$ to $`Q^2=8`$ GeV<sup>2</sup> provides an important confirmation of the applicability of leading twist QCD to this process. The magnitude of $`Q^2F_{\gamma \pi }(Q^2)`$ is remarkably consistent with the predicted form, assuming the asymptotic distribution amplitude and including the LO QCD radiative correction with $`\alpha _V(e^{3/2}Q)/\pi 0.12`$. One could allow for some broadening of the distribution amplitude with a corresponding increase in the value of $`\alpha _V`$ at small scales. Radyushkin , Ong , and Kroll  have also noted that the scaling and normalization of the photon-to-pion transition form factor tends to favor the asymptotic form for the pion distribution amplitude and rules out broader distributions such as the two-humped form suggested by QCD sum rules .

The two-photon annihilation process $`\gamma ^{}\gamma `$ hadrons, which is measurable in single-tagged $`e^+e^{}e^+e^{}\mathrm{hadrons}`$ events, provides a semi-local probe of $`C=+`$ hadron systems $`\pi ^0,\eta ^0,\eta ^{},\eta _c,\pi ^+\pi ^{}`$, etc. The $`\gamma ^{}\gamma \pi ^+\pi ^{}`$ hadron pair process is related to virtual Compton scattering on a pion target by crossing. The leading twist amplitude is sensitive to the $`1/x1/(1x)`$ moment of the two-pion distribution amplitude coupled to two valence quarks .

Two-photon reactions, $`\gamma \gamma H\overline{H}`$ at large s = $`(k_1+k_2)^2`$ and fixed $`\theta _{\mathrm{cm}}`$, provide a particularly important laboratory for testing QCD since these cross-channel “Compton” processes are the simplest calculable large-angle exclusive hadronic scattering reactions. The helicity structure, and often even the absolute normalization can be rigorously computed for each two-photon channel . In the case of meson pairs, dimensional counting predicts that for large $`s`$, $`s^4d\sigma /dt(\gamma \gamma M\overline{M}`$ scales at fixed $`t/s`$ or $`\theta _{\mathrm{c}.\mathrm{m}.}`$ up to factors of $`\mathrm{ln}s/\mathrm{\Lambda }^2`$. The angular dependence of the $`\gamma \gamma H\overline{H}`$ amplitudes can be used to determine the shape of the process-independent distribution amplitudes, $`\varphi _H(x,Q)`$. An important feature of the $`\gamma \gamma M\overline{M}`$ amplitude for meson pairs is that the contributions of Landshoff pitch singularities are power-law suppressed at the Born level – even before taking into account Sudakov form factor suppression. There are also no anomalous contributions from the $`x1`$ endpoint integration region. Thus, as in the calculation of the meson form factors, each fixed-angle helicity amplitude can be written to leading order in $`1/Q`$ in the factorized form $`[Q^2=p_T^2=tu/s;\stackrel{~}{Q}_x=\mathrm{min}(xQ,(lx)Q)]`$:

$$_{\gamma \gamma M\overline{M}}=_0^1dx_0^1dy\varphi _{\overline{M}}(y,\stackrel{~}{Q}_y)T_H(x,y,s,\theta _{\mathrm{c}.\mathrm{m}.}\varphi _M(x,\stackrel{~}{Q}_x),$$
(32)
where $`T_H`$ is the hard-scattering amplitude $`\gamma \gamma (q\overline{q})(q\overline{q})`$ for the production of the valence quarks collinear with each meson, and $`\varphi _M(x,\stackrel{~}{Q})`$ is the amplitude for finding the valence $`q`$ and $`\overline{q}`$ with light-front fractions of the meson’s momentum, integrated over transverse momenta $`k_{}<\stackrel{~}{Q}.`$ The contribution of non-valence Fock states are power-law suppressed. Furthermore, the helicity-selection rules  of perturbative QCD predict that vector mesons are produced with opposite helicities to leading order in $`1/Q`$ and all orders in $`\alpha _s`$. The dependence in $`x`$ and $`y`$ of several terms in $`T_{\lambda ,\lambda ^{}}`$ is quite similar to that appearing in the meson’s electromagnetic form factor. Thus much of the dependence on $`\varphi _M(x,Q)`$ can be eliminated by expressing it in terms of the meson form factor. In fact, the ratio of the $`\gamma \gamma \pi ^+\pi ^{}`$ and $`e^+e^{}\mu ^+\mu ^{}`$ amplitudes at large $`s`$ and fixed $`\theta _{CM}`$ is nearly insensitive to the running coupling and the shape of the pion distribution amplitude:

$$\frac{\frac{d\sigma }{dt}(\gamma \gamma \pi ^+\pi ^{})}{\frac{d\sigma }{dt}(\gamma \gamma \mu ^+\mu ^{})}\frac{4|F_\pi (s)|^2}{1\mathrm{cos}^2\theta _{\mathrm{c}.\mathrm{m}.}}.$$
(33)
The comparison of the PQCD prediction for the sum of $`\pi ^+\pi ^{}`$ plus $`K^+K^{}`$ channels with recent CLEO data  is shown in Fig. 3. The CLEO data for charged pion and kaon pairs show a clear transition to the scaling and angular distribution predicted by PQCD  for $`W=\sqrt{(}s_{\gamma \gamma }>2`$ GeV. See Fig. 3. It is clearly important to measure the magnitude and angular dependence of the two-photon production of neutral pions and $`\rho ^+\rho ^{}`$ cross sections in view of the strong sensitivity of these channels to the shape of meson distribution amplitudes. QCD also predicts that the production cross section for charged $`\rho `$-pairs (with any helicity) is much larger that for that of neutral $`\rho `$ pairs, particularly at large $`\theta _{\mathrm{c}.\mathrm{m}.}`$ angles. Similar predictions are possible for other helicity-zero mesons. The cross sections for Compton scattering on protons and the crossed reaction $`\gamma \gamma p\overline{p}`$ at high momentum transfer have also been evaluated , providing important tests of the proton distribution amplitude.

It is particularly compelling to see a transition in angular dependence between the low energy chiral and PQCD regimes. The success of leading-twist perturbative QCD scaling for exclusive processes at presently experimentally accessible momentum transfer can be understood if the effective coupling $`\alpha _V(Q^{})`$ is approximately constant at the relatively small scales $`Q^{}`$ relevant to the hard scattering amplitudes . The evolution of the quark distribution amplitudes In the low-$`Q^{}`$ domain at also needs to be minimal. Sudakov suppression of the endpoint contributions is also strengthened if the coupling is frozen because of the exponentiation of a double logarithmic series.

A debate has continued  on whether processes such as the pion and proton form factors and elastic Compton scattering $`\gamma p\gamma p`$ might be dominated by higher-twist mechanisms until very large momentum transfer. If one assumes that the light-front wavefunction of the pion has the form $`\psi _{\mathrm{soft}}(x,k_{})=A\mathrm{exp}(b\frac{k_{}^2}{x(1x)})`$, then the Feynman endpoint contribution to the overlap integral at small $`k_{}`$ and $`x1`$ will dominate the form factor compared to the hard-scattering contribution until very large $`Q^2`$. However, this ansatz for $`\psi _{\mathrm{soft}}(x,k_{})`$ has no suppression at $`k_{}=0`$ for any $`x`$; i.e., the wavefunction in the hadron rest frame does not fall-off at all for $`k_{}=0`$ and $`k_z\mathrm{}`$. Thus such wavefunctions do not represent well soft QCD contributions. Endpoint contributions are also suppressed by the QCD Sudakov form factor, reflecting the fact that a near-on-shell quark must radiate if it absorbs large momentum. One can show  that the leading power dependence of the two-particle light-front Fock wavefunction in the endpoint region is $`1x`$, giving a meson structure function which falls as $`(1x)^2`$ and thus by duality a non-leading contribution to the meson form factor $`F(Q^2)1/Q^3`$. Thus the dominant contribution to meson form factors comes from the hard-scattering regime.

Radyushkin  has argued that the Compton amplitude is dominated by soft end-point contributions of the proton wavefunctions where the two photons both interact on a quark line carrying nearly all of the proton’s momentum. This description appears to agree with the Compton data at least at forward angles where $`t<10`$ GeV<sup>2</sup>. From this viewpoint, the dominance of the factorizable PQCD leading twist contributions requires momentum transfers much higher than those currently available. However, the endpoint model cannot explain the empirical success of the perturbative QCD fixed $`\theta _{c.m.}`$ scaling $`s^7d\sigma /dt(\gamma p\pi ^+n)\mathrm{const}`$ at relatively low momentum transfer in pion photoproduction .

Clearly much more experimental input on hadron wavefunctions is needed, particularly from measurements of two-photon exclusive reactions into meson and baryon pairs at the high luminosity $`B`$ factories. For example, the ratio $`\frac{d\sigma }{dt}(\gamma \gamma \pi ^0\pi ^0)/\frac{d\sigma }{dt}(\gamma \gamma \pi ^+\pi ^{})`$ is particularly sensitive to the shape of pion distribution amplitude. Baryon pair production in two-photon reactions at threshold may reveal physics associated with the soliton structure of baryons in QCD . In addition, fixed target experiments can provide much more information on fundamental QCD processes such as deeply virtual Compton scattering and large angle Compton scattering.

## 7 Self-Resolved Diffractive Reactions and Light Cone Wavefunctions

Diffractive multi-jet production in heavy nuclei provides a novel way to measure the shape of the LC Fock state wavefunctions and test color transparency. For example, consider the reaction  $`\pi A\mathrm{Jet}_1+\mathrm{Jet}_2+A^{}`$ at high energy where the nucleus $`A^{}`$ is left intact in its ground state. The transverse momenta of the jets balance so that $`\stackrel{}{k}_i+\stackrel{}{k}_2=\stackrel{}{q}_{}<R_{}^{1}{}_{A}{}^{}.`$ The light-front longitudinal momentum fractions also need to add to $`x_1+x_21`$ so that $`\mathrm{\Delta }p_L<R_A^1`$. The process can then occur coherently in the nucleus.

Because of color transparency, the valence wavefunction of the pion with small impact separation, will penetrate the nucleus with minimal interactions, diffracting into jet pairs . The $`x_1=x`$, $`x_2=1x`$ dependence of the di-jet distributions will thus reflect the shape of the pion valence light-front wavefunction in $`x`$; similarly, the $`\stackrel{}{k}_1\stackrel{}{k}_2`$ relative transverse momenta of the jets gives key information on the derivative of the underlying shape of the valence pion wavefunction . The diffractive nuclear amplitude extrapolated to $`t=0`$ should be linear in nuclear number $`A`$ if color transparency is correct. The integrated diffractive rate should then scale as $`A^2/R_A^2A^{4/3}`$. Preliminary results on a diffractive dissociation experiment of this type E791 at Fermilab using 500 GeV incident pions on nuclear targets  appear to be consistent with color transparency. The measured longitudinal momentum distribution of the jets  is consistent with a pion light-cone wavefunction of the pion with the shape of the asymptotic distribution amplitude, $`\varphi _\pi ^{\mathrm{asympt}}(x)=\sqrt{3}f_\pi x(1x)`$. Data from CLEO  for the $`\gamma \gamma ^{}\pi ^0`$ transition form factor also favor a form for the pion distribution amplitude close to the asymptotic solution to the perturbative QCD evolution equation .

The diffractive dissociation of a hadron or nucleus can also occur via the Coulomb dissociation of a beam particle on an electron beam (e.g. at HERA or eRHIC) or on the strong Coulomb field of a heavy nucleus (e.g. at RHIC or nuclear collisions at the LHC) . The amplitude for Coulomb exchange at small momentum transfer is proportional to the first derivative $`_ie_i\frac{}{\stackrel{}{k}_{Ti}}\psi `$ of the light-front wavefunction, summed over the charged constituents. The Coulomb exchange reactions fall off less fast at high transverse momentum compared to pomeron exchange reactions since the light-front wavefunction is effective differentiated twice in two-gluon exchange reactions.

It will also be interesting to study diffractive tri-jet production using proton beams $`pA\mathrm{Jet}_1+\mathrm{Jet}_2+\mathrm{Jet}_3+A^{}`$ to determine the fundamental shape of the 3-quark structure of the valence light-front wavefunction of the nucleon at small transverse separation . For example, consider the Coulomb dissociation of a high energy proton at HERA. The proton can dissociate into three jets corresponding to the three-quark structure of the valence light-front wavefunction. We can demand that the produced hadrons all fall outside an opening angle $`\theta `$ in the proton’s fragmentation region. Effectively all of the light-front momentum $`_jx_j1`$ of the proton’s fragments will thus be produced outside an “exclusion cone”. This then limits the invariant mass of the contributing Fock state $`_n^2>\mathrm{\Lambda }^2=P^{+2}\mathrm{sin}^2\theta /4`$ from below, so that perturbative QCD counting rules can predict the fall-off in the jet system invariant mass $``$. At large invariant mass one expects the three-quark valence Fock state of the proton to dominate. The segmentation of the forward detector in azimuthal angle $`\varphi `$ can be used to identify structure and correlations associated with the three-quark light-front wavefunction . An interesting possibility is that the distribution amplitude of the $`\mathrm{\Delta }(1232)`$ for $`J_z=1/2,3/2`$ is close to the asymptotic form $`x_1x_2x_3`$, but that the proton distribution amplitude is more complex. This ansatz can also be motivated by assuming a quark-diquark structure of the baryon wavefunctions. The differences in shapes of the distribution amplitudes could explain why the $`p\mathrm{\Delta }`$ transition form factor appears to fall faster at large $`Q^2`$ than the elastic $`pp`$ and the other $`pN^{}`$ transition form factors . One can use also measure the dijet structure of real and virtual photons beams $`\gamma ^{}A\mathrm{Jet}_1+\mathrm{Jet}_2+A^{}`$ to measure the shape of the light-front wavefunction for transversely-polarized and longitudinally-polarized virtual photons. Such experiments will open up a direct window on the amplitude structure of hadrons at short distances. The light-front formalism is also applicable to the description of nuclei in terms of their nucleonic and mesonic degrees of freedom . Self-resolving diffractive jet reactions in high energy electron-nucleus collisions and hadron-nucleus collisions at moderate momentum transfers can thus be used to resolve the light-front wavefunctions of nuclei.

## 8 Higher Fock States and the Intrinsic Sea

One can identify two contributions to the heavy quark sea, the “extrinsic” contributions which correspond to ordinary gluon splitting, and the “intrinsic” sea which is multi-connected via gluons to the valence quarks. The leading $`1/m_Q^2`$ contributions to the intrinsic sea of the proton in the heavy quark expansion are proton matrix elements of the operator  $`\eta ^\mu \eta ^\nu G_{\alpha \mu }G_{\beta \nu }G^{\alpha \beta }`$ which in light-cone gauge $`\eta ^\mu A_\mu =A^+=0`$ corresponds to three or four gluon exchange between the heavy-quark loop and the proton constituents in the forward virtual Compton amplitude. The intrinsic sea is thus sensitive to the hadronic bound-state structure . The maximal contribution of the intrinsic heavy quark occurs at $`x_Qm_Q/_im_{}`$ where $`m_{}=\sqrt{m^2+k_{}^2}`$; i.e. at large $`x_Q`$, since this minimizes the invariant mass $`_n^2`$. The measurements of the charm structure function by the EMC experiment are consistent with intrinsic charm at large $`x`$ in the nucleon with a probability of order $`0.6\pm 0.3\%`$  which is consistent with the recent estimates based on instanton fluctuations .

Chang and Hou  have recently discussed the consequences of intrinsic charm in heavy quark states such as the $`B`$, $`\mathrm{\Lambda }_B`$, and $`\mathrm{{\rm Y}}`$, such as an anomalous momentum distribution for $`BJ/\psi X`$. The characteristic momenta characterizing the $`B`$ meson is most likely higher by a factor of 2 compared to the momentum scale of light mesons, This effect is analogous to the higher momentum scale of muonium $`\mu ^+e^{}`$ versus that of positronium $`e^+e^{}`$ in atomic physics because of the larger reduced mass. Thus one can expect a higher probability for intrinsic charm in heavy hadrons compared to light hadrons.

One can also distinguish “intrinsic gluons”  which are associated with multi-quark interactions and extrinsic gluon contributions associated with quark substructure. One can also use this framework to isolate the physics of the anomaly contribution to the Ellis-Jaffe sum rule . Thus neither gluons nor sea quarks are solely generated by DGLAP evolution, and one cannot define a resolution scale $`Q_0`$ where the sea or gluon degrees of freedom can be neglected.

It is usually assumed that a heavy quarkonium state such as the $`J/\psi `$ always decays to light hadrons via the annihilation of its heavy quark constituents to gluons. However, as Karliner and I  have shown, the transition $`J/\psi \rho \pi `$ can also occur by the rearrangement of the $`c\overline{c}`$ from the $`J/\psi `$ into the $`|q\overline{q}c\overline{c}`$ intrinsic charm Fock state of the $`\rho `$ or $`\pi `$. On the other hand, the overlap rearrangement integral in the decay $`\psi ^{}\rho \pi `$ will be suppressed since the intrinsic charm Fock state radial wavefunction of the light hadrons will evidently not have nodes in its radial wavefunction. This observation provides a natural explanation of the long-standing puzzle  why the $`J/\psi `$ decays prominently to two-body pseudoscalar-vector final states, breaking hadron helicity conservation , whereas the $`\psi ^{}`$ does not.

The higher Fock state of the proton $`|uuds\overline{s}`$ should resemble a $`|K\mathrm{\Lambda }`$ intermediate state, since this minimizes its invariant mass $``$. In such a state, the strange quark has a higher mean momentum fraction $`x`$ than the $`\overline{s}`$ . Similarly, the helicity of the intrinsic strange quark in this configuration will be anti-aligned with the helicity of the nucleon . This $`Q\overline{Q}`$ asymmetry is a striking feature of the intrinsic heavy-quark sea.

## 9 Non-Perturbative Solutions of Light-Front Quantized QCD

Is there any hope of computing light-front wavefunctions from first principles? The solution of the light-front Hamiltonian equation $`H_{LC}^{QCD}|\mathrm{\Psi }=M^2|\mathrm{\Psi }`$ is an eigenvalue problem which in principle determines the masses squared of the entire bound and continuum spectrum of QCD. If one introduces periodic or anti-periodic boundary conditions, the eigenvalue problem is reduced to the diagonalization of a discrete Hermitian matrix representation of $`H_{LC}^{QCD}.`$ The light-front momenta satisfy $`x^+=\frac{2\pi }{L}n_i`$ and $`P^+=\frac{2\pi }{L}K`$, where $`_in_i=K.`$ The number of quanta in the contributing Fock states is restricted by the choice of harmonic resolution. A cutoff on the invariant mass of the Fock states truncates the size of the matrix representation in the transverse momenta. This is the essence of the DLCQ method , which has now become a standard tool for solving both the spectrum and light-front wavefunctions of one-space one-time theories – virtually any $`1+1`$ quantum field theory, including “reduced QCD” (which has both quark and gluonic degrees of freedom) can be completely solved using DLCQ . The method yields not only the bound-state and continuum spectrum, but also the light-front wavefunction for each eigensolution .

In the case of theories in 3+1 dimensions, Hiller, McCartor, and I  have recently shown that the use of covariant Pauli-Villars regularization with DLCQ allows one to obtain the spectrum and light-front wavefunctions of simplified theories, such as (3+1) Yukawa theory. Dalley et al. have shown how one can use DLCQ in one space-one time, with a transverse lattice to solve mesonic and gluonic states in $`3+1`$ QCD . The spectrum obtained for gluonium states is in remarkable agreement with lattice gauge theory results, but with a huge reduction of numerical effort. Hiller and I  have shown how one can use DLCQ to compute the electron magnetic moment in QED without resort to perturbation theory.

One can also formulate DLCQ so that supersymmetry is exactly preserved in the discrete approximation, thus combining the power of DLCQ with the beauty of supersymmetry . The “SDLCQ” method has been applied to several interesting supersymmetric theories, to the analysis of zero modes, vacuum degeneracy, massless states, mass gaps, and theories in higher dimensions, and even tests of the Maldacena conjecture . Broken supersymmetry is interesting in DLCQ, since it may serve as a method for regulating non-Abelian theories .

There are also many possibilities for obtaining approximate solutions of light-front wavefunctions in QCD. QCD sum rules, lattice gauge theory moments, and QCD inspired models such as the bag model, chiral theories, provide important constraints. Guides to the exact behavior of LC wavefunctions in QCD can also be obtained from analytic or DLCQ solutions to toy models such as “reduced” $`QCD(1+1).`$ The light-front and many-body Schrödinger theory formalisms must match In the nonrelativistic limit.

It would be interesting to see if light-front wavefunctions can incorporate chiral constraints such as soliton (Skyrmion) behavior for baryons and other consequences of the chiral limit in the soft momentum regime. Solvable theories such as $`QCD(1+1)`$ are also useful for understanding such phenomena. It has been shown that the anomaly contribution for the $`\pi ^0\gamma \gamma `$ decay amplitude is satisfied by the light-front Fock formalism in the limit where the mass of the pion is light compared to its size .

## 10 Non-Perturbative Calculations of the Pion Distribution Amplitude

The distribution amplitude $`\varphi (x,\stackrel{~}{Q})`$ can be computed from the integral over transverse momenta of the renormalized hadron valence wavefunction in the light-cone gauge at fixed light-cone time :

$$\varphi (x,\stackrel{~}{Q})=d^2\stackrel{}{k_{}}\theta \left(\stackrel{~}{Q}^2\frac{\stackrel{}{k_{}}^2}{x(1x)}\right)\psi ^{(\stackrel{~}{Q})}(x,\stackrel{}{k_{}}),$$
(34)
where a global cutoff in invariant mass is identified with the resolution $`\stackrel{~}{Q}`$. The distribution amplitude $`\varphi (x,\stackrel{~}{Q})`$ is boost and gauge invariant and evolves in $`\mathrm{ln}\stackrel{~}{Q}`$ through an evolution equation . Since it is formed from the same product of operators as the non-singlet structure function, the anomalous dimensions controlling $`\varphi (x,Q)`$ dependence in the ultraviolet $`\mathrm{log}Q`$ scale are the same as those which appear in the DGLAP evolution of structure functions . The decay $`\pi \mu \nu `$ normalizes the wave function at the origin: $`a_0/6=_0^1𝑑x\varphi (x,Q)=f_\pi /(2\sqrt{3}).`$ One can also compute the distribution amplitude from the gauge invariant Bethe-Salpeter wavefunction at equal light-cone time. This also allows contact with both QCD sum rules and lattice gauge theory; for example, moments of the pion distribution amplitudes have been computed in lattice gauge theory .

Dalley  has recently calculated the pion distribution amplitude from QCD using a combination of the discretized DLCQ method for the $`x^{}`$ and $`x^+`$ light-cone coordinates with the transverse lattice method  in the transverse directions, A finite lattice spacing $`a`$ can be used by choosing the parameters of the effective theory in a region of renormalization group stability to respect the required gauge, Poincaré, chiral, and continuum symmetries. The overall normalization gives $`f_\pi =101`$ MeV compared with the experimental value of $`93`$ MeV. Figure 4 (a) compares the resulting DLCQ/transverse lattice pion wavefunction with the best fit to the diffractive di-jet data (see the next section) after corrections for hadronization and experimental acceptance . The theoretical curve is somewhat broader than the experimental result. However, there are experimental uncertainties from hadronization and theoretical errors introduced from finite DLCQ resolution, using a nearly massless pion, ambiguities in setting the factorization scale $`Q^2`$, as well as errors in the evolution of the distribution amplitude from 1 to $`10\mathrm{GeV}^2`$. Instanton models also predict a pion distribution amplitude close to the asymptotic form . In contrast, recent lattice results from Del Debbio et al.  predict a much narrower shape for the pion distribution amplitude than the distribution predicted by the transverse lattice. A new result for the proton distribution amplitude treating nucleons as chiral solitons has recently been derived by Diakonov and Petrov . Dyson-Schwinger models  of hadronic Bethe-Salpeter wavefunctions can also be used to predict light-cone wavefunctions and hadron distribution amplitudes by integrating over the relative $`k^{}`$ momentum. There is also the possibility of deriving Bethe-Salpeter wavefunctions within light-cone gauge quantized QCD  in order to properly match to the light-cone gauge Fock state decomposition.

## 11 Conclusions

In these lectures I have shown how the universal, process-independent and frame-independent light-front Fock-state wavefunctions encode the properties of a hadron in terms of its fundamental quark and gluon degrees of freedom. Knowledge of such wavefunctions will be critical for progress in understanding exclusive $`B`$ decays.

I have shown how, given the proton’s light-front wavefunctions, one can compute not only the moments of the quark and gluon distributions measured in deep inelastic lepton-proton scattering, but also the multi-parton correlations which control the distribution of particles in the proton fragmentation region and dynamical higher twist effects. Light-front wavefunctions also provide a systematic framework for evaluating exclusive hadronic matrix elements, including time-like heavy hadron decay amplitudes, form factors, and the generalized form factors that appear in deeply virtual Compton scattering. The light-front Hamiltonian formalism also provides a physical factorization scheme for separating hard and soft contributions in both exclusive and inclusive hard processes.

The leading-twist QCD predictions for exclusive two-photon processes such as the photon-to-pion transition form factor and $`\gamma \gamma `$ hadron pairs are based on rigorous factorization theorems. The recent data from the CLEO collaboration on $`F_{\gamma \pi }(Q^2)`$ and the sum of $`\gamma \gamma \pi ^+\pi ^{}`$ and $`\gamma \gamma K^+K^{}`$ channels are in excellent agreement with the QCD predictions. It is particularly compelling to see a transition in angular dependence between the low energy chiral and PQCD regimes. The success of leading-twist perturbative QCD scaling for exclusive processes at presently experimentally accessible momentum transfer can be understood if the effective coupling $`\alpha _V(Q^{})`$ is approximately constant at the relatively small scales $`Q^{}`$ relevant to the hard scattering amplitudes . The evolution of the quark distribution amplitudes In the low-$`Q^{}`$ domain at also needs to be minimal. Sudakov suppression of the endpoint contributions is also strengthened if the coupling is frozen because of the exponentiation of a double logarithmic series.

One of the formidable challenges in QCD is the calculation of non-perturbative wavefunctions of hadrons from first principles. The recent calculation of the pion distribution amplitude by Dalley  using light-cone and transverse lattice methods is particularly encouraging. The predicted form of $`\varphi _\pi (x,Q)`$ is somewhat broader than but not inconsistent with the asymptotic form favored by the measured normalization of $`Q^2F_{\gamma \pi ^0}(Q^2)`$ and the pion wavefunction inferred from diffractive di-jet production.

Clearly much more experimental input on hadron wavefunctions is needed, particularly from measurements of two-photon exclusive reactions into meson and baryon pairs at the high luminosity $`B`$ factories. For example, the ratio $`\frac{d\sigma }{dt}(\gamma \gamma \pi ^0\pi ^0)/\frac{d\sigma }{dt}(\gamma \gamma \pi ^+\pi ^{})`$ is particularly sensitive to the shape of pion distribution amplitude. Baryon pair production in two-photon reactions at threshold may reveal physics associated with the soliton structure of baryons in QCD . In addition, fixed target experiments can provide much more information on fundamental QCD processes such as deeply virtual Compton scattering and large angle Compton scattering.

A remarkable new type of diffractive jet production reaction, “self-resolving diffractive interactions” can provide direct empirical information on the light-front wavefunctions of hadrons. The recent E791 experiment at Fermilab has not only determined the main features of the pion wavefunction, but has also confirmed color transparency, a fundamental test of the gauge properties of QCD. Analogous reaction involving nuclear projectiles can resolve the light-front wavefunctions of nuclei in terms of their nucleon and mesonic degrees of freedom. It is also possible to measure the light-front wavefunctions of atoms through high energy Coulomb dissociation.

There has been notable progress in computing light-front wavefunctions directly from the QCD light-front Hamiltonian, using DLCQ and transverse lattice methods. Even without full non-perturbative solutions of QCD, one can envision a program to construct the light-front wavefunctions using measured moments constraints from QCD sum rules, lattice gauge theory, and data from hard exclusive and inclusive processes. One can also be guided by theoretical constraints from perturbation theory which dictate the asymptotic form of the wavefunctions at large invariant mass, $`x1`$, and high $`k_{}`$. One can also use ladder relations which connect Fock states of different particle number; perturbatively-motivated numerator spin structures; conformal symmetry, guidance from toy models such as “reduced” $`QCD(1+1)`$; and the correspondence to Abelian theory for $`N_C0`$, as well as many-body Schrödinger theory in the nonrelativistic domain.

## Acknowledgments

Work supported by the Department of Energy under contract number DE-AC03-76SF00515. I wish to thank the organizers of this meeting, Hsiang-nan Li and Wei-Min Zhang, for their outstanding hospitality in Taiwan. Much of this work is based on collaborations, particularly with Markus Diehl, Paul Hoyer, Dae Sung Hwang, Peter Lepage, Bo-Qiang Ma, Hans Christian Pauli, Johan Rathsman, Ivan Schmidt, and Prem Srivastava.