Search

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., Comment: Updates to the list of authors; Preprint number changed from theory to experiment; Updates to sections 4 and 6, including additional figures

We discuss the prospects for improving the precision on the hadronic corrections to the anomalous magnetic moment of the muon, and the plans of the Muon $g-2$ Theory Initiative to update the Standard Model prediction., Comment: Contribution to the US Community Study on the Future of Particle Physics (Snowmass 2021)

We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant $\alpha$ and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including $\mathcal{O}(\alpha^5)$ with negligible numerical uncertainty. The electroweak contribution is suppressed by $(m_\mu/M_W)^2$ and only shows up at the level of the seventh significant digit. It has been evaluated up to two loops and is known to better than one percent. Hadronic contributions are the most difficult to calculate and are responsible for almost all of the theoretical uncertainty. The leading hadronic contribution appears at $\mathcal{O}(\alpha^2)$ and is due to hadronic vacuum polarization, whereas at $\mathcal{O}(\alpha^3)$ the hadronic light-by-light scattering contribution appears. Given the low characteristic scale of this observable, these contributions have to be calculated with nonperturbative methods, in particular, dispersion relations and the lattice approach to QCD. The largest part of this review is dedicated to a detailed account of recent efforts to improve the calculation of these two contributions with either a data-driven, dispersive approach, or a first-principle, lattice-QCD approach. The final result reads $a_\mu^\text{SM}=116\,591\,810(43)\times 10^{-11}$ and is smaller than the Brookhaven measurement by 3.7$\sigma$. The experimental uncertainty will soon be reduced by up to a factor four by the new experiment currently running at Fermilab, and also by the future J-PARC experiment. This and the prospects to further reduce the theoretical uncertainty in the near future-which are also discussed here-make this quantity one of the most promising places to look for evidence of new physics., Comment: 196 pages, 103 figures, version published in Phys. Rept., bib files for the citation references are available from: https://muon-gm2-theory.illinois.edu

In light of recent experimental results, we revisit the dispersive analysis of the $\omega \to 3\pi$ decay amplitude and of the $\omega\pi^0$ transition form factor. Within the framework of the Khuri-Treiman equations, we show that the $\omega \to 3\pi$ Dalitz-plot parameters obtained with a once-subtracted amplitude are in agreement with the latest experimental determination by BESIII. Furthermore, we show that at low energies the $\omega\pi^0$ transition form factor obtained from our determination of the $\omega \to 3\pi$ amplitude is consistent with the data from MAMI and NA60 experiments.

Khuri-Treiman equations have proven to be a useful theoretical tool in the analysis of 3-body decays, specially into the $3\pi$ final state. In this work we present in full detail the necessary generalization of the formalism to study the decays of particles with arbitrary spin, parity, and charge conjugation. To this extent, we find it most convenient to work with helicity amplitudes instead of the so-called invariant amplitudes, specially when dealing with the unitarity relations. The isobar expansions in the three possible ($s$-, $t$-, and $u$-) final channels are related with the appropriate crossing matrices. We pay special attention to the kinematical singularities and constraints of the helicity amplitudes, showing that these can be derived by means of the crossing matrix.

We use Regge phenomenology to study the structure of the poles of the $N^*$ and $\Delta^*$ spectrum. We employ the available pole extractions from partial wave analysis of meson scattering and photoproduction data. We assess the importance of the imaginary part of the poles (widths) to obtain a consistent determination of the parameters of the Regge trajectory. We compare the several pole extractions and show how Regge phenomenology can be used to gain insight in the internal structure of baryons. We find that the majority of the states in the parent Regge trajectories are compatible with a mostly compact three-quark state picture.

The Dalitz decays eta -> e^+e^-g and omega -> pi^0 e^+e^- have been measured in the g p -> eta p and g p -> omega p reactions, respectively, with the A2 tagged-photon facility at the Mainz Microtron, MAMI. The value obtained for the slope parameter of the electromagnetic transition form factor of eta, Lambda^{-2}_eta=(1.97+/-0.11_tot) GeV^{-2}, is in good agreement with previous measurements of the eta -> e^+e^-g and eta -> mu^+mu^-g decays. The uncertainty obtained in the value of Lambda^{-2}_eta is lower than in previous results based on the eta -> e^+e^-g decay. The value obtained for the omega slope parameter, Lambda^{-2}_omega_pi^0 = (1.99+/-0.21_tot) GeV^{-2}, is somewhat lower than previous measurements based on omega -> pi^0 mu^+mu^-, but the results for the omega transition form factor are in better agreement with theoretical calculations, compared to earlier experiments., Comment: 19 figures, 2 tables

In view of the recent high-statistic KLOE data for the $\eta \to \pi^+ \pi^- \pi^0$ decay, a new determination of the quark mass double ratio has been done. Our approach relies on a dispersive model that takes into account rescattering effects between three pions via subenergy unitarity. The latter is essential to reproduce the Dalitz plot distribution. A simultaneous description of the KLOE and WASA-at-COSY data is achieved in terms of just two real parameters. From a global fit, we determine $Q=21.6 \pm 1.1$. The predicted slope parameter for the neutral channel $\alpha=-0.025\pm 0.004$ is in reasonable agreement with the PDG average value., Comment: 7 pages, 2 figures

We present a unitary multichannel model for $\bar{K}N$ scattering in the resonance region that fulfills unitarity. It has the correct analytical properties for the amplitudes once they are extended to the complex-$s$ plane and the partial waves have the right threshold behavior. To determine the parameters of the model, we have fitted single-energy partial waves up to $J=7/2$ and up to 2.15 GeV of energy in the center-of-mass reference frame obtaining the poles of the $\Lambda^*$ and $\Sigma^*$ resonances, which are compared to previous analyses. We provide the most comprehensive picture of the $S=-1$ hyperon spectrum to date. Important differences are found between the available analyses making the gathering of further experimental information on $\bar{K}N$ scattering mandatory to make progress in the assessment of the hyperon spectrum., Comment: Joint Physics Analysis Center, 28 pages, 19 figures, 3 tables

We compare the low-energy partial wave analyses $\pi N$ scattering with a high-energy data via finite energy sum rules. We construct a new set of amplitudes by matching the imaginary part from the low-energy analysis with the high-energy, Regge parametrization and reconstruct the real parts using dispersion relations., Comment: Joint Physics Analysis Center ; 15 pages, 12 figures

We apply the generalized Veneziano model ($B_5$ model) in the double-Regge exchange limit to the $\gamma p \to K^+ K^- p$ reaction. Four different cases defined by the possible combinations of the signature factors of leading Regge exchanges (($K^*,a_2/f_2$), ($K^*,\rho/\omega$), ($K_2^*,a_2/f_2$), and ($K_2^*,\rho/\omega$)) have been simulated through the Monte Carlo method. Suitable event candidates for the double-Regge exchange high-energy limit were selected employing \textit{Van Hove} plots as a better alternative to kinematical cuts in the $K^+ K^- p$ Dalitz plot. In this way we predict and analyze the double-Regge contribution to the $K^+ K^- p$ Dalitz plot, which constitutes one of the major backgrounds in the search for strangeonia, hybrids and exotics using $\gamma p\rightarrow K^+K^-p$ reaction. We expect that data currently under analysis, and that to come in the future, will allow verification of the double-Regge behavior and a better assessment of this component of the amplitude., Comment: 12 pages, 10 figures

In this work, we present different representations of Khuri-Treiman equation, the advantage and disadvantage of each representations are discussed. Using a toy model for scattering amplitude, we study the sensitivity of solution of KT equation to the left hand cut of this toy model and to the different approximate methods. At last, we briefly discuss Watson's theorem when three particles in final states are involved., Comment: EPJA publication version

The decays $\omega/\phi \rightarrow 3\pi$ are considered in the dispersive framework that is based on the isobar decomposition and sub-energy unitarity. The inelastic contributions are parametrized by the power series in a suitably chosen conformal variable that properly accounts for the analytic properties of the amplitude. The Dalitz plot distributions and integrated decay widths are presented. Our results indicate that the final state interactions may be sizable. As a further application of the formalism we also compute the electromagnetic transition form factors of $\omega/\phi \rightarrow \pi^0\gamma^*$., Comment: 14 pages, 9 figures

A new measurement of the rare, doubly radiative decay eta->pi^0 gamma gamma was conducted with the Crystal Ball and TAPS multiphoton spectrometers together with the photon tagging facility at the Mainz Microtron MAMI. New data on the dependence of the partial decay width, Gamma(eta->pi^0 gamma gamma), on the two-photon invariant mass squared, m^2(gamma gamma), as well as a new, more precise value for the decay width, Gamma(eta->pi^0 gamma gamma) = (0.33+/-0.03_tot) eV, are based on analysis of 1.2 x 10^3 eta->pi^0 gamma gamma decays from a total of 6 x 10^7 eta mesons produced in the gamma p -> eta p reaction. The present results for dGamma(eta->pi^0 gamma gamma)/dm^2(gamma gamma) are in good agreement with previous measurements and recent theoretical calculations for this dependence., Comment: 10 pages, 8 figures, 1 table

We study the reactions $\gamma\gamma\rightarrow \pi^0\pi^0$, $\pi^+\pi^-$, $K^0\bar{K}^0$, $K^+K^-$, $\eta \eta$ and $\pi^0\eta$ based on a chiral Lagrangian with dynamical light vector mesons as formulated within the hadrogenesis conjecture. At present our chiral Lagrangian contains 5 unknown parameters that are relevant for the photon fusion reactions. They parameterize the strength of interaction terms involving two vector meson fields. These parameters are fitted to photon fusion data $\gamma\gamma\rightarrow \pi^0\pi^0$, $\pi^+\pi^-, \pi^0\eta$ and to the decay $\eta\rightarrow\pi^0\gamma\gamma$. In order to derive gauge invariant reaction amplitudes in the resonance region constraints from micro-causality and exact coupled-channel unitarity are used. Our results are in good agreement with the existing experimental data from threshold up to about 0.9 GeV for the two-pion final states. The $a_0$ meson in the $\pi^0\eta$ channel is dynamically generated and an accurate reproduction of the $\gamma\gamma\rightarrow \pi^0\eta$ data is achieved up to 1.2 GeV. Based on our parameter sets we predict the $\gamma\gamma\rightarrow $ $K^0\bar{K}^0$, $K^+K^-$, $\eta \eta$ cross sections., Comment: 17 pages, 10 figures, all assumptions used in our method are clarified

A theoretical study of Goldstone boson scattering based on the chiral Lagrangian with vector meson fields is presented. In application of a recently developed novel approach we extrapolate subthreshold partial-wave amplitudes into the physical region. The constraints set by micro-causality and coupled-channel unitarity are kept rigourously. It is shown that already the leading order subthreshold amplitudes lead to s- and p-wave $\pi\pi$ and $\pi K$ phase shifts are in agreement with the experimental data up to about 1.2 GeV., Comment: 4 pages, 4 figures, contribution to MESON 2012 workshop, Krakow, Poland, 31 May - 5 June 2012

We present a study of Goldstone boson scattering based on the flavor SU(3) chiral Lagrangian formulated with vector mesons in the tensor field representation. A coupled-channel channel computation is confronted with the empirical s- and p-wave phase shifts, where good agreement with the data set is obtained up to about 1.2 GeV. There are two relevant free parameters only, the chiral limit value of the pion decay constant and the coupling constant characterizing the decay of the rho meson into a pair of pions. We apply a recently suggested approach that implements constraints from micro- causality and coupled-channel unitarity. Generalized potentials are obtained from the chiral Lagrangian and are expanded in terms of suitably constructed conformal variables. The partial-wave scattering amplitudes are defined as solutions of non-linear integral equations that are solved by means of an N/D ansatz., Comment: 15 pages, 8 figures, typos corrected, accepted for publication in Physics Letters B

Dynamics of hadro-quarkonium system is formulated, based on the channel coupling of a light hadron (h) and heavy quarkonium (Q\bar{Q}) to intermediate open-flavor heavy-light mesons (Q\bar{q}, \bar{Q}q). The resulting effective interaction is defined by overlap integrals of meson wavefunctions and (hq\bar q) coupling, where h is pi, rho, omega, phi, without fitting parameters. Equations for hadro-quarkonium amplitudes and resonance positions are written explicitly, and numerically calculated for the special case of pi Upsilon(nS) (n=1,2,3). It is also shown, that the recently observed by Belle two peaks Z_b(10610) and Z_b(10650) are in agreement with the proposed theory. It is demonstrated, that theory predicts peaks at the BB*, B*B* thresholds in all available pi Upsilon(nS) channels. Analytic nature of these peaks is investigated, and shown to be due to a common multichannel resonance poles close to the BB*, B*B* thresholds. The general mechanism of these hadro-quarkonium resonances does not assume any molecular or four-quark (tetraquark) dynamics., Comment: 18 pages, 10 figures, minor changes, version to be published in Phys. Rev. D

We present a pedagogical case study how to combine micro-causality and unitarity based on a perturbative approach. The method we advocate constructs an analytic extrapolation of partial-wave scattering amplitudes that is constrained by the unitarity condition. Suitably constructed conformal mappings help to arrive at a systematic approximation of the scattering amplitude. The technique is illustrated at hand of a Yukawa interaction. The typical case of a superposition of strong short-range and weak long-range forces is investigated., Comment: 12 pages, 12 figures

Hyperfine splittings (HFS) are calculated within the Field Correlator Method, taking into account relativistic corrections. The HFS in bottomonium and the $B_q$ (q=n,s) mesons are shown to be in full agreement with experiment if a universal coupling $\alpha_{HF}=0.310$ is taken in perturbative spin-spin potential. It gives $M(B^*)-M(B)=45.7(3)$ MeV, $M(B_s^*)-M(B_s)=46.7(3)$ MeV ($n_f=4$), while in bottomonium $\Delta_{HF}(b\bar b)=M(\Upsilon(9460))-M(\eta_b(1S))=63.4$ MeV for $n_f=4$ and 71.1 MeV for $n_f=5$ are obtained; just latter agrees with recent BaBar data. For unobserved excited states we predict $M(\Upsilon(2S))-M(\eta_b(2S))=36(2)$ MeV, $M(\Upsilon(3S))-M(\eta(3S))=28(2)$ MeV, and also $M(B_c^*)=6334(4)$ MeV, $M(B_c(2S))=6868(4)$ MeV, $M(B_c^*(2S))=6905(4)$ MeV. The mass splittings between $D(2^3S_1)-D(2^1S_0)$, $D_s(2^3S_1)-D_s(2^1S_0)$ are predicted to be $\sim 70$ MeV, which are significantly smaller than in several other studies., Comment: 13 pages

The dynamical mechanism of channel coupling with the decay channels is applied to the case of coupled charmonium - $DD^*$ states with $J^{PC}=1^{++}$. A pole analysis is done and the $DD^*$ production cross section is calculated in qualitative agreement with experiment. The sharp peak at the $D_0D^*_0$ threshold and flat background are shown to be due to Breit-Wigner resonance, shifted by channel coupling from the original position of 3954 MeV for the $2^3P_1$, $Q\bar Q$ state. A similar analysis, applied to the $n=2$, $^3P_2$, $^1P_1$, $^3P_0$, allows us to associate the first one with the observed $Z(3930)$ J=2 and explains the destiny of $^3P_0$., Comment: 5 pages, 4 figures. Accepted for publication in Phys. Rev. Lett

A universal description of the hyperfine splittings (HFS) in bottomonium and the $B_q (q=n,s,c)$ mesons is obtained with a universal strong coupling constant $\alpha_s(\mu)=0.305(2)$ in a spin-spin potential. Other characteristics are calculated within the Field Correlator Method, taking the freezing value of the strong coupling independent of $n_f$. The HFS $M(B^*)- M(B)=45.3(3)$ MeV, $M(B_s^*) - M(B_s)=46.5(3)$ MeV are obtained in full agreement with experiment both for $n_f=3$ and $n_f=4$. In bottomonium, $M(\Upsilon(9460))- M(\eta_b)=70.0(4)$ MeV for $n_f=5$ agrees with the BaBar data, while a smaller HFS, equal to 64(1) MeV, is obtained for $n_f=4$. We predict HFS $M(\Upsilon(2S))-M(\eta_b(2S))=36(1)$ MeV, $M(\Upsilon(3S))- M(\eta(3S))=27(1)$ MeV, and $M(B_c^*) - M(B_c)= 57.5(10)$ MeV, which gives $M(B_c^*)=6334(1)$ MeV, $M(B_c(2 {}^1S_0))=6865(5)$ MeV, and $M(B_c^*(2S {}^3S_1))=6901(5)$ MeV., Comment: 5 pages revtex4

The mechanism of channel coupling via decay products is used to study energy shifts, level mixing as well as the possibility of new near-threshold resonances in $c\bar c, b\bar b$ systems. The Weinberg eigenvalue method is formulated in the multichannel problems, which allows to describe coupled-channel resonances and wave functions in a unitary way, and to predict new states due to channel coupling. Realistic wave functions for all single-channel states and decay matrix elements computed earlier are exploited, and no new fitting parameters are involved. Examples of level shifts, widths and mixings are presented; the dynamical origin of X(3872) and the destiny of the single-channel $2^3P_1(c\bar c)$ state are clarified. As a result a sharp and narrow peak in the state with quantum numbers $J^{PC}=1^{++}$ is found at 3.872 GeV, while the single-channel resonance originally around 3.940 GeV, becomes increasingly broad and disappears with growing coupling to open channels., Comment: 27 pages, 13 figures The article is accepted for publication in Physical Review D

The dielectron widths of $\Upsilon(nS) (n=1,...,7)$ and vector decay constants are calculated using the Relativistic String Hamiltonian with a universal interaction. For $\Upsilon(nS) (n=1,2,3)$ the dielectron widths and their ratios are obtained in full agreement with the latest CLEO data. For $\Upsilon(10580)$ and $\Upsilon(11020)$ a good agreement with experiment is reached only if the 4S--3D mixing (with a mixing angle $\theta=27^\circ\pm 4^\circ$) and 6S--5D mixing (with $\theta=40^\circ\pm 5^\circ$) are taken into account. The possibility to observe higher "mixed $D$-wave" resonances, $\tilde\Upsilon(n {}^3D_1)$ with $n=3,4,5$ is discussed. In particular, $\tilde\Upsilon(\approx 11120)$, originating from the pure $5 {}^3D_1$ state, can acquire a rather large dielectron width, $\sim 130$ eV, so that this resonance may become manifest in the $e^+e^-$ experiments. On the contrary, the widths of pure $D$-wave states are very small, $\Gamma_{ee}(n{}^3 D_1) \leq 2$ eV., Comment: 13 pages, no figures

The possible description of Z(4430) as a pseudoresonance structure in $\pi \psi'$ reaction, is considered. The analysis is performed with single-scattering contribution to $\pi \psi'$ elastic scattering via $D^*D_1(2420)$ intermediate energy., Comment: 3 pages, 4 figures

The possibility to observe new bottomonium states with $J^{PC}= 1^{--}$ in the region $10.7-11.1$ GeV is discussed. The analysis of the di-electron widths shows that the $(n+1)^3S_1$ and $n^3D_1$ states ($n\geq 3$) may be mixed with a rather large mixing angle, $\theta\sim30^\circ$ and this effect provides the correct values of $\Gamma_{ee}(\Upsilon(10580))$ and $\Gamma_{ee}(\Upsilon(11020))$. On the other hand, the $S-D$ mixing gives rise to an increase by two orders of magnitude of the di-electron widths of the mixed $\tilde\Upsilon(n^3D_1$) resonances ($n=3,4,5$), which originate from pure $D-$wave states. The value $\Gamma_{ee}(\tilde\Upsilon(3D))=0.095^{+0.028}_{-0.025}$ keV is obtained, being only $\sim 3$ times smaller than the di-electron width of $\Upsilon(10580)$, while $\Gamma_{ee}(\tilde\Upsilon(5D))\sim 135$ eV appears to be close to $\Gamma_{ee}(\Upsilon(11020))$ and therefore this resonance may become manifest in the $e^+e^-$ experiments. The mass differences between $M(nD)$ and $M((n+1)S) (n=4,5)$ are shown to be rather small, $50\pm 10$ MeV., Comment: 5 pages, no figures

The di-electron widths of $\psi(4040)$, $\psi(4160)$, and $\psi(4415)$, and their ratios are shown to be in good agreement with experiment, if in all cases the $S-D$ mixing with a large mixing angle $\theta\approx 34^\circ$ is taken. Arguments are presented why continuum states give small contributions to the wave functions at the origin. We find that the Y(4360) resonance, considered as a pure $3 {}^3D_1$ state, would have very small di-electron width, $\Gamma_{ee}(Y(4360))=0.060$ keV. On the contrary, for large mixing between the $4 {}^3S_1$ and $3 {}^3D_1$ states with the mixing angle $\theta=34.8^\circ$, $\Gamma_{ee}(\psi(4415))=0.57$ keV coincides with the experimental number, while a second physical resonance, probably Y(4360), has also a rather large $\Gamma_{ee} (Y(\sim 4400))=0.61$ keV. For the higher resonance Y(4660), considered as a pure $5 {}^3S_1$ state, we predict the di-electron width $\Gamma_{ee}(Y(4660))=0.70$ keV, but it becomes significantly smaller, namely 0.31 keV, if the mixing angle between the $5 {}^3S_1$ and $4 {}^3D_1$ states $\theta=34^\circ$. The mass and di-electron width of the $6 {}^3S_1$ charmonium state are calculated., Comment: 19 pages, no figures

The vector and pseudoscalar decay constants are calculated in the framework of the Field Correlator Method. Di-electron widths: $\Gamma_{ee}(J/\psi)=5.41$ keV, $\Gamma_{ee}(\psi'(3686))=2.47$ keV, $\Gamma_{ee}(\psi''(3770))=0.248$ keV, in good agreement with experiment, are obtained with the same coupling, $\alpha_s=0.165$, in QCD radiative corrections. We show that the larger $\alpha_s=0.191\pm 0.004$ is needed to reach agreement with experiment for $\Gamma_{\gamma\gamma}(\eta_c)=7.22$ keV, $\Gamma_{\gamma\gamma} (\chi(^3P_0))=3.3$ keV, $\Gamma_{\gamma\gamma}(\chi(^3P_2))= 0.54$ keV, and also for $\Gamma(J/\psi\to 3g)=59.5$ keV, $\Gamma(J/\psi\to \gamma 2g)=5.7$ keV. Meanwhile even larger $\alpha_s=0.238$ gives rise to good description of $\Gamma(\psi'\to 3g)=52.7$ keV, $\Gamma(\psi'\to \gamma 2g)= 3.5$ keV, and provides correct ratio of the branching fractions: $\frac{\mathcal{B}(J/\psi\to light hadrons)}{\mathcal{B}(\psi'\to light hadrons)}=0.24.$, Comment: 8 pages, no figures