---
title: 'Lavrentiev''s Gap: Approximation & Regularity'
url: https://www.emergentmind.com/topics/lavrentiev-s-gap
type: topic
---

# Lavrentiev's Gap: Approximation & Regularity

Lavrentiev's gap is the strict inequality between the infimum of a variational integral over a natural “large” admissible class and the infimum over a denser but more regular subclass. In the standard abstract form, if \(\mathcal G:X\to[0,\infty]\) and \(Y\subset X\) is dense in the topology of \(X\), the phenomenon occurs when
\[
\inf_{u\in X}\mathcal G(u)\;<\;\inf_{u\in Y}\mathcal G(u),
\]
equivalently when some finite-energy map in \(X\) cannot be approximated by smoother competitors without increasing the energy. First pointed out by M.A. Lavrentiev in 1927, the gap has become a central obstruction in the calculus of variations, affecting approximation theory, regularity, numerical discretization, and the structure of admissible classes in problems with singular or nonuniform growth [2305.04726].

## 1. Definitions and admissible classes

The precise form of Lavrentiev's gap depends on the ambient variational setting. In one-dimensional scalar problems, a common formulation compares \(W^{1,1}\) or \(AC\) against \(C^1\) or \(W^{1,\infty}\). For the autonomous problem
\[
\min F(y)=\int_0^1L(y(t),y'(t))\,dt,\qquad y(0)=0,\quad y\in W^{1,1}([0,1],\mathbb R),
\]
one sets
\[
I_1=\inf_{y\in W^{1,1},\,y(0)=0}F(y),\qquad
I_2=\inf_{y\in C^1,\,y(0)=0}F(y),
\]
and a strict Lavrentiev gap occurs when \(I_1<I_2\) [2209.03820]. Equivalent formulations appear throughout the literature: \(AC\) versus \(C^1\), \(W^{1,1}\) versus \(W^{1,\infty}\), Sobolev–Orlicz classes versus closures of smooth maps, or regular invertible deformation classes versus their weak closures [1610.03111], [2010.03264], [2603.22873].

The functional-analytic content is always the same. The gap is not merely a failure of topological density, since the smaller class is often dense in the ambient topology. Rather, it is a failure of **energy density**: there is no sequence of smoother admissible maps whose energies converge to the infimum attained or approximated in the larger class. This is why the phenomenon directly obstructs mollification arguments, Euler–Lagrange approximation, and conforming discretizations.

In generalized-growth settings the definition is expressed relative to the natural energy space. For instance, for the double-phase functional
\[
\mathcal F[u]=\int_\Omega \bigl(|\nabla u|^p+a(x)|\nabla u|^q\bigr)\,dx,
\]
one asks whether
\[
\inf_{v\in u_0+W(\Omega)}\mathcal F[v]
\;=\;
\inf_{w\in u_0+C_c^\infty(\Omega)}\mathcal F[w],
\]
where \(W(\Omega)\) is the corresponding Musielak–Orlicz–Sobolev class [2303.05877]. For manifold-valued maps, the same question becomes the modular or strong density of smooth maps in \(W^{1,\varphi}(M,N)\) [2512.18447]. In nonlinear elasticity, the gap can arise between a regular invertible class and its weak \(H^1\)-closure, showing that even the admissible class itself may fail to be weakly closed [2603.22873].

## 2. Canonical one-dimensional examples

The classical prototype is Manià’s problem,
\[
J(v)=\int_0^1 [v'(x)]^6\,[v(x)^3-x]^2\,dx,\qquad v(0)=0,\ v(1)=1.
\]
The singular map \(u(x)=x^{1/3}\) satisfies \(J(u)=0\), hence it realizes the \(W^{1,1}\)-infimum. However, \(u'(x)=\tfrac13 x^{-2/3}\) is unbounded near \(0\), and any Lipschitz competitor has strictly positive energy; a numerical study quoting Ferriero states that for any sequence of Lipschitz trajectories \(y_n\to x^{1/3}\) a.e., one has \(I[y_n]\to\infty\) [1610.03111], [1003.0934]. This example became the standard model for how a singular minimizer can be invisible to regular approximants.

A crucial refinement is that endpoint conditions matter. In Mariconda’s one-dimensional non-autonomous theory, the same Manià integrand yields a true gap for the two-endpoint problem \(y(0)=0\), \(y(1)=1\), but no gap for the one-endpoint problem with only \(y(1)=1\). The truncations
\[
y_h(s)=
\begin{cases}
h^{-1/3},&0\le s\le 1/h,\\
s^{1/3},&s>1/h,
\end{cases}
\]
satisfy \(y_h(1)=1\), converge to \(y(s)=s^{1/3}\) in \(W^{1,1}\), and satisfy \(F(y_h)=0\) for all \(h\) [2201.07222]. This sharply separates the one-endpoint and two-endpoint geometries.

A second instructive one-dimensional example is due to Cerf–Mariconda. Define
\[
L(y,v)=
\begin{cases}
\bigl(v^2-\frac1{4y^2}\bigr)^2v^2,& y\neq 0,\\[2mm]
1,& y=0,
\end{cases}
\]
and let \(y_*(t)=\sqrt t\). Then \(F(y_*)=0\), so \(y_*\) is a \(W^{1,1}\)-minimizer with \(I_1=0\), while any Lipschitz \(y\) with \(y(0)=0\) satisfies \(F(y)\ge 1\). Hence \(I_1=0<I_2=1\) [2209.03820]. The mechanism is explicit: along the singular minimizer the term \(v^2-\frac1{4y^2}\) vanishes exactly, while regular competitors cannot reproduce that singular profile near \(t=0\).

These examples established two enduring facts. First, polynomial or even smooth dependence on the gradient does not preclude a gap. Second, the phenomenon is often driven by a singular minimizer whose derivative blows up on a small set while the integrand is engineered to vanish exactly along that singular trajectory.

## 3. One-dimensional non-occurrence theorems

One-dimensional scalar problems also admit some of the sharpest positive results on the **absence** of the Lavrentiev phenomenon. For autonomous functionals with one endpoint condition,
\[
\min F(y)=\int_0^1L(y(t),y'(t))\,dt,\qquad y(0)=0,
\]
Cerf–Mariconda introduced **Condition (R)**: there exist locally Lipschitz functions \(\rho^-,\rho^+:\mathbb R\to\mathbb R\) such that \(\rho^-(z)<0<\rho^+(z)\) for all \(z\), and for every bounded interval \(J\subset\mathbb R\),
\[
\sup_{z\in J}L(z,\rho^-(z))<+\infty,\qquad
\sup_{z\in J}L(z,\rho^+(z))<+\infty.
\]
If \(y\in W^{1,p}([0,1])\), \(F(y)<+\infty\), and \(L\) satisfies Condition (R) on \(y([0,1])\), then there exists a sequence of Lipschitz functions \(y_k\) with \(y_k(0)=0\), \(y_k\to y\) in \(W^{1,p}\), and \(F(y_k)\to F(y)\). In particular, no gap occurs at \(y\); if Condition (R) holds on every bounded interval, then
\[
\inf_{W^{1,1},\,y(0)=0}F(y)
=
\inf_{\operatorname{Lip},\,y(0)=0}F(y).
\]
Condition (R) is strictly weaker than the classical rectangle-boundedness condition of Alberti–Serra Cassano, because it requires boundedness only along two graphs \(v=\rho^\pm(z)\), not on full state-velocity rectangles [2209.03820].

Mariconda’s non-autonomous theory reaches a related conclusion under a different structural package. Writing
\[
L(s,y,v)=\Lambda(s,y,v)\Psi(s,y),
\]
the hypotheses include radial convexity in \(v\), a time-Lipschitz condition (S) on \(\Lambda\), boundedness of \(\Psi\) along the image of \(y\), continuity of \(s\mapsto \Psi(s,z)\), and local boundedness of \(\Lambda\) near the graph of \(y\). Under these assumptions, one obtains a constructive sequence of Lipschitz reparametrizations \(y_h\) with the same one-endpoint data, \(y_h\to y\) in \(W^{1,p}\), and \(F(y_h)\to F(y)\); for two endpoints one needs additional hypotheses such as positivity of \(\Psi\) and either real-valuedness of \(\Lambda\) or blow-up of \(\Lambda\) near \(\partial\operatorname{Dom}(\Lambda)\) [2201.07222].

The underlying method is reparametrization rather than direct smoothing. High-derivative sets are slowed down, and compensating intervals are inserted so that the final map remains Lipschitz and preserves the endpoint constraints. In the autonomous one-endpoint case, Cerf–Mariconda refine this by replacing problematic slopes on bad intervals with the bounded velocities \(\rho^\pm\), then reconnecting the endpoint through short differential-equation arcs \(z'=\rho^\pm\) [2209.03820]. These constructions show that, at least in one dimension, the absence of a gap can often be reduced to a geometric control of the integrand along selected velocity profiles.

## 4. Nonuniform growth, sharp thresholds, and density scales

In higher dimensions and nonstandard-growth settings, Lavrentiev’s gap is tightly linked to the density of smooth maps in Musielak–Orlicz-type energy spaces. For double-phase energies
\[
\mathcal F[u]=\int_\Omega\bigl(|\nabla u|^p+a(x)|\nabla u|^q\bigr)\,dx,
\]
the classical criterion assumes \(a\in C^{0,\alpha}(\Omega)\), \(0\le a\le \Lambda\), and \(p<q\le p+\alpha\), which guarantees density of \(C_c^\infty(\Omega)\) in the natural space and hence absence of the gap. A sharp extension is the scale \(SP^\gamma(\Omega)\), defined by
\[
a\in SP^\gamma(\Omega)
\iff
\exists C>0\ \forall x,y\in\Omega:\ a(x)\le C\bigl(a(y)+|x-y|^\gamma\bigr).
\]
If \(a\in SP^\gamma(\Omega)\) and \(q\le p+\gamma\), then no Lavrentiev gap arises; if competitors are \(C^{0,\theta}\), the weaker bound \(q\le p+\gamma/(1-\theta)\) suffices. Conversely, under \(1<p<n<n+\gamma<q\), there exist \(a\in SP^\gamma(B_1)\) and boundary data \(u_0\) such that the gap occurs [2303.05877].

The smoothness of the weight can be pushed further. For \(a\in C^{k,\alpha}(\overline\Omega)\), the sharp threshold becomes
\[
q\le p+(k+\alpha)\max(1,p/N).
\]
If this holds, then
\[
\inf_{u\in u_0+W_0^{1,1}(\Omega)}I[u]
=
\inf_{u\in u_0+C_c^\infty(\Omega)}I[u],
\]
while if
\[
q>p+(k+\alpha)\max(1,p/N),
\]
there exist \(\Omega\), \(a\in C^{k,\alpha}(\overline\Omega)\), and \(u_0\) with a genuine gap. In particular, if \(a\in C^\infty(\overline\Omega)\), no additional restrictions on \(p\) and \(q\) are required [2509.06567].

Generalized Orlicz and borderline double-phase models exhibit a complementary dichotomy. In the two-dimensional log-weighted model
\[
\Phi_{p,\alpha,\beta}(x,t)
=
\frac1p\,t^p\log^{-\beta}(e+t)
+
a(x)\frac1p\,t^p\log^\alpha(e+t),
\]
with checkerboard coefficient \(a(x)\), one has \(H^{1,\Phi}=W^{1,\Phi}\) and hence no gap if \(\min\{\alpha,\beta\}\le 1\), whereas \(H^{1,\Phi}\subsetneq W^{1,\Phi}\) and a genuine gap occurs if \(\alpha>1\) and \(\beta>1\) [2010.03264]. In anisotropic Musielak–Orlicz settings, Borowski–Chlebicka–Miasojedow prove absence of the gap under a local balance condition
\[
M_B^+(\xi)\le M_B^-(C\xi)+1
\quad\text{for }|\xi|\le c\,r^{\gamma-1},
\]
recovering in particular the sharp double-phase bound \(q\le p+\alpha\) in the isotropic \(C^{0,\alpha}\)-case and its coordinatewise analogues in multi-phase anisotropic models [2210.15217].

Parabolic double-phase problems admit an analogous density theory. If \(a\in C^{\alpha,\alpha/2}(\Omega_T)\) and
\[
q\le p+\frac{p\alpha}{n+2},
\]
then every finite-energy map in the natural parabolic class admits smooth approximants with convergence in \(L^p(I;W^{1,p}(B))\) and convergence of the stationary energy \(\mathcal P\). The threshold improves to
\[
q\le p+\max\Bigl\{\alpha,\frac{p\alpha}{n+2}\Bigr\}
\]
for bounded solutions, and to
\[
q\le p+\max\Bigl\{\frac{s\alpha}{n+s},\frac{p\alpha}{n+2}\Bigr\}
\]
under \(C(0,T;L^s)\)-regularity [2603.14235].

A broader scalar theorem shows that even in non-autonomous, non-convex problems, the gap can be discarded under an anti-jump balance condition \((H^{\rm conv})\), while dropping boundedness or convexity in the second variable and any \(\Delta_2\)-type assumption in the last variable [2410.14995]. This suggests that the decisive issue is not regularity of the integrand in isolation, but a quantitative balance between its spatial oscillation and its growth with respect to the gradient.

## 5. Singular sets, fractal barriers, and topological obstructions

A recurring misconception is that the absence of a Lavrentiev gap forces the singular set of a minimizer to be small. Gratwick disproved this sharply: given any closed Lebesgue-null set \(E\subset[a,b]\) and an arbitrary smooth strictly convex superlinearity \(w\), there exist a smooth strictly convex Lagrangian \(L\in C^\infty(\mathbb R^3)\) and a unique minimizer \(u\in AC(a,b)\) whose singular set
\[
\Sigma(u)=\{x\in[a,b]:u'(x)=\pm\infty\}
\]
coincides exactly with \(E\), while still admitting smooth competitors \(u_k\in C^\infty([a,b])\) with \(u_k\to u\) uniformly and \(I[u_k]\to I[u]\). Thus non-occurrence of the Lavrentiev phenomenon does not imply that the singular set is small [1503.06694].

Another common belief is that the gap is tied to an exponent crossing the ambient dimension. Fractal constructions show this is not generally correct. Balci–Diening–Surnachev replace a point singularity by a fractal contact set \(S\) of prescribed Hausdorff dimension and construct gaps for variable-exponent, double-phase, and weighted \(p\)-energies, thereby showing that the dimensional threshold is not essential in the classical sense [1906.04639]. In a related direction, variational problems with differential forms admit a unifying “separating pair” construction \((u,A)\) that violates Stokes’ theorem on a negligible set and produces non-density in double-phase, borderline double-phase, and variable-exponent models [2305.04726].

For manifold-valued maps, the issue is entangled with topology. Under growth conditions expressed by integral criteria on \(\varphi^-_M\), or under \(k\)-connectedness of the target manifold \(N\) with upper-growth exponent \(\gamma\le k+1\), smooth maps are strongly or weakly dense in \(W^{1,\varphi}(M,N)\), and hence no Lavrentiev gap occurs. But when the local comparability condition fails, there are vectorial double-phase counterexamples with smooth boundary data and a genuine gap between \(W^{1,\varphi}_{u_0}\) and smooth \(N\)-valued competitors [2512.18447]. The scalar density problem and the topological obstruction problem are therefore distinct but structurally compatible facets of the same phenomenon.

## 6. Numerical analysis and discrete approximation

Lavrentiev’s gap is a decisive obstruction for conforming numerical schemes because such schemes typically minimize over spaces contained in \(W^{1,\infty}\) or in smooth finite-dimensional trial classes. When a gap is present, the discrete minimum converges, at best, to the infimum over the regular subclass rather than to the true Sobolev infimum. This failure is already visible in Manià’s problem: piecewise-linear approximations and standard finite-difference discretizations drive the discrete integral away from zero as the mesh is refined, matching the theoretical obstruction [1003.0934].

A direct remedy is the enhanced finite element method of Schnake–Feng. On a quasi-uniform mesh \(T_h\), one defines the componentwise cut-off
\[
[\chi_h^\alpha(s)]_i=
\begin{cases}
s_i,& |s_i|\le h^{-\alpha},\\
\operatorname{sign}(s_i)h^{-\alpha},& |s_i|>h^{-\alpha},
\end{cases}
\]
and minimizes the modified energy
\[
E_h^{\rm enh}(v_h)=\int_\Omega f\bigl(x,v_h(x),\chi_h^\alpha(\nabla v_h(x))\bigr)\,dx.
\]
For Manià’s example, where the gradient growth exponent is \(p=6\), any \(\alpha<1/6\) works; in numerical experiments, the enhanced method approaches \(x^{1/3}\) well, while the standard FEM diverges [1610.03111].

The one-dimensional theory for Manià’s problem was later placed on a \(\Gamma\)-convergence foundation. Feng–Siktar consider the cut-off functional
\[
J_h(v)=\int_0^1\bigl[X_{h,\alpha}(v'(x))\bigr]^6\,(v(x)^3-x)^2\,dx
\]
on the piecewise affine space \(X_h\), where
\[
X_{h,\alpha}(t)=\operatorname{sgn}(t)\min\{|t|,h^{-\alpha}\}.
\]
If
\[
\alpha<\min\Bigl\{\frac16+\frac1{6p'},\frac1{5p'}\Bigr\},
\qquad p'=\frac p{p-1},
\]
then \(J_h\) \(T\)-converges to \(J\) with respect to the strong \(W^{1,p}(0,1)\)-topology, and the discrete minimizers converge strongly in \(W^{1,p}(0,1)\) to the unique singular minimizer \(u_*(x)=x^{1/3}\). The proof uses a “layered” interpolant \(K_h\) with strong \(W^{1,p}\)-stability and approximation properties [2410.06434].

A different strategy is nonconformity. For convex non-autonomous integrands with non-standard growth, the Crouzeix–Raviart finite element scheme minimizes over a broken affine space \(V_h^{CR}\) with piecewise gradient \(\nabla_h\) and meshwise quadrature
\[
\mathcal F_h(v_h)=\sum_{T\in\mathcal T_h}\int_T[\phi(x_T,\nabla v_h)-f\cdot v_h]\,dx.
\]
Under either quadrature assumptions \((A1)\)–\((A3)\) or the simpler condition \((B)\), one has
\[
\lim_{h\to0}\mathcal F_h(u_h)=\mathcal F(u)=\min_W\mathcal F(v),
\]
together with strong \(L^{p_-}\)-convergence of a subsequence and weak convergence of broken gradients. In the piecewise-constant exponent benchmark, conforming FEM converges to \(\inf_H\mathcal F\), while the Crouzeix–Raviart scheme converges to \(\min_W\mathcal F\) [2106.06837]. The numerical message is unambiguous: when the gap is genuine, conformity can be a liability rather than an advantage.

## 7. Nonlinear elasticity and current directions

The scope of Lavrentiev’s gap now extends beyond scalar and generalized-growth models into three-dimensional nonlinear elasticity. For the neo-Hookean energy
\[
E[u]=\int_\Omega W(Du(x))\,dx,\qquad
W(F)=|F|^2+H(|\det F|),
\]
with \(H\) convex, \(\lim_{t\to\infty}H(t)/t=+\infty\), and \(\lim_{s\to0^+}H(s)=+\infty\), a recent result constructs boundary data for which the infimum over the regular invertible class \(A^r\) is strictly larger than the infimum over its sequential weak \(H^1\)-closure \(\overline{A^r}\). More precisely, for any \(\lambda\in(0,2\pi)\), there exists \(\delta\in(0,1]\) such that
\[
\inf_{u\in A^r_\delta}E[u]
\;\ge\;
\inf_{u\in \overline{A^r_\delta}}E[u]+\lambda.
\]
The mechanism is a dipole-type singularity: a deformation creates self-contact across a sphere \(\Gamma\), and the inverse deformation develops a jump whose singular gradient has total variation \(2\pi\). The gap therefore reflects inverse-regularity failure rather than classical cavitation [2603.22873].

This development repositions the phenomenon in a broader landscape. In scalar theory, modern results increasingly identify structural “anti-jump” or “balance” conditions under which the gap disappears [2410.14995]. In nonuniform ellipticity, the sharpness of Hölder, \(SP^\gamma\), \(C^{k,\alpha}\), and anisotropic balance thresholds has become a central theme [2303.05877], [2509.06567], [2210.15217]. In geometry and topology, density results for manifold-valued maps coexist with counterexamples based on fractal barriers and differential forms [2512.18447], [2305.04726].

A plausible implication is that Lavrentiev’s gap is best understood not as an isolated pathology, but as a precise indicator of mismatch between the natural energy space and the approximation class imposed by analysis, topology, or computation. In one dimension this mismatch can sometimes be neutralized by reparametrization; in double-phase and Musielak–Orlicz problems it is governed by sharp continuity scales; in numerical analysis it dictates whether conforming approximations are intrinsically blind; and in nonlinear elasticity it can encode physically meaningful singular mechanisms that are invisible to regular invertible competitors.

Source: https://www.emergentmind.com/topics/lavrentiev-s-gap