Inverse Exponential Potentials: Quantum & Cosmology

• Inverse exponential potentials are defined by an exponential decay with respect to the inverse variable, leading to sharply confining wells or plateaus and enabling exact analytical solutions in both quantum and cosmological settings.
• They facilitate rigorous spectral analysis and resonance quantization using modified Bessel, confluent hypergeometric, and Heun functions, which are crucial for modeling scattering properties and bound-state spectra.
• In inflationary cosmology, these potentials generate flat slow-roll conditions compatible with observational data and support graceful exit mechanisms through additional steep exponential terms.

The term "inverse exponential potential" designates a distinct class of potentials characterized by decaying or singular behavior as a function of a spatial or field variable. Across quantum mechanics and cosmological model-building, such potentials have been formulated as $V(r) = \lambda e^{-r}$ in Schrödinger theory and %%%%1%%%% in inflationary cosmology. Their defining trait is an exponential dependence on the inverse of the relevant variable, leading to sharply confining wells or plateaus in the associated dynamical equations. Rigorous analytical methods yield exact or conditionally exact solutions, underpinning rich spectra of resonances, scattering properties, and cosmological predictions.

1. Quantum Inverse Exponential Potentials: Formulation and Analytic Solutions

For a particle of mass $m$ subject to $V(r) = \lambda e^{-r}$ with $\lambda > 0$, the radial Schrödinger equation is cast as

$$-\frac{\hbar^2}{2m} \frac{d^2}{dr^2}\psi(r) + \lambda e^{-r} \psi(r) = E \psi(r),$$

subject to $\psi(0)=0$ and an outgoing wave condition at $r\to\infty$. By introducing dimensionless variables ($\epsilon\equiv 2mE/\hbar^2$, $\Lambda\equiv 2m\lambda/\hbar^2$), the equation reduces to

$$-\psi''(r) + \Lambda e^{-r} \psi(r) = \epsilon \psi(r),$$

which, via variable transformations $t = 2\sqrt{\Lambda} e^{-r/2}$ and $\mu = \sqrt{-4\epsilon}$, becomes the standard modified Bessel equation

$t^2 \psi_{tt} + t \psi_t - (t^2 + \mu^2)\psi = 0.$

The general solution in terms of modified Bessel functions is

$\psi(r) = A I_\mu (t) + B I_{-\mu} (t),$

where $I_\pm$ are the modified Bessel functions of the first kind. The Siegert (outgoing-wave) boundary condition at $r\to\infty$ selects $A\ne0$, $B=0$, while boundary regularity at $r=0$ imposes the transcendental quantization condition $I_\mu (2\sqrt{\Lambda}) = 0$. The energy eigenvalues are $E_n = (\hbar^2/2m)(-\mu_n^2/4)$, where the (typically complex) $\mu_n$ solve the Bessel function zero condition (Garcia, 2023).

An equivalent representation via confluent hypergeometric functions arises from the substitution $y = e^{-r}$, recasting the problem in terms of $U(a,b,z)$ (Tricomi's function) and leading to an equivalent quantization equation $U(a,b,2\sqrt{\Lambda})=0$.

2. Spectral Structure, Resonances, and Associated Families

The zeros $\mu_n$ of $I_\mu(z)$ populate the $\mathrm{Re}\,\mu<0$ half-plane and occur as either negative-real "virtual" states or complex-conjugate resonance pairs. Under complex $r$-plane dilation ($r\mapsto r e^{i\theta}$), the potential interpolates between the inverse-exponential "barrier" $V(r)=\Lambda e^{-r}$ and the "wall" $V_w(r)=\Lambda e^{r}$. The quantization condition generalizes to a family of transcendental equations: $$\sin(\pi\nu)^{-1} \left[ e^{m\pi i\nu} I_{-\nu}(2\sqrt{\Lambda}) - e^{-m\pi i\nu} I_\nu(2\sqrt{\Lambda}) \right] = 0,\quad m=0, \pm1, \pm2, \ldots$$ where for $m=0$ one obtains the wall quantization $K_\nu(2\sqrt{\Lambda})=0$ (with $K$ the modified Bessel function of the second kind), and for $|m|=1$ the original barrier $\mu_n$ sequence (Garcia, 2023). In the limit $\Lambda\rightarrow\infty$, two solution classes emerge:

• Barrier-type resonances: Roots tending to the $\mu_n$ defined by $I_\mu$ zeros, exhibiting exponential convergence as shown via steepest-descent (Debye) analysis.
• Rational wall-type values: Other roots "freeze" at rational values $k/m$ determined by $\sin(m\pi\nu)=0$ with $\sin(\pi\nu)\neq0$.

The Riccati–Padé method (RPM) efficiently reproduces both classes of resonances numerically, offering exponential convergence in the sequence of Hankel determinant roots as a function of series truncation order (Garcia, 2023).

3. Generalizations and Conditionally Exactly Solvable Exponential-Type Potentials

Further generalization yields conditionally exactly solvable potentials of the form

$$V_\pm(x;m) = \frac{m^2}{e^x-1} \pm \frac{m}{2} \frac{e^x}{(e^x - 1)^{3/2}},$$

for $x>0$ and $m>0$, with ground-state-absence due to the non-integrability of the zero mode. These partner potentials are multiplicative shape-invariant under supersymmetric quantum mechanics: $V_+(x,m) = V_-(x,-m).$ Their continuum eigenfunctions are exact linear combinations of two hypergeometric functions, and the exact scattering amplitudes $S_\pm(\omega)$ are derived in closed form (in terms of products and quotients of Gamma functions and hypergeometric parameters). For $V_+$, the potential is a monotonically decaying repulsive exponential wall; for $V_-$, the potential profile combines a singular well and finite barrier, again supporting only the continuum (Lopez-Ortega, 2016).

4. Singular Inverse–Exponential Potentials and Bound-State Spectra

Potentials with short-range exponential decay but inverse square-root singularity at the origin take the form

$$V(x) = -V_0 \frac{e^{-x/\sigma}}{\sqrt{1-e^{-x/\sigma}}},\quad V_0>0,\;\sigma>0,\;x>0,$$

with the behaviors $V(x)\sim -V_0 \sqrt{\sigma/x}$ as $x\to0^+$ and $V(x)\sim -V_0 e^{-x/\sigma}$ as $x\to\infty$. The corresponding Schrödinger equation reduces via $z=e^{-x/\sigma}$ and ansatz $\psi(x) = (z+1)^{\alpha_1} (z-1)^{\alpha_2} u(z)$ to the general Heun equation, with bound-state quantization given exactly in terms of ${}_3F_2$ (Clausen function) roots (Ishkhanyan, 2018). Only a finite number of bound states exist, determined semiclassically by $$n_\mathrm{max} \approx \left\lfloor 2\sqrt{m\sigma^2 V_0 /\hbar^2} - 1 \right\rfloor$$.

5. Inverse Exponential Potentials in Inflationary Cosmology

In cosmology, the “inverse exponential” (IExp) potential takes the form

$$V(\phi) = V_0 \exp(-\alpha M_p/\phi),\quad \phi>0,\;\alpha>0,$$

where $V_0$ sets the energy scale and $M_p = (8\pi G)^{-1/2}$ is the reduced Planck mass (Hossain, 1 Feb 2026). This potential motivates a minimal and theoretically robust monotonic slow-roll inflationary scenario. Key properties include:

• Tracker parameterization: $\lambda(\phi)=\alpha/\phi^2$, $\Gamma(\phi)=1-2\phi/\alpha$, yielding flat slope and strong concavity for large $\phi$.
• Slow-roll parameters: $\epsilon_V = \alpha^2 M_p^4/(2\phi^4)$, $$\eta_V = \alpha^2 M_p^4/\phi^4 - 2\alpha M_p^3/\phi^3$$; inflation ends when $\epsilon_V(\phi_\mathrm{end})=1$, with $\phi_\mathrm{end}=[\alpha M_p^2/\sqrt{2}]^{1/2}$.
• Number of e-foldings: $N_* \simeq (\phi_*^3-\phi_\mathrm{end}^3)/(3M_p^3)$, leading to $\phi_* \simeq [3\alpha N_*]^{1/3} M_p$ for $N_*\gg1$.
• Observational predictions: For $N_*\approx50$–$60$ and typical $\alpha\sim\mathcal{O}(0.1$–$10)$, predicts $n_s \approx 0.97$, $r \lesssim 0.01$–$0.03$, fully compatible with the 1$\sigma$ region of the combined SPA+BK+DESI2 CMB/BAO datasets.

To achieve a graceful exit from inflation, an additional steep exponential term $e^{-\beta\phi/M_p}$ with $\beta\gg1$ is introduced, generating a global minimum and enabling reheating by post-inflationary field oscillation. This yields maximum reheating temperatures $T_\mathrm{rh}\sim10^{12}$–$10^{13}$ GeV, contingent on the post-inflationary effective equation of state (Hossain, 1 Feb 2026).

6. Physical Properties, Methodological Implications, and Exact Solvability

Inverse exponential potentials, whether singular or regular, exhibit quintessential "short-range" characteristics, exponential decay, or singular confinement at the origin. In quantum contexts, exact or conditionally exactly solvable analytic structures arise, linked to modified Bessel or hypergeometric function solutions. The spectrum may be discrete, continuous, or resonant depending on boundary conditions and parameter regimes. Supersymmetry-based methods, complex variable continuation, and the Riccati–Padé approach provide diverse analytic and numerical access to resonance structure and spectral quantization (Garcia, 2023, Lopez-Ortega, 2016, Ishkhanyan, 2018).

Table: Core Formulations of Inverse Exponential Potentials in Quantum and Cosmological Contexts

Physical context	Potential Formulation	Solution Structure
Schrödinger (barrier)	$V(r)=\lambda e^{-r}$	Bessel/Hypergeometric, transcendental
Schrödinger (well)	$V(x)=-V_0 e^{-x/\sigma}/\sqrt{1-e^{-x/\sigma}}$	Heun, Gauss ${}_2F_1$, finite spectrum
SUSY partners	$$V_\pm(x;m)=\frac{m^2}{e^x-1}\pm \frac{m}{2}\frac{e^x}{(e^x-1)^{3/2}}$$	Two-term hypergeometric, scattering
Inflationary potential	$V(\phi)=V_0 \exp(-\alpha M_p/\phi)$	Slow-roll inflation, graceful exit

7. Significance and Outlook

The inverse exponential potential encapsulates a unified analytic framework underlying resonance phenomena in quantum theory and slow-roll inflation in early-universe cosmology. Its analytic tractability—via Bessel, confluent hypergeometric, and Heun-class reductions—unlocks rigorous access to spectral properties and dynamical predictions. In cosmology, its predictive agreement with SPA+BK+DESI2 $n_s$–$r$ constraints underscores its observational viability, while quantum versions remain central in the theoretical investigation of exactly and conditionally exactly solvable models (Garcia, 2023, Hossain, 1 Feb 2026, Lopez-Ortega, 2016, Ishkhanyan, 2018).