This document presents the initial scientific case for upgrading the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) to 22 GeV. It is the result of a community effort, incorporating insights from a series of workshops conducted between March 2022 and April 2023. With a track record of over 25 years in delivering the world's most intense and precise multi-GeV electron beams, CEBAF's potential for a higher energy upgrade presents a unique opportunity for an innovative nuclear physics program, which seamlessly integrates a rich historical background with a promising future. The proposed physics program encompass a diverse range of investigations centered around the nonperturbative dynamics inherent in hadron structure and the exploration of strongly interacting systems. It builds upon the exceptional capabilities of CEBAF in high-luminosity operations, the availability of existing or planned Hall equipment, and recent advancements in accelerator technology. The proposed program cover various scientific topics, including Hadron Spectroscopy, Partonic Structure and Spin, Hadronization and Transverse Momentum, Spatial Structure, Mechanical Properties, Form Factors and Emergent Hadron Mass, Hadron-Quark Transition, and Nuclear Dynamics at Extreme Conditions, as well as QCD Confinement and Fundamental Symmetries. Each topic highlights the key measurements achievable at a 22 GeV CEBAF accelerator. Furthermore, this document outlines the significant physics outcomes and unique aspects of these programs that distinguish them from other existing or planned facilities. In summary, this document provides an exciting rationale for the energy upgrade of CEBAF to 22 GeV, outlining the transformative scientific potential that lies within reach, and the remarkable opportunities it offers for advancing our understanding of hadron physics and related fundamental phenomena.

We perform the first simultaneous extraction of parton collinear and transverse degrees of freedom from low-energy fixed-target Drell-Yan data in order to compare the transverse momentum dependent (TMD) parton distribution functions (PDFs) of the pion and proton. We demonstrate that the transverse separation of the quark field encoded in TMDs of the pion is more than $4 \sigma$ smaller than that of the proton. Additionally, we find the transverse separation of the quark field decreases as its longitudinal momentum fraction decreases. In studying the nuclear modification of TMDs, we find clear evidence for a transverse EMC effect. We comment on possible explanations for these intriguing behaviors, which call for a deeper examination of tomography in a variety of strongly interacting quark-gluon systems.

We extract pion parton distribution functions (PDFs) in a Monte Carlo global QCD analysis of experimental data together with reduced Ioffe time pseudo-distributions and matrix elements of current-current correlators generated from lattice QCD. By including both experimental and lattice QCD data, our analysis rigorously quantifies both the uncertainties of the pion PDFs and systematic effects intrinsic to the lattice QCD observables. The reduced Ioffe time pseudo-distributions significantly decrease the uncertainties on the PDFs, while the current-current correlators are limited by the systematic effects associated with the lattice. Consistent with recent phenomenological determinations, the behavior of the valence quark distribution of the pion at large momentum fraction is found to be similar to(1 - x)(beta eff) with beta(eff) approximate to 1.0-1.2.

This report describes the physics case, the resulting detector requirements, and the evolving detector concepts for the experimental program at the Electron-Ion Collider (EIC). The EIC will be a powerful new high-luminosity facility in the United States with the capability to collide high-energy electron beams with high-energy proton and ion beams, providing access to those regions in the nucleon and nuclei where their structure is dominated by gluons. Moreover, polarized beams in the EIC will give unprecedented access to the spatial and spin structure of the proton, neutron, and light ions. The studies leading to this document were commissioned and organized by the EIC User Group with the objective of advancing the state and detail of the physics program and developing detector concepts that meet the emerging requirements in preparation for the realization of the EIC. The effort aims to provide the basis for further development of concepts for experimental equipment best suited for the science needs, including the importance of two complementary detectors and interaction regions. This report consists of three volumes. Volume I is an executive summary of our findings and developed concepts. In Volume II we describe studies of a wide range of physics measurements and the emerging requirements on detector acceptance and performance. Volume III discusses general-purpose detector concepts and the underlying technologies to meet the physics requirements. These considerations will form the basis for a world-class experimental program that aims to increase our understanding of the fundamental structure of all visible matter

A valuable approach to the analysis of hadron physics observables is provided by QCD's equations-of-motion; namely, the Dyson-Schwinger equations. Drawing from a diverse collection of predictions, we revisit: $\gamma \gamma^* \to$ neutral pseudoscalar transition form factors, their corresponding valence-quark distribution amplitudes and a recent result on the pion distribution functions.

We perform the first global QCD analysis of pion valence, sea quark, and gluon distributions within a Bayesian Monte Carlo framework with threshold resummation on Drell-Yan cross sections at next-to-leading log accuracy. Exploring various treatments of resummation, we find that the large-x asymptotics of the valence quark distribution ∼(1-x)^{β_{v}} can differ significantly, with β_{v} ranging from ≈1 to >2.5 at the input scale. Regardless of the specific implementation, however, the resummation induced redistribution of the momentum between valence quarks and gluons boosts the total momentum carried by gluons to ≈40%, increasing the gluon contribution to the pion mass to ≈40  MeV.

We perform a new Monte Carlo QCD analysis of pion parton distribution functions, including, for the first time, transverse momentum dependent pion-nucleus Drell-Yan cross sections together with p(T)-integrated Drell-Yan and leading neutron electroproduction data from HERA. We assess the sensitivity of the Monte Carlo fits to kinematic cuts, factorization scale, and parametrization choice, and we discuss the impact of the various datasets on the pion's quark and gluon distributions. This study provides the necessary step towards the simultaneous analysis of collinear and transverse momentum dependent pion distributions and the determination of the pion's three-dimensional structure.

We perform the first global QCD analysis of parton distribution functions (PDFs) in the pion, combining πA Drell-Yan data with leading neutron electroproduction from HERA within a Monte Carlo approach based on nested sampling. Inclusion of the HERA data allows the pion PDFs to be determined down to much lower values of x, with relatively weak model dependence from uncertainties in the chiral splitting function. The combined analysis reveals that gluons carry a significantly higher pion momentum fraction, ∼30%, than that inferred from Drell-Yan data alone, with sea quarks carrying a somewhat smaller fraction, ∼15%, at the input scale. Within the same effective theory framework, the chiral splitting function and pion PDFs can be used to describe the d[over ¯]-u[over ¯] asymmetry in the proton.