Screened Potential Model Overview

• Screened Potential Model is an effective framework that replaces bare long-range interactions with forms that weaken via exponential, saturating, or threshold-dependent behavior.
• It is applied across fields such as atomic physics, heavy-quark spectroscopy, molecular dynamics, and plasma physics to simulate medium effects like vacuum polarization and string breaking.
• Implementations reveal varied spectral consequences, where some models preserve discrete bound states while others exhibit critical parameters beyond which states vanish.

A screened potential model is an effective interaction framework in which an unscreened long-range potential is replaced by a screened form that weakens with distance, saturates at large separation, or changes across threshold regions. In the literature surveyed here, screened potentials appear as Yukawa-type factors multiplying Coulombic terms, saturating confinement terms of the form $\lambda(1-e^{-\mu r})/\mu$, threshold-dependent piecewise Cornell-like interactions, and Yukawa-screened exchange kernels in electronic-structure theory (Osobonye et al., 2020, Bokade et al., 6 Jan 2025, Gonzalez, 2014, Tran et al., 2011). The screening parameter is correspondingly model dependent: it may control the range suppression of an inverse-distance force, the flattening scale of confinement, the onset of threshold-induced saturation, or the decay length of an exchange kernel.

1. Formal definitions and recurring parameterizations

Several mathematically distinct constructions are all described as screened potential models because each replaces an unscreened interaction by one with reduced long-range strength. In atomic and molecular Schrödinger problems, the screened interaction is often a Coulomb term multiplied by an exponential factor, as in the Yukawa potential $-e^{-\delta r}/r$, the exponential cosine screened Coulomb potential $-e^{-\delta r}\cos(\delta r)/r$, or generalized screened Kratzer-Hellmann forms (Roy, 2019, Osobonye et al., 2020). In heavy-quark spectroscopy, the confining part is frequently screened through a saturating term such as

$$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$

or through threshold-dependent flattening of a Cornell potential (Bokade et al., 6 Jan 2025, Gonzalez, 2014). In solid-state hybrid DFT, the bare Coulomb kernel is replaced by the Yukawa interaction

$$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$

thereby removing the long-range Hartree-Fock tail (Tran et al., 2011).

Context	Representative screened form	Screening role
Atomic or molecular bound states	$-e^{-\delta r}/r$	exponential range suppression
Heavy-quark confinement	$\lambda(1-e^{-\mu r})/\mu$	linear at short distance, saturating at large distance
Threshold-based quark models	piecewise Cornell-like potential with constants beyond crossing radii	threshold-induced flattening
Screened exchange in solids	$$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$	suppression of long-range HF exchange

Two structural motifs recur. First, short-distance behavior is often retained: screened confinement models reduce to linear confinement for $r\ll 1/\mu$, and Yukawa-type kernels reduce to Coulombic behavior at small separation. Second, large-distance behavior is modified: the interaction either decays exponentially, saturates to a constant, or is replaced by threshold plateaus (0903.5506, Hao et al., 2022, Tran et al., 2011). This suggests that “screening” is not a single mathematical prescription but a family of effective reductions of long-range interaction strength.

2. Spectral consequences in screened Coulomb problems

The spectral effect of screening depends on the asymptotic structure of the chosen model. A notable case is the screened Coulomb potential

$V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$

for which contradictory claims had appeared regarding the disappearance of bound states at critical screening values. The 2023 analysis of this problem concludes that such critical values $-e^{-\delta r}/r$0 are incorrect: the long-range tail remains attractive, with $-e^{-\delta r}/r$1, and the effective potential does not become repulsive for any value of $-e^{-\delta r}/r$2 and $-e^{-\delta r}/r$3 (Fernández, 2023). The Hellmann-Feynman theorem gives

$-e^{-\delta r}/r$4

so energies move upward with increasing $-e^{-\delta r}/r$5, but this monotonic increase does not imply loss of the discrete spectrum. The same work derives a small-$-e^{-\delta r}/r$6 approximation,

$-e^{-\delta r}/r$7

and a large-$-e^{-\delta r}/r$8 asymptotic law for $-e^{-\delta r}/r$9 states,

$-e^{-\delta r}\cos(\delta r)/r$0

consistent with the pseudospectral calculations of Xu et al. and explicitly rejecting the interpretation advanced by Stachura and Hancock (Fernández, 2023).

By contrast, other screened Coulomb families do exhibit genuine critical screening. For the Hulthén, Yukawa, and exponential cosine screened Coulomb potentials, the critical parameter $-e^{-\delta r}\cos(\delta r)/r$1 is defined by $-e^{-\delta r}\cos(\delta r)/r$2, and states cease to be bound for $-e^{-\delta r}\cos(\delta r)/r$3 (Roy, 2019). Using the generalized pseudospectral method, accurate $-e^{-\delta r}\cos(\delta r)/r$4 values were reported for all hydrogenic states with $-e^{-\delta r}\cos(\delta r)/r$5. Representative results include the exact Hulthén $-e^{-\delta r}\cos(\delta r)/r$6-state threshold $-e^{-\delta r}\cos(\delta r)/r$7, the Yukawa value $-e^{-\delta r}\cos(\delta r)/r$8, and the exponential cosine screened Coulomb value $-e^{-\delta r}\cos(\delta r)/r$9 (Roy, 2019).

These results establish an important distinction. In one class of screened Coulomb models, screening changes level positions without terminating the bound spectrum; in another, screening drives states to threshold and defines finite critical parameters. The difference is controlled by the detailed long-range structure of the screened interaction rather than by the generic presence of a screening parameter.

3. Heavy-quark, meson, and baryon spectroscopy

In hadron spectroscopy, screened potential models are used primarily to encode unquenching, string breaking, or meson-meson threshold effects within an effective valence description. One influential construction is the Generalized Screened Potential Model, developed for bottomonium and later applied to charmonium. It starts from a Cornell potential,

$$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$0

but replaces the single global interaction by threshold-dependent potential branches. Below threshold the potential is Cornell-like; at and beyond threshold-crossing radii it flattens to constants determined by open-flavor meson-meson masses. In this formulation, screening is explicitly tied to meson-meson configurations such as $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$1, but the model is presented as an effective spectral model rather than an explicit coupled-channel calculation (Gonzalez, 2014). For bottomonium, this produces a denser spectrum above the first open-flavor threshold: in the $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$2 sector between $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$3 and $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$4 MeV the GSPM yields three states, $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$5, $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$6, and $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$7 MeV, whereas the Cornell model yields only one state at $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$8 MeV (Gonzalez, 2014). In charmonium, the same threshold mechanism was used to assign or predict extra $$V_S(r)=\lambda\left(\frac{1-e^{-\mu r}}{\mu}\right)+V_0,$$9 states, including $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$0 as a $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$1 threshold-induced state and additional states near $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$2 and $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$3 MeV (González, 2015).

A second major line replaces the linear confining term by a saturating screened interaction. In charmonium and bottomonium this is typically written as

$$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$4

which behaves linearly at small $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$5 and flattens at large $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$6 (Bokade et al., 6 Jan 2025). The same mechanism appears in nonrelativistic and relativistic screened potential models for charmonium, bottomonium, bottom mesons, $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$7 mesons, and beauty baryons (M et al., 10 Mar 2025, Bokade et al., 2024, Feng et al., 2022, Hao et al., 2022, Kaushal et al., 20 May 2025). A consistent phenomenological conclusion is that screening lowers higher excited masses relative to unscreened linear potentials. In the 2009 screened charmonium analysis, for example, the $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$8, $$\frac{e^{-\lambda|\mathbf r-\mathbf r'|}}{|\mathbf r-\mathbf r'|},$$9, $-e^{-\delta r}/r$0, and $-e^{-\delta r}/r$1 levels were all substantially lower than in the unscreened comparison model, and this compression was used to motivate assignments such as $-e^{-\delta r}/r$2, $-e^{-\delta r}/r$3, and $-e^{-\delta r}/r$4 (0903.5506).

Relativistic screened models extend this logic by combining screened confinement with spinless Salpeter kinematics, perturbative spin-dependent forces, and explicit $-e^{-\delta r}/r$5-$-e^{-\delta r}/r$6 mixing analyses. In bottomonium, a relativistic screened potential model interprets $-e^{-\delta r}/r$7 as a $-e^{-\delta r}/r$8-$-e^{-\delta r}/r$9 mixed state, $\lambda(1-e^{-\mu r})/\mu$0 as $\lambda(1-e^{-\mu r})/\mu$1-$\lambda(1-e^{-\mu r})/\mu$2, $\lambda(1-e^{-\mu r})/\mu$3 as a pure $\lambda(1-e^{-\mu r})/\mu$4, and $\lambda(1-e^{-\mu r})/\mu$5 and $\lambda(1-e^{-\mu r})/\mu$6 as $\lambda(1-e^{-\mu r})/\mu$7-$\lambda(1-e^{-\mu r})/\mu$8 mixed states (Bokade et al., 6 Jan 2025). Related nonrelativistic work with screened confinement, $\lambda(1-e^{-\mu r})/\mu$9 corrections, and the Matrix-Numerov method treats bottomonium, $$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$0 diquarks, and $$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$1 baryons within a unified potential framework (Kaushal et al., 20 May 2025). In light-baryon spectroscopy, a screened hypercentral constituent quark model uses

$$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$2

finds that screening is noticeable in light systems at higher mass scale, and reports that hyperfine splitting is smaller than in the linear-potential case (Menapara et al., 2023).

For open-bottom mesons, screened nonrelativistic quark models combined with the $$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$3 decay model have been used to assign $$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$4, $$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$5, $$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$6, and $$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$7, and likewise a series of recently observed $$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$8 structures, by requiring simultaneous consistency of masses and strong widths (Feng et al., 2022, Hao et al., 2022). Across these applications, screening functions as an effective proxy for vacuum polarization, pair creation, or open-channel coupling, while preserving the practical solvability of a two-body or reduced few-body bound-state problem.

4. Molecular, thermodynamic, and time-dependent screened models

In molecular and related one-body quantum problems, screened potential models often combine inverse-distance terms with exponential attenuation and are solved approximately by algebraic reduction methods. One example is the screened Kratzer potential

$$e^{-\lambda|\mathbf r-\mathbf r'|}/|\mathbf r-\mathbf r'|$$9

used for diatomic molecules, where $r\ll 1/\mu$0 is the screening parameter and the radial Schrödinger equation is treated with the Greene-Aldrich approximation and factorization methods (Amadi et al., 2020). The resulting energy spectrum is then used to build partition functions and thermodynamic observables for HCl, LiH, and H$r\ll 1/\mu$1 under modified Dirac delta and uniform superstatistical distributions. The reported trends include decreasing partition function with increasing $r\ll 1/\mu$2, entropy changes that depend on the deformation parameters, and increasing mean energy and heat capacity with $r\ll 1/\mu$3 (Amadi et al., 2020).

A more general construction is the screened Kratzer-Hellmann potential,

$r\ll 1/\mu$4

which reduces to Hellmann, screened Kratzer, Yukawa, and Coulomb limits under appropriate parameter choices (Osobonye et al., 2020). Using the Nikiforov-Uvarov method, this model yields approximate bound-state energies, normalized eigenfunctions in terms of Jacobi polynomials, rotational-vibrational partition functions, and information-theoretic quantities. The numerical analysis reports that increasing $r\ll 1/\mu$5 raises the energy eigenvalues, that Shannon entropy results satisfy the Bialynicki-Birula–Mycielski inequality, and that the Fisher information results satisfy the stated Stam/Cramér-Rao-type bound (Osobonye et al., 2020).

Screened Kratzer structures have also been used for heavy quarkonia and for the time-like analogue of quantum mechanics. A series-expansion solution of a Kratzer plus screened Coulomb model was applied to charmonium and bottomonium mass spectra, with special limits reducing to pure Kratzer, Coulomb, and screened Coulomb cases (Inyang et al., 2021). In the Feinberg-Horodecki framework, a time-dependent screened Kratzer-Hellmann potential generates quantized momentum eigenvalues rather than energy eigenvalues; special cases again include the screened Kratzer, Hellmann, screened Coulomb, and Coulomb potentials (Farout et al., 2020). A plausible implication is that screened-potential methodology is valued not only for reproducing spectroscopy but also for preserving analytically tractable deformations of standard solvable problems.

5. Plasmas, white dwarfs, and screened exchange in solids

In plasma physics, a screened potential model describes the effective interaction generated by the dielectric response of an electron gas around a test ion. For a classical ion in a weakly correlated quantum plasma, the benchmark screened interaction is obtained from the static RPA dielectric function, while collisions can be incorporated through the Mermin dielectric function (Moldabekov et al., 2015). Within this setting, the Yukawa potential appears as the lowest-order long-wavelength approximation, but the principal comparison in the literature is among linearized quantum-hydrodynamic screened potentials. The Shukla-Eliasson potential is found to be qualitatively different from the full RPA result and to predict an attractive minimum absent from the benchmark potential, whereas the Stanton-Murillo potential at any temperature and the Akbari-Moghanjoughi potential at $r\ll 1/\mu$6 are significantly more accurate (Moldabekov et al., 2015). The same analysis identifies the correct fermionic Bohm prefactor as $r\ll 1/\mu$7 at $r\ll 1/\mu$8.

In white-dwarf asteroseismology, the screened Coulomb potential is not used as a macroscopic stellar potential but as a microscopic ingredient in diffusion coefficients. For the DBV star PG 0112+104, WDEC model grids built from MESA-derived starting models were evolved with pure-Coulomb and screened-Coulomb diffusion physics; the screened version changed the C/O-He interface and Brunt-Väisälä frequency near $r\ll 1/\mu$9 and reduced the best-fit $V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$0 from $V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$1 s to $V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$2 s, a $V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$3 improvement (Chen, 2017). For the DAV star HS 0507+0434B, replacing the pure Coulomb interaction by the screened form

$V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$4

reduced $V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$5 from $V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$6 s to $V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$7 s, a $V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$8 improvement, and was interpreted as altering the C/He and He/H transition zones relevant for mode trapping (Chen, 2020).

In all-electron electronic-structure theory, screened potentials enter through exchange rather than through an external binding potential. The WIEN2k implementation of screened hybrid functionals replaces the Hartree-Fock kernel by the Yukawa interaction

$V(r)=-r^{-1}e^{-C/r}, \qquad C>0,$9

leading to the YS-PBE0 functional (Tran et al., 2011). Screening removes the long-range HF singularity, improves numerical convergence, and yields transition energies and structural properties intermediate between PBE and unscreened PBE0. In Cu$-e^{-\delta r}/r$00O, YS-PBE0 gives a gap of $-e^{-\delta r}/r$01 eV and an electric-field gradient of $-e^{-\delta r}/r$02 V/m$-e^{-\delta r}/r$03, substantially improving over semilocal approximations (Tran et al., 2011).

6. Conceptual distinctions, controversies, and limitations

The expression “screened potential model” therefore denotes a methodological family rather than a unique theory. In some contexts screening means exponential suppression of a Coulomb tail; in others it means flattening of confinement through string breaking, threshold-dependent piecewise potentials, screened microscopic diffusion, or removal of long-range exact exchange (Roy, 2019, Gonzalez, 2014, Chen, 2017, Tran et al., 2011). A common misconception is that screening always implies either a repulsive core or the eventual disappearance of all bound states. The screened Coulomb analysis of $-e^{-\delta r}/r$04 shows that this is false for that model, whereas the hydrogenic Yukawa, Hulthén, and exponential cosine screened Coulomb systems do possess genuine critical screening parameters (Fernández, 2023, Roy, 2019).

A second recurring issue concerns what is being modeled effectively. In the GSPM and related quarkonium constructions, screening is explicitly described as an effective spectral encoding of meson-meson thresholds rather than a full coupled-channel treatment (Gonzalez, 2014). In white-dwarf applications, the screened potential alters diffusion and composition interfaces rather than the bulk stellar structure (Chen, 2017). In screened hybrid DFT, the target is not confinement or binding spectra but exchange nonlocality in a periodic medium (Tran et al., 2011). This suggests that the term is best understood operationally: a screened potential model replaces a bare interaction by a reduced-range or saturated effective interaction appropriate to a specified physical regime.

Limitations are correspondingly model specific. Threshold-based hadron models omit explicit channel widths and dynamical mixing (González, 2015). Higher bottomonium states in screened quark models remain sensitive to neglected $-e^{-\delta r}/r$05-$-e^{-\delta r}/r$06 mixing and coupled-channel effects above open-flavor thresholds (Kaushal et al., 20 May 2025). Nonrelativistic screened charmonium models reproduce masses more successfully than decay widths (M et al., 10 Mar 2025). In plasma screening, simplified hydrodynamic potentials can fail qualitatively when benchmarked against kinetic-theory results (Moldabekov et al., 2015). Even so, across these disparate domains the screened potential model remains a durable organizing principle: retain the successful short-distance structure of the unscreened interaction, encode medium or threshold effects through screening, and use the resulting effective potential to obtain tractable spectra, wave functions, or response properties.