$S-P-D$ Mixing in Vector Quarkonia from the Salpeter Equation with Optimized Wave Function Representations

This paper proposes a novel mechanism based on the instantaneous Bethe-Salpeter (Salpeter) equation for investigating wave function mixing in vector mesons such as $ψ(3770)$. Conventional theories typically treat $ψ(3770)$ as a $2S-1D$ mixed state; however, considering only tensor forces or relativistic corrections alone often leads to mixing angles that are too small and inconsistent with experimental data. Phenomenological $2S-1D$ mixing requires experimental data as input to determine the mixing angles, resulting in limited theoretical studies on states like $Υ(1D, 2D)$ in the absence of experimental data. To more accurately describe $S-D$ mixing and its relativistic effects, this paper systematically compares eight possible relativistic wave function representations ($φ_1$ to $φ_8$) by solving the Salpeter equation and calculates the mass spectra and dileptonic decay widths of charmonium and bottomonium. The study finds that the wave function representation $φ_2$ can simultaneously reproduce the experimental data of both charmonium and bottomonium well. Further analysis reveals that, in addition to $S-D$ mixing, the wave functions of vector mesons contain a non-negligible $P$-wave component, meaning they are $S-P-D$ mixed states. We predict the mixing angles for bottomonium $Υ(1D)$ and $Υ(2D)$ to be $(1.78^{+0.32}{-0.25})^\circ$ and $(5.44^{+1.10}{-0.76})^\circ$, with dileptonic decay widths of $2.29^{+0.86}{-0.69}$ eV and $10.5^{+4.2}{-3.1}$ eV, respectively.

💡 Research Summary

The authors address a long‑standing discrepancy in the description of vector quarkonia, especially the ψ(3770) state, whose measured electronic width is far larger than that predicted for a pure D‑wave. Conventional non‑relativistic potential models treat ψ(3770) as a 2S–1D mixed state, attributing the mixing either to tensor forces or to coupled‑channel effects. However, such approaches either yield mixing angles that are too small (≈5–10°) or require experimental input to fix the angle, limiting predictive power for states without data, such as the bottomonium D‑wave states Υ(1D) and Υ(2D).

To overcome these limitations, the paper adopts the instantaneous Bethe–Salpeter (Salpeter) equation, which treats the bound‑state wave function relativistically while keeping the interaction kernel non‑relativistic (a modified Cornell potential). In the instantaneous approximation the most general 1⁻⁻ wave function contains eight independent Dirac structures. The authors construct eight distinct wave‑function ansätze (φ₁–φ₈), each incorporating different subsets of S‑wave, P‑wave, and D‑wave components. By inserting each ansatz into the Salpeter equation they obtain a set of coupled integral equations for the four independent radial functions f₃, f₄, f₅, and f₆. The constraints ϕ⁺⁻ = ϕ⁻⁺ = 0 reduce the number of independent equations to four, matching the number of unknown radial functions.

For each φₖ they solve the equations numerically, compute the mass eigenvalues of the low‑lying charmonium and bottomonium states, and evaluate the electronic decay constant F_V, which determines the dileptonic width Γ(V→e⁺e⁻). Comparison with experimental masses and widths shows that most ansätze reproduce the spectra within ≈10 MeV, but only φ₂ simultaneously yields the correct electronic widths for both S‑dominant (J/ψ, ψ(2S)) and D‑dominant (ψ(3770)) states. φ₂ has the structure

 ϕ₂ = (ε·q) f₁ + /q f₃ + /P/q f₄ + (M f₅ + /P f₆) /ε,

where the (ε·q) term carries a P‑wave contribution, the /q and /P/q terms contain both S‑ and D‑wave pieces, and the (M f₅ + /P f₆) /ε term represents the pure S‑wave component. The presence of both (ε·q) and /ε indicates that the physical vector meson is an S‑P‑D mixed state, not merely an S‑D mixture.

Applying the φ₂ framework to bottomonium, the authors fit the parameters (quark masses, string tension, Λ_QCD, etc.) to the known Υ(1S), Υ(2S), and Υ(3S) masses, then predict the masses of Υ(1D) and Υ(2D) (≈10 154 MeV and ≈10 454 MeV) in good agreement with existing measurements where available. More importantly, they predict the mixing angles for these D‑wave states: θ(Υ(1D)) = (1.78⁺⁰·³²₋₀·₂⁵)° and θ(Υ(2D)) = (5.44⁺¹·¹₀₋₀·₇₆)°, both significantly smaller than the angles often assumed in phenomenological S‑D mixing models. Corresponding dileptonic widths are Γ(Υ(1D)→e⁺e⁻) = 2.29⁺⁰·⁸⁶₋₀·⁶⁹ eV and Γ(Υ(2D)→e⁺e⁻) = 10.5⁺⁴·²₋₃·₁ eV, providing concrete targets for future high‑luminosity e⁺e⁻ colliders.

The paper thus demonstrates that a fully relativistic treatment of the bound‑state wave function, combined with a systematic exploration of its Dirac structure, resolves the ψ(3770) width puzzle and predicts previously inaccessible bottomonium properties. The identification of a non‑negligible P‑wave component (S‑P‑D mixing) is a novel insight, suggesting that relativistic corrections and tensor forces must be treated on an equal footing. The methodology offers a predictive framework for other heavy‑quark systems where experimental data are scarce, and it underscores the importance of choosing an appropriate wave‑function representation when solving the Salpeter equation.