Free Boundary & Obstacle Problems

• Free boundary/obstacle problems are mathematical models where solutions satisfy PDEs with additional constraints, leading to unknown interfaces that separate contact and non-contact regions.
• They are applied in finance, phase transitions, and materials science, using variational methods and nonlocal operators to model complex phenomena.
• Advanced techniques such as blow-up analysis, boundary Harnack principles, and geometric measure theory are key to classifying regular versus singular free boundaries.

Free boundary and obstacle problems are central to the analysis of partial differential equations (PDEs) in both classical and modern frameworks, featuring unknown interfaces or domains to be determined as part of the solution. These problems arise in diverse contexts such as materials science, finance (e.g., American options), phase transitions, and inverse problems. The rigorous mathematical theory requires detailed analysis of the regularity, structure, and classification of both solutions and their free boundaries, with tools spanning variational methods, viscosity solution theory, blow-up and monotonicity analysis, and geometric measure theory.

1. Formulations and Classifications

Classical and Variational Obstacle Problems

The canonical obstacle problem seeks $u: \Omega \to \mathbb{R}$ above a prescribed obstacle $\psi$, minimizing the energy

$$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$

in the convex set $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$, with $f$ typically in $L^2(\Omega)$. The solution $u$ satisfies the variational inequality and complementary conditions: $$\begin{aligned} &u \geq \psi, \quad (\Delta u - f)(u - \psi) = 0 \ \text{a.e.}, \ & \Delta u \leq f \ \text{in } \Omega, \quad \Delta u = f \ \text{in } \{u > \psi\}. \end{aligned}$$ The free boundary $\Gamma = \partial\{u > \psi\}$ splits the domain into contact and non-contact regions.

Generalizations

Obstacle-type problems extend in many directions:

• Integro-differential/Nonlocal: Replace $\Delta$ with operators of the form

$\psi$0

where $\psi$1 satisfies structural ellipticity and symmetry conditions (Caffarelli et al., 2016).

• Fully Nonlinear and Degenerate: Equations like $\psi$2 for convex, uniformly elliptic $\psi$3, and degeneracy parameter $\psi$4 (Silva et al., 2019).
• Non-stationary/Parabolic: Obstacle problems for $\psi$5, modeling time-evolving phase boundaries and financial instruments (Barrios et al., 2016, Kukuljan, 2022).
• Thin (Signorini) and Higher Order: Boundary or codimension-one obstacle constraints, including those for the bi-Laplacian (Danielli et al., 14 Sep 2025), or with variable/anisotropic structure (Andreucci et al., 2024).

2. Free Boundary Regularity and Structure

Dichotomy and Blow-Up Analysis

At a free boundary point, blow-up techniques are crucial for classification:

Context	Blow-up Types	Regular Points	Singular Points
Classical obstacle ($\psi$6)	Quadratics	$\psi$7	$\psi$8 ($\psi$9, $$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$0) (Figalli et al., 2017)
Fractional Laplacian ($$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$1)	$$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$2, quadratics	Homogeneity $$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$3	Homogeneity $$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$4 (Barrios et al., 2015, Caffarelli et al., 2016)
Parabolic, $$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$5	$$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$6	Nondegenerate expansion	Flat (degenerate) behavior (Barrios et al., 2016)
Degenerate fully nonlinear/gradient	$$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$7 separation	$$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$8 at FB	Porosity, lower regularity (Silva et al., 2019)

Regular points admit nontrivial, low-order homogeneous blow-ups (half-space or sectoral), and the free boundary locally aligns with the zero set of these profiles. Singular points correspond to higher-order homogeneous solutions with lower-dimensional structure in the free boundary.

Regularity and Stratification

The celebrated Caffarelli dichotomy shows that the regular set is relatively open in the free boundary, forming a $$J[u] = \int_\Omega \frac{1}{2} |\nabla u|^2 - f u \, dx$$9 (or $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$0, depending on context) hypersurface (Figalli et al., 2017, Figalli, 2018). Singular points are contained in a union of lower-dimensional $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$1 manifolds, with improved $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$2 or $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$3 structure in specific cases and dimensions. Finer stratification accounts for the dimension of the kernel of the blow-up polynomial's Hessian: $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$4 and each $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$5 is locally $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$6-rectifiable (Figalli et al., 2017). In nonlocal and thin obstacle problems, analogous stratifications hold—with the regular/degenerate dichotomy governed by the homogeneity and growth of the blow-up profiles (Barrios et al., 2015, Danielli et al., 14 Sep 2025, Andreucci et al., 2024).

3. Advanced Operator Classes and Nonlocal Problems

Integro-Differential and Fully Nonlinear Contexts

For fractional and general integro-differential operators of order $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$7 with ellipticity and symmetry (e.g., $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$8), purely nonlocal blow-up and boundary Harnack methods yield:

• Near regular points: $$K_\psi = \{ v \in H^1(\Omega) : v|_{\partial\Omega} = g,\, v \geq \psi \ \text{a.e.}\}$$9, with $f$0 and the free boundary $f$1 (Caffarelli et al., 2016).
• For fully nonlinear (convex) nonlocal operators, similar expansions and regularity hold provided interior $f$2 and boundary $f$3 regularity is available (Caffarelli et al., 2016).

Key innovations are the avoidance of monotonicity formulas and reduction to local PDEs; instead, boundary Harnack, blow-up classification, and nonlocal barrier arguments are fundamental.

Degenerate Cases and Singular Forcing

For operators of the form $f$4, if $f$5, obstacle constraints can induce higher regularity at free boundary points than in the unconstrained equation. Solutions are locally $f$6 or fully $f$7 at $f$8, exhibiting quadratic separation from the obstacle, even when interior regularity may be limited (Silva et al., 2019).

For degenerate right-hand sides (e.g., $f$9 in the Laplacian obstacle), the blow-ups are cubic rather than quadratic, and only Lipschitz regularity (with corners) for the free boundary may hold at the degenerate point, in contrast to full $L^2(\Omega)$0 regularity in the classical case (Liu, 11 Feb 2026).

4. Methodologies and Analytical Tools

Technique	Key Applications
Blow-up/Scaling Analysis	Classify free boundary points, derive expansions, stratify
Monotonicity Formulas	Weiss, Almgren, Monneau: ensure uniqueness, rate of decay
Boundary Harnack Principle	Proportional derivatives, regularity of level-sets (Caffarelli et al., 2016)
Comparison Principles	Nondegeneracy estimates, barrier construction
Penalization/Viscosity	Existence and uniqueness for fully nonlinear/degenerate PDE

Notably, in sign-changing, nonhomogeneous, or mixed operator settings, nonstandard monotonicity and geometric compactness arguments are often substituted for classical variation techniques (Salo et al., 27 Jun 2025, Araújo et al., 8 Jan 2025).

5. Geometric Features, Rectifiability, and Applications

Rectifiability and Geometric Measure

Recent advances establish $L^2(\Omega)$1-rectifiability and finite Minkowski content for free boundaries under minimal regularity assumptions—even with merely Lipschitz or Sobolev coefficients (Andreucci et al., 2024). The stratified singular sets described above connect with deep concepts in geometric measure theory (De Giorgi structure theorem, BV functions).

Application Paradigms

• Mathematical Finance: Obstacle problems describe American option pricing, both in local and nonlocal (jump) diffusion models, tying regularity of the free boundary to the early exercise region (Barrios et al., 2016).
• Inverse Problems: In scattering, free-boundary methods expose geometric rigidity (C$L^2(\Omega)$2 regularity) of non-scattering obstacles under harmonic data, crucial for uniqueness and stability in inverse scattering (Salo et al., 27 Jun 2025).
• Biology/Hele-Shaw-type Models: Propagation with obstacles and interacting fronts are governed by geometric PDEs with coupled obstacle constraints, exhibiting regularity up to the interface (Bevilacqua et al., 3 Jul 2025, Araújo et al., 8 Jan 2025).
• Boundary and Thin Obstacle: New multiscale analytic tools (e.g., boundary Harnack for C$L^2(\Omega)$3 moving domains) enable iterative improvements of regularity for parabolic nonlocal problems (Kukuljan, 2022) and higher-order bi-Laplacian thin obstacle problems (Danielli et al., 14 Sep 2025).

6. Open Directions and Unresolved Issues

Major active directions include:

• Fine regularity and singularity classification in fully nonlinear, degenerate, or nonlocal settings beyond translation-invariant operators.
• Parabolic analogues for subcritical regimes and with drift or lower-order perturbations.
• Free boundary regularity for anisotropic, heterogeneous, and quasilinear obstacle problems; quantitative homogenization under random media (Aleksanyan et al., 2021).
• Analyticity, higher smoothness, or structural stratification of singular sets under minimal assumptions (Figalli et al., 2017, Andreucci et al., 2024).
• New phenomena at boundaries and in "no-sign" problems with sign-changing data, and in problems with mixed or coupled operators (Araújo et al., 8 Jan 2025, Salo et al., 27 Jun 2025).

7. Summary Table: Key Results by Problem Type

Problem Type	Regularity at Regular Points	Free Boundary Regularity	Singular Set Structure	Notes
$L^2(\Omega)$4, classical	$L^2(\Omega)$5	Analytic (C$L^2(\Omega)$6)	Stratified by C$L^2(\Omega)$7 manifolds, lower dimension	(Figalli et al., 2017, Figalli, 2018)
$L^2(\Omega)$8 (fractional)	$L^2(\Omega)$9	$u$0 for $u$1	Stratified C$u$2 or C$u$3 manifolds	(Caffarelli et al., 2016, Barrios et al., 2015)
Parabolic nonlocal ($u$4)	$u$5	$u$6, upgrades to $u$7	As above, for space-time slices	(Kukuljan, 2022, Barrios et al., 2016)
Fully nonlinear, degenerate	$u$8, or $u$9 if $$\begin{aligned} &u \geq \psi, \quad (\Delta u - f)(u - \psi) = 0 \ \text{a.e.}, \ & \Delta u \leq f \ \text{in } \Omega, \quad \Delta u = f \ \text{in } \{u > \psi\}. \end{aligned}$$0	Free boundary $$\begin{aligned} &u \geq \psi, \quad (\Delta u - f)(u - \psi) = 0 \ \text{a.e.}, \ & \Delta u \leq f \ \text{in } \Omega, \quad \Delta u = f \ \text{in } \{u > \psi\}. \end{aligned}$$1 at regular points	Porous singular part, stratified	(Silva et al., 2019)
Thin obstacle with variable coefficients	$$\begin{aligned} &u \geq \psi, \quad (\Delta u - f)(u - \psi) = 0 \ \text{a.e.}, \ & \Delta u \leq f \ \text{in } \Omega, \quad \Delta u = f \ \text{in } \{u > \psi\}. \end{aligned}$$2, stratified by contact order	$$\begin{aligned} &u \geq \psi, \quad (\Delta u - f)(u - \psi) = 0 \ \text{a.e.}, \ & \Delta u \leq f \ \text{in } \Omega, \quad \Delta u = f \ \text{in } \{u > \psi\}. \end{aligned}$$3 manifold at reg. pts.	Hausdorff dimension stratification	(Andreucci et al., 2024, Danielli et al., 14 Sep 2025)
Constraint maps (vectorial obstacles)	$$\begin{aligned} &u \geq \psi, \quad (\Delta u - f)(u - \psi) = 0 \ \text{a.e.}, \ & \Delta u \leq f \ \text{in } \Omega, \quad \Delta u = f \ \text{in } \{u > \psi\}. \end{aligned}$$4, $$\begin{aligned} &u \geq \psi, \quad (\Delta u - f)(u - \psi) = 0 \ \text{a.e.}, \ & \Delta u \leq f \ \text{in } \Omega, \quad \Delta u = f \ \text{in } \{u > \psi\}. \end{aligned}$$5 locally	$$\begin{aligned} &u \geq \psi, \quad (\Delta u - f)(u - \psi) = 0 \ \text{a.e.}, \ & \Delta u \leq f \ \text{in } \Omega, \quad \Delta u = f \ \text{in } \{u > \psi\}. \end{aligned}$$6 near regular points	Discrete/stratified singular points (2D/ND)	(Figalli et al., 2023)

---

The modern theory of free boundary and obstacle problems thus integrates fine analytic tools, expansive operator classes (including nonlocality, nonlinearity, degeneracy), and geometric measure analysis to yield a rich, nuanced regularity and stratification framework. The interplay of blow-up techniques, monotonicity, and geometric compactness underpins both fundamental results and ongoing research directions.