
Parton distribution functions (PDFs) are nonperturbative quantities
describing the relation between a hadron and quarks and gluons within it. We
propose to extract PDFs from QCD global analysis of "data" generated by lattice
QCD calculations of good "lattice cross sections", which are basically
singlehadron matrix elements that are lattice QCD calculable and perturbative
QCD factorizable into the PDFs. To demonstrate the existence of good "lattice
cross sections", we take quasiquark distribution introduced by Ji [1] as a
case study to show that it could be factorized into the PDFs to all orders in
perturbation theory if it can be multiplicatively renormalized. We calculate
the factorized coefficients at the nexttoleading order in $\alpha_s$.

In the framework of quantum chromodynamics (QCD), parton distribution
functions (PDFs) quantify how the momentum and spin of a hadron are divided
among its quark and gluon constituents. Two main approaches exist to determine
PDFs. The first approach, based on QCD factorization theorems, realizes a QCD
analysis of a suitable set of hardscattering measurements, often using a
variety of hadronic observables. The second approach, based on firstprinciple
operator definitions of PDFs, uses lattice QCD to compute directly some
PDFrelated quantities, such as their moments. Motivated by recent progress in
both approaches, in this document we present an overview of latticeQCD and
globalanalysis techniques used to determine unpolarized and polarized proton
PDFs and their moments. We provide benchmark numbers to validate present and
future latticeQCD calculations and we illustrate how they could be used to
reduce the PDF uncertainties in current unpolarized and polarized global
analyses. This document represents a first step towards establishing a common
language between the two communities, to foster dialogue and to further improve
our knowledge of PDFs.

Following our previous proposal [1], we construct a class of good "lattice
cross sections" (LCSs), from which we could study partonic structure of hadrons
from ab initio lattice QCD calculations. These good LCSs, on the one hand, can
be calculated directly in lattice QCD, and on the other hand, can be factorized
into parton distribution functions (PDFs) with calculable coefficients, in the
same way as QCD factorization for factorizable hadronic cross sections. PDFs
could be extracted from QCD global analysis of the lattice QCD generated data
of LCSs. We also show that proposed functions for lattice QCD calculation of
PDFs in the literature are special cases of these good LCSs.

Heavy quarkonium ($\Upsilon$) production at low transverse momentum
($P_\perp$) in highenergy hadronic collisions is revisited from the point of
view of transverse momentum dependent (TMD) framework. We perform resummation
of double logarithmic correction associated with initialstate soft gluon
shower for $b\bar b$ production by employing CollinsSoperSterman (CSS)
formalism. We show that the CSS formalism provides a nice description of
$\Upsilon$ production data in p+$\rm\bar{p}$ collisions at Tevatron and p+p
collisions at the LHC.

We study the production and polarization of $J/\psi$ mesons within a jet in
protonproton collisions at the LHC. We define the $J/\psi$jet fragmentation
function as a ratio of differential jet cross sections with and without the
reconstructed $J/\psi$ in the jet. We demonstrate that this is a very useful
observable to help explore the $J/\psi$ production mechanism, and to
differentiate between different NRQCD global fits based on inclusive $J/\psi$
cross sections. Furthermore, we propose to measure the polarization of $J/\psi$
mesons inside the jet, which can provide even more stringent constraints for
the heavy quarkonium production mechanism.

Quasiparton distribution functions have received a lot of attentions in both
perturbative QCD and lattice QCD communities in recent years because they not
only carry good information on the parton distribution functions, but also
could be evaluated by lattice QCD simulations. However, unlike the parton
distribution functions, the quasiparton distribution functions have
perturbative ultraviolet power divergences because they are not defined by
twist2 operators. In this paper, we identify all sources of ultraviolet
divergences for the quasiparton distribution functions in coordinatespace,
and demonstrate that power divergences, as well as all logarithmic divergences
can be renormalized multiplicatively to all orders in QCD perturbation theory.

In recent years, the quasi parton distribution has been introduced for
extracting the parton distribution functions from lattice QCD simulations. The
quasi and standard distribution share the same perturbative collinear
singularity and the renormalized quasi distribution can be factorized into the
standard distribution with a perturbative matching factor. The quasi parton
distribution is known to have powerlaw UV divergences, which do not exist in
the standard distribution. We discuss in this talk the nonperturbative
renormalization scheme for the power divergence. We also demonstrate the
perturbative matching of the quasi quark distribution between continuum and
lattice at the oneloop.

Polarization measurements provide a detailed method to test the Standard
Model and to search for new physics. Most previous studies depend on
preselected coordinates, which blurs the significance of the results. The
construction of two rotationinvariant observables in vector boson decay into a
fermion pair has been proved to be a big success. In this work, we show that
there are more rotationinvariant observables and provide a general recipe to
find all of them in an arbitrary decay process. Taking spin1/2 and spin1
particle decay processes as examples, we calculate the explicit expressions of
all rotationinvariant observables, which can serve as a robust test of the
detector acceptance and help the analysis of experimental data.

A completely new strategy to calculate parton distribution functions on the
lattice has recently been proposed. In this method, lattice calculable
observables, called quasi distributions, are related to normal distributions.
The quasi distributions are known to contain powerlaw UV divergences arise
from a Wilson line in the nonlocal operator, while the normal distributions
only have logatithmic UV divergences. We propose possible method to subtract
the power divegence to make the matching of the quasi with the normal
distributions welldefined. We also demonstrate the matching of the quasi quark
distribution between continuum and lattice implementing the power divergence
subtraction. The matching calculations are carried out by oneloop
perturbation.

We present the nuclear transverse momentum broadening for DrellYan lepton
pair production in p+A collisions at nexttoleading order (NLO) in powers of
strong coupling constant $\alpha_s$. We verify that the transverse momentum
weighted differential cross section in NLO perturbative QCD (pQCD) at twist4
can be factorized into the convolution of parton distribution function of an
active parton in the projectile proton, a twist4 multiparton correlation
function of the target nucleus, and perturbatively calculable hard coefficient
function. We identify a QCD evolution equation for such a twist4 nuclear
gluonquark correlation function, and verify its universality  its process
independence. This finding demonstrates the prediction power of pQCD
factorization approach for studying parton multiple scattering in nuclear
medium.

In this talk, we present the QCD factorization formula for heavy quarkonium
production at large $p_T$ with factorized leadingpower and nexttoleading
power contributions in the $1/p_T$ expansion. We show that the leading order
analytical calculations in this QCD factorization approach can reproduce
effectively the full nexttoleading order numerical results derived using
nonrelativistic QCD (NRQCD) factorization formalism. We demonstrate that the
nexttoleading power contributions are crucial to the description of the
channels that are the most relevant for the rate as well as polarization of
$J/\psi$ production at current collider energies.

Study of the polarized heavy quarkonium production in recently proposed QCD
factorization formalism requires knowledge of a large number of input
fragmentation functions (FFs) from a single parton or a heavy quarkantiquark
pair to a polarized heavy quarkonium. In this paper, we calculate these FFs at
the input scale in terms of nonrelativistic QCD (NRQCD) factorization. We
derive all relevant polarized NRQCD longdistance matrix elements based
symmetries and propose a selfconsistent scheme to define them in arbitrary $d$
dimensions. We then calculate polarized input FFs contributed from all $S$wave
and $P$wave NRQCD intermediate states. With our calculation of the polarized
input FFs, and the partonic hard part available in literatures, the QCD
factorization formalism is ready to be applied to polarized heavy quarkonium
production.

We introduce a modified factorization formalism in quantum chromodynamics for
hadronic production of $W$ and $Z$ bosons at large transverse momentum $p_T$.
When $p_T$ is much larger than the invariant mass $Q$ of the vector boson, this
new factorization formalism systematically resums the large fragmentation
logarithms, $\alpha_s^m\ln^m(p_T^2/Q^2)$, to all orders in the strong coupling
$\alpha_s$. Using our modified factorization formalism, we calculate the
nexttoleading order (NLO) predictions for $W$ and $Z$ boson production at
high $p_T$ at the CERN Large Hadron Collider and at a future 100 TeV
protonproton collider. Our NLO results are about $5\%$ larger in
normalization, and they show improved convergence and moderate reduction of the
scale variation compared to the NLO predictions derived in a conventional
fixedorder perturbative expansion.

We briefly summarize the important and critical roles that PHENIX Experiment
has played in determining the proton's internal structure in terms of quarks
and gluons, and their dynamics. Some pioneering measurements by PHENIX
Experiment on the motion and polarization of quarks and gluons, as well as
their correlations inside a fast moving proton are presented. Some future
opportunities and potentials of PHENIX Experiment are also discussed.

In this talk, we review a QCD factorization based approach to extract parton
distribution and correlation functions from lattice QCD calculation of single
hadron matrix elements of quarkgluon operators. We argue that although the
lattice QCD calculations are done in the Euclidean space, the nonperturbative
collinear behavior of the matrix elements are the same as that in the Minkowski
space, and could be systematically factorized into parton distribution
functions with infrared safe matching coefficients. The matching coefficients
can be calculated perturbatively by applying the factorization formalism on to
asymptotic partonic states.

In this talk, I demonstrate that the proposed ElectronIon Collider (EIC)
will be an ideal and unique future facility to address many overarching
questions about QCD and strong interaction physics at one place. The EIC will
be the world's first polarized electronproton (and light ion), as well as the
first electronnucleus collider at flexible collision energies. With its high
luminosity and beam polarization, the EIC distinguishes itself from HERA and
the other fixed target electronhadron facilities around the world. The EIC is
capable of taking us to the next QCD frontier to explore the glue that binds us
all.

We calculate the ${\cal O}(\alpha_s^3)$ shortdistance, QCD
collinearfactorized coefficient functions for all partonic channels that
include the production of a heavy quark pair at short distances. This provides
the first power correction to the collinearfactorized inclusive hadronic
production of heavy quarkonia at large transverse momentum, $p_T$, including
the full leadingorder perturbative contributions to the production of heavy
quark pairs in all color and spin states employed in NRQCD treatments of this
process. We discuss the role of the first power correction in the production
rates and the polarizations of heavy quarkonia in high energy hadronic
collisions. The consistency of QCD collinear factorization and nonrelativistic
QCD factorization applied to heavy quarkonium production is also discussed.

We present a factorization formalism for inclusive production of heavy
quarkonia of large transverse momentum, $p_T$ at collider energies, including
both leading power (LP) and nexttoleading power (NLP) behavior in $p_T$. We
demonstrate that both LP and NLP contributions can be factorized in terms of
perturbatively calculable shortdistance partonic coefficient functions and
universal nonperturbative fragmentation functions, and derive the evolution
equations that are implied by the factorization. We identify projection
operators for all channels of the factorized LP and NLP infrared safe
shortdistance partonic hard parts, and corresponding operator definitions of
fragmentation functions. For the NLP, we focus on the contributions involving
the production of a heavy quark pair, a necessary condition for producing a
heavy quarkonium. We evaluate the first nontrivial order of evolution kernels
for all relevant fragmentation functions, and discuss the role of NLP
contributions.

A polarized $ep/eA$ collider (ElectronIon Collider, or EIC), with polarized
proton and lightion beams and unpolarized heavyion beams with a variable
centerofmass energy $\sqrt{s} \sim 20$ to $\sim100$~GeV (upgradable to
$\sim 150$ GeV) and a luminosity up to $\sim 10^{34} \, \textrm{cm}^{2}
\textrm{s}^{1}$, would be uniquely suited to address several outstanding
questions of Quantum Chromodynamics, and thereby lead to new qualitative and
quantitative information on the microscopic structure of hadrons and nuclei.
During this meeting at Jefferson Lab we addressed recent theoretical and
experimental developments in the spin and the threedimensional structure of
the nucleon (sea quark and gluon spatial distributions, orbital motion,
polarization, and their correlations). This minireview contains a short
update on progress in these areas since the EIC White
paper~\cite{Accardi:2012qut}.

We show that when the factorized cross section for heavy quarkonium
production includes nexttoleading power (NLP) contributions associated with
the production of the heavy quark pair at short distances, it naturally
reproduces all high $p_T$ results calculated in nonrelativistic QCD (NRQCD)
factorization. This extended formalism requires fragmentation functions for
heavy quark pairs, as well as for light partons. When these fragmentation
functions are themselves calculated using NRQCD, we find that two of the four
leading NRQCD production channels, ${^3\hspace{0.6mm}S_{1}^{[1]}}$ and
${^1\hspace{0.6mm}S_{0}^{[8]}}$, are dominated by the NLP contributions for a
very wide $p_T$ range. The large nexttoleading order corrections of NRQCD are
absorbed into the leading order of the first power correction. The impact of
this finding on the heavy quarkonium production and its polarization is
discussed.

Difference between the $\bar u$ and $\bar d$ sea quark distributions in the
proton was first observed in the violation of the Gottfried sum rule in
deepinelastic scattering (DIS) experiments. The parton momentum fraction $x$
dependence of this difference has been measured over the region $0.02 < x <
0.35$ from DrellYan and semiinclusive DIS experiments. The DrellYan data
suggested a possible signchange for $\bar d(x)\bar u(x)$ near $x \sim 0.3$,
which has not yet been explained by existing theoretical models. We present an
independent evidence for the $\bar d(x)\bar u(x)$ signchange at $x \sim 0.3$
from an analysis of the DIS data. We further discuss the $x$dependence of
$\bar d(x)\bar u(x)$ in the context of meson cloud model and the lattice QCD
formulation.

We study single jet production in electronnucleus collisions e^ + N_A > J
+ X, using the 1jettiness (\tau_1) global event shape. It inclusively
quantifies the pattern of radiation in the final state, gives enhanced
sensitivity to soft radiation at wide angles from the nuclear beam and
finalstate jet, and facilitates the resummation of large Sudakov logarithms
associated with the veto on additional jets. Through their effect on the
observed pattern of radiation, 1jettiness can be a useful probe of nuclear
PDFs and power corrections from dynamical effects in the nuclear medium. This
formalism allows for the standard jet shape analysis while simultaneously
providing sensitivity to soft radiation at wide angles from the jet. We use a
factorization framework for crosssections differential in $\tau_1$ and the
transverse momentum (P_{J_T}) and rapidity (y) of the jet, in the region
\tau_1<< P_{J_T}. The restriction $\tau_1 << P_{J_T}$ allows only soft
radiation between the nuclear beam and jet directions, thereby acting as a veto
on additional jets. This region is also insensitive to the details of the jet
algorithm, allowing for better theoretical control over resummation, while
providing enhanced sensitivity to nuclear medium effects. We give numerical
results at leading twist, with resummation at the nexttonexttoleading
logarithmic (NNLL) level of accuracy, for a variety of nuclear targets. Such
studies would be ideal for the EIC and the LHeC proposals for a future
electronion collider, where a range of nuclear targets are planned.

We highlight the progress, current status, and open challenges of QCDdriven
physics, in theory and in experiment. We discuss how the strong interaction is
intimately connected to a broad sweep of physical problems, in settings ranging
from astrophysics and cosmology to stronglycoupled, complex systems in
particle and condensedmatter physics, as well as to searches for physics
beyond the Standard Model. We also discuss how success in describing the strong
interaction impacts other fields, and, in turn, how such subjects can impact
studies of the strong interaction. In the course of the work we offer a
perspective on the many research streams which flow into and out of QCD, as
well as a vision for future developments.

Recently, a new perturbative QCD factorization formalism for heavy quarkonium
production at a large transverse momentum was proposed. Phenomenological
application of this new approach relies on our knowledge of a large number of
universal fragmentation functions (FFs) at an input factorization scale
$\mu_0\gtrsim 2m_Q$ with heavy quark mass $m_Q$, which are nonperturbative, and
in principle, should be extracted from data. With heavy quark mass $m_Q\gg
\Lambda_{\rm QCD}$, we calculate these input FFs in terms of nonrelativistic
QCD (NRQCD) factorization. In a companion paper, we calculated these input FFs
in NRQCD factorization approach including all contributions from $S$wave NRQCD
$Q\bar{Q}$states. In this paper, we calculate contributions to the heavy
quarkpair FFs from all $P$wave NRQCD $Q\bar{Q}$states.

A QCD factorization formalism was recently proposed for evaluating heavy
quarkonium production at large $p_T$ at collider energies. With systematically
calculated shortdistance partonic hard parts and evolution kernels of
fragmentation functions (FFs), the predictive power of this factorization
approach relies on our knowledge of a large number of universal FFs at an input
factorization scale $\mu_0\gtrsim 2m_Q$ with heavy quark mass $m_Q$. With the
large heavy quark mass, the relative motion of the heavy quark pair inside a
heavy quarkonium is effectively nonrelativistic. We evaluate these universal
input FFs using nonrelativistic QCD (NRQCD) factorization, and express the
large number of FFs in terms of a few universal NRQCD longdistance matrix
elements (LDMEs) with perturbatively calculated coefficients. We derive
complete contributions to the singleparton FFs at both ${\cal O}(\alpha_s)$
and ${\cal O}(\alpha_s^2)$, and the heavy quark pair FFs at ${\cal
O}(\alpha_s)$. We present detailed derivation for all contributions involving
LDMEs of $S$wave NRQCD $Q\bar{Q}$states ($P$wave contributions in a
companion paper). Our results bridge the gap between the QCD factorization
formalism and its phenomenological applications.