Light-Cone Harmonic Oscillator Model

• The Light-Cone Harmonic Oscillator Model is a light-front framework that employs harmonic oscillator wavefunctions to represent meson longitudinal and transverse momentum distributions.
• It integrates spin structures via Melosh–Wigner rotations and couples a Gaussian transverse profile with parameterized longitudinal corrections for precise endpoint control.
• Recent improvements calibrate the model using QCD sum-rule moments and lattice inputs, enhancing its predictions for exclusive heavy-to-light transition form factors.

Searching arXiv for recent and foundational papers on the light-cone harmonic oscillator model and closely related formulations. First, I’ll look for the 2024 paper on the improved LCHO model for the $\phi$-meson LCDA. Searching arXiv for: (Hu et al., 2024) light-cone harmonic oscillator phi meson The Light-Cone Harmonic Oscillator Model designates a class of constructions in which harmonic-oscillator wavefunctions are formulated in light-cone or light-front variables. In hadron physics, it is used to build meson light-cone wavefunctions that encode both longitudinal momentum-fraction dependence and transverse-momentum structure, typically within the Brodsky–Huang–Lepage prescription and with Melosh–Wigner spin rotations, after which leading-twist light-cone distribution amplitudes are obtained by integrating over transverse momenta up to a factorization scale. In a distinct relativistic usage, the same label has also been attached to a Lorentz-invariant oscillator derived from a quaternionic non-relativistic oscillator and written in terms of the light-cone interval. These usages share harmonic-oscillator language but correspond to different theoretical objects and applications (Hu et al., 2024, Zhong et al., 2021, Arbab, 2017).

1. Light-cone wavefunction framework in QCD

In the QCD literature represented by the improved $\phi$-meson and pion studies, the Light-Cone Harmonic Oscillator Model is a light-cone wavefunction ansatz for valence hadronic structure. The central object is a meson light-cone wavefunction, written as a product of a spin wavefunction and a momentum-space wavefunction. For the longitudinally polarized $\phi$ meson,

$$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$

with a longitudinal spin factor derived from Melosh–Wigner rotation and a spatial part of BHL harmonic-oscillator form,

$$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$

The model therefore couples a Gaussian transverse profile to an $x$-dependent width through the invariant mass and explicitly retains constituent $s$-quark mass effects (Hu et al., 2024).

The pion formulation has the same structural logic. Its full light-cone wavefunction is written as

$$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$

where the spin-space wavefunction is again generated by Wigner–Melosh rotation and the spatial part is

$$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$

In both cases the LCDA is obtained by integrating the LCWF over transverse momentum up to $\mu$, so the LCHO model is not merely a transverse ansatz; it is a combined model for spin, longitudinal profile, and confinement-motivated $\phi$0 dependence (Zhong et al., 2021).

Normalization conditions are an essential part of the construction. For the $\phi$1 meson, the wavefunction is constrained by

$\phi$2

with $\phi$3. For the pion, analogous constraints come from $\phi$4 and $\phi$5 (Hu et al., 2024, Zhong et al., 2021).

2. Longitudinal profile, endpoint control, and model improvement

A defining development of the improved LCHO literature is the replacement of overly rigid longitudinal profiles by compact parameterizations that control endpoint behavior and shape. In the $\phi$6-meson study, two improved forms are introduced: $\phi$7 and

$\phi$8

Model I is described as Gegenbauer-like, whereas Model II uses a multiplicative power of $\phi$9 to control end-point behavior and better accommodate evolution toward the asymptotic $\phi$0 form, with a small $\phi$1 modulation refining the shape (Hu et al., 2024).

The pion study develops the same idea in a slightly different notation. Earlier choices employed simple Gegenbauer-inspired forms, but the improved proposal introduces

$\phi$2

These forms are symmetric under $\phi$3, positive on $\phi$4, and intended to tune endpoint behavior and overall width without adding many Gegenbauer terms. The stated motivation is to avoid relying on poorly determined higher Gegenbauer moments while still using higher moment information (Zhong et al., 2021).

For the $\phi$5-meson transverse twist-2 LCDA, the longitudinal correction is implemented more minimally through

$\phi$6

while the closed-form LCDA contains a difference of error functions. The paper emphasizes two consequences: approach to the asymptotic shape $\phi$7 as $\phi$8, and exponential endpoint suppression as $\phi$9, which mitigates endpoint singularities in LCSR convolutions (Yang et al., 1 Nov 2025).

A recurrent misconception is that the LCHO model is only a Gaussian in $$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$0. The improved literature shows that the decisive phenomenological freedom lies equally in the longitudinal profile $$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$1, especially in the treatment of endpoint suppression, central peaking, and compatibility with ERBL evolution (Hu et al., 2024, Zhong et al., 2021).

3. Moment constraints, sum rules, and ERBL evolution

The improved LCHO program fixes model parameters by matching moments of the LCDA rather than by ad hoc shape selection. The standard moments are

$$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$2

with even $$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$3 for the isospin-symmetric channels discussed here. In the $$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$4-meson analysis, the first ten $$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$5-moments are derived from QCD sum rules under background field theory. The correlator is

$$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$6

with interpolating currents

$$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$7

The OPE is performed in deep Euclidean space and includes perturbative and vacuum-condensate contributions up to dimension six, with the $$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$8 pole isolated and the continuum modeled by quark–hadron duality at $$\Psi_{2;\phi}^{\parallel}(x,\mathbf{k}_\perp)=\chi_{2;\phi}(x,\mathbf{k}_\perp)\,\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp),$$9, followed by a Borel transform (Hu et al., 2024).

The numerical $$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$0 moments at $$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$1 are

$$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$2

$$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$3

A stabilized Borel window is chosen by capping continuum contributions for $$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$4 at $$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$5, $$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$6, $$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$7, $$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$8, $$\psi_{2;\phi}^{\rm R}(x,\mathbf{k}_\perp)=A_{2;\phi}^{\parallel}\,\varphi_{2;\phi}^{\parallel}(x)\, \exp\!\left[-\frac{\mathbf{k}_\perp^2+m_s^2}{8\,\beta_{2;\phi}^2\,x(1-x)}\right].$$9, respectively, and dimension-six terms below $x$0. Matching these moments yields, at $x$1, Model I parameters

$x$2

and Model II parameters

$x$3

At the LCSR working scale $x$4, Model II becomes

$x$5

with goodness-of-fit $x$6 (Hu et al., 2024).

The pion paper follows the same philosophy but highlights a specific normalization problem: the sum rule for the $x$7 moment $x$8 cannot be normalized. It therefore proposes a “more reasonable” extraction,

$x$9

to reduce correlated systematics. The resulting moments at $s$0 are

$s$1

$s$2

with corresponding Gegenbauer moments

$s$3

For the improved LCHO-IV fit with $s$4,

$s$5

with $s$6 (Zhong et al., 2021).

Scale dependence is treated through ERBL evolution. For the longitudinal vector case,

$s$7

with

$s$8

The $s$9 transverse study uses the analogous chiral-odd evolution,

$$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$0

and fixes its LCHO parameters from normalization, the average transverse momentum $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$1, and the Dimou–Lyon–Zwicky input $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$2 (Yang et al., 1 Nov 2025).

4. The $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$3-meson implementation and the single-peak LCDA

The 2024 improved $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$4-meson study is a paradigmatic realization of the LCHO program. It starts from a longitudinally polarized $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$5 LCWF with explicit Melosh–Wigner spin structure for a vector meson, $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$6-quark mass effects in both spin and spatial parts, and a refined $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$7-dependence with controlled endpoint behavior. The LCDA is reconstructed by integrating the LCWF over transverse momentum up to the factorization scale. The resulting $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$8 tends to be a single-peak behavior, which is stated to be consistent with the latest Lattice QCD result, with the earlier QCD SR double-hump LCDA disfavored. The paper also notes consistency with lattice-QCD in the large-momentum effective theory framework and with DSE/Bethe–Salpeter models (Hu et al., 2024).

The same LCDA is then inserted into a light-cone sum-rule analysis of $$\Psi_{2;\pi}(x,k_\perp)=\sum_{\lambda_1,\lambda_2}\chi^{\lambda_1\lambda_2}_{2;\pi}(x,k_\perp)\,\Psi^R_{2;\pi}(x,k_\perp),$$9 transition form factors. Using the correlation function with current $$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$0 and pseudoscalar interpolator $$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$1, the correlator is decomposed into invariant amplitudes, matched between the OPE expanded near $$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$2 up to twist-4 LCDAs and the hadronic representation with the $$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$3 pole and continuum. The adopted numerical setup is

$$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$4

with Borel windows

$$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$5

At the large recoil point,

$$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$6

and

$$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$7

The paper states that these ratios are in good agreement with BESIII, PDG and LQCD within uncertainties (Hu et al., 2024).

For the full kinematic region, the form factors are extrapolated with the simplified series expansion

$$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$8

with residuals

$$\Psi^R_{2;\pi}(x,k_\perp)=A_{2;\pi}\,\phi_{2;\pi}(x)\, \exp\!\left[-\frac{k_\perp^2+m_q^2}{8\beta_{2;\pi}^2x(1-x)}\right].$$9

The integrated semileptonic widths are

$\mu$0

$\mu$1

leading to

$\mu$2

These branching fractions are described as showing good agreement with BESIII, CLEO, and BABAR, and the paper further quotes integrated asymmetry and polarization observables such as

$\mu$3

Within this application, the LCHO model functions as a parameter-controlled nonperturbative input for exclusive heavy-to-light phenomenology (Hu et al., 2024).

5. Generalizations to the pion and the $\mu$4 meson

The improved pion LCHO model shows how the same framework extends beyond vector mesons and beyond a pure Gegenbauer truncation. After fitting the BFTSR moments up to $\mu$5 order, the resulting twist-2 pion DA is described as “more closely to the AdS/QCD and lattice result, but is narrower than that by Dyson–Schwinger equation.” At $\mu$6 the paper quotes

$\mu$7

and the inverse moment

$\mu$8

The improved DA is then applied to the pion–photon transition form factor and to the $\mu$9 form factor, with the prediction

$\phi$00

at $\phi$01, stated to be consistent with other theory (Zhong et al., 2021).

The 2025 $\phi$02-meson study extends the model to transverse twist-2 vector LCDAs and illustrates a different but related strategy: using a right-handed chiral current so that the transverse twist-2 LCDA $\phi$03 dominates the LCSR for $\phi$04. The LCWF is written in BHL form with constituent light-quark mass $\phi$05 and Gaussian transverse dependence,

$\phi$06

supplemented by the longitudinal correction $\phi$07. The central parameter sets are

$\phi$08

at $\phi$09, and

$\phi$10

at $\phi$11 (Yang et al., 1 Nov 2025).

The resulting $\phi$12 form factors at $\phi$13 are

$\phi$14

with

$\phi$15

After simplified $\phi$16-series extrapolation, the paper quotes

$\phi$17

and branching fractions

$\phi$18

in good agreement with BESIII and CLEO. These results show that the LCHO framework is not confined to a single channel or polarization sector; it has been adapted to longitudinal and transverse twist-2 amplitudes and to both light and heavy-to-light observables (Yang et al., 1 Nov 2025).

6. The relativistic oscillator formulation and the light-cone interval

A separate line of work uses “Light-Cone Harmonic Oscillator Model” in a different sense. In “On relativistic harmonic oscillator,” a relativistic quantum harmonic oscillator in $\phi$19 dimensions is derived from a quaternionic non-relativistic harmonic oscillator. The quaternionic equation yields the Klein–Gordon wave equation with a covariant mass,

$\phi$20

where

$\phi$21

so the oscillator potential is

$\phi$22

The dependence is through the Lorentz scalar $\phi$23, and the paper presents Lorentz-invariant solutions with Gaussian envelopes in that variable (Arbab, 2017).

After separation of variables, the total energy spectrum is

$\phi$24

rather than the non-relativistic $\phi$25D value $\phi$26. The stated reason is that the time coordinate contributes

$\phi$27

to the total energy, while each spatial coordinate contributes $\phi$28. The model also derives a quantized effective mass spectrum,

$\phi$29

and imposes the condition $\phi$30 for physical states (Arbab, 2017).

Its light-cone content is formulated through the “quantum interval”

$\phi$31

with null, timelike, and spacelike regimes determined by the refractive index $\phi$32. In $\phi$33 dimensions, with

$\phi$34

this becomes

$\phi$35

This formulation is conceptually distinct from the hadronic BHL/Melosh–Wigner LCHO model. The former is a Lorentz-invariant oscillator with a covariant mass and medium-dependent interval, whereas the latter is a phenomenological light-cone wavefunction model for hadron structure (Arbab, 2017).

7. Limitations, uncertainties, and current interpretive status

Across the hadronic literature, the LCHO model is explicitly uncertainty-limited by its nonperturbative inputs and by the rigidity of the harmonic-oscillator ansatz. The $\phi$36-meson study identifies dominant sources from the QCD sum-rule inputs for $\phi$37-moments, including $\phi$38, $\phi$39, $\phi$40, $\phi$41, four-quark condensates, their RG evolution, the continuum threshold $\phi$42, and Borel-window choices. On the LCSR side, the uncertainties arise from $\phi$43, $\phi$44, $\phi$45, $\phi$46, the threshold $\phi$47, Borel parameters, and twist-3/4 $\phi$48 LCDAs. The paper also notes that changing the constituent $\phi$49 from $\phi$50 to $\phi$51 alters shapes mildly, and that these uncertainties propagate to the LCDA peak sharpness, endpoint suppression, and especially to $\phi$52 and $\phi$53 (Hu et al., 2024).

The pion analysis makes analogous points. Its BFTSR systematics are controlled by choosing Borel windows with dimension-6 contributions below $\phi$54 and continuum contributions at $\phi$55–$\phi$56 depending on $\phi$57, but the model remains sensitive to the constituent mass $\phi$58 and the oscillator scale $\phi$59. As $\phi$60 increases from $\phi$61 to $\phi$62, $\phi$63 decreases from about $\phi$64 to about $\phi$65 and the DA becomes flatter; the preferred range is stated as $\phi$66. The study also reports that including a Jacobi factor in the spatial wavefunction has small numerical impact on $\phi$67 (Zhong et al., 2021).

The $\phi$68 transverse analysis identifies its main uncertainties as the LCDA parameters $\phi$69 through $\phi$70 and $\phi$71, the decay constant $\phi$72, the charm mass $\phi$73, the sum-rule choices for $\phi$74 and $\phi$75, the factorization scale, LO ERBL evolution, and the truncation of the Gegenbauer series to the first non-trivial term. The quoted uncertainties on $\phi$76 and $\phi$77 are at the $\phi$78–$\phi$79 level and the branching-fraction errors are about $\phi$80–$\phi$81 (Yang et al., 1 Nov 2025).

The principal limitation shared by the hadronic LCHO constructions is the Gaussian transverse profile itself. The $\phi$82-meson paper explicitly notes that more flexible transverse-momentum dependence, such as power-law tails, could be explored. It also notes that higher-twist contributions are included only up to twist-4, and that the $\phi$83-moments are computed with truncations at dimension-six condensates. At the same time, the same paper states that the framework and parameter-fixing strategy can be extended to other vector mesons such as $\phi$84, $\phi$85, and $\phi$86, and that transverse $\phi$87 LCDAs can be treated analogously. This suggests that the Light-Cone Harmonic Oscillator Model is best understood not as a single fixed formula but as a family of light-cone ansätze whose defining features are harmonic transverse structure, explicit spin treatment, and parameter determination from moment constraints and exclusive-process phenomenology (Hu et al., 2024).