Modulus of continuity for solutions of non-local heat equations

Basic Information

• Paper ID: 2507.19023
• Title: Modulus of continuity for solutions of non-local heat equations
• Authors: Ben Andrews and Sophie Chen
• Classification: math.AP (Analysis of PDEs)
• Publication Date: November 14, 2025 (arXiv v2)
• Paper Link: https://arxiv.org/abs/2507.19023

Abstract

This paper generalizes the modulus of continuity method for parabolic equations (such as the method used to prove the fundamental gap conjecture) to non-local heat equations. Specifically, the authors prove that if a solution to a non-local heat equation has an initial modulus of continuity satisfying simple criteria, then this modulus is preserved at all subsequent times. The research covers three cases: non-local heat equations on ℝⁿ, domain heat equations on one-dimensional bounded intervals, and non-local heat equations with nonlinear terms. In extending the one-dimensional results, the authors discover a counterexample showing that an analogue of the non-local Payne-Weinberger inequality will depend on multiple geometric properties of bounded convex domains, rather than merely on the diameter.

Research Background and Motivation

Problem Background

1. Success of Classical Theory: The modulus of continuity method has achieved remarkable success in classical parabolic equation theory, particularly playing a key role in proving important results such as the Fundamental Gap Conjecture.
2. Rise of Non-local Operators: Non-local operators such as the fractional Laplacian have widespread applications in physics, probability theory, and geometric analysis, making the corresponding non-local heat equations important objects of study.
3. Theoretical Gap: Although the modulus of continuity theory for classical heat equations is mature, the corresponding theory for non-local heat equations has not been systematically established.

Research Significance

1. Theoretical Completeness: Generalizing classical methods to the non-local setting enriches the theoretical framework of partial differential equations.
2. Potential Applications: Provides technical tools for proving important results such as non-local versions of the Payne-Weinberger inequality.
3. Geometric Insights: Reveals through counterexamples the essential differences between non-local and classical problems regarding geometric dependence.

Limitations of Existing Methods

Classical modulus of continuity methods rely on local properties (such as gradient estimates) and cannot be directly applied to non-local operators, requiring the development of new techniques to handle long-range interactions.

Core Contributions

1. Modulus of Continuity Theorem on ℝⁿ (Theorem 2.3): Proves that on the full space ℝⁿ, the modulus of continuity of solutions to non-local heat equations is preserved during evolution, provided the initial modulus of continuity satisfies specific criteria.
2. Results for One-Dimensional Domain Heat Equations (Theorem 3.3): Establishes similar modulus of continuity preservation properties for domain heat equations on compact intervals I = -D/2, D/2.
3. Construction of High-Dimensional Counterexample (Section 3.2.1): Constructs an ingenious counterexample proving that the second eigenvalue of the domain fractional Laplacian cannot depend solely on the domain's diameter, contrasting with the classical Payne-Weinberger inequality.
4. Nonlinear Extension (Theorem 4.2): Extends the method to non-local heat equations with gradient-dependent nonlinear terms, demonstrating the flexibility of the technique.
5. Coupling-by-Reflection Technique: Develops the "coupling-by-reflection" technique to cleverly handle non-local integral terms, which is the core technical innovation of the proof.

Detailed Methodology

Problem Formulation

Consider the non-local linear operator: $L_\Omega u(x) := \int_\Omega \rho(z-x)(u(z)-u(x))dz, \quad x \in \Omega$

where:

• $\Omega \subset \mathbb{R}^n$ is an open set
• $\rho$ is a non-negative, non-increasing, rotationally symmetric kernel function
• $u$ belongs to an appropriate function space

Objective: Prove that if $u(x,t)$ satisfies the non-local heat equation $u_t(x,t) = L_\Omega u(x,t)$ and the initial data $u_0$ has modulus of continuity $\phi(\cdot, 0)$, then $\phi(\cdot, t)$ is a modulus of continuity for $u(\cdot, t)$ at all times $t \geq 0$, i.e., $|u(y,t) - u(x,t)| \leq 2\phi\left(\frac{|y-x|}{2}, t\right)$

Core Technical Framework

1. Auxiliary Function Construction

Define the auxiliary function: $Z_\epsilon(x,y,t) = u(y,t) - u(x,t) - 2\phi\left(\frac{|y-x|}{2}, t\right) - \epsilon e^{Ct}E(x,y)$

where $E(x,y)$ is a regularization term (such as $\Psi(x) + \Psi(y)$), and $C > K$ is a carefully chosen constant.

Key Idea: Prove that for all $\epsilon > 0$, $Z_\epsilon < 0$.

2. Proof by Contradiction Structure

Assume there exists $\epsilon_0 > 0$ such that $Z_{\epsilon_0}$ attains a non-negative maximum of 0 at some point $(x_0, y_0, t_0)$. Then:

• Time derivative condition: $\frac{\partial Z_\epsilon}{\partial t}\big|_{(x_0,y_0,t_0)} \geq 0$
• Spatial maximum condition: $\nabla Z_\epsilon = 0$

3. Coupling-by-Reflection Technique

This is the core technical innovation of the paper. Choose coordinates so that $e_1 = \frac{y-x}{|y-x|}$, and define the reflection map: $Rr = r - 2\langle r, e_1 \rangle e_1$

Key Observation: Due to the rotational symmetry of the kernel $\rho$, we have $\rho(r) = \rho(Rr)$.

Decompose the integration domain into three disjoint sets:

• $S = \{r : \langle r, e_1 \rangle = \frac{|x-y|}{2}\}$ (middle hyperplane)
• $L = \{r : \langle r, e_1 \rangle < \frac{|x-y|}{2}\}$ (left half-space)
• $R = \{r : \langle r, e_1 \rangle > \frac{|x-y|}{2}\}$ (right half-space)

Through variable substitution and reflection symmetry, the two integrals can be merged onto the same domain, thereby effectively applying the maximum principle.

4. Dimension Reduction Technique

By integrating out the $(n-1)$ dimensions, reduce the $n$-dimensional problem to a one-dimensional problem: $\tilde{\rho}(w) = \int_{\mathbb{R}^{n-1}} \rho(w,p)dp$

Require the odd extension $\tilde{\phi}$ of the modulus of continuity to satisfy the one-dimensional non-local heat equation: $\tilde{\phi}_t(r,t) = \int_{\mathbb{R}} \tilde{\rho}(w)(\tilde{\phi}(r+w,t) - \tilde{\phi}(r,t))dw$

Technical Innovations

1. Clever Choice of Regularization Function: For the fractional Laplacian $(−\Delta)^s$, choose $\Psi(x) = (1+|x|^2)^{\alpha/2}$, where $0 < \alpha < \min(2s, 2)$, ensuring: $K := \sup_{x \in \mathbb{R}^n} L_{\mathbb{R}^n}\Psi(x) < \infty$
2. Detailed Analysis of Boundary Cases: In the one-dimensional case (Theorem 3.3), systematically handle various geometric configurations near interval endpoints through translation and symmetrization techniques.
3. Control of Nonlinear Terms: In Theorem 4.2, successfully control the contribution of the nonlinear term $q(|Du|)$ through Lipschitz conditions and gradient growth bounds: $|D\Psi(x)| \leq K_\nabla \Psi(x)$

Main Results in Detail

Theorem 2.3 (Non-local Modulus of Continuity on ℝⁿ)

Statement: Let $u: \mathbb{R}^n \times [0,\infty) \to \mathbb{R}$ satisfy the non-local heat equation with $\|u\|_{L^\infty} \leq 1$. If $\phi: [0,\infty) \times [0,\infty) \to [0,\infty)$ satisfies:

• (a) $\phi' > 0$ and $\lim_{r \to 0^+}\phi(r,t) = 0$
• (b) $\phi(\cdot, 0)$ is a modulus of continuity for $u(\cdot, 0)$
• (c) The odd extension $\tilde{\phi}$ satisfies the one-dimensional non-local heat equation
• (d) There exists a regularization function $\Psi$ such that $\sup L_{\mathbb{R}^n}\Psi < \infty$

then $\phi(\cdot, t)$ is a modulus of continuity for $u(\cdot, t)$ at all times.

Proof Outline:

1. At the maximum point $(x_0, y_0, t_0)$, the time derivative yields: $0 \leq \int_{\mathbb{R}^n} \rho(r)(u(y+Rr) - u(x+r))dr - 2\phi_t(s) - CE(\Psi(x)+\Psi(y))$
2. After applying coupling-by-reflection, the integral term becomes: $2\int_{\mathbb{R}^n} \rho(r)(\tilde{\phi}(s-r_1) - \tilde{\phi}(s))dr = 2\tilde{\phi}_t(s)$
3. The regularization term satisfies: $E(L_{\mathbb{R}^n}\Psi(y) + L_{\mathbb{R}^n}\Psi(x)) \leq 2EK$
4. Since $C > K$ and $\Psi \geq 1$, we obtain a contradiction: $0 \leq 2EK - CE(\Psi(x)+\Psi(y)) \leq 2E(K-C) < 0$

Theorem 3.3 (Modulus of Continuity on One-Dimensional Compact Intervals)

For domain heat equations on $I = [-D/2, D/2]$, the proof structure is similar but requires careful handling of boundary cases.

Key Technical Challenges:

• When $|x-y| \geq 2\delta$ (non-overlapping): directly apply reflection technique
• When $|x-y| < 2\delta$ (overlapping): decompose the integration domain into $L_1 = [y-x-\delta, \frac{y-x}{2}]$ and $L_2 = [-\delta, \frac{y-x}{2}]$, utilizing $\chi_{L_2} - \chi_{L_1} \geq 0$ to apply the maximum principle

Boundary Treatment:

• Symmetric configuration ($x = -y$): adjust integration interval to $[-\delta, \frac{D}{2}-y]$
• Asymmetric configuration: translate so the midpoint is at the origin, separating symmetric and asymmetric parts

Counterexample Construction (Section 3.2.1)

Consider the rectangular domain $\Omega = [0,L] \times [-\epsilon, \epsilon]$, where $0 < L \ll \epsilon$, and take the trial function $u(x,y) = y$.

Calculation:

1. Zero mean condition: $\int_\Omega u = 0$
2. Energy estimate: $E_{reg}[u,u] = \int_\Omega \int_\Omega \frac{(y_1-y_2)^2}{|p-q|^{2+2s}}dpdq \leq 4L^2\epsilon^2$
3. $L^2$ norm: $\|u\|_{L^2}^2 = \frac{2L\epsilon^3}{3}$

Conclusion: $\lambda_2 \leq \frac{E_{reg}[u,u]}{\|u\|_{L^2}^2} \leq \frac{CL}{\epsilon} \to 0^+, \quad L \to 0^+$

This shows that as the rectangular domain degenerates to the line segment $[-\epsilon, \epsilon]$, the second eigenvalue approaches 0, contradicting a lower bound depending solely on diameter.

Geometric Insight: This counterexample reveals the sensitivity of non-local operators to domain shape (rather than merely diameter), which is an essential distinction from the classical Payne-Weinberger inequality.

Theorem 4.2 (Nonlinear Case)

Extend to the equation: $u_t = L_{\mathbb{R}^n}u + q(|Du|)$

where $q$ satisfies:

• (i) $q$ is non-decreasing with $q(0) = 0$
• (ii) $q$ is $C^1$ and Lipschitz: $|q'(p)| \leq c$

Key Estimate: At the maximum point, $T_q = q(|Du(y)|) - q(|Du(x)|) \leq c||Du(y)| - |Du(x)||$

Using the spatial derivative condition $Du(x) = \phi' e_1 - ED\Psi(x)$, we obtain: $T_q \leq cK_\nabla E(\Psi(x) + \Psi(y))$

Choosing $C > K + cK_\nabla$ yields a contradiction.

Technical Details and Proof Strategy

Application of Maximum Principle

At the maximum point $(x_0, y_0, t_0)$, we have: $Z_\epsilon(x_0, y_0, t_0) = 0, \quad Z_\epsilon \leq 0 \text{ for } t \leq t_0$

This gives two key conditions:

1. Time Derivative: $\frac{\partial Z_\epsilon}{\partial t} \geq 0$
2. Spatial Maximum: For all $(x,y)$, $Z_\epsilon(x,y,t_0) \leq Z_\epsilon(x_0,y_0,t_0)$

Fine Treatment of Integral Terms

In the inequality $0 \leq \int_L (\rho(r) - \rho(g(r)))(u(y+Rr) - u(x+r))dr + \cdots$

the key observation is:

• Since $r \in L$ implies $g(r) \in R$
• $\rho$ is non-increasing, so $\rho(r) - \rho(g(r)) \geq 0$
• The maximum principle can be safely applied

Clever Use of Symmetry

By choosing coordinates so that $e_1 = \frac{y-x}{|y-x|}$, the reflection transformation has a concise form and preserves the kernel function: $\rho(r) = \rho(|r|) = \rho(|Rr|) = \rho(Rr)$

This is the key symmetry that makes the entire proof possible.

Related Work

Classical Modulus of Continuity Theory

This paper's methods originate from a series of works by Andrews and Clutterbuck:

• AC10: Proof of the fundamental gap conjecture using modulus of continuity methods
• AC09, AC13a, AC13b: Gradient estimates and modulus of continuity for quasilinear parabolic equations

Non-local Operator Theory

• Spectral theory of fractional Laplacian
• Domain fractional operators and their applications
• Research on non-local Neumann boundary conditions

Payne-Weinberstein Inequality

The classical result states that the second eigenvalue of the Neumann Laplacian on convex domains satisfies: $\lambda_2 \geq \frac{\pi^2}{D^2}$ depending only on the diameter $D$. This paper's counterexample shows that the non-local case does not possess this property.

Conclusions and Discussion

Main Conclusions

1. Complete Theory on Full Space: On ℝⁿ, the modulus of continuity preservation for non-local heat equations is completely proved, providing a foundational tool for non-local parabolic equation theory.
2. Success of One-Dimensional Theory: On compact intervals, despite the presence of boundaries, the modulus of continuity property holds, with proof techniques requiring careful analysis of various geometric configurations.
3. Barriers to High-Dimensional Generalization: Through an ingenious counterexample, the paper reveals the essential difficulty of non-local problems in higher dimensions—geometric dependence transcends the single invariant of diameter.
4. Flexibility of the Method: The nonlinear extension demonstrates the adaptability of the technical framework.

Limitations

1. Unresolved High-Dimensional Domain Problem: Conjecture 3.20 (modulus of continuity theorem on high-dimensional bounded convex domains) remains open, with the counterexample suggesting direct generalization is impossible.
2. Restrictions on Kernel Functions: The method relies on rotational symmetry and monotonicity of the kernel; more general kernels may require new techniques.
3. Existence of Regularization Function: Condition (d) requires the existence of $\Psi$ satisfying specific growth conditions, which may not hold for general operators.
4. Optimality Questions: Are the conditions on the modulus of continuity optimal? Do weaker sufficient conditions exist?

Future Directions

1. New Methods for High-Dimensional Domain Problems: Develop theoretical frameworks not relying on the single invariant of diameter, potentially involving more complex geometric quantities such as domain shape and inradius.
2. Handling Non-symmetric Kernels: Explore the possibility of removing the rotational symmetry assumption.
3. Applications to Other Non-local Problems:
   • Non-local mean curvature flow
   • Fractional Yamabe problem
   • Non-local phase transition models
4. Numerical Verification: Conduct numerical simulations of counterexamples and boundary cases for additional insights.
5. Optimal Transport Connections: Explore connections with coupling methods in optimal transport theory.

In-Depth Evaluation

Strengths

1. Strong Technical Innovation:
   • The coupling-by-reflection technique cleverly handles non-local integral terms, providing a new tool for treating long-range interactions
   • The construction of regularization functions is ingenious, balancing control at infinity with boundedness of operator action
   • Boundary treatment in the one-dimensional case demonstrates superior technical mastery
2. Significant Theoretical Contributions:
   • First systematic generalization of the modulus of continuity method to non-local heat equations
   • The counterexample construction provides profound geometric insights, revealing essential features of non-local problems
   • Points new directions for research on non-local Payne-Weinberger inequalities
3. Clear Proof Structure:
   • Unified proof-by-contradiction framework, easy to understand and generalize
   • Clear motivation and technical details for each step
   • Figures (such as Figures 1, 2, 3) aid understanding of geometric configurations
4. Completeness of Results:
   • Covers full space, bounded domains, and nonlinear cases
   • Presents both positive results and counterexamples, showing the complete picture

Weaknesses

1. Absence of High-Dimensional Results:
   • Conjecture 3.20 remains unproved, leaving an important theoretical gap
   • While the counterexample shows direct generalization is impossible, no alternative approach is proposed
2. Technical Nature of Conditions:
   • The modulus of continuity must satisfy a one-dimensional non-local heat equation (condition (c)), which may be difficult to verify in practice
   • The existence condition for regularization functions is strong, limiting applicability
3. Insufficient Application Demonstration:
   • No concrete examples of modulus of continuity functions (e.g., for power-law kernels)
   • Lacks quantitative comparisons with classical results
4. Absence of Numerical Verification:
   • The counterexample has only theoretical analysis, lacking numerical computational support
   • No numerical experiments verifying the practical effectiveness of the theorems

Impact

1. Theoretical Significance:
   • Provides foundational tools for non-local parabolic equation theory
   • Reveals essential differences between non-local and classical problems
   • May inspire research on other non-local problems
2. Methodological Contributions:
   • The coupling-by-reflection technique is applicable to other non-local problems
   • The proof-by-contradiction framework has universal applicability
3. Potential Applications:
   • Regularity theory for non-local diffusion processes
   • Numerical analysis of fractional equations
   • Continuity estimates for stochastic processes
4. Research Direction Guidance:
   • Clearly indicates that high-dimensional problems require new geometric invariants
   • Provides important clues for research on non-local Payne-Weinberger inequalities

Applicable Scenarios

1. Theoretical Research:
   • Regularity theory for non-local parabolic equations
   • Spectral theory of non-local operators
   • Non-local problems in geometric analysis
2. Concrete Applications:
   • Analysis of solution properties for fractional diffusion equations
   • Regularity of non-local phase transition models
   • Theoretical foundations for anomalous diffusion processes
3. Method Generalization:
   • Generalizable to other non-local evolution equations
   • Potentially applicable to non-local elliptic problems
   • Provides theoretical support for numerical methods

Summary of Technical Highlights

1. Geometric Insight of Coupling-by-Reflection: Through reflection symmetry, unify the treatment of two integrals, which is the core of the entire proof.
2. Cleverness of Dimension Reduction: Reduce high-dimensional problems to one dimension while preserving the essential structure.
3. Dual Role of Regularization Term: Both controls behavior at infinity and provides the crucial negative term in the proof by contradiction.
4. Complete Analysis of Boundary Cases: Systematically handle various geometric configurations, demonstrating technical maturity.
5. Wisdom in Counterexample Construction: Simple yet profound, revealing the essential difficulty of the problem.

Key References

• AC10 B. Andrews and J. Clutterbuck, "Proof of the fundamental gap conjecture", 2010
• AC13b B. Andrews and J. Clutterbuck, "Sharp modulus of continuity for parabolic equations on manifolds and lower bounds for the first eigenvalue", Anal. PDE, 2013

Overall Assessment: This is an excellent paper with strong technical content and significant theoretical contributions. The authors successfully generalize the classical modulus of continuity method to non-local heat equations, develop innovative technical tools (particularly the coupling-by-reflection technique), and reveal essential features of non-local problems through ingenious counterexamples. Although the high-dimensional bounded domain case remains incompletely resolved, the paper establishes a solid foundation for research in this direction and clearly indicates future research paths. This work has important value for both non-local partial differential equation theory and geometric analysis.