[Paper Review] Regularity theory for nonlinear integral operators

Luis Caffarelli, Chi Hin Chan+1 more|arXiv (Cornell University)|Mar 8, 2010

Numerical methods in inverse problems17 references16 citations

TL;DR

This paper establishes Hölder regularity for solutions to a class of nonlocal parabolic equations arising in nonlinear integral variational problems, using a De Giorgi-type iteration method adapted to nonlocal operators. It proves that solutions with initial data in $L^2$ become instantaneously bounded and Hölder continuous in space-time for any positive time, with the regularity quantified by a Hölder exponent depending on the kernel parameters and the nonlinearity's convexity structure.

ABSTRACT

This article is dedicated to the proof of the existence of classical solutions for a class of non-linear integral variational problems. Those problems are involved in nonlocal image and signal processing.

Motivation & Objective

• To develop a regularity theory for nonlocal evolution equations of variational type with measurable kernels.
• To establish the existence and Hölder continuity of classical solutions to a class of nonlinear integral equations arising in image and signal processing.
• To extend De Giorgi's method to nonlocal operators by analyzing the Hölder regularity of first derivatives of solutions.
• To prove that solutions to the time-dependent nonlocal equation become instantaneously regular in space and time for any positive time.
• To provide a rigorous foundation for the regularity of solutions to nonlocal variational problems with convex, $C^2$ nonlinearities and singular kernels.

Proposed method

• Adapts the De Giorgi-Nash-Moser iteration scheme to nonlocal operators by studying the Hölder regularity of spatial derivatives of solutions.
• Transforms the original equation into a nonlocal evolution equation for the derivative $w = D_e \theta$, with a new kernel $K(t,x,y) = \phi''(\theta(y)-\theta(x))K(y-x)$, which satisfies the same structural conditions as the original kernel.
• Uses difference quotients $D_e^h \theta$ as approximations to derivatives to make the iteration argument rigorous in the absence of classical differentiability.
• Applies a sequence of rescaling and normalization steps to iteratively improve oscillation bounds on the solution, relying on measure-theoretic estimates on super-level sets.
• Employs a weighted cutoff function $\psi_{\varepsilon,\lambda}$ to control the growth of the solution and derive decay estimates in oscillation over nested cylinders.
• Establishes the final Hölder continuity by constructing a sequence of rescaled functions and applying a refined oscillation lemma to obtain exponential decay of oscillation in dyadic cylinders.

Experimental results

Research questions

• RQ1 Can the De Giorgi method be extended to nonlocal equations of the form $\partial_t \theta = \int \phi'(\theta(y)-\theta(x))K(y-x)\,dy$ with singular kernels?
• RQ2 What regularity properties, such as Hölder continuity, can be established for solutions to such nonlocal variational equations?
• RQ3 Does the solution become instantaneously regular in space and time, even if the initial data is only in $L^2$?
• RQ4 How does the Hölder exponent of the solution depend on the kernel's singularity strength $s$ and the convexity of $\phi$?
• RQ5 Can the nonlocal analog of the De Giorgi oscillation lemma be formulated and proven using measure-theoretic and iterative techniques?

Key findings

• Solutions to the nonlocal parabolic equation (2.3) with initial data in $H^1(\mathbb{R}^N)$ exist globally and are classical in the $L^2$ sense.
• For any $t_0 > 0$, the spatial gradient $\nabla_x \theta$ is Hölder continuous on $(t_0, \infty) \times \mathbb{R}^N$, with the regularity independent of the initial data's smoothness beyond $H^1$.
• The Hölder exponent $\alpha$ of the gradient is positive and depends only on the parameters $s$, $\Lambda$, and the structure of $\phi''$, with $\alpha = \frac{\ln(1 - \lambda^*/4)}{\ln(K^s)}$ for some $K < 1$.
• The proof relies on an iterative oscillation decay argument applied to the derivative $w = D_e \theta$, showing that oscillation decays geometrically over dyadic cylinders.
• The method establishes that solutions become bounded and Hölder continuous in any compact subset of $\mathbb{R}^N \times (0,\infty)$, even if the initial data is only in $L^2$.
• The result confirms the instantaneous smoothing effect of the nonlocal operator, analogous to the parabolic regularity theory in the local case.