On positive solutions of Lane-Emden equations on the integer lattice graphs

Basic Information

• Paper ID: 2510.08947
• Title: On positive solutions of Lane-Emden equations on the integer lattice graphs
• Authors: Huyuan Chen (Fudan University), Bobo Hua (Fudan University), Feng Zhou (East China Normal University)
• Classification: math.AP (Analysis of PDEs)
• Publication Date: October 10, 2025
• Paper Link: https://arxiv.org/abs/2510.08947v1

Abstract

This paper investigates the existence and nonexistence of positive solutions to the Lane-Emden equation $-\Delta u = Q|u|^{p-2}u$ on the $d$-dimensional integer lattice $\mathbb{Z}^d$ as well as on half-space and quadrant domains, where $d \geq 2$, $p > 0$, and $Q$ is a Hardy-type positive potential function satisfying $Q(x) \sim (1+|x|)^{-\alpha}$ with $\alpha \in [0,+\infty]$. The authors identify the Sobolev supercritical region in the parameter space $(\alpha, p)$, where existence of positive solutions is established via variational methods. Conversely, in the Serrin subcritical region, nonexistence is proved through iterative analysis of decay behavior at infinity. Notably, in the full space and half-space domains, an intermediate region exists between the Sobolev critical line and the Serrin critical line where the existence of positive solutions remains an open problem, whereas no such intermediate region exists in the quadrant domain.

Research Background and Motivation

Problem Background

The Lane-Emden equation is a classical model for semilinear elliptic differential equations, originating from astrophysics to describe the structure of self-gravitating, spherically symmetric polytropic fluids in hydrostatic equilibrium. The standard form is: $-\Delta_{\mathbb{R}^d}u = |u|^{p-2}u \quad \text{in } \mathbb{R}^d$

Research Motivation

1. Extension from continuous to discrete: While Lane-Emden equations in the continuous setting have been extensively studied, research on discrete lattice graphs remains relatively limited.
2. Identification of critical exponents: Determining the precise critical exponents for existence of positive solutions on different domains (full space, half-space, quadrant).
3. Effect of potential functions: Investigating the influence of Hardy-type potential functions $Q(x)$ on solution existence.
4. Geometric dependence: Exploring how different geometric domains affect critical exponents.

Limitations of Existing Research

• Gu-Huang-Sun 15 established nonexistence for $d \geq 3$ and $p \leq \frac{d}{d-2}+1$
• Hua-Li 16 proved existence for $p > \frac{2d}{d-2}$
• The existence problem in the range $\frac{d}{d-2}+1 < p \leq \frac{2d}{d-2}$ remains open

Core Contributions

1. Complete critical theory: Establishes a complete critical exponent theory for three distinct domains (full space $\mathbb{Z}^d$, half-space $\mathbb{Z}^d_+$, quadrant $\mathbb{Z}^d_*$).
2. Unified analytical framework: Provides a unified framework for handling different domains through integral equation methods and variational techniques.
3. Precise critical exponents: Identifies domain-geometry-dependent Serrin and Sobolev exponents:
   • Full space: Serrin exponent $1+\frac{d-\alpha}{d-2}$, Sobolev exponent $1+\frac{d-\alpha}{d-2}+\frac{2-\alpha}{d-2}$
   • Half-space: Serrin exponent $1+\frac{d-\alpha}{d-1}$, Sobolev exponent $1+\frac{d-\alpha}{d-1}+\frac{1-\alpha}{d-1}$
   • Quadrant: Serrin exponent $1+\frac{d-\alpha}{d}$, Sobolev exponent $1+\frac{d-\alpha}{d}-\frac{\alpha}{d}$
4. Identification of open problems: Clearly delineates intermediate regions in full space and half-space where positive solution existence remains open.

Methodology Details

Problem Formulation

Investigation of positive solution existence for the Lane-Emden equation on three distinct domains:

1. Full space: $-\Delta u = Q|u|^{p-2}u$ in $\mathbb{Z}^d$
2. Half-space: $\{-\Delta u = Q|u|^{p-2}u$ in $\mathbb{Z}^d_+$, $u = 0$ on $\partial\mathbb{Z}^d_+\}$
3. Quadrant: $\{-\Delta u = Q|u|^{p-2}u$ in $\mathbb{Z}^d_*$, $u = 0$ on $\partial\mathbb{Z}^d_*\}$

Core Methodological Framework

1. Integral Equation Transformation

Converting the original differential equation into an integral equation: $u = \Phi_{d,\beta} * (Q|u|^{p-2}u)$ where $\Phi_{d,\beta}$ is the fundamental solution for the corresponding domain.

2. Variable Substitution Technique

Introducing the substitution $v = Q^{1/p'}|u|^{p-2}u$, transforming the equation into: $|v|^{p'-2}v = Q^{1/p}\Phi_{d,\beta} * (Q^{1/p}v)$ which possesses a variational structure.

3. Energy Functional

Defining the energy functional: $J_0(v) = \frac{1}{p'}\int_{\mathbb{Z}^d}|v|^{p'}dx - \frac{1}{2}\int_{\mathbb{Z}^d}vK_{p,\beta}(v)dx$ where $K_{p,\beta}(v) = Q^{1/p}\Phi_{d,\beta} * (Q^{1/p}v)$ is the Birman-Schwinger operator.

Technical Innovations

1. Precise Fundamental Solution Estimates

Constructing precise fundamental solution estimates for each domain:

• Full space: $\Phi_d(x,y) \leq c_1(1+|x-y|)^{2-d}$
• Half-space: $\Phi_{d,+}(x,y) \leq c_1(1+|x-y|)^{1-d}$ and $\Phi_{d,+}(x,y) \geq \frac{1}{c}x_1(1+|x-y|)^{-d}$
• Quadrant: $\Phi_{d,*}(x,y) \leq c_1(1+|x-y|)^{-d}$ and $\Phi_{d,*}(x,y) \geq \frac{1}{c}x_1x_2(1+|x-y|)^{-d-2}$

2. Non-vanishing Lemma

Proving the crucial non-vanishing result (Lemma 3.1), ensuring compactness of Palais-Smale sequences.

3. Iterative Nonexistence Proof

Systematically proving nonexistence in subcritical cases through constructing test functions $w_\tau(x) = |x|^{-\tau}$ and iterative sequences $\{\tau_j\}$.

Experimental Setup

Theoretical Verification Framework

This is primarily theoretical research verified through rigorous mathematical proofs, including:

1. Existence Proof Strategy

• Sobolev supercritical case: Using Mountain Pass theorem
• Critical case: Using eigenvalue theory
• Sublinear case: Using upper and lower solution methods

2. Nonexistence Proof Strategy

• Constructing appropriate test functions
• Using integral divergence arguments
• Iteratively improving lower bound estimates

3. Fundamental Solution Construction

• Half-space: $\Phi_{d,+}(x,y) = \Phi_d(x-y) - \Phi_d(x-y^*)$
• Quadrant: $\Phi_{d,*}(x,y) = \frac{1}{2}(2\Phi_d(x-y) - \Phi_d(x-y^*) - \Phi_d(x-y^\#))$

Main Results

Theorem 1.1 (Full space $\mathbb{Z}^d$)

1. Existence: When $\alpha \in [0,+\infty]$ and $p \in (2,+\infty) \cap (2^*_{1,\alpha},+\infty)$, the problem admits at least one positive solution.
2. Sublinear existence: When $\alpha > 2$ and $p \in [1+\frac{d-\alpha}{d-2}, 2) \cap (1,2)$, the problem admits a unique positive solution.
3. Nonexistence: When $\alpha \in (-\infty,d)$ and $p \in (1, 1+\frac{d-\alpha}{d-2})$ or $p = 1+\frac{d-\alpha}{d-2} > 2$, the problem has no positive solution.

Theorem 1.2 (Half-space $\mathbb{Z}^d_+$)

Similar results with modified critical exponents:

• Serrin exponent: $1+\frac{d-\alpha}{d-1}$
• Sobolev exponent: $2^*_{1/2,\alpha} = \frac{2(d-\alpha)}{d-1}$

Theorem 1.3 (Quadrant $\mathbb{Z}^d_*$)

• Serrin exponent: $1+\frac{d-\alpha}{d}$
• Note: No intermediate open region exists in the quadrant domain

Key Findings

1. Geometric Dependence

Critical exponents strongly depend on domain geometry:

• Dimension reduction leads to lower critical exponents
• The number of boundary conditions affects critical behavior

2. Open Regions

Identifying open regions in full space and half-space: $\{(\alpha,p) : \alpha \in (0,2), 1+\frac{d-\alpha}{d-2} < p \leq 1+\frac{d+2-2\alpha}{d-2}\}$

3. Decay Estimates

When $Q \geq C > 0$, all positive solutions satisfy $\lim_{|x|\to\infty} u(x) = 0$.

Related Work

Continuous Case Studies

• Pohožaev identity: Nonexistence for $p \in (1, \frac{2d}{d-2})$
• Aubin-Talenti bubble solutions: Exact solutions in the critical case $p = \frac{2d}{d-2}$
• Variational methods: Existence for $p > \frac{2d}{d-2}$

Discrete Case Studies

• Gu-Huang-Sun 15: Nonexistence for $d \geq 3$, $p \leq \frac{d}{d-2}+1$
• Hua-Li 16: Existence for $p > \frac{2d}{d-2}$
• Elliptic equations on graphs: Semilinear elliptic problems on general graphs

Potential Function Research

• Ni 26: Scalar curvature equations in conformal geometry
• Bianchi et al. 3: Existence with radially symmetric decreasing potentials
• Cao-Peng 5: Polynomial decay solutions with vanishing potentials at infinity

Conclusions and Discussion

Main Conclusions

1. Complete classification: Provides complete classification of the parameter space $(\alpha,p)$ for three domains.
2. Geometric effects: Demonstrates that domain geometry fundamentally affects critical exponents.
3. Open problems: Identifies precise regions requiring further investigation.

Limitations

1. Intermediate regions: Intermediate regions in full space and half-space remain open.
2. General potential functions: Results primarily target Hardy-type potentials $Q(x) \sim (1+|x|)^{-\alpha}$.
3. Higher-dimensional generalizations: Some techniques may be difficult to extend to more general graph structures.

Future Directions

1. Resolution of intermediate regions: New techniques needed to address regions between critical lines.
2. General graphs: Extension to more general graph structures.
3. Numerical verification: Numerical experiments to validate theoretical predictions.

In-Depth Evaluation

Strengths

1. Theoretical completeness: Provides a complete theoretical framework for Lane-Emden equations on three domains.
2. Technical innovation: Integral equation methods and precise fundamental solution estimates are important contributions.
3. Geometric insights: Reveals profound effects of domain geometry on critical behavior.
4. Open problem identification: Clearly delineates directions requiring further research.

Weaknesses

1. Computational complexity: Fundamental solution construction and estimation involve intricate calculations.
2. Technical sophistication: Some proof techniques may be difficult to generalize.
3. Practical applications: Limited discussion of practical application backgrounds.

Impact

1. Theoretical contribution: Establishes foundations for discrete Lane-Emden equation theory.
2. Methodological value: Integral equation methods applicable to other discrete elliptic problems.
3. Interdisciplinary impact: Connects analysis, geometry, and combinatorial mathematics.

Applicable Scenarios

1. Mathematical physics: Discretized physical models
2. Network science: Diffusion processes on large networks
3. Numerical analysis: Discrete studies of continuous problems

References

The paper cites 45 important references covering classical and cutting-edge work in elliptic equation theory, variational methods, and graph analysis.

Overall Assessment: This is a high-quality theoretical mathematics paper making significant contributions to discrete Lane-Emden equation theory. With high technical depth and complete results, it establishes a solid foundation for further development in this field.