
Forward photoproduction of $J/\psi$ can be used to extract Generalized Parton
Distributions(GPD's) of gluons. We analyze the process at twist3 level and
study relevant classifications of twist3 gluon GPD's. At leading power or
twist2 level the produced $J/\psi$ is transversely polarized. We find that at
twist3 the produced $J/\psi$ is longitudinally polarized. Our study shows that
in high energy limit the twist3 amplitude is only suppressed by the inverse
power of the heavy quark mass relatively to the twist2 amplitude. This
indicates that the power correction to the crosssection of unpolarized
$J/\psi$ can have a sizeable effect. We have also derived the amplitude of the
production of $h_c$ at twist3. The production of other quarkonia has been
briefly discussed.

We study single spin asymmetries at oneloop accuracy in semiinclusive DIS
with a transversely polarized hadron in the initial state. Two measurable spin
observables are predicted in the framework of QCD collinear factorization. One
of the spin observables is the Sivers weighted asymmetry, another one is the
Collins weighted asymmetry. The prediction takes a form of convolutions of
perturbative coefficient functions and nonperturbative functions, which are
twist2 transversity distributions, twist3 parton distributions, twist2 and
twist3 parton fragmentation functions. These nonperturbative functions can be
extracted from measurements of the spin observables and provide valuable
information of the inner structure of hadrons. The measurements can be done in
current COMPASS and JLab experiments and in future experiments of EIC. The
perturbative coefficient functions are calculated at oneloop level. There are
collinear divergences in the calculation involving chiralityeven and
chiralityodd twist operators. We find that all collinear divergences can be
correctly subtracted so that the final results are finite.

We derive the complete set of evolutions of chiralityodd twist3
fragmentation functions at oneloop level. There are totally nine real twist3
fragmentation functions, among which seven are independent. The
renormalizationscale dependence of the nine functions has an important
implication for studies of single transversespin asymmetries. We find that the
evolutions of the three complex fragmentation functions defined by
quarkgluonquark operator are mixed with themselves. There is no mixing with
the fragmentation functions defined only with bilinear quark field operators.
In the large$N_c$ limit the evolutions of the three complex fragmentation
functions are simplified and reduced to six homogenous equations.

We study transversemomentumdependent factorization at twist3 for DrellYan
processes. The factorization can be derived straightforwardly at leading order
of $\alpha_s$. But at this order we find that lightcone singularities already
exist and effects of soft gluons are not correctly factorized. We regularize
the singularities with gauge links off the lightcone and introduce a soft
factor to factorize the effects of soft gluons. Interestingly, the soft factor
must be included in the definition of subtracted TMD parton distributions to
correctly factorize the effects of soft gluons. We derive the CollinsSoper
equation for one of twist3 TMD parton distributions. The equation can be
useful for resummation of large logarithms terms appearing in the corresponding
structure function in collinear factorization. However, the derived equation is
nonhomogeneous. This will make the resummation complicated.

We study single spin asymmetries at oneloop accuracy in DrellYan processes
in which one of the initial hadrons is transversely polarized. The
spindependent part of differential crosssections can be factorized with
various hadronic matrix elements of twist2 and twist3 operators. These
operators can be of even and oddchirality. In this work, the studied
observables of asymmetries are differential crosssections with different
weights. These weights are selected so that the observables are spindependent
and their virtual corrections are completely determined by quark form factor.
In the calculations of oneloop corrections we meet collinear divergences in
the contributions involving chiralityodd and chiralityeven operators. We find
that all of the divergences can be correctly subtracted. Therefore, our results
give an explicit example of QCD factorization with twist3 operators,
especially, QCD factorization with chiralityodd twist3 operators.

Inclusive production of a quarkonium $\eta_{c,b}$ in hadron collisions at low
transverse momentum can be used to extract various
TransverseMomentumDependent(TMD) gluon distributions of hadrons, provided the
TMD factorization for the process holds. The factorization involving
unpolarized TMD gluon distributions of unpolarized hadrons has been examined
with onshell gluons at oneloop level. In this work we study the factorization
at oneloop level with diagram approach in the most general case, where all TMD
gluon distributions at leading twist are involved. We find that the
factorization holds and the perturbative effects are represented by one
perturbative coefficient. Since the initial gluons from hadrons are offshell
in general, there exists the socalled superleading region found recently. We
find that the contributions from this region can come from individual diagrams
at oneloop level, but they are cancelled in the sum. Our factorized result for
the differential crosssection is explicitly gaugeinvariant.

We study semiinclusive DIS with a transversely polarized target in the
approach of collinear factorization. The effects related with the transverse
polarization are at twist3. We derive the complete result of twist3
contributions to the relevant hadronic tensor at leading order of $\alpha_s$,
and construct correspondingly experimental observables. Measuring these
observables will help to extract the twist2 transversity distribution, twist3
distributions and twist3 fragmentation functions of the produced unpolarized
hadron. A detailed comparison with the approach of
transversemomentumdependent factorization is made.

The collinear factorization at twist3 for DrellYan processes is studied
with the motivation to solve the discrepancy in literature about the single
spin asymmetry in the lepton angular distribution, and to show how QCD gauge
invariance is realized in the hadronic tensor. The obtained result here agrees
with our early result derived with a totally different approach. In addition to
the asymmetry we can construct another two observables to identify the spin
effect. We show that the gauge invariance of different contributions in the
hadronic tensor is made in different ways by summing the effects of gluon
exchanges. More interestingly is that we can show that the virtual correction
to one structure function of the hadronic tensor, hence to some weighted SSA
observables, is completely determined by the quark form factor. This will
simplify the calculation of higher order corrections. The corresponding result
in semiinclusive DIS is also given for the comparison with DrellYan
processes.

Quarkonium production at low transverse momentum in hadron collisions can be
used to extract TransverseMomentumDependent(TMD) gluon distribution
functions, if TMD factorization holds there. We show that TMD factorization for
the case of Pwave quarkonium with $J^{PC}=0^{++}, 2^{++}$ holds at oneloop
level, but is violated beyond oneloop level. TMD factorization for other
Pwave quarkonium is also violated already at oneloop.

We study the oneloop correction in TransverseMomentumDependent(TMD)
factorization for DrellYan processes at small transverse momentum of the
lepton pair. We adopt the socalled subtractive approach, in which one can
systematically construct contributions for subtracting longdistance effects
represented by diagrams. The perturbative parts are obtained after the
subtraction. We find that the perturbative coefficients of all structure
functions in TMD factorization at leading twist are the same. The perturbative
parts can also be studied with scattering of partons instead of hadrons. In
this way, the factorization of many structure functions can only be examined by
studying the scattering of multiparton states, where there are many diagrams.
These diagrams have no similarities to those treated in the subtractive
approach. As an example, we use existing results of one structure function
responsible for SingleSpinAsymmetry, to show that these diagrams in the
scattering of multiparton states are equivalent to those treated in the
subtractive approach after using Ward identity.

We study DIS with a spin1 hadron or a nucleus polarized in arbitrary
direction. The differential crosssection in this case will have an azimuthal
dependence. We derive the dependence by including contributions from twist2
and twist3 QCD operators. The twist3 contribution is computed at treelevel.
A spin1 hadron or nucleus can have nonzero tensor polarization. We find that
all structure functions related to the tensorpolarization of the initial
hadron can be extracted by study the azimuthal dependence. A nonzero result of
the structure functions of a nucleus from experiment will indicate nontrivial
inner structure of the nucleus.

We show in detail that the $k_T$factorization for exclusive processes is
gaugedependent and inconsistent.

We have proposed to use an effective theory to describe interactions of an
$N\bar N$system. The effective theory can be constructed in analogy to the
existing effective theory for an $NN$system. In this work we study the
nexttoleading order correction to $N\bar N$ scattering near the threshold in
the effective theory. We find that the experimental data can be well described
with the effective theory.

We use multiparton states to examine the leading order collinear
factorization of single transversespin asymmetries in DrellYan processes.
Twist3 operators are involved in the factorization. We find that the socalled
softgluonpole contribution in the factorization must exist in order to make
the factorization correct. This contribution comes from the corresponding
crosssection at oneloop, while the hardpole contribution in the
factorization comes from the crosssection at treelevel. Although the two
contributions come from results at different orders, their perturbative
coefficient functions in the factorization are at the same order. This is in
contrast to factorizations only involving twist2 operators. The
softgluonpole contribution found in this work is in agreement with that
derived in a different way. For the hardpole contributions we find an extra
contribution from an extra parton process contributing to the asymmetries. We
also solve a part of discrepancy in evolutions of the twist3 operator. The
method presented here for analyzing the factorization can be generalized to
other processes and can be easily used for studying factorizations at higher
orders, because the involved calculations are of standard scattering
amplitudes.

The observed enhancement of $p\bar p$production near the threshold in
radiative decays of $J/\psi$ and $e^+e^$annihilations can be explained with
final state interactions among the produced $N\bar N$ system, where the
enhancement is essentially determined by $N\bar N$ elastic scattering
amplitudes. We propose to use an effective theory for interactions in a $N\bar
N$ system near its threshold. The effective theory is similar to the wellknown
one for interactions in a $NN$ system but with distinctions. It is interesting
to note that in the effective theory some corrections to scattering amplitudes
at treelevel can systematically be summed into a simple form. These
corrections are from rescattering processes. With these corrected amplitudes we
are able to describe the enhancement near the threshold in radiative decays of
$J/\psi$ and $e^+e^$annihilations, and the $p\bar p$ elastic scattering near
the threshold.

Glauber gluons in DrellYan processes are soft gluons with the transverse
momenta much larger than their momentum components along the directions of
initial hadrons. Their existence has been a serious challenge in proving the
factorization of DrellYan processes. The recently proposed soft collinear
effect theory of QCD can provide a transparent way to show factorizations for a
class of processes, but it does not address the effect of glauber gluons. In
this letter we first confirm the existence of glauber gluons through an
example. We then add glauber gluons into the effective theory and study their
interaction with other particles. In the framework of the effective theory with
glauber gluons we are able to show that the effects of glauber gluons in
DrellYan processes are canceled and the factorization holds in the existence
of glauber gluons. Our work completes the proof or argument of factorization of
DrellYan process in the framework of the soft collinear effective theory.

In the collinear factorization of the form factor for the transition
$\gamma^* \pi^0 \to \pi$ the hard part contains double log terms. These terms
will spoil the perturbative expansion of the hard part. A simple exponentiation
for resummation leads to divergent results. We study the resummation of these
double log's. We make an analysis to show the origin of the double log's. With
the understanding of the origin one can introduce soft factors and Nonstandard
LightCone Wave Functions(NLCWF) to derive a factorized form for the form
factor, where the hard part does not contain the double log's. There is a
perturbative relation between NLCWF and the standard LightCone Wave
Function(LCWF). Beside the renormalization scale $\mu$ the introduced NLCWF's
and soft factors have extra scales to characterize the double log's. Using the
evolutions of the extra scales and the relation we can do the resummation of
the double log's perturbatively in sense that LCWF's are the only
nonpertubative objects in the resumed formula. Our results with some models of
LCWF show that there is a significant difference between numerical predictions
with or without the resummation.

Two factorization approaches have been proposed for single transverse spin
asymmetries. One is the collinear factorization, another is the
transversemomentumdependent factorization. They have been previously derived
in a formal way by using diagram expansion at hadron level. If the two
factorizations hold or can be proven, they should also hold when we replace
hadrons with parton states. We examine these two factorizations at parton level
with massless partons. It is nontrivial to generate these asymmetries at parton
level with massless partons because the asymmetries require helicityflip and
nonzero absorptive parts in scattering amplitudes. By constructing suitable
parton states with massless partons we derive the two factorizations for the
asymmetry in DrellYan processes. It is found from our results that the
collinear factorization derived at parton level is not the same as that derived
at hadron level. Our results with massless partons confirm those derived with
single massive parton state in our previous works.

Using partonic states we verify the transverse momentum dependent
factorization for single transverse spin asymmetries in semi inclusive DIS at
leading but nontrivial order of alpha_s. The factorization has been previously
derived in a formal way by using diagram expansion at hadron level. We find
with our partonic results that two relevant structure functions satisfy the
factorization. Our results also satisfy the collinear factorization but with
the perturbative coefficient different than that derived formally.

A new method is proposed to calculate wave functions in $k_T$factorization
in \cite{LiMi} as a comment about our paper \cite{FMW}. We point out that the
results obtained with the method are in conflict with the translation
invariance and depend on the chosen contours for loopintegrals. Therefore, the
method is in principle unacceptable and the results with the method cannot be
correct.

This paper has been withdrawn

In the $k_T$factorization for exclusive processes, the nontrivial
$k_T$dependence of perturbative coefficients, or hard parts, is obtained by
taking offshell partons. This brings up the question of whether the
$k_T$factorization is gauge invariant. We study the $k_T$factorization for
the case $\pi \gamma^* \to \gamma$ at oneloop in a general covariant gauge.
Our results show that the hard part contains a lightcone singularity that is
absent in the Feynman gauge, which indicates that the $k_T$factorization is
{\it not} gauge invariant. These divergent contributions come from the
$k_T$dependent wave function of $\pi$ and are not related to a special
process. Because of this fact the $k_T$factorization for any process is not
gauge invariant and is violated. Our study also indicates that the
$k_T$factorization used widely for exclusive Bdecays is not gauge invariant
and is violated.

The k_Tfactorization has been widely used for exclusive decays of Bmesons.
In this factorization the pertubative coefficients are extracted from
scattering of offshell partons. Because the offshellness of partons the
extracted pertubative coefficients in general are not gaugeinvariant. We show
that these perturbative coefficients contain gaugedependent singularities
beyond treelevel. This leads to that the k_Tfactorization is gaugedependent
and violated in the general covariant gauge. This is not consistent with
expectations of physics. In the factorization there exist two different wave
functions for a Bmeson in general. We show that the two wave functions satisfy
different evolution equations whose solutions are used to resum large log's.
Based on this fact, the rseummation of large log's at leading log approximation
in the widely used k_Tfactorization with only one wave function is
inconsistent.

Single transversespin asymmetries have been studied intensively both in
experiment and theory. Theoretically, two factorization approaches have been
proposed. One is by using transversemomentumdependent factorization and the
asymmetry comes from the so called Sivers function. Another is by using
collinear factorization where the nonperturbative effect is parameterized by a
twist3 hadronic matrix element. However, the factorized formulas for the
asymmetries in the two approaches are derived at hadron level formally by
diagram expansion, where one works with various parton density matrices of
hadrons. If the two factorizations hold, they should also hold at parton level.
We examine this for DrellYan processes by replacing hadrons with partons. By
calculating the asymmetry, Sivers function and the twist3 matrix element at
nontrivial leading order of $\alpha_s$, we find that we can reproduce the
result of the transversemomentumdependent factorization. But we can only
verify the result of the collinear factorization partly. Two formally derived
relations between Sivers function and the twist3 matrix element are also
examined with negative results.

We study the observed enhancement of a $p\bar p$ system near the threshold in
the process $J/\psi \to \gamma p\bar p$ and $e^+ e^ \to p\bar p$. From early
studies the enhancement can be explained by final state interactions, which are
in general taken into account with some potential models. In this work we offer
a simple approach within quantum field theory to explain the observed
enhancement. We point out that among different final state interactions the
rescattering in a $N\bar N$ system though exchange of $\pi$ is the most
important. The effects of the rescattering is completely fixed by the
wellknown coupling $g_{\pi NN}$. Our results show that the enhancement in
$J/\psi \to \gamma p\bar p$ and $e^+ e^ \to p\bar p$ can be well described
with the rescattering effects.