This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Nucleon resonances and transition form factors

V. D. Burkert Thomas Jefferson National Accelerator Facility, Newport News, VA, USA

.1 Abstract

This is a contribution to the review ”50 Years of Quantum Chromodynamics” edited by F. Gross and E. Klempt, to be published in Journal EPJC. This contribution reviews the nucleon resonance transition form factors determined from meson electro-production experiments at electron accelerator facilities, i.e. this contribution focuses on ”space-like” transition form factors and amplitudes. Comparisons are made when available to LQCD and to approaches with traceable links to strong QCD and to advanced quark model calculations.

.2 Introduction and Formalism

Meson photoproduction has become an essential tool in the search for new excited light-quark baryon states. As discussed in the previous section, many new excited states have been discovered thanks to high precision photoproduction data in different final states Anisovich et al. (2017), and are now included in recent editions of the Review of Particle Physics (RPP) Workman (2022). The exploration of the internal structure of excited states and the effective degrees of freedom contributing to s-channel resonance excitation requires the use of electron beams, which is the subject of this contribution, where the virtuality (Q2Q^{2}) of the exchanged photon can be varied to pierce through the peripheral meson cloud and probe the quark core and its spatial structure. Electroproduction can thus say something about if a resonance is generated through short distance photon interaction with the small quark core, or through interaction with a more extended hadronic system.

The experimental exploration of resonance transition form factors reaches over 60 years with many review articles describing this history. Here we refer to a few recent ones Stoler (1993); Burkert and Lee (2004); Aznauryan and Burkert (2012a); Aznauryan et al. (2013). A review of recent electroproduction experiments in hadron physics and their interpretation within modern approaches of strong interaction physics can be found in Ref. Brodsky et al. (2020).

Electroproduction of final states with pseudoscalar mesons (e.g. NπN\pi, pηp\eta, KΛK\Lambda) have been employed at Jefferson Laboratory mostly with the CEBAF Large Acceptance Spectrometer (CLAS) operating at an instantaneous luminosity of 103410^{34}sec-1cm-2. In Hall A and Hall C, pairs of individual well-shielded focusing magnetic spectrometers are employed with more specialized aims and limited acceptance, but operating at much higher luminosity. This experimental program led to new insights into the scale dependence of effective degrees of freedom, e.g. meson-baryon, constituent quarks, and dressed quark contributions. Several excited states, shown in Fig. 1 assigned to their primary SU(6)O(3)SU(6)\otimes O(3) supermultiplets, have been studied this way, mostly with CLAS in Hall B.

Refer to caption
Figure 1: Excited states of the proton that have been studied in electroproduction to determine their resonance transition amplitudes or form factors. States highlighted in red are discussed in this subsection. Graphics from Ref. Burkert (2020).

Most of the resonance couplings have been extracted from single pseudoscalar meson production. In electroproduction, there are 6 complex helicity amplitudes, requiring a minimum of 11 independent measurements for a complete 111With the exception of an overall phase that cannot be determined model-independent determination of the amplitudes. In addition, measurements of isospin amplitudes require additional measurements. Following this, the complex amplitudes would need to be subjected to analyses of their phase motions to determine resonance masses on the (real) energy axis, or poles in the (complex) energy plane. Fortunately, in the lower mass range a variety of constraints can be applied to limit the number of unknowns when fitting the cross section data. These include the masses of states quite well known from hadronic processes or from meson photoproduction. Also, the number of possible angular momenta is limited to lπ2l_{\pi}\leq 2 in the examples discussed in the following. Additional constraints come from the Watson theorem Watson (1954) that relates the electromagnetic phases to the hadronic ones, and the use of dispersion relations, assuming the imaginary parts of the amplitude are given by the resonance contribution, and the real parts determined through dispersion integrals and additional pole terms. Other approaches use unitary isobar models that parameterize all known resonances and background terms, and unitarize the full amplitudes in a K-matrix procedure. In the following, we show results based on both approaches, where additional systematic uncertainties have been derived from the differences in the two procedures.

The availability of electron accelerators with the possibility of generating high beam currents at CEBAF at Jefferson Lab in the US and MAMI at Mainz in Germany, has enabled precise studies of the internal structure of excited states in the NN^{*} and the Δ\Delta^{*} sectors employing s-channel resonance excitations in large ranges of photon virtuality Q2Q^{2}. This has enabled probing the degrees of freedom relevant in the resonance excitation as a function of the distance scale probed. This is the topic we will elucidate in the following sections and the relevance to (approximations to) QCD. First we briefly discuss the formalism needed for a quantitative analysis of the single pseudoscalar meson electroproduction.

Refer to caption
Figure 2: The kinematics of single π+\pi^{+} electro-production off protons in the laboratory system.

.2.1 Formalism in the analysis of electroproduction of single pseudoscalar mesons

The simplest process used in the study of resonance transition amplitudes is single pion or kaon production, e.g. epeπ+nep\to e\pi^{+}n. Single π+\pi^{+} and π0\pi^{0} production are most suitable for the study of the lower-mass excited states as they couple dominantly to the excited states with masses up to 1.7\approx 1.7 GeV. It may then be useful to describe in more detail the analysis techniques and the formalism used. The unpolarized differential cross section for single pseudoscalar meson production can be written in the one-photon exchange approximation as:

dσdEfdΩedΩπ\displaystyle\frac{d\sigma}{dE_{f}d\Omega_{e}d\Omega_{\pi}} =\displaystyle= ΓdσdΩπ,\displaystyle\Gamma{\frac{d\sigma}{d\Omega_{\pi}}}\,, (1)

where Γ\Gamma is the virtual photon flux,

Γ=αem2π2Q2(W2M2)Ef2MEi11ϵ,\displaystyle\Gamma=\frac{\alpha_{em}}{2\pi^{2}Q^{2}}\frac{(W^{2}-M^{2})E_{f}}{2ME_{i}}\frac{1}{1-\epsilon}, (2)

where MM is the proton mass, WW the mass of the hadronic final state, ϵ\epsilon is the photon polarization parameter, Q2Q^{2} the photon virtuality, EiE_{i} and EfE_{f} represent the initial and the final electron energies, respectively. Moreover,

ϵ=[1+2(1+ν2Q2)tan2θe2]1\displaystyle\epsilon=\left[1+2\left(1+\frac{\nu^{2}}{Q^{2}}\right)\tan^{2}\frac{\theta_{e}}{2}\right]^{-1} (3)

and

dσdΩπ\displaystyle\frac{d\sigma}{d\Omega_{\pi}} =\displaystyle= σT+ϵσL+ϵσTTcos2ϕπ\displaystyle\sigma_{T}+\epsilon\sigma_{L}+\epsilon\sigma_{TT}\cos 2\phi_{\pi}
+2ϵ(1+ϵ)σLTcosϕπ.\displaystyle\qquad\qquad\qquad+\sqrt{2\epsilon(1+\epsilon)}\sigma_{LT}\cos{\phi}_{\pi}\,.

The kinematics for single π+\pi^{+} production is shown in Fig. 2.

The observables of the process γvpπN\gamma_{v}p\to\pi N^{\prime} can be expressed in terms of six parity-conserving helicity amplitudes  Aznauryan and Burkert (2012a); Walker (1969); Berends et al. (1967) :

Hi=λπ;λN|T|λγν;λp,\displaystyle{H_{i}=\left<\lambda_{\pi};\lambda_{N}|T|\lambda_{\gamma_{\nu}};\lambda_{p}\right>,} (4)

where λ\lambda denotes the helicity of the respective particle, λπ=0\lambda_{\pi}=0, λN=±12\lambda_{N}=\pm\frac{1}{2}, λγv=±1,0\lambda_{\gamma_{v}}=\pm 1,0, and λp=±12\lambda_{p}=\pm\frac{1}{2}, and HiH_{i} are complex functions of Q2Q^{2}, WW, and θπ\theta^{*}_{\pi}.

.2.2 Multipoles and partial wave decompositions

The response functions in (1) are given by:

σT\displaystyle\sigma_{T} =\displaystyle= pπW2KM(|H1|2+|H2|2+|H3|2+|H4|2),\displaystyle\frac{\vec{\it p}_{\pi}W}{2KM}(|H_{1}|^{2}+|H_{2}|^{2}+|H_{3}|^{2}+|H_{4}|^{2}), (5)
σL\displaystyle\sigma_{L} =\displaystyle= pπW2KM(|H5|2+|H6|2),\displaystyle\frac{\vec{\it p}_{\pi}W}{2KM}(|H_{5}|^{2}+|H_{6}|^{2}), (6)
σTT\displaystyle\sigma_{TT} =\displaystyle= pπW2KMRe(H2H3H1H4),\displaystyle\frac{\vec{\it p}_{\pi}W}{2KM}Re(H_{2}H_{3}^{*}-H_{1}H_{4}^{*}), (7)
σLT\displaystyle\sigma_{LT} =\displaystyle= pπW2KMRe[H5(H1H4)+H6(H2+H3)],\displaystyle\frac{\vec{\it p}_{\pi}W}{2KM}Re[H_{5}^{*}(H_{1}-H_{4})+H_{6}^{*}(H_{2}+H_{3})]\,, (8)

where pπ\vec{\it p}_{\pi} is the pion 3-momentum in the hadronic center-of-mass system, and KK is the equivalent real photon lab energy needed to generate a state with mass WW:

K=W2M22M.\displaystyle K=\frac{W^{2}-M^{2}}{2M}\,. (9)

The helicity amplitudes Hi,i=1H_{i},i=1–6, can be expanded into Legendre polynomials:

H1\displaystyle H_{1} =\displaystyle= 12sinθcosθ2l=1(Bl+B(l+1))(Pl′′Pl+1′′)\displaystyle\frac{1}{\sqrt{2}}\sin\theta\cos{\frac{\theta}{2}}\sum_{l=1}^{\infty}(B_{l+}-B_{(l+1)-})(P^{\prime\prime}_{l}-P^{\prime\prime}_{l+1})
H2\displaystyle H_{2} =\displaystyle= 2cosθ2l=1(Al+A(l+1))(PlPl+1)\displaystyle\sqrt{2}\cos{\frac{\theta}{2}}\sum_{l=1}^{\infty}(A_{l+}-A_{(l+1)-})(P^{\prime}_{l}-P^{\prime}_{l+1})
H3\displaystyle H_{3} =\displaystyle= 12sinθsinθ2l=1(Bl++B(l+1))(Pl′′+Pl+1′′)\displaystyle\frac{1}{\sqrt{2}}\sin\theta\sin{\frac{\theta}{2}}\sum_{l=1}^{\infty}(B_{l+}+B_{(l+1)-})(P^{\prime\prime}_{l}+P^{\prime\prime}_{l+1})
H4\displaystyle H_{4} =\displaystyle= 2sinθ2l=1(Al++A(l+1))(Pl+Pl+1)\displaystyle\sqrt{2}\sin{\frac{\theta}{2}}\sum_{l=1}^{\infty}(A_{l+}+A_{(l+1)-})(P^{\prime}_{l}+P^{\prime}_{l+1})
H5\displaystyle H_{5} =\displaystyle= 2cosθ2l=1(Cl+C(l+1))(PlPl+1)\displaystyle\sqrt{2}\cos{\frac{\theta}{2}}\sum_{l=1}^{\infty}(C_{l+}-C_{(l+1)-})(P^{\prime}_{l}-P^{\prime}_{l+1})
H6\displaystyle H_{6} =\displaystyle= 2sinθ2l=1(Cl++C(l+1))(Pl+Pl+1),\displaystyle\sqrt{2}\sin{\frac{\theta}{2}}\sum_{l=1}^{\infty}(C_{l+}+C_{(l+1)-})(P^{\prime}_{l}+P^{\prime}_{l+1})\,, (10)

where the Al+A_{l+} and Bl+B_{l+} etc., are the transverse partial wave helicity elements for λγp=12\lambda_{\gamma p}=\frac{1}{2} and λγp=32\lambda_{\gamma p}=\frac{3}{2}, and C±C_{\pm} the longitudinal partial wave helicity elements. In the subscript, l+l+ and (l+1)(l+1)- define the π\pi orbital angular momenta, and the sign ±\pm is related to the total angular momentum J=lπ±12J=l_{\pi}\pm\frac{1}{2}. In the analysis of data on the NΔ(1232)N\Delta(1232) transition, linear combinations of partial wave helicity elements are expressed in terms of electromagnetic multipoles:

Ml+\displaystyle M_{l+} =\displaystyle= 12(l+1)[2Al+(l+2)Bl+]\displaystyle\frac{1}{2(l+1)}[2A_{l+}-(l+2)B_{l+}] (11)
El+\displaystyle E_{l+} =\displaystyle= 12(l+1)(2Al++lBl+)\displaystyle\frac{1}{2(l+1)}(2A_{l+}+lB_{l+}) (12)
Ml+1,\displaystyle M_{l+1,-} =\displaystyle= 12(l+1)(2Al+1,lBl+1,)\displaystyle\frac{1}{2(l+1)}(2A_{l+1,-}-lB_{l+1,-}) (13)
El+1,\displaystyle E_{l+1,-} =\displaystyle= 12(l+1)[2Al+1,+(l+2)Bl+1,]\displaystyle\frac{1}{2(l+1)}[-2A_{l+1,-}+(l+2)B_{l+1,-}] (14)
Sl+\displaystyle S_{l+} =\displaystyle= 1l+1Q2Q2Cl+\displaystyle\frac{1}{l+1}\sqrt{\frac{{\vec{\it Q^{*}}^{2}}}{Q^{2}}}C_{l+} (15)
Sl+1,\displaystyle S_{l+1,-} =\displaystyle= 1l+1Q2Q2Cl+1,,\displaystyle\frac{1}{l+1}\sqrt{\frac{{\vec{\it Q^{*}}^{2}}}{Q^{2}}}C_{l+1,-}\,, (16)

where Q\vec{\it Q}^{*} is the photon 3-momentum in the hadronic rest frame. The electromagnetic multipoles are often used to describe the transition from the nucleon ground state to the Δ(1232)\Delta(1232), which is dominantly described as a magnetic dipole transition M1+M_{1+}. The electromagnetic multipoles as well as the partial wave helicity elements are complex quantities and contain both non-resonant and resonant contributions. In order to compare the results to model predictions and LQCD, an additional analysis must be performed to separate the resonant parts A^±\hat{\rm A}_{\pm}, B^±\hat{\rm B}_{\pm}, etc., from the non-resonant parts of the amplitudes. In a final step, the known hadronic properties of a given resonance can be used to determine photocoupling helicity amplitudes that characterize the electromagnetic vertex:

A^l±\displaystyle\hat{A}_{l\pm} =\displaystyle= FCπNIA1/2,\displaystyle\mp FC^{I}_{\pi N}A_{1/2}, (17)
B^l±\displaystyle\hat{B}_{l\pm} =\displaystyle= ±F16(2j1)(2j+3)CπNIA3/2,\displaystyle\pm F\sqrt{\frac{16}{(2j-1)(2j+3)}}C^{I}_{\pi N}A_{3/2}, (18)
S^l±\displaystyle\hat{S}_{l\pm} =\displaystyle= F222J+1CπNIS1/2,\displaystyle-F\frac{2\sqrt{2}}{2J+1}C^{I}_{\pi N}S_{1/2}, (19)
F\displaystyle F =\displaystyle= 1(2j+1)πKpπΓπΓ2\displaystyle\sqrt{\frac{1}{(2j+1)}\pi}\frac{K}{p_{\pi}}\frac{\Gamma_{\pi}}{\Gamma^{2}}

where the CπNIC^{I}_{\pi N} are isospin coefficients. The total transverse absorption cross section for the transition into a specific resonance is given by:

σT=2MWRΓ(A1/22+A3/22).\displaystyle\sigma_{T}=\frac{2M}{W_{R}\Gamma}(A^{2}_{1/2}+A^{2}_{3/2}). (20)

Experiments in the region of the Δ(1232)32+\Delta(1232)\frac{3}{2}^{+} resonance often determine the electric quadrupole ratio REMR_{EM}

REM=Im(E1+)Im(M1+)R_{EM}=\frac{Im(E_{1+})}{Im(M_{1+})} (21)

and the scalar quadrupole ratio RSMR_{SM}

RSM=Im(S1+)Im(M1+)R_{SM}=\frac{Im(S_{1+})}{Im(M_{1+})} (22)

where E1+E_{1+}, S1+S_{1+}, and M1+M_{1+} are the electromagnetic transition multipoles at the mass of the Δ(1232)32+\Delta(1232)\frac{3}{2}^{+} resonance.

It is worth noting that the electric and the scalar transition amplitudes are connected at the so-called pseudo-threshold Qpt2=(WM)2Q^{2}_{pt}=-(W-M)^{2} through the Siegert theorem Siegert (1937), whose impact on the extraction of nucleon resonance transition amplitudes from pion electro-production is discussed in Tiator (2016); Ramalho (2016).

.2.3 Resonance analysis tools

A model-independent determination of the amplitudes contributing to the electro-excitation of resonances in single pseudoscalar pion production epeNπep\to e^{\prime}N\pi (see kinematics of single pion production in Fig. 2) requires a large number of independent measurements at each value of the electron kinematics WW, Q2Q^{2}, the hadronic cms angle cosθπ\cos{\theta^{\pi}}, and the azimuthal angle ϕπ\phi^{\pi} describing the angle between the electron scattering plane and the hadronic decay plane. Such a measurement requires full exclusivity of the final state and employing both polarized electron beams and the measurements of the nucleon recoil polarization.

Refer to caption
Figure 3: JLab/Hall A data for epepπ0\vec{e}p\to e\vec{p}\pi^{0} response functions at W=1.232 GeV and Q2=1.0Q^{2}=1.0 Kelly et al. (2007). Notations refer to transverse (t), normal (n) and longitudinal (l) components of the proton recoil polarization. The curves correspond to results obtained using SAID (short dashed), MAID (dashed-dotted), and the dynamical models DMT Kamalov et al. (2001) (dotted), and SL Sato and Lee (2001) (long-dashed/green). The other curves correspond to Legendre and multipole fits performed by the authors.

Such measurements would in general require full 4π4\pi coverage for the hadronic final state. The only measurement that could claim to be complete was carried out at Jefferson Lab in Hall A Kelly et al. (2007) employing a limited kinematics centered at resonance for epepπ0\vec{e}p\to e^{\prime}\vec{p}\pi^{0} at W=1.232W=1.232 GeV, and Q21Q^{2}\approx 1 GeV2. Figure 5 shows the 16 response functions extracted from this measurement. The results of this measurement in terms of the magnetic NΔN\Delta transition form factor and the quadrupole ratios are included in Fig. 4 among other data. They coincide very well with results of other experiments Aznauryan et al. (2009); Ungaro et al. (2006); Joo et al. (2002); Frolov et al. (1999) using different analysis techniques that may be also applied to broader kinematic conditions, especially higher mass resonances. Details of the latter are discussed in  Aznauryan and Burkert (2012a); Tiator et al. (2011). We briefly summarize them here:

Refer to caption   Refer to caption Refer to caption
Figure 4: The NΔ(1232)N\Delta(1232) transition amplitudes. Left: The magnetic NΔ\Delta transition form factor normalized to the dipole form factor and compared with the Light-Front Relativistic Quark Model (LFRQM)Aznauryan and Burkert (2012b, 2016) with running quark mass, and with results using the Dyson-Schwinger Equation Segovia et al. (2014). Both predictions are close to the data at high Q2Q^{2}. At Q2<3Q^{2}<3GeV2 meson-baryon contributions are significant. Middle: The electric (top) and scalar (bottom) quadrupole/magnetic-dipole ratios REMR_{EM} and RSMR_{SM}. Right: REMR_{EM} and RSMR_{SM} from Lattice QCD Alexandrou et al. (2008); Behrndt and Cvetic (2005) compared to data in the low Q2Q^{2} domain.
  • Dispersion Relations have been employed in two ways: One is based on fixed-t dispersion relations for the invariant amplitudes and was successfully used throughout the nucleon resonance region. Another way is based on DR for the mulipole amplitudes of the Δ(1232)\Delta(1232) resonance, and allows getting functional forms of these amplitudes with one free parameter for each of them. It was employed for the analysis of the more recent data.

  • The Unitary Isobar Model (UIM) was developed in  Drechsel et al. (1999) from the effective Lagrangian approach for pion photoproduction Peccei (1969). Background contributions from t-channel ρ\rho and ω\omega exchanges are introduced and the overall amplitude is unitarized in a K-matrix approximation.

  • Dynamical Models have been developed, as SAID from pion photoproduction data Arndt et al. (2002), the Sato-Lee model was developed in Sato and Lee (1996). Its essential feature is the consistent description of πN\pi N scattering and the pion electroproduction from nucleons. It was utilized in the study of Δ(1232)\Delta(1232) excitations in the epepπ0ep\to ep\pi^{0} channel Sato and Lee (2001). The Dubna-Mainz-Taipei model Kamalov and Yang (1999) builds unitarity via direct inclusion of the πN\pi N final state in the T-matrix of photo- and electroproduction.

.3 Light-quark baryon resonance electroproduction

In order to learn from the meson electroproduction data about the internal spin and spatial electromagnetic structure, it is essential to have advanced models available with links to the fundamentals of QCDQCD.

Refer to caption
Figure 5: Sample of results of an analysis by the JLab group of the Legendre moments of epeπ+n\vec{e}p\to e\pi^{+}n structure functions in comparison with experimental data Park et al. (2008) at Q2=2.44Q^{2}=2.44 GeV2. The solid (dashed) curves correspond to results obtained using the DR (UIM) approach.

While most of the analyses have focused on single pseudoscalar meson production, such as γvpNπ,pη,KΛ\gamma_{v}p\to N\pi,p\eta,K\Lambda, KΣK\Sigma, more recent work included the pπ+πp\pi^{+}\pi^{-} final state both in real photoproduction Golovatch et al. (2019) as well as in electroproduction Mokeev et al. (2020). The 2-pion final state has more sensitivity to excited NN^{*} and Δ\Delta^{*} states in the mass range above 1.6 GeV, with several states dominantly coupling to NππN\pi\pi final states, enabling the study of their electromagnetic transition form factors in the future.

.3.1 The NΔ(1232)32+N\Delta(1232)\frac{3}{2}^{+} transition

The Δ++\Delta^{++} isobar was first observed 70 years ago in Enrico Fermi’s experiment that used a π+\pi^{+} meson beam scattered off the protons in a hydrogen target Anderson et al. (1952). The cross section showed a sharp rise above threshold towards a mass near 1200 MeV. While the energy of the pion beam was not high enough to see the maximum and the fall-off following the peak, a strong indication of the first baryon resonance was observed. It took 12 more years and the development of the underlying symmetry in the quark model before a microscopic explanation of this observation could emerge. There was, however, a problem; while the existence of the Δ+,0,\Delta^{+,0,-} could be explained within the model, the existence of the Δ(1232)++\Delta(1232)^{++}, which within the quark model would correspond to a state |uuu|u{\uparrow}u{\uparrow}u{\uparrow}\rangle, was forbidden as it would have an overall symmetric wave function. It took the introduction of para Fermi statistics Greenberg (1964) what later became ”color”, to make the overall wave function anti-symmetric. In this way the Δ++(1232)\Delta^{++}(1232) resonance may be considered a harbinger of the development of QCD.

The nucleon to Δ(1232)32+\Delta(1232)\frac{3}{2}^{+} transition is now well measured in a large range of Q2Q^{2} Ungaro et al. (2006); Joo et al. (2002); Frolov et al. (1999). At the real photon point, it is explained by a dominant magnetic transition from the nucleon ground state to the Δ(1232)\Delta(1232) excited state. Additional contributions are related to small D-wave components in both the nucleon and the Δ(1232)\Delta(1232) wave functions leading to electric quadrupole and scalar quadrupole transitions. These remain in the few % ranges at small Q2Q^{2}. The magnetic transition is to 65%\approx 65\% given by a simple spin flip of one of the valence quarks as seen in Fig. 4. The remaining 35% of the magnetic dipole strength is attributed to meson-baryon contributions. There exist extensive theoretical reviews of the NΔ(1232)N\Delta(1232) transition in the lower Q2Q^{2} range are available Pascalutsa and Vanderhaeghen (2006), and more recent reviews that cover the full Q2Q^{2} covered by data Aznauryan and Burkert (2012a); Tiator et al. (2011).

The electric quadrupole ratio REMR_{EM} was found as:

REM0.02.\displaystyle R_{EM}\approx-0.02. (23)

There has been a longstanding prediction of asymptotic pQCD, that REM+1R_{EM}\to+1 at Q2Q^{2}\to\infty. Results on the magnetic transition form factor GMn,Ash\rm G_{Mn,Ash}, defined in the Ash convention Ash et al. (1967), and on the quadrupole transition ratios are shown in Fig. 4.

Refer to captionRefer to caption Refer to caption
Figure 6: Helicity transition amplitudes for the proton to Roper N(1440)12+N(1440)\frac{1}{2}^{+} excitation in units of 103GeV1/210^{-3}\rm{GeV}^{-1/2} compared to various model calculations; see text. Left: Transverse A1/2A_{1/2} amplitude. Middle: Scalar S1/2S_{1/2} amplitude. Right: Helicity amplitudes of the Roper resonance at low Q2Q^{2}. Data are compared to calculations within Effective Field Theory Bauer et al. (2014) shown in solid black lines. The other broken lines are parts of the full calculations. The data are from Aznauryan et al. (2009); Mokeev et al. (2012); Drechsel et al. (2007). The open red circle at Q20.1Q^{2}\approx 0.1GeV2 is the result of an analysis of epepπ0ep\to ep\pi^{0} data from Štajner et al. (2017).

GMn,Ash\rm G_{Mn,Ash} is shown normalized to the dipole form factor, but indicates a much faster Q2Q^{2} fall-off compared to that. In comparison to the advanced LFRQM with momentum-dependent constituent quark mass, and with the Dyson-Schwinger Equation (DSE-QCD) results, there is good agreement at the high-Q2Q^{2} end of the data. The discrepancy at small Q2=0Q^{2}=0 is likely due to the meson-baryon contributions at small Q2Q^{2}, which are not modeled in either of the calculations.

The quadrupole ratio REMR_{EM} shows no sign of departing significantly from its value at Q2=0Q^{2}=0, even at the highest Q26.5Q^{2}\approx 6.5 GeV2. Both calculations barely depart from REM=0R_{EM}=0, and remain near zero at all Q2>2Q^{2}>2 GeV2. This indicates that the negative constant value shown by the data is likely due to meson-baryon contributions that are not included in the theoretical models. For the scalar quadrupole ratio RSMR_{SM} the asymptotic prediction in holographic QCD (hQCD)  Grigoryan et al. (2009) is:

RSM=ImS1+ImM1+1,atQ2,\displaystyle R_{SM}=\frac{ImS_{1+}}{ImM_{1+}}\to-1,~{}~{}{\rm at~{}Q^{2}\to\infty}, (24)

while REMR_{EM} in hQCD is predicted to approach +1 asymptotically. The RSMR_{SM} data show indeed a strong trend towards increasing negative values at larger Q2Q^{2}, semi-quantitatively described by both calculations at Q2<4Q^{2}<4 GeV2. The Dyson-Schwinger equation approach predicts a flattening of RSMR_{SM} at Q2>4Q^{2}>4 GeV2, while the Light Front relativistic Quark Model predicts a near constant negative slope of RSM(Q2)R_{SM}(Q^{2}) also at higher Q2Q^{2}.

.3.2 The Roper resonance N(1440)12+N(1440)\frac{1}{2}^{+}

The Roper resonance, discovered in 1964 Roper (1964) in a phase shift analysis of elastic πN\pi N scattering data, has been differently interpreted for half a century. In the non-relativistic quark model (nrQM), the state is the first radial excitation of the nucleon ground state with a mass expected around 1750 MeV, much higher than the measured Breit-Wigner mass of 1440\approx 1440 MeV. This discrepancy is now understood as the consequence of a dynamical coupled channel effect that shifts the mass below the mass of the N(1535)1/2N(1535){1/2}^{-} state, the negative-parity partner of the nucleon Suzuki et al. (2010). Another problem with the quark model was the sign of the transition form factor A1/2(Q2=0)A_{1/2}(Q^{2}=0), predicted in the nrQM as large and positive, while experimental analyses showed a negative value.

These discrepancies resulted in different interpretations of the state that could only be resolved with electroproduction data from CLAS at Jefferson Lab, the development of continuous QCD approximations in the Dyson-Schwinger equation approach Segovia et al. (2015) and Light Front Relativistic QM with momentum-dependent quark masses Aznauryan and Burkert (2012b) shown in Fig. 6, and Lattice data Mathur et al. (2005); Lin and Cohen (2012). A recent review of the history and current status of the Roper resonance, is presented in a colloquium-style article published in Review of Modern Physics Burkert and Roberts (2019).

Descriptions of the baryon resonance transitions form factors, including the Roper resonance N(1440)12+N(1440)\frac{1}{2}^{+}, have also been carried out within holographic models de Teramond and Brodsky (2012); Ramalho and Melnikov (2018). In the range Q2<0.6Q^{2}<0.6 GeV2, calculations based on meson-baryon degrees of freedom and effective field theory Bauer et al. (2014) have been successfully performed, as may be seen in Fig. 6. Earlier model descriptions, such as the Isgur-Karl model that describe the nucleon as a system of 3 constituent quarks in a confining potential and a one-gluon exchange contribution leading to a magnetic hyperfine splitting of states Isgur and Karl (1978, 1977), and the relativized version of Capstick Capstick and Isgur (1986) have popularized the model that became the basis for many further developments and variations, e.g. the light front relativistic quark model, and the hypercentral quark model Giannini et al. (2003). Other models were developed in parallel. The cloudy bag model Thomas et al. (1981) describes the nucleon as a bag of 3 constituent quarks surrounded by a cloud of pions. It has been mostly applied to nucleon resonance excitations in real photoproduction, Q2=0Q^{2}=0 Bermuth et al. (1988); Thomas et al. (1981), with some success in the description of the Δ(1232)32+\Delta(1232)\frac{3}{2}^{+} and the Roper resonance transitions.

There is agreement with the data at Q2>1.5Q^{2}>1.5 GeV2 for these two states, while the meson-baryon contributions for the Δ(1232)\Delta(1232) are more extended, and agreement with the quark based calculations is reached at Q2>4Q^{2}>4 GeV2. The calculations deviate significantly from the data at lower Q2Q^{2}, which indicates the presence of non-quark core effects. For the Roper resonance such contributions have been described successfully in dynamical meson-baryon models Obukhovsky et al. (2011) and in effective field theory Bauer et al. (2014). Calculations on the Lattice for the N-Roper transition form factors F1pRF_{1}^{pR} and F2pRF_{2}^{pR}, which are combinations of the transition amplitudes A1/2A_{1/2} and S1/2S_{1/2}, have been carried out with dynamical quarks Lin and Cohen (2012). The results agree well with the data in the range Q2<1.0Q^{2}<1.0 GeV2, where data and calculations overlap Fig. 7.

New electroproduction data on the Roper Štajner et al. (2017) and on several higher mass states have been obtained in the 2-pion channel, specifically in epepπ+πep\to e^{\prime}p\pi^{+}\pi^{-} Mokeev et al. (2016).

Refer to caption
Figure 7: Dirac and Pauli transition form factors F1F_{1} and F2F_{2} for the proton to N(1440)1/2+N(1440)1/2^{+} transition compared to Lattice QCD calculations Lin and Cohen (2012) with pion masses (in GeV): 0.39 (red squares), 0.45 (orange triangles), and 0.875 (green circles) on the Nf=2+1N_{f}=2+1 anisotropic lattices, compared to CLAS results (black circles). The F1F_{1} and F2F_{2} form factors are linear combinations of the A1/2A_{1/2} and S1/2S_{1/2} amplitudes.

The mass of the Roper state has been computed on the Lattice and extrapolated to the physical pion mass, showing good agreement with the physical value measured with a Breit-Wigner parametrization. It should be noted that the Roper mass measured at the pole in the complex plane is significantly different from the value obtained using the BW ansatz.

Refer to caption
Figure 8: The transverse helicity transition amplitudes of N(1520)32N(1520)\frac{3}{2}^{-} versus Q2Q^{2}, compared to the LFRQM, A1/2A_{1/2} (left), A3/2A_{3/2} (middle). The shaded area indicates the contribution from non-quark contributions as estimated from the difference of the measured data points and the LF RQM contribution, likely due to hadronic contributions. Right: Helicity asymmetry AhelA_{hel}, as defined in Eq. 26. Graphics from Ref. Aznauryan and Burkert (2012a)

Supported by an extensive amount of single pion electroproduction data, covering the full phase space in the pion polar and azimuthal center-of-mass angles, and accompanied by several theoretical modeling, we can summarize our current understanding of the N(1440)12+N(1440)\frac{1}{2}^{+} state as follows:

  • The Roper resonance is, at heart, the first radial excitation of the nucleon.

  • It consists of a well-defined dressed-quark core, which plays a role in determining the system’s properties at all length scales, but exerts a dominant influence on probes with Q2>mN2Q^{2}>m_{N}^{2}, where mNm_{N} is the nucleon mass.

  • The core is augmented by a meson cloud, which both reduces the Roper’s core mass by 20%\approx 20\%, thereby solving the mass problem that was such a puzzle in constituent quark model treatments, and, at low Q2Q^{2}, contributes an amount to the electroproduction transition form factors that is comparable in magnitude with that of the dressed quark core, but vanishes rapidly as Q2Q^{2} is increased beyond mN2m_{N}^{2}.

As stated in the conclusions of  Burkert and Roberts (2019): ”The fifty years of experience with the Roper resonance have delivered lessons that cannot be emphasized too strongly. Namely, in attempting to predict and explain the QCD spectrum, one must fully consider the impact of meson-baryon final state interactions and the coupling between channels and states that they generate, and look beyond merely locating the poles in the S-matrix, which themselves reveal little structural information, to also consider the Q2Q^{2} dependencies of the residues, which serve as a penetrating scale-dependent probe of resonance composition.”

.3.3 The helicity structure of the N(1520)32N(1520)\frac{3}{2}^{-}

The N(1520)32N(1520)\frac{3}{2}^{-} state corresponds to the lowest excited nucleon resonance with JP=32J^{P}=\frac{3}{2}^{-}. Its helicity structure is defined by the Q2Q^{2} dependence of the two transverse transition amplitudes A1/2A_{1/2} and A3/2A_{3/2}. They are both shown in Fig. 8. A particularly interesting feature of this state is that at the real photon point, A3/2A_{3/2} is strongly dominant, while A1/2A_{1/2} is very small. However, at high Q2Q^{2}, A1/2A_{1/2} is becoming dominant, while A3/2A_{3/2} drops rapidly. This behavior is qualitatively consistent with the expectation of asymptotic QCD, which predicts the transition helicity amplitudes to behave like:

A1/2aQ3,A3/2bQ5.\displaystyle\hskip 71.13188ptA_{1/2}\propto\frac{a}{Q^{3}},A_{3/2}\propto\frac{b}{Q^{5}}\,. (25)

The helicity asymmetry

Ahel=A1/22A3/22A1/22+A3/22,\displaystyle\hskip 71.13188ptA_{hel}=\frac{A_{1/2}^{2}-A_{3/2}^{2}}{A_{1/2}^{2}+A_{3/2}^{2}}, (26)

shown in Fig. 8, illustrates this rapid change in the helicity structure of the γvpN(1520)3/2\gamma_{v}pN(1520){3/2}^{-} transition. At Q2>2Q^{2}>2~{}GeV2, A1/2A_{1/2} fully dominates the process.

Refer to caption
Figure 9: Right: The transverse transition helicity amplitude A1/2A_{1/2} versus Q2Q^{2}. At Q2>2Q^{2}>2GeV2 the data are well described by the light-cone sum rules LCSR Braun et al. (2009). The light front relativistic quark model (LFRQM) Aznauryan and Burkert (2017) describes that data at Q2>1Q^{2}>1GeV2.

.3.4 Transition Form Factors of N(1535)12N(1535)\frac{1}{2}^{-} - A state with a hard quark core.

This state is the parity partner state to the ground state nucleon, with the same spin 1/2 but with opposite parity, its quark content requires an orbital L=1 excitation in the transition from the proton. In the SU(6)O(3)SU(6)\otimes O(3) symmetry scheme, the state is a member of the [70,1][70,1^{-}] super multiplet. This state couples equally to NπN\pi and to NηN\eta final state. It has therefore be probed using both decay channels epepηep\to ep\eta and epeNπ+,0ep\to eN\pi^{+,0}. Because of isospin I=1/2I=1/2 for nucleon states, the coupling to the charged π+n\pi^{+}n channel is preferred over π0p\pi^{0}p owing to the Clebsch-Gordon coefficients.

The A1/2A_{1/2} helicity amplitude for the γpN(1535)12\gamma pN(1535)\frac{1}{2}^{-} resonance excitation shown in Fig. 9 represents the largest range in Q2Q^{2} of all nucleon states for which resonance transition form factors have been measured as part of the broad experimental program at JLab.

For this state, as well as for the N(1440)12+N(1440)\frac{1}{2}^{+} state, advanced relativistic quark model calculations Aznauryan and Burkert (2015a), DSE-QCD calculations Segovia et al. (2015) and Light Cone sum rule results Anikin et al. (2015) are available, employing QCD-based modeling of the excitation of the quark core for the first time.

Refer to caption Refer to caption Refer to caption
Figure 10: The transverse amplitudes of the proton to N(1675)52N(1675)\frac{5}{2}^{-} transition compared to the LF RQM Aznauryan and Burkert (2015b), hypercentral QM Santopinto and Giannini (2012), and contributions from meson-baryon (MB) coupled channel dynamics Julia-Diaz et al. (2008). Left: A1/2A_{1/2}, Middle:A3/2A_{3/2}. Both quark models predict very small amplitudes for the proton, while the meson-baryon contributions estimate is large and is close to the data. Right: A1/2A_{1/2} for neutron target (only photoproduction data available) compared to the LFRQM and hCQM. Both quark models predict large amplitudes for neutrons, more than factor 10 compared to protons at Q2=0Q^{2}=0. Assuming similar meson-baryon contributions as in the proton case with opposite sign could quantitatively explain the single measured value at the photon point.

The transverse transition amplitude A1/2A_{1/2} of N(1535)12N(1535)\frac{1}{2}^{-} is a prime example of the power of meson electroproduction to unravel the internal structure of the resonance transition. In another section of the Volume ”50 Years of Quantum Chromodynamics”  Burkert et al. (2022), the nature of this state is discussed as a possible example of a dynamically generated resonance. The electroproduction data shown here reveal structural aspects of the state and its nature that require a different interpretation. The transition form factor A1/2A_{1/2} of the state, shown in Fig. 9, is quantitatively reproduced over a large range in Q2Q^{2} by two alternative approaches, the LFRQM and the LCSR. Both calculations are based on the assumptions of the presence of a 3-quark core. Notice that there is deviation from the quark calculations at Q2<12Q^{2}<1-2 GeV2, highlighted as the shaded area in Fig. 9, which may be assigned to the presence of non-quark contributions. Attempts to compute the transition form factors within strictly dynamical models have not succeeded in explaining the available data Jido et al. (2008). The discrepancy could be resolved if the character of the probe, meson (pion) in the case of hadron interaction and short range photon interaction in the case of electroproduction, probe different parts of the resonance’s spatial structure: peripheral in case of meson scattering and short distance behavior in electro-production. The peripheral meson scattering and low Q2Q^{2} meson photo-production reveal the dynamical features of the state, whereas high Q2Q^{2} electroproduction reveals the structure of the quark core.

.3.5 The N(1675)52N(1675)\frac{5}{2}^{-} resonance reveals the meson-baryon contributions

In previous discussions we have concluded that meson-baryon degrees of freedom provide significant strength to the resonance excitation in the low-Q2Q^{2} domain where quark based approaches LF RQM, DSE/QCD, and LC SR calculations fail to reproduce the transition amplitudes quantitatively. Our conclusion rests, in part, with this assumption. But, how can we be certain of the validity of this conclusion?

The N(1675)52N(1675)\frac{5}{2}^{-} resonance allows us to test this assumption, quantitatively. Figure 10 shows our current knowledge of the transverse helicity amplitudes A1/2(Q2)A_{1/2}(Q^{2}) and A3/2(Q2)A_{3/2}(Q^{2}), for proton target compared to RQM Aznauryan and Burkert (2017) and hypercentral CQM Santopinto and Giannini (2012) calculations. The specific quark transition for a JP=52J^{P}=\frac{5}{2}^{-} state belonging to the [SU(6)O(3)]=[70,1][SU(6)\otimes O(3)]=[70,1^{-}] supermultiplet configuration, in non-relativistic approximation prohibits the transition from the proton in a single quark transition. This suppression, known as the Moorhouse selection rule Moorhouse (1966), is valid for the transverse transition amplitudes A1/2A_{1/2} and A3/2A_{3/2} at all Q2Q^{2}. It should be noted that this selection rule does apply to the transition from a proton target, it does not apply to the transition from the neutron, which is consistent with the data. Modern quark models that go beyond single quark transitions, confirm quantitatively the suppression resulting in very small amplitudes from protons but large ones from neutrons. Furthermore, a direct computation of the hadronic contribution to the transition from protons confirms this (Fig. 10). The measured helicity amplitudes off the protons are almost exclusively due to meson-baryon contributions as the dynamical coupled channel (DCC) calculation indicates (dashed line). The close correlation of the DCC calculation and the measured data for the case when quark contributions are nearly absent, supports the phenomenological description of the helicity amplitudes in terms of a 3-quark core that dominate at high Q2Q^{2} and meson-baryon contributions that can make important contributions at lower Q2Q^{2}.

.3.6 Resonance lightfront transition charge densities.

Knowledge of the helicity amplitudes in a large Q2Q^{2} allows for the determination of the transition charge densities on the light cone in transverse impact parameter space (bx,byb_{x},b_{y}Tiator and Vanderhaeghen (2009). The relations between the helicity transition amplitudes and the Dirac and Pauli resonance transition form factors are given by:

A1/2=eQK(4MNM)1/2{F1NN+F2NN}\displaystyle A_{1/2}=e\frac{Q_{-}}{\sqrt{K}(4M_{N}M^{*})^{1/2}}\{F_{1}^{NN^{*}}+F_{2}^{NN^{*}}\}~{}~{}~{}~{}~{} (27)
S1/2=eQK(4MNM)1/2(Q+Q2M)(M+MN)Q2\displaystyle S_{1/2}=e\frac{Q_{-}}{\sqrt{K}(4M_{N}M^{*})^{1/2}}\left(\frac{Q_{+}Q_{-}}{2M^{*}}\right)\frac{(M^{*}+M_{N})}{Q^{2}}
×{F1NNQ2(M+MN)2F2NN},\displaystyle\times\{F_{1}^{NN^{*}}-\frac{Q^{2}}{(M^{*}+M_{N})^{2}}F_{2}^{NN^{*}}\}\,,~{}~{}~{}~{}~{} (28)

where MM^{*} is the mass of the excited state NN^{*}, K=M2MN22MK=\frac{{M^{*}}^{2}-M^{2}_{N}}{2M^{*}} is the equivalent photon energy, Q+Q_{+} and QQ_{-} are short hands for Q±M±MN)2+Q2Q_{\pm}\equiv\sqrt{M^{*}\pm M_{N})^{2}+Q^{2}}. The charge and magnetic lightfront transition densities ρ0NN\rho_{0}^{NN^{*}} and ρTNN\rho_{T}^{NN^{*}}, respectively, are given as:

ρ0NN(b)=0dQ2πJ0(bQ)F1NN(Q2)\displaystyle\rho_{0}^{NN^{*}}(\vec{b})=\int_{0}^{\infty}\frac{dQ}{2\pi}J_{0}(bQ)F_{1}^{NN^{*}}(Q^{2}) (29)
ρTNN(b)=ρ0NN(b)+sin(ϕbϕs)×\displaystyle\rho_{T}^{NN^{*}}(\vec{b})=\rho_{0}^{NN^{*}}(\vec{b})+\sin(\phi_{b}-\phi_{s})\times
0dQ2πQ2(M+MN)J1(bQ)F2NN(Q2).\displaystyle\int_{0}^{\infty}\frac{dQ}{2\pi}\frac{Q^{2}}{(M^{*}+M_{N})}J_{1}(bQ)F_{2}^{NN^{*}}(Q^{2})\,. (30)

Refer to caption

Figure 11: Left panels: N(1440)N(1440), top: projection of charge densities on byb_{y}, bottom: transition charge densities when the proton is spin polarized along bxb_{x}. Right panels: same for N(1535)N(1535). Note that the densities are scaled with b2b^{2} to emphasize the outer wings. Color code:negative charge is blue, positive charge is red. Note that all scales are the same for ease of comparison Burkert (2018b). Graphics credit: F.X. Girod.

Similar transverse charge transition densities can be defined for JP=32+J^{P}=\frac{3}{2}^{+} states such as the Δ(1232)32+\Delta(1232)\frac{3}{2}^{+}. This has been studied in Carlson and Vanderhaeghen (2008) and is shown in Fig. 12.

A comparison of N(1440)12+N(1440)\frac{1}{2}^{+} and N(1535)12N(1535)\frac{1}{2}^{-} is shown in Figure 11. There are clear differences in the charge transition densities between the two states. The Roper state has a softer positive core and a wider negative outer cloud than N(1535)12N(1535)\frac{1}{2}^{-} and develops a larger shift in byb_{y} when the proton is polarized along the bxb_{x} axis.

.3.7 Single Quark Transition Model

Many of the exited states for which there is information about the transition form factors available have been assigned as members of the [SU(6),LP][SU(6),L^{P}] = [70,1][70,1^{-}] super multiplet of the [SU(6)O(3)][SU(6)\otimes O(3)] symmetry group. In a model, where only single quark transitions to the excited states are considered Hey and Weyers (1974); Cottingham and Dunbar (1979); Burkert et al. (2003), only 3 of the amplitudes need to be known to determine the remaining 16 transverse helicity amplitudes for all states in [70,1][70,1^{-}] including on neutrons. However, the picture is now more complicated due to the strong admixture of meson-baryon components to the single quark transition especially in the lower Q2Q^{2} range. This requires a model to separate the single quark contributions from the hadronic part before projections for other states can be made Ramalho (2014).

.4 Higher mass baryons and hybrid baryons

The existence of baryons containing significant active gluonic components in the wave function has been predicted some decade ago Dudek and Edwards (2012) employing Lattice QCD simulations. The lowest such ”hybrid” state is expected to be a JP=12+J^{P}=\frac{1}{2}^{+} nucleon state. LQCD projects a mass of 1.3 GeV above the nucleon mass, i.e. approximately 2.2-2.3 GeV, and several other states should appear close by in JPJ^{P} = 12+\frac{1}{2}^{+} and JPJ^{P} = 32+\frac{3}{2}^{+}, as seen in Fig. 13.

How do we identify these states? Hybrid baryons have same spin-parity as other ordinary baryons. In contrast to hybrid mesons, there are no hybrid baryons with ”exotic” quantum numbers. One possibility is to search for more states than the quark model predicts in some mass range. The other possibility is to study the transition form factors of excited states. Hybrid states may be identified as states with a different Q2Q^{2} behavior than what is expected from a 3-quark state. The sensitivity Li et al. (1992) is demonstrated for the Roper resonance that projected a very rapid drop of the A1/2(Q2)A_{1/2}(Q^{2}) with Q2Q^{2}, and S1/2(Q2)0S_{1/2}(Q^{2})\sim 0 prediction. Both are incompatible with what we know today about the Roper resonance. Precision electroproduction data in the mass range above 2 GeV will be needed to test high mass states for their potential hybrid character, e.g. from experiments at CLAS12 Lanza and D’Angelo (2021).

Refer to caption

Figure 12: Quark transverse transition charge density corresponding to the pΔ+p\to\Delta^{+} transition. Light color indicates positive charge, dark color indicates negative charge. Top: p and Δ\Delta are unpolarized. Middle: p and Δ\Delta are polarized along bxb_{x} axis generating an electric dipole along the byb_{y} axis. Bottom: Quadrupole contribution to transition density. Graphics adapted from Carlson and Vanderhaeghen (2008).

.5 Conclusions and Outlook

In this contribution we have focused on more recent results of nucleon resonance transition amplitudes and their interpretation within LQCD and within most advanced approaches, e.g. in light front relativistic quark models and approaches with traceable links to first principle QCD such as Dyson-Schwinger Equations Roberts (2008) and light cone sum rules Braun et al. (2009). These calculations describe the transition form factors at Q22Q^{2}\geq 2 GeV2, while at lower Q2Q^{2} values hadronic degrees of freedom must be included and could even dominate contributions of the quark core.

For the lowest mass states, Δ(1232)32+\Delta(1232)\frac{3}{2}^{+} and the Roper N(1440)12+N(1440)\frac{1}{2}^{+}, LQCD calculations have been carried out that are consistent with the data within large uncertainties. These calculations are about one decade old, and new data, with higher precision in more extended kinematic range have been added to the database that warrant new Lattice calculations at the physical pion mass to be carried out.

Over the past decade, eight baryon states in the mass range from 1.85 to 2.15 GeV have been either discovered or evidence for the existence of states has been significantly strengthened. Some of these states are in the mass range and have JPCJ^{PC} quantum numbers that could have significant contributions of gluonic components. Such ”hybrid” states are in fact predicted in LQCD Dudek and Edwards (2012). These states appear with the same quantum numbers as ordinary quark excitations, and can only be isolated from ordinary states due to the Q2Q^{2} dependence of their helicity amplitudes Li et al. (1992), which is expected to be quite different from ordinary 3-quark excitation. The study of hybrid baryon excitations then requires new electroproduction data especially at low Q2Q^{2} Lanza and D’Angelo (2021) with different final states and with masses above 2 GeV.

Refer to caption

Figure 13: Projections of excited baryons with dominant gluonic components (marked in blue shades) in LQCD with 400MeV pions. The lowest hybrid baryon is projected with mass 1.3 GeV above the nucleon mass. The 1/2+1/2^{+} and 3/2+3/2^{+} states are clustered in a narrow mass range of about 200 MeV.

Despite the very significant progress made in recent years to further establish the light-quark baryon spectrum and explore the internal structure of excited states and the relationship to QCD Carman et al. (2020); Brodsky et al. (2020), much remains to be done. A vast amount of precision data already collected needs to be included in the multi-channel analysis frameworks, and polarization data are still to be analyzed. There are approved proposals to study resonance excitation at much higher Q2Q^{2} and with higher precision at Jefferson Lab with CLAS12 Tian et al. (2022); Burkert (2018a), which may begin to reveal the transition to the bare quark core contributions at short distances.

A new avenue of experimental research has recently been opened up with the first data-based extraction of a gravitational property of the proton, its internal pressure distribution, which is represented by the gravitational form factor Dq(t)D^{q}(t). It is one of the form factors of the QCD matrix element of the energy-momentum tensor, its internal pressure and shear stress distribution in space Burkert et al. (2018, 2021). These properties, as well as the distribution of mass and angular momentum, and are accessible directly in gravitational interaction, which is highly impractical. However the relevant gravitational form factor Dq(t)D^{q}(t) for the ground state nucleon can be accessed indirectly through the process of deeply virtual Compton scattering and in time-like Compton scattering Ji (1997); Chatagnon et al. (2021). Both processes, having a J=1J=1 photon in the initial state as well as in the final state, contain components of J=2J=2 that couple to the proton through a tensor interaction, as gravity does Polyakov and Schweitzer (2018).

Mechanical properties of resonance transitions have recently been explored for the N(1535)12N(938)N(1535)\frac{1}{2}^{-}\to N(938) gravitational transition form factors calculations in Özdem and Azizi (2020) and in Polyakov and Tandogan (2020). In order to access these novel gravitational transition form factors experimentally, the nucleon to resonance transition generalized parton distributions must be studied. The framework for studying the NN(1535)N\to N(1535) transition GPDs, which would enable experimental access to these mechanical properties, remains to be developed. The required effort is quite worthwhile as a new avenue of hadron physics has opened up that remains to be fully explored.

References