Higher Radial & Orbital Excitations

• Higher radial and orbital excitations are quantum states with extra nodes or increased angular momentum, serving as benchmarks for nonperturbative QCD and atomic structure tests.
• Advanced methodologies like variational GEVP and structured light excitation enable precise analysis of these states in atomic, hadronic, and multiquark systems.
• Experimental observations in heavy-light meson spectra and exotic tetraquark decays critically validate theoretical models and drive progress in understanding complex quantum dynamics.

Higher radial and orbital excitations refer to quantum states of composite systems—such as atoms, hadrons, or multiquark states—in which the constituent particles exhibit additional nodes (radial excitations) or higher orbital angular momentum (orbital excitations) relative to the ground state. These excitations are manifest across atomic, hadronic, and multiquark spectra, and their characterization is pivotal for understanding the structure, spectra, and interaction mechanisms governing complex quantum systems. Experimental identification and theoretical modeling of such excitations provide stringent tests of nonperturbative QCD, atomic structure theory, and potential models.

1. Definitions and Spectroscopic Framework

Quantum states are classified by their radial quantum number $n$ and orbital angular momentum $l$ (or $L$ in hadronic/atomic notation). States with $n > 1$ but unchanged $l$ are termed higher radial excitations; those with $l > 0$ (i.e., $P$-wave, $D$-wave, etc.) are orbital excitations. In multi-body systems, further complexity arises due to possible excitations in relative coordinates ("$\rho$-mode," "$\lambda$-mode") and color-spin coupling, as relevant for tetraquarks and complex hadrons (Wang et al., 2021).

In the atomic context, Rydberg states refer to highly excited configurations with large $l$0; in hadronic physics, higher radial and orbital excitations are systematically mapped via quantum numbers $l$1 for mesons and baryons. For exotic multiquark systems, the assignment involves $l$2, internal color structure, and sometimes exotic quantum numbers (e.g., $l$3, $l$4) forbidden to simple $l$5 configurations.

2. Theory and Model Descriptions

Atomic Systems and Rydberg States

Atomic Rydberg states are accurately modeled using single-particle Schrödinger equations with effective potentials incorporating core polarization and quantum defect corrections. The effective potential of Marinescu et al. includes a nonlocal, core-polarized term and reproduces both inner-core electronic behavior and the asymptotic hydrogenic tail (Sanayei et al., 2015).

Radial and orbital excitations obey

$l$6

where $l$7 is the principal quantum number, $l$8 is the quantum defect, and $l$9 is the fine-structure correction. Uniform WKB-type approximations (Langer form near the outer turning point; Fock ansatz near $L$0) allow for analytic expressions for radial wavefunctions across the full range of $L$1, including for large $L$2 (Sanayei et al., 2015).

A subtle "anomaly" for $L$3 in Rb, Cs arises due to a tiny classically allowed region inside the core, causing three turning points in the effective potential and invalidating standard two-point WKB treatments. Numerical solutions restore agreement with experiment (Sanayei et al., 2015).

Hadronic Spectroscopy

Hadronic models, notably the relativized Godfrey–Isgur (GI) quark model with modified ("screened") linear confinement, predict mass spectra and decay properties of higher radial and orbital excitations in systems such as open-charm ($L$4) mesons (Song et al., 2015). The effective Hamiltonian comprises a central potential (Coulomb + screened linear), smeared spin–spin, and Breit–Fermi corrections. Radial excitations correspond to changes in $L$5, and orbital excitations increase $L$6.

For heavy-light mesons in the static limit, interpolating operators such as $L$7 project onto $L$8-wave ($L$9) configurations, while $n > 1$0, with a covariant derivative, access $n > 1$1-wave ($n > 1$2) and multi-hadron thresholds (Blossier et al., 2016). The Generalized Eigenvalue Problem (GEVP) approach is used to extract the spectrum from a matrix of two-point correlators with various Gaussian smearings and Dirac structures.

Multiquark States and Tetraquarks

Fully-charmed tetraquarks ($n > 1$3) exhibit a spectrum of S-wave radial excitations ($n > 1$4) and $n > 1$5-wave orbital excitations. The nonrelativistic four-quark model incorporates one-gluon-exchange (OGE) and scalar linear confinement, matched to the charmonium spectrum. In this framework, two types of $n > 1$6 orbital excitations are distinguished: $n > 1$7-mode (inside one diquark or antidiquark cluster) and $n > 1$8-mode (between clusters), each with possible color configurations $n > 1$9 or $l$0 (Wang et al., 2021).

An essential finding is that the $l$1 configuration is typically more deeply bound than $l$2 for the lowest states, due to stronger inter-cluster attraction.

3. Experimental Spectra and Assignments

Charmed Meson Family

The modified GI model reproduces the observed ordering and approximate mass values of the higher $l$3 meson excitations. For example, $l$4 is identified as $l$5, and $l$6 as mixed $l$7–$l$8 (mostly $l$9) (Song et al., 2015). Orbital excitations such as $l > 0$0 correspond to $l > 0$1 assignments, with mixing angles reflecting heavy-quark symmetry.

Predicted masses for yet-unobserved higher excitations (2P, 2D, 1F, etc.) cluster in the 3.05–3.15 GeV region, with widths and dominant decay modes calculable in the $l > 0$2 quark pair creation model.

Selected Higher Excitations in the $l > 0$3-meson Family

State	Theoretical $l > 0$4 (GeV)	Assignment/Comment
$l > 0$5	2.534	$l > 0$6
$l > 0$7	2.593	$l > 0$8
$l > 0$9	2.779	$P$0
$P$1	2.976	$P$2
$P$3	3.015	$P$4
$P$5	3.053	Not yet observed

Tetraquarks and Exotic Signals

The $P$6 spectrum calculated with the parameters fixed to reproduce charmonium benchmarks results in S-wave ground states at 6405–6547 MeV depending on $P$7, with first radial excitations at 6867–6917 MeV. Low-lying $P$8-wave states ($P$9-modes in $D$0) appear at 6589–6608 MeV, with the mass gap between $D$1 and $D$2 notably smaller than between $D$3 and $D$4 states (Wang et al., 2021).

The observed $D$5 state in $D$6 at LHCb coincides with these predictions, suggesting a possible identification as either a $D$7 ($D$8, $D$9) or $\rho$0-wave ($\rho$1, $\rho$2) tetraquark (Wang et al., 2021).

4. Methodologies for Excitation Extraction

Lattice QCD Techniques

In heavy-light mesons, radial and orbital excitations are isolated via matrices of two-point correlators—constructed from interpolators with different smearings and Dirac structures—and the variational GEVP method. The use of derivative-based interpolators, such as $\rho$3, couples strongly to $\rho$4-wave and two-hadron states, making the inclusion or exclusion of these operators a diagnostic for the nature of extracted excitations (Blossier et al., 2016). Charge-density and axial-density spatial distributions further serve as diagnostics; nodes in axial densities confirm the one-particle nature of the radial excitation, while drifting integrated vector charge signals multihadron contamination.

Numerical extractions yield, for the $\rho$5 system, a first radial–ground splitting of approximately $\rho$6 MeV and an orbital–ground splitting of $\rho$7 GeV (E5 ensemble, $\rho$8 fm).

Atomic Excitation with Structured Light

Laguerre–Gauss beams, which carry photon orbital angular momentum (OAM) $\rho$9, can induce atomic transitions to higher-$\lambda$0 Rydberg states inaccessible by conventional dipole selection rules (Rodrigues et al., 2015). The selection rules are generalized: $\lambda$1 allowing, for example, excitation from $\lambda$2 to $\lambda$3 and beyond in a single-photon process. The radial matrix elements scale as $\lambda$4 with $\lambda$5, $\lambda$6, etc., such that for $\lambda$7, $\lambda$8 in $\lambda$9 is $l$00, respectively. Rates for such transitions are substantial (5–20% of usual dipole rates for $l$01) and can be fully calculated with the tabulated coupling coefficients (Rodrigues et al., 2015).

5. Decay Properties and Phenomenological Implications

Higher radial and orbital excitations frequently decay via OZI-allowed strong processes. The $l$02 (quark pair creation, QPC) model gives partial and total widths, predicting, for instance, the main decay channels and branching ratios of $l$03 to be consistent with experiment, although total widths tend to be underpredicted—a known issue of the model (Song et al., 2015). For Rydberg states in atoms, properties such as the Fermi-contact hyperfine splitting constants scale precisely with excitation number: $l$04 (Sanayei et al., 2015).

Exotic multiquarks—such as tetraquark $l$05 or $l$06 states—are predicted to have distinctive decay modes ($l$07 for $l$08) and mass gaps (P-wave–S-wave) that facilitate identification through channels inaccessible to ordinary mesons (Wang et al., 2021).

6. Challenges and Future Prospects

Extraction of higher excitations is complicated by proximity to multi-hadron thresholds and the potential for operator bases to have suppressed overlap with certain states, leading to "missing" levels in variational analyses (Blossier et al., 2016). Optimal disentanglement requires a judicious basis of interpolators, thorough analysis of density distributions, and, where relevant, the inclusion of explicit two-hadron operators.

For heavy-quark and tetraquark systems, future studies will benefit from lighter pion mass ensembles in lattice QCD, explicit inclusion of scattering operators, and systematic continuum/chiral extrapolation of energy splittings. In atomic physics, further exploration of OAM-carrying photon interactions with complex atoms or ions may extend the reach of high-$l$09 excitations beyond current limits.

A plausible implication is that identification of narrow, parity-exotic states in the predicted mass bands by future experiments will provide critical evidence for the nature and dynamics of higher radial and orbital excitations, both in conventional and exotic hadrons.