Pith. sign in

REVIEW 3 major objections 4 minor 7 references

Asymptotics of nonlinear Robin energies

T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper derives the first-order asymptotic expansions of nonlinear Robin energies as the boundary-penalty parameter tends to zero and to infinity, showing exactly how the energy approaches or diverges from the Neumann and Dirichlet…

desk verdict Nonlinear Robin asymptotics are genuinely new and the proofs are solid; the Neumann constant has a sign-of-exponent typo and the Dirichlet expansion is over-advertised without its regularity condition. read the letter →

arxiv 2506.06914 v1 pith:XVFBPOP4 submitted 2025-06-07 math.AP math.SP

Giuseppe Buttazzo , Roberto Ognibene This is my paper
classification math.APmath.SP MSC 35J2035J2535J6635J9249J45
keywords Robinboundaryconditionp-LaplacianasymptoticanalysisnonlinearellipticPDEsDirichletlimitNeumannvariationalmethodsenergyexpansion
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper proves first-order asymptotic expansions for a family of nonlinear Robin energies on a bounded Lipschitz domain as the boundary-penalty parameter $\alpha$ tends to $0^+$ and $+\infty$. In the Dirichlet limit, the minimum equals the Dirichlet energy plus a correction of order $\alpha^{-1/(q-1)}$ whose coefficient is an explicit boundary integral of a power of the normal derivative of the limiting $p$-harmonic solution. In the Neumann limit, if the source has zero mean, the energy approaches the Neumann energy linearly in $\alpha$; if the source has nonzero mean, the energy diverges to $-\infty$ at rate $\alpha^{-1/(q-1)}$, with a coefficient set by the mean of $f$ and the boundary measure. A sympathetic reader should care because these results turn the intuitive interpolation between Neumann and Dirichlet regimes into quantitative statements, obtained by elementary variational estimates rather than spectral theory.

What carries the argument

The central object is the energy functional $E_\alpha(u) = \frac1p\int_\Omega|\nabla u|^p dx + \frac{\alpha}{q}\int_{\partial\Omega}|u|^q dS - \int_\Omega f u dx$, whose unique minimizer $u_\alpha$ satisfies the weak nonlinear Robin condition $|\nabla u_\alpha|^{p-2}\partial_\nu u_\alpha + \alpha|u_\alpha|^{q-2}u_\alpha = 0$ on the boundary. The mechanism is a rescaling $u = u_\infty + \alpha^{-\gamma}v$ near the Dirichlet limit, followed by convexity to drop the nonnegative gradient excess and isolate a boundary minimization; the scalar inequality $\min_{t\in\mathbb{R}}(\frac{a}{q}|t|^q - bt) = -\frac{q-1}{q}|b|^{q/(q-1)}a^{-1/(q-1)}$ produces the explicit constants and the exponent $\gamma=1/(q-1)$ in both main expansions.

What would settle it

On the unit ball $\mathbb{B}_d$ in $\mathbb{R}^d$ with $p=q=2$ and $f\equiv 1$, the Dirichlet solution is $u_\infty=(1-|x|^2)/(2d)$ and $\partial_\nu u_\infty=-1/d$. The claimed expansion predicts $\alpha(E_\alpha-E_\infty)\to -\frac{1}{2}\int_{\partial\mathbb{B}_d}(\partial_\nu u_\infty)^2 dS = -\frac{\mathrm{Vol}(\mathbb{B}_d)}{2d}$. Because the linear Robin problem is exactly solvable by radial functions, the minimizer being $u_\alpha=u_\infty+(\alpha d)^{-1}$, one can compute $E_\alpha-E_\infty$ directly and verify this constant; any other limit would falsify Theorem 3.1.

Watch

Extended reading notes

Core claim

On a bounded Lipschitz domain $\Omega\subset\mathbb{R}^d$, the paper studies $E_\alpha = \inf\{ \frac1p\int_\Omega|\nabla u|^p dx + \frac{\alpha}{q}\int_{\partial\Omega}|u|^q dS - \int_\Omega f u dx : u\in W^{1,p}(\Omega)\}$ and proves, with $\gamma=1/(q-1)$, that as $\alpha\to+\infty$, $$E_\$\alpha$ = E_\infty - \frac{q-1}{q} \$alpha^{{-\gamma}}$\int_{\partial\$\Omega$}|\partial_\nu u_\infty|^{\frac{(p-1)q}{q-1}} dS + o(\$alpha^{{-\gamma}}$),$$ provided the normal derivative of the Dirichlet solution $u_\infty$ has the needed trace regularity. As $\alpha\to 0^+$, the behavior splits: when $\int_\Omega f dx=0$, $E_\alpha = E_0 + \frac{\alpha}{q}\int_{\partial\Omega}|u_0|^q dS + o(\alpha)$; when $\int_\Omega f dx\neq 0$, $E_\alpha = -\alpha^{-\gamma}\frac{q-1}{q}\bigl|\int_\Omega f dx\bigr|^{q/(q-1)} H^{d-1}(\partial\Omega)^{-1/(q-1)} + o(\alpha^{-\gamma})$. The two-sided estimates preceding these expansions also identify the sharp remainder bounds.

Load-bearing premise

For the Dirichlet-limit expansion to hold with a finite leading constant, a certain power of the normal derivative of the limiting Dirichlet solution must be regular enough on the boundary to be extended into the domain; this assumption is not automatic for $p$-harmonic functions on Lipschitz domains, and if it fails the two-sided bounds in Theorem 3.1 break down.

Editorial extensions

If this is right

  • The Robin minimizers converge to the Dirichlet solution, and the energy gap has leading order $\alpha^{-1/(q-1)}$ with a coefficient fixed by the boundary trace of the Dirichlet solution's normal derivative.
  • In the zero-mean-source case, the approach to the Neumann energy is linear in $\alpha$, with slope $\frac1q\int_{\partial\Omega}|u_0|^q dS$.
  • In the nonzero-mean case, the leading divergence is $-\frac{q-1}{q}\alpha^{-1/(q-1)}|\int_\Omega f|^{q/(q-1)} H^{d-1}(\partial\Omega)^{-1/(q-1)}$, so only the mean of $f$ and the boundary measure matter at leading order.
  • The identity $\frac{d}{d\alpha}E_\alpha = \frac1q\int_{\partial\Omega}|u_\alpha|^q dS$ holds for every $\alpha>0$, giving a variational link between the energy slope and the boundary $L^q$ norm of the minimizer.
  • For $p=q=2$ the Dirichlet-limit formula recovers the linear Robin-to-Dirichlet asymptotics previously studied, so the nonlinear result contains the linear case.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The exponent $1/(q-1)$ gives a natural coupling scaling: if a physical interface penalty is scaled as $\alpha\sim\varepsilon^{-(q-1)}$, the Dirichlet-limit correction stays of order one as $\varepsilon\to0$; the paper does not discuss this scaling interpretation.
  • The same rescaling suggests that the minimizers themselves converge to $u_\infty$ in $W^{1,p}$ at the rate $\alpha^{-1/(q-1)}$, up to the trace-regularity obstructions; the paper only states convergence of the energies.
  • Because the incompatible Neumann case involves only $|\int_\Omega f|$ and $H^{d-1}(\partial\Omega)$, one could optimize the divergence rate over domains of fixed volume by varying the boundary measure; this shape question is not treated in the paper.
  • A natural test of the regularity threshold is to construct a $p$-harmonic Dirichlet solution on a Lipschitz polyhedron whose normal derivative has borderline trace regularity; depending on the outcome, the leading constant may need reinterpretation.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies the asymptotic behavior of the minimum E_α of the nonlinear Robin energy J_α(u) = (1/p)∫_Ω|∇u|^p dx + (α/q)∫_{∂Ω}|u|^q dS − ∫_Ω f u dx on a bounded Lipschitz domain Ω, with p,q>1 and f∈(W^{1,p}(Ω))*. It first establishes basic properties: concavity and local Lipschitz continuity of α↦E_α, existence and uniqueness of the minimizer u_α, and differentiability with dE_α/dα = (1/q)∫_{∂Ω}|u_α|^q dS. The main results are first-order expansions: as α→+∞, E_α approaches the Dirichlet energy with rate α^{-1/(q-1)} under a boundary regularity assumption on ∂_ν u_∞; as α→0+, under the compatibility condition ∫_Ω f=0, E_α approaches the Neumann energy linearly in α, while if ∫_Ω f≠0, E_α diverges like α^{-1/(q-1)}. The proofs rely on convexity, explicit minimization of boundary terms, and approximation by test functions. A central displayed formula in the divergent Neumann case has an incorrect exponent on H^{d-1}(∂Ω) that must be corrected.

Significance. If corrected, the paper provides a natural and interesting nonlinear counterpart to the known linear (p=q=2) Robin asymptotics. Its strengths are the self-contained variational method, the explicit form of the leading constants, and the sharp distinction between the compatible and incompatible Neumann regimes. The proofs of Theorems 3.1, 4.1 and 4.2 are coherent and the convexity/testing arguments are sound. However, the advertised Dirichlet-limit expansion is conditional on a boundary regularity hypothesis on the conormal derivative of the Dirichlet solution, and this hypothesis is not automatic on Lipschitz domains; the abstract and introduction do not state this caveat. Additionally, the leading constant in Theorem 4.2 and Eq. (1.4) is printed with the wrong sign on the boundary-measure exponent.

major comments (3)
  1. [§4.2, Theorem 4.2 and §1, Eq. (1.4)] The exponent of H^{d-1}(∂Ω) is wrong. Testing E_α with the constant function t = sign(∫_Ω f)(|∫_Ω f|/(αH^{d-1}(∂Ω)))^{1/(q-1)} gives min_t [ (αH^{d-1}(∂Ω)/q)|t|^q − t∫_Ω f ] = −(q−1)/q |∫_Ω f|^{q/(q−1)} (αH^{d-1}(∂Ω))^{−1/(q−1)}. Thus the coefficient contains H^{d-1}(∂Ω)^{−1/(q−1)}, and not H^{d-1}(∂Ω)^{1/(q−1)} as stated. The correction is needed in the two-sided bounds, in the 'In particular' conclusion of Theorem 4.2, and in Eq. (1.4). This is a load-bearing formula: the explicit dependence of the divergent constant on the domain is reversed.
  2. [§3, Theorem 3.1, and §1, Abstract] The Dirichlet-limit expansion is presented in the abstract and introduction without qualification, but Theorem 3.1 proves the sharp expansion only under the extra hypothesis that either |∂_ν u_∞|^{(p−q)/(q−1)}∂_ν u_∞ ∈ W^{1−1/p,p}(∂Ω) or, for the final conclusion, ∂_ν u_∞ ∈ L^{(p−1)q/(q−1)}(∂Ω). This regularity is not automatic for p-harmonic functions on generic Lipschitz domains; even in the linear case p=q=2 the conormal derivative of the Dirichlet solution need not be square-integrable on a Lipschitz boundary. Since the explicit constant in the leading term is the paper's central quantitative claim, the authors should either state the expansion as conditional in the abstract and introduction or add a precise discussion of conditions under which the hypothesis is satisfied.
  3. [§4.1, Theorem 4.1 and Abstract] The same presentation issue arises in the Neumann-limit compatible case: the abstract says the energy 'linearly approaches' the Neumann energy, but Theorem 4.1 is proved under the assumption u_0 ∈ L^q(∂Ω), which is not automatic for a W^{1,p} trace when q is large. The theorem itself states the assumption, but the advertised claim should be qualified accordingly, or the paper should explain when the assumption holds.
minor comments (4)
  1. [§3, proof of Theorem 3.1] In the sentence 'the fact that ρ_α → 0 as α → 0+ is trivial', the limit should be α → +∞, consistent with the section.
  2. [§4.2, proof of Theorem 4.2] The paragraph beginning 'For what concerns the upper bound, we let v_u = ...' actually proves the lower bound; the label should be corrected.
  3. [§3, Theorem 3.1] The expression |∂_ν u_∞|^{(p−q)/(q−1)}∂_ν u_∞ is ambiguous when p<q and ∂_ν u_∞=0, since a negative power of zero occurs. It should be defined by continuity as sign(∂_ν u_∞)|∂_ν u_∞|^{(p−1)/(q−1)}.
  4. [§1, Introduction] The heuristic divergence estimate after Eq. (1.2) appears misprinted: the proposed constant test function does not yield the stated bound E_α ≤ H^{d-1}(∂Ω)/q |∫_Ω f|^q − α^{-1/q}, nor does it show divergence to −∞. The correct constant test giving divergence is t = sign(∫_Ω f)(|∫_Ω f|/(αH^{d-1}(∂Ω)))^{1/(q−1)}, which leads to the bound in Theorem 4.2.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the Robin-energy expansions follow from exact variational identities with explicitly computed coefficients, and the sole self-citation is purely contextual.

full rationale

The derivation chain is self-contained. Theorem 3.1 obtains the Dirichlet-limit expansion from the exact change of variables u ↦ u∞ + α^{−γ}u with γ = 1/(q−1), rewriting Eα = E∞ + α^{−γ}Rα (Eq. 3.1); the bulk term is controlled by convexity of t ↦ |t|^p and the boundary term is minimized pointwise via Remark 2.6, giving the coefficient −(q−1)/q ∫∂Ω |∂νu∞|^{(p−1)q/(q−1)} dS. This constant is computed from the problem data (f, Ω, u∞), not fitted to the asymptotics. Theorem 4.1 uses the exact sandwich estimate (4.1) together with weak convergence of traces, and Theorem 4.2 splits the energy into the zero-mean fluctuation problem Kf and a one-dimensional optimization over constants whose minimum is given explicitly by Remark 2.6. The only extra hypotheses are boundary regularity of ∂νu∞ (Theorem 3.1) and u0 ∈ L^q(∂Ω) (Theorem 4.1); these are genuine structural conditions for Lipschitz domains, not assumptions that already contain the target expansions, so the advertised claims may be narrower than stated but are not circular. The sole self-citation, [Ogn24], appears only in the introduction's literature overview and is never used in any proof, so no load-bearing argument reduces to it. No fitted parameters, no renaming of known results, and no imported uniqueness theorems are present. The authors' explicit statement 'this paper is self-contained: all results are proved relying only on elementary variational techniques' is borne out by the proofs.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claims rest only on standard Sobolev space machinery and structural assumptions on the data (f ∈ (W^{1,p})*, Lipschitz domain, regularity of the Dirichlet normal derivative). No constants are fitted to data and no new entities are postulated.

assumptions (4)
  • standard math Standard trace theorems and compact embedding of W^{1,p}(Ω) into L^p(Ω) and L^q(∂Ω) under appropriate range conditions.
    Used in Lemma 2.3 and throughout for compactness and coercivity of the energy functional.
  • standard math Direct method of calculus of variations: strict convexity implies existence and uniqueness of minimizers.
    Invoked in Lemma 2.4 for existence and uniqueness of uα and in the existence claims for the limit problems.
  • standard math Convexity inequality |t+v|^q ≥ |t|^q + q|t|^{q-2}t v and differentiability of convex L^p-norms.
    Used in the lower bound of Theorem 4.2 and in the convexity estimates of Theorem 3.1.
  • standard math Integration by parts / Green's identity for the p-Laplacian with boundary trace.
    Used in Theorem 3.1 to pass from the Euler-Lagrange equation to the energy identity involving the normal derivative of the Dirichlet solution.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Asymptotics of nonlinear Robin energies." pith.science (2026). https://pith.science/paper/XVFBPOP4

@misc{pith2026250606914,
  author       = {Pith},
  title        = {Pith review of: Asymptotics of nonlinear Robin energies},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XVFBPOP4}},
  note         = {Machine review of arXiv:2506.06914}
}
abstract

This paper investigates the asymptotic behavior of a class of nonlinear variational problems with Robin-type boundary conditions on a bounded Lipschitz domain. The energy functional contains a bulk term (the $p$-norm of the gradient), a boundary term (the $q$-norm of the trace) scaled by a parameter $\alpha>0$, and a linear source term. By variational methods, we derive first-order expansions of the minimum as $\alpha\to 0^+$ (Neumann limit) and as $\alpha\to+\infty$ (Dirichlet limit). In the Dirichlet limit, the energy converges to the one of Dirichlet problem with a power-type quantified rate (depending only on $q$), while the Neumann limit exhibits a dichotomy: under a compatibility condition, the energy linearly approaches the one of Neumann problem, otherwise, it diverges as a power of $\alpha$ depending only on $q$.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

7 extracted references · 6 canonical work pages

  1. [1]

    Alber and A

    H.-D. Alber and A. G. Ramm , Asymptotics of the solution to Robin problem, J. Math. Anal. Appl. 349 no. 1 (2009), 156--164 (English). doi:10.1016/j.jmaa.2008.08.038

    work page doi:10.1016/j.jmaa.2008.08.038 2009

  2. [2]

    Belgacem , H

    F. Belgacem , H. BelHadjAli , A. BenAmor , and A. Thabet , Robin Laplacian in the large coupling limit: convergence and spectral asymptotic, Ann. Sc. Norm. Super. Pisa, Cl. Sci. (5) 18 no. 2 (2018), 565--591 (English). doi:10.2422/2036-2145.201601_008

    work page arXiv 2018

  3. [3]

    Bisterzo and G

    A. Bisterzo and G. Siclari , Quantitative spectral stability for operators with compact resolvent, Calculus of Variations and Partial Differential Equations 64 no. 6 (2025). doi:10.1007/s00526-025-03020-7

    work page doi:10.1007/s00526-025-03020-7 2025

  4. [4]

    Brasche and M

    J. Brasche and M. Demuth , Dynkin's formula and large coupling convergence, J. Funct. Anal. 219 no. 1 (2005), 34--69 (English). doi:10.1016/j.jfa.2004.06.007

    work page doi:10.1016/j.jfa.2004.06.007 2005

  5. [5]

    Kova r \' k and K

    H. Kova r \' k and K. Pankrashkin , On the \(p\) - Laplacian with Robin boundary conditions and boundary trace theorems, Calc. Var. Partial Differ. Equ. 56 no. 2 (2017), 29 (English), Id/No 49. doi:10.1007/s00526-017-1138-4

    doi:10.1007/s00526-017-1138-4 2017

  6. [6]

    On asymptotics of Robin eigenvalues in the Dirichlet limit

    R. Ognibene , On asymptotics of R obin eigenvalues in the D irichlet limit, Preprint (2024). Available at https://arxiv.org/abs/2407.19505

    work page Pith review arXiv 2024

  7. [7]

    Pankrashkin , An eigenvalue estimate for a Robin \(p\) -laplacian in \(C^1\) domains, Proc

    K. Pankrashkin , An eigenvalue estimate for a Robin \(p\) -laplacian in \(C^1\) domains, Proc. Am. Math. Soc. 148 no. 10 (2020), 4471--4477 (English). doi:10.1090/proc/15116

    work page doi:10.1090/proc/15116 2020

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.