---
title: Singular Lane–Emden–Fowler Equations
url: https://www.emergentmind.com/topics/singular-lane-emden-fowler-equation
type: topic
---

# Singular Lane–Emden–Fowler Equations

Searching arXiv for recent papers on singular Lane–Emden–Fowler equations and related numerical/analytic methods.
The singular Lane–Emden–Fowler equation denotes a family of nonlinear differential equations in which singular behavior is generated either by a coefficient that blows up at a distinguished point, typically \(x=0\), or by boundary data that force the unknown to vanish where the nonlinearity contains a negative power of the solution. In the material surveyed here, two principal settings appear. The first is the semilinear elliptic PDE
\[
-\Delta u=f(X)\,u^{-\gamma},
\]
posed in a bounded Lipschitz domain with Dirichlet data, where the singularity is induced by the boundary condition \(u=0\) on part of \(\partial\Omega\) [2503.16095]. The second is the classical ODE family
\[
y''+\frac{\alpha}{x}y'+Q(x)\,y=f(x,y),
\]
or closely related Emden–Fowler forms, where the term \(\alpha/x\) is singular at the origin [2410.05409]. Across both settings, the subject combines local singular analysis, existence and uniqueness theory, asymptotics, boundary regularity, and a wide range of numerical methodologies.

## 1. Problem classes and sources of singularity

In the elliptic PDE setting studied by Guo, Li, and Zhang, one works in a bounded Lipschitz domain \(\Omega\subset\mathbb R^n\), \(n\ge2\), with a positive Hölder-continuous weight \(f\) satisfying \(0<\lambda\le f\le\Lambda\), and seeks a classical solution
\[
u>0\ \text{in }\Omega,\qquad -\Delta u=f(X)\,u^{-\gamma},\qquad u=0\ \text{on }\Gamma,
\]
together with arbitrary continuous data on \(\partial\Omega\setminus\Gamma\). The singularity arises precisely because \(u\to0\) on \(\Gamma\), making the right-hand side unbounded [2503.16095].

In the ODE literature, the singular Lane–Emden–Fowler equation typically refers to a second-order equation with a radial term \(\alpha/x\). A representative form is
\[
\frac{d^2y}{dx^2}+\frac{\alpha}{x}\frac{dy}{dx}+Q(x)\,y=f(x,y),
\qquad a\le x\le b,
\]
with singular initial conditions at \(x=0\) [2410.05409]. Closely related formulations include
\[
u''+\frac{N-1}{x}u'=h(x,u)
\]
for singular initial-value or boundary-value problems [2301.01041], and
\[
y''+\frac{k}{x}y'+f(x,y)=0
\]
for Emden–Fowler type IVPs [1604.05344]. In each case, the term \(k/x\) or \((N-1)/x\) becomes unbounded as \(x\to0\), so regularity conditions such as \(y'(0)=0\) or \(u'(0)=0\) are imposed to ensure bounded behavior.

A separate class of singularity appears in perturbed boundary-value problems with singular endpoints at both \(x=0\) and \(x=1\). One example is
\[
p(x)\,u''+q(x)\,u'+r(x)\,u+\delta\,u^p=0,
\]
where \(q\) and \(r\) contain poles at \(x=0\) and \(p(x)\) vanishes at \(x=1\) [1810.01410]. Another variant places the singular coefficient at the upper endpoint \(t=1\), in the time-singular form
\[
\frac{d^2x}{dt^2}=-\frac{a(t,x(t))}{1-t}\frac{dx}{dt}+g(t,x(t)),
\qquad 0\le t<1,
\]
so that \((1-t)^{-1}\) becomes singular as \(t\to1^-\) [2405.08736].

These formulations are related by the common feature that a formally lower-order term, or the nonlinearity itself, becomes dominant near a singular point or boundary. This suggests that “singular Lane–Emden–Fowler equation” is best viewed as a structural class rather than a single canonical equation.

## 2. Boundary regularity and cone-frequency classification for the PDE

For the semilinear elliptic problem in a Lipschitz domain, the principal geometric device is the limiting cone at a boundary point. If \(\Sigma\subset\partial B_1\) is a spherical open set with Lipschitz boundary, its cone is
\[
\mathrm{Cone}_\Sigma=\{r\theta:r>0,\theta\in\Sigma\}.
\]
Let \(\lambda_\Sigma\) be the first Dirichlet eigenvalue on \(\Sigma\), and define the minimal frequency \(\phi_\Sigma>0\) by
\[
\phi_\Sigma(n+\phi_\Sigma-2)=\lambda_\Sigma.
\]
The pair \((\Sigma,\gamma)\) is called sub-critical if \(2/(1+\gamma)<\phi_\Sigma\), critical if \(2/(1+\gamma)=\phi_\Sigma\), and super-critical if \(2/(1+\gamma)>\phi_\Sigma\) [2503.16095].

This classification governs the boundary growth rate. If \(u\) solves \(-\Delta u=f\,u^{-\gamma}\) in a Lipschitz cylinder \(GC_{3R}\), then in \(GC_R\) there is a constant \(C=C(n,L,\gamma,\lambda,\Lambda)\) such that:

| Regime | Upper growth rate |
|---|---|
| Sub-critical | \(u(X)\le C\,R^{-2/(1+\gamma)}\|u\|_{L^\infty}|X/R|^{2/(1+\gamma)}\) |
| Critical | \(u(X)\le C\,R^{-2/(1+\gamma)}\|u\|_{L^\infty}|X/R|^{\phi_\Sigma}\ln(2R/|X|)\) |
| Super-critical | \(u(X)\le C\,R^{-2/(1+\gamma)}\|u\|_{L^\infty}|X/R|^{\phi_\Sigma}\) |

Moreover, if \(\Omega=\mathrm{Cone}_\Sigma\) near \(0\), these rates are optimal [2503.16095].

Along the inward normal ray \(X=t e_n\), the corresponding lower bounds are also regime-dependent:
\[
u(t e_n)\gtrsim t^{2/(1+\gamma)} \quad\text{in the sub-critical case},
\]
\[
u(t e_n)\gtrsim t^{\phi_\Sigma} \quad\text{in the super-critical case},
\]
and
\[
u(t e_n)\gtrsim t^{\phi_\Sigma}(\ln(1/t))^{\phi_\Sigma/2}
\quad\text{in the critical case}
\]
[2503.16095]. The critical logarithmic correction is a distinctive feature of the borderline frequency balance.

The significance of this framework is that boundary regularity is not described by a single universal exponent. Instead, the interaction between the singular power \(u^{-\gamma}\) and the local cone geometry determines the asymptotic profile. A plausible implication is that non-smooth boundary geometry is not merely a technical complication but part of the principal asymptotic mechanism.

## 3. Well-posedness, comparison principles, and analytic continuation

For the Lipschitz-domain PDE, the paper establishes well-posedness together with boundary growth and boundary Harnack estimates [2503.16095]. Several auxiliary results support this theory. A maximum-principle lemma states that if \(u\) and \(v\) are respectively super- and sub-solutions of \(-\Delta w=f\,w^{-\gamma}\) and \(u\ge v\) on \(\partial\Omega\), then \(u\ge v\) in \(\Omega\). A non-degeneracy lemma states that any super-solution \(u\) in \(B_r\) satisfies \(u(0)\gtrsim r^{2/(1+\gamma)}\), preventing excessively flat behavior near the boundary [2503.16095].

For singular ODEs, existence and uniqueness are often obtained by reformulating the equation as an equivalent regular problem. Seiler and Seiß rewrite
\[
u''+\frac{N-1}{x}u'=h(x,u)
\]
as a three-dimensional autonomous system in \((x,u,u')\):
\[
\frac{dx}{ds}=x,\qquad \frac{du}{ds}=x\,u',\qquad \frac{du'}{ds}=x\,h(x,u)-(N-1)\,u'.
\]
The point \(\rho=(0,u_0,0)\) is a stationary point, and the unique solution of the singular initial-value problem is exactly the integral curve lying on the one-dimensional unstable manifold of \(\rho\) [2301.01041]. Under the hypotheses \(g'(y)>0\) and \(\partial_{u'}f(y,u_0,u_1)<0\) in the generalized formulation \(g(x)u''-f(x,u,u')=0\), they state that there is a unique smooth solution through the singular point whose prolonged graph is that unstable manifold [2301.01041].

In the perturbed endpoint-singular problem on \([0,1]\), local analytic solutions are constructed separately at \(x=0\) and \(x=1\) using convergent Frobenius-type series, and a Lyapunov function
\[
H(x)=\frac{p(x)}{2}u'(x)^2+\frac{\delta}{p+1}u(x)^{p+1}+\frac{r}{2}u(x)^2
\]
is used to show \(H'<0\) on \([0,1]\), ruling out movable blow-up in \((0,1)\) [1810.01410]. Under the stated matching and Lyapunov-type conditions, the paper concludes that there exists at least one nontrivial solution analytic at both singular endpoints.

These approaches illustrate three distinct well-posedness paradigms: comparison and barrier methods for singular semilinear PDEs, dynamical-systems desingularization for ODEs at \(x=0\), and endpoint power-series matching plus Lyapunov control for problems singular at both ends.

## 4. Boundary Harnack theory and the failure of the classical linear picture

A central contribution in the Lipschitz-domain PDE theory is the boundary Harnack principle for the singular equation. If \(u,v\ge0\) both solve
\[
-\Delta w=f\,w^{-\gamma}
\]
in \(GC_{3R}\) with zero boundary data on \(\Gamma\), then in \(GC_R\),
\[
C^{-1}R^{2/(1+\gamma)}\|v\|_\infty^{-1}\le \frac{u}{v}
\le C\,\|u\|_\infty R^{-2/(1+\gamma)},
\]
with \(C=C(n,L,\gamma,\lambda,\Lambda)\) [2503.16095]. The paper identifies this as a generalization of the classical Kemper estimate for harmonic functions and states that, to the authors’ knowledge, it is the first Kemper-type estimate for singular semilinear equations [2503.16095].

The semilinear setting differs sharply from the linear one. In particular, the boundedness of \(u/v\) does not imply continuity up to \(\Gamma\), and this failure is formalized in Theorem 1.6 of the paper [2503.16095]. Under sub-critical or critical hypotheses, one can show \(u/v\to1\) at the boundary, while under additional smoothness or convexity of \(\Gamma\) one recovers full continuity [2503.16095].

A key technical input is the lack of a suitable upper barrier in the singular semilinear case. To overcome this, the paper introduces iterative upper-barrier constructions and a Campanato-type iteration. One lemma states that if \(0\le u\le v\) in \(B_1\), with \(-\Delta v\ge0\) and \(-\Delta u\le-t\,u\) for \(0<t\ll1\), then \(u\le(1-c\,t)\,v\) in \(B_{1/2}\); this underlies the iteration proving \(u/v\to1\) in the sub-critical and critical regimes [2503.16095].

The resulting theory corrects a common overextension of linear intuition: for singular semilinear equations, the ratio of two positive solutions need not inherit the boundary continuity properties familiar from harmonic analysis. This distinction is one of the main conceptual differences between the singular Lane–Emden–Fowler PDE and linear elliptic theory.

## 5. Asymptotic constructions, barriers, and auxiliary functions

The proof of sharp boundary growth in cones is built around a rescaled difference
\[
V(X)=2^{2/(1+\gamma)}U(X/2)-U(X),
\]
where \(U\) solves the constant-weight problem \(-\Delta U=U^{-\gamma}\). This function is nonnegative, subharmonic in the cone, and vanishes on its boundary [2503.16095]. By comparison with the homogeneous harmonic profile
\[
H_\Sigma(X)=|X|^{\phi_\Sigma}E_\Sigma(X/|X|),
\]
one obtains
\[
V(X)\le C\,H_\Sigma(X)
\]
[2503.16095].

The barrier iteration is then encoded in the dyadic quantities
\[
r_k=2^{-k},\qquad
A_k=r_k^{-2/(1+\gamma)}\max_{|X|\le r_k}U(X),
\]
for which
\[
A_{k+1}-A_k\lesssim r_k^{\phi_\Sigma-2/(1+\gamma)}
=2^{(2/(1+\gamma)-\phi_\Sigma)k}.
\]
A geometric or log-geometric summation yields the three growth regimes [2503.16095]. The auxiliary function \(V\) thus measures how far \(U\) is from exact homogeneity of degree \(2/(1+\gamma)\), and its decay controls the dyadic increments of the barrier sequence.

In endpoint-singular ODE problems, analogous asymptotic information is encoded in local power-series expansions. At \(x=0\), one sets
\[
u(x)=\sum_{k=0}^\infty h_k\,x^k,\qquad h_0=c,\ h_1=0,
\]
while at \(x=1\), with \(y=1-x\),
\[
u(y)=\sum_{k=0}^\infty d_k\,y^k,\qquad d_0=b,
\]
and derives recurrences whose convergence is guaranteed by a normal-form reduction and an analytic existence theorem [1810.01410]. In the time-singular model, the classical Lane–Emden cases \(n=0,1,5\) are given explicitly as
\[
x(t)=1-\frac{t^2}{6},\qquad
x(t)=\frac{\sin t}{t},\qquad
x(t)=\left(1+\frac{t^2}{3}\right)^{-1/2},
\]
and for other integer indices a Frobenius-type power series about \(t=0\) is used [2405.08736].

These constructions indicate that singular Lane–Emden–Fowler problems are often controlled not by a direct closed-form solution but by carefully chosen comparison objects: harmonic profiles in cones, dyadic barrier sequences, unstable manifolds, or endpoint Frobenius expansions.

## 6. Numerical methods and computational formulations

The numerical literature represented here is diverse, but it is organized around a shared objective: regularize or bypass the singular structure while preserving the nonlinear dynamics.

One family uses basis expansions and collocation. The orthonormal polynomial wavelet method considers
\[
t\,y''(t)+k\,y'(t)+t\,f\bigl(t,y(t),t^k y'(t)\bigr)=0,\qquad 0<t\le1,
\]
with \(y'(0)=\alpha\) and \(a\,y(1)+b\,y'(1)=\beta\), expands the unknown in wavelet bases generated by orthonormal polynomials, and enforces the ODE at midpoint collocation nodes [1911.02004]. Nonlinearity is handled either by Newton’s quasilinearization or by Newton–Raphson on the discrete system [1911.02004]. The paper states that, as the resolution is increased, the computed solutions converge to exact or known solutions, and reports errors of order \(10^{-8}\)–\(10^{-14}\) in a Chandrasekhar stellar-structure test and maximum errors \(<10^{-10}\) in a thermal explosion problem [1911.02004].

A second family uses rational spectral or reproducing-kernel representations. The Rational Chebyshev of Second Kind collocation method solves
\[
y''(x)+\frac{\alpha}{x}y'(x)+f(x)\,g(y(x))=h(x),\qquad x>0,
\]
on \([0,\infty)\) by mapping the semi-infinite interval to \([-1,1]\) and using the ansatz
\[
y_N(x)=A+Bx+x^2\sum_{k=0}^N a_k\,R_k(x),
\]
so that the center conditions are built in [1508.07240]. The reproducing-kernel Hilbert-space method rewrites the singular IVP as \(L\,u=F(x,u)\), constructs a reproducing kernel \(R_x(y)\), and represents the solution by a convergent series in an orthonormal system derived from kernel functions [1704.06830]. The latter paper states that the \(n\)-term approximation converges uniformly on \([a,T]\) and reports maximum absolute errors below \(10^{-6}\) in its examples [1704.06830].

A third family is based on dynamical systems and shooting. For
\[
u''+\frac{N-1}{x}u'=h(x,u),
\]
Seiler and Seiß compute the unstable manifold of the associated autonomous vector field and combine this with a shooting procedure for boundary-value problems [2301.01041]. The paper states that no special series expansions at \(x=0\) are needed, and reports that default Maple RKF4/5 tolerances \(\tau_{\rm rel}=10^{-6}\), \(\tau_{\rm abs}=10^{-7}\) gave 6–7 digits of accuracy on test problems [2301.01041].

A fourth family uses neural networks. The shifted Legendre neural network method builds the trial solution
\[
y_{\rm trial}(x;\Theta)=y_0+y_1(x-a)+(x-a)^2N(x;\Theta),
\]
so that \(y(0)=y_0\) and \(y'(0)=y_1\) are satisfied automatically and the singular \((\alpha/x)\) term does not cause numerical blow-up [2410.05409]. For the standard \(m=5\) Lane–Emden case, the paper reports max absolute error \(\lesssim5\times10^{-3}\), while for an inhomogeneous Lane–Emden–Fowler problem it reports max error \(O(10^{-11})\) [2410.05409]. Physics-informed neural networks were also benchmarked for second-, third-, and fourth-order singular ODEs, using both soft constraints in the loss and hard-constraint trial solutions such as \(\hat y(x)=1+x^2\mathcal{NN}(x)\) or \(\hat y(x)=x^3\mathcal{NN}(x)\) [2307.07302]. The hard-constraint formulation guarantees exact boundary satisfaction, and the paper reports errors ranging from \(\sim10^{-6}\) to \(\sim10^{-4}\) across the benchmark suite [2307.07302].

A fifth line of work emphasizes high-order iterative discretization. For the singular boundary-value problem
\[
y''(x)+\frac{\alpha}{x}y'(x)=f(x,y(x)),\qquad 0<x\le1,\qquad y'(0)=0,\qquad \mu y(1)+\sigma y'(1)=A,
\]
Dang et al. derive a continuous fixed-point iteration using the Green’s function and then discretize it by a corrected trapezoidal, Euler–Maclaurin-based rule [2508.19729]. They present three discrete iterative schemes for three specific cases and state that the resulting methods achieve eighth-order accuracy and convergence [2508.19729]. In one benchmark, the reported errors decrease from \(3.87\times10^{-10}\) at \(N=8\) to \(8.89\times10^{-18}\) at \(N=64\), with observed orders \(8.73\), \(9.30\), and \(7.34\) on successive refinements [2508.19729].

Taken together, these methods show no single dominant numerical paradigm. Collocation, spectral rational bases, reproducing kernels, unstable-manifold shooting, perturbation iteration, and neural-network residual minimization each target a different structural aspect of the singularity.

## 7. Applications, interpretations, and recurrent themes

The classical Lane–Emden equation appears in astrophysics through stellar structure and polytropic gas spheres. The time-singular study states that the radial Lane–Emden equation models phenomena in astrophysics such as stellar structure and is governed by polytropics with applications in isothermal gas spheres [2405.08736]. In the same source, the dimensionless variable \(x(t)\) is interpreted as proportional to the density profile of a self-gravitating polytropic gas sphere, with the polytropic index
\[
n=\frac{1}{\gamma-1}
\]
determined by the adiabatic exponent [2405.08736].

Other model problems represented in the numerical literature include Chandrasekhar’s stellar structure equation, thermal explosion models, shallow membrane cap equations, human-head heat conduction, Michaelis–Menten kinetics, oxygen-uptake problems, and the Thomas–Fermi equation after the Majorana transform [1911.02004; 2301.01041]. These examples show that the singular Lane–Emden–Fowler framework extends well beyond polytropes.

A recurrent theme across the analytic and computational work is that the singularity is often geometrically or structurally removable, but not negligible. In ODEs it may be “blown up” into a regular equilibrium of an autonomous system [2301.01041]; in PDEs it may be encoded by cone geometry and frequency [2503.16095]; in numerical schemes it is absorbed into the ansatz, Green’s function, or choice of collocation points [1911.02004; 2410.05409]. This suggests that successful analysis of singular Lane–Emden–Fowler equations depends less on suppressing the singularity than on reformulating it in a structure-compatible way.

A further recurrent theme is that singular problems do not automatically inherit the regularity properties of their nonsingular analogues. The failure of continuity of \(u/v\) in the semilinear boundary Harnack problem [2503.16095] is one example; the need for trial solutions that encode exact endpoint behavior in neural and collocation methods is another [2410.05409; 2307.07302]. In this sense, the singular Lane–Emden–Fowler equation functions as a meeting point of asymptotic analysis, nonlinear elliptic theory, singular dynamical systems, and high-accuracy numerical computation.

Source: https://www.emergentmind.com/topics/singular-lane-emden-fowler-equation