Pith. sign in

REVIEW 2 major objections 6 minor 1 cited by

Multi-bubble solutions for the Brezis-Nirenberg problem in four dimensions

T0 review · 2 major / 6 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 The four-dimensional Brezis-Nirenberg problem admits positive multi-bubble solutions.

desk verdict Solid construction of positive multi-bubbles in 4D with a repairable gap in the degree step and some minor presentational issues. read the letter →

arxiv 2505.24387 v1 pith:TOLOISFB submitted 2025-05-30 math.AP

Angela Pistoia , Giuseppe Mario Rago , Giusi Vaira This is my paper
classification math.AP MSC 35B4435B3335J25
keywords Brezis-Nirenbergproblemmulti-bubblesolutionsblow-upanalysisfourdimensionsLyapunov-SchmidtreductionRobinfunctionGreen'sstablecriticalset
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

The paper proves that the four-dimensional Brezis-Nirenberg problem, which had resisted all previous multi-bubble constructions, admits positive solutions that blow up at any prescribed number $k$ of interior points, provided a spectral condition holds. The condition is that the smallest eigenvalue $\Lambda_1$ of a matrix built from the Robin function and the Green's function has a stable critical set on which it is strictly positive. This fills a long-standing gap: the same rate condition was already known to be necessary from the fine multibubble asymptotic analysis of the problem, and the paper shows it is sufficient under these hypotheses. The proof uses a Lyapunov-Schmidt reduction whose decisive step is solving a non-variational reduced problem by identifying its vector field with $\nabla\Lambda_1$ and applying a Brouwer-degree argument. Explicit dumbbell domains realize the condition for arbitrarily many bubbles, and thin annuli realize it for two or four bubbles.

What carries the argument

The central object is the $k\times k$ matrix $M(\xi)$ whose diagonal entries are the Robin function $\tau_\Omega(\xi_i)$ and whose off-diagonal entries are $-G(\xi_i,\xi_j)$, where $G$ is the Dirichlet Green's function. Its smallest eigenvalue $\Lambda_1(\xi)$ is simple with a strictly positive eigenvector by the Perron-Frobenius theorem, and it selects the exponential concentration rates: the ansatz sets $\delta_1 = e^{-8\pi^2\lambda/\varepsilon}$ and $\delta_i = d_i\delta_1$, and the reduced problem forces $\lambda=\Lambda_1(\xi)$. The load-bearing step is showing that the reduced vector field $F_3(\xi,d(\xi),\Lambda_1(\xi))$ equals $\nabla\Lambda_1(\xi)$; together with a degree-preserving lemma, a stable critical set of $\Lambda_1$ with positive values gives a non-zero Brouwer degree for the full reduced system. The Lyapunov-Schmidt reduction then turns this zero of the reduced problem into an actual solution of the PDE.

What would settle it

Numerically solve the reduced system (2.24) in a dumbbell with three wells satisfying condition (4.34): if no solution exists near the predicted set with $\lambda=\Lambda_1(\xi)>0$ even though $\Lambda_1$ has a stable positive critical set, the degree argument would be invalid. Alternatively, evaluate the explicit series for $\Lambda_1(r)$ in the thin annulus for $k=3$; a minimum $\le 0$ for some $\rho$ sufficiently close to $1$ would disprove the positivity conjecture that the annulus examples rely on.

Watch

Extended reading notes

Core claim

On a bounded smooth domain in $\mathbb{R}^4$, for any integer $k\ge 1$, if the smallest eigenvalue $\Lambda_1(\xi)$ of the matrix $M(\xi)$ with diagonal entries $\tau_\Omega(\xi_i)$ and off-diagonal entries $-G(\xi_i,\xi_j)$ has a stable critical set $K$ with $\Lambda_1>0$ on $K$, then there exists a family of positive solutions to $-\Delta u = u^3 + \varepsilon u$ in $\Omega$ with $u=0$ on $\partial\Omega$ that blows up at $k$ distinct points $\xi^0=(\xi_1^0,\dots,\xi_k^0)\in K$. The concentration rates are exponentially small in $\varepsilon$ and satisfy $\varepsilon \log \delta_{i,\varepsilon}^{-1}\to \Lambda_1(\xi^0)$ as $\varepsilon\to 0$. This is the four-dimensional counterpart of the previously established necessary condition from the fine multibubble asymptotic analysis, and it closes the open question of existence of positive multi-bubble solutions in dimension four. The construction is illustrated by a dumbbell domain with $k$ components and by a thin annulus with two or four symmetric blow-up points.

Load-bearing premise

The entire construction rests on the hypothesis that the smallest eigenvalue $\Lambda_1$ of the bubble-interaction matrix has a stable critical set where it is strictly positive; the paper verifies this only in the dumbbell and in the two- and four-point annulus cases, and if $\Lambda_1\le 0$ at the would-be blow-up locations the exponential ansatz cannot produce the claimed rates.

Editorial extensions

If this is right

  • The necessary rate condition $\varepsilon \log \delta_{i,\varepsilon}^{-1}\to \Lambda_1(\xi^0)$ from the fine multibubble analysis becomes sufficient under the stable-critical-set hypothesis, completing the asymptotic classification for four-dimensional multi-bubble blow-up.
  • In a dumbbell obtained by joining $k$ subdomains by thin necks, positive solutions concentrate at $k$ points, with $k$ arbitrary, as soon as the Robin function of one well is strictly smaller than the others.
  • In a thin annulus, positive solutions concentrate at two or four symmetric points for sufficiently small thickness; the authors conjecture the number of peaks grows as the thickness decreases.
  • The reduced problem is solved without a variational structure by identifying the reduced vector field with $\nabla\Lambda_1$ and using Brouwer degree, a mechanism that may carry over to other critical problems with exponentially small interaction scales.

Reading between the lines

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

  • If $\Lambda_1>0$ is the only real obstruction, the method suggests that any domain whose Green's-function matrix has a positive smallest eigenvalue with a stable critical set should carry multi-bubble solutions; testing thin annuli for $k=3,5,\dots$ with the explicit series could turn the paper's conjecture into a theorem.
  • The same eigenvector-based reduction may apply to sign-changing bubbles or to slightly supercritical problems, where the reduced equations typically lack variational structure.
  • Because the concentration rates are exponentially small in $1/\varepsilon$, these solutions are invisible to polynomial-in-$\varepsilon$ expansions; any numerical detection would need to resolve exponentially thin scales.
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

2 major / 6 minor

Summary. The paper constructs multi-bubble positive solutions for the four-dimensional Brezis-Nirenberg problem −Δu = u^3 + εu in Ω with u = 0 on ∂Ω. The main result (Theorem 1.2) states that if K is a stable critical set of the smallest eigenvalue Λ1(ξ) of the matrix M(ξ) built from the Robin function and Green function, and if Λ1 > 0 on K, then for all small ε there exists a family of positive solutions blowing up at k points ξ^0_i ∈ K with concentration rates δ_{i,ε} satisfying ε log δ_{i,ε}^{-1} → Λ1(ξ^0) as ε → 0. The proof uses a Lyapunov-Schmidt reduction with an ansatz of k Aubin-Talenti bubbles whose rates are exponentially small in 1/ε, leading to a 5k-dimensional reduced system (2.24). The reduced system is solved by degree theory: the authors reduce the problem to a stable critical set of Λ1 via a key lemma. Concrete examples are provided for dumbbell domains (arbitrary k) and for thin annuli (k = 2, 4).

Significance. If correct, this is the first construction of positive multi-bubble solutions for the Brezis-Nirenberg problem in dimension four, and it shows that the necessary concentration-rate condition identified by König and Laurain is sufficient under a natural spectral hypothesis. The proof is based on well-documented estimates (Appendix A) and derives the exponential rate formula from the reduced equations rather than imposing it. The paper also gives concrete domains where the hypothesis Λ1 > 0 is verified. The main weakness is a localized but load-bearing error in the degree argument at Eq. (3.31), which appears repairable. The paper is otherwise coherent and the computations are detailed.

major comments (2)
  1. [Section 3.1, Eq. (3.31)] The claimed identity φ(ξ) = ∇Λ1(ξ) is false. The computation displayed immediately above (3.31) yields (1 + Σ_{j=2}^k d_j^2) ∂_{ξ_1}Λ1 = (M̃^1(1,d)^T)_1 and (1 + Σ_{j=2}^k d_j^2) ∂_{ξ_i}Λ1 = d_i (M̃^i(1,d)^T)_i for i = 2, ..., k. Since F3 is defined as the vector of the (M̃^ℓ(1,d)^T)_i components, it follows that F3_i = c_i ∂_{ξ_i}Λ1 with c_1 = 1 + Σ d_j^2 and c_i = (1 + Σ d_j^2)/d_i for i ≥ 2. These coefficients are positive but not all equal, so F3(ξ,d(ξ),Λ1(ξ)) is a positive diagonal multiple of ∇Λ1, not ∇Λ1 itself. Consequently the inference "by (3.31) deg(φ,Θ,0) ≠ 0" is not justified as written. The gap is repairable because multiplication by the positive diagonal matrix diag(c_i) is an orientation-preserving homeomorphism and preserves the Brouwer degree, and the zero sets coincide; nevertheless, the proof must be corrected.
  2. [Section 3.1, Eqs. (3.27)-(3.28)] The existence and uniqueness statement for d(ξ) ∈ (0,+∞)^{k-1} is not fully proved. The text solves F1 = F2 = 0 and derives the eigenvector relation M(ξ)(1,d)^T = Λ1(ξ)(1,d)^T, but it never explicitly verifies that the vector d defined by (3.27) has positive entries. Positivity is needed for the domain of the reduced problem and for the sign of the diagonal coefficients in the corrected form of (3.31). The missing step is immediate: by (1.3) and the Perron-Frobenius theorem, the eigenvector e(ξ) has positive entries and first component 1, so (1,d(ξ)) must equal e(ξ); this should be stated explicitly.
minor comments (6)
  1. [Eq. (1.3)] The eigenvector e(ξ) is said to belong to R^4; it should be R^k, since M(ξ) is a k×k matrix.
  2. [Section 4.1, proof of Proposition 4.1] The reference "Definition ??" should be "Definition 1.1", and the set denoted D in the proof ("that is (ξ1,...,ξk) ∈ D") should be Θ.
  3. [Section 4.2, Eq. (4.40)] The expression "min_{r∈(ρ,1)} Λ1(ρ)" should read "min_{r∈(ρ,1)} Λ1(r)".
  4. [Abstract and Introduction] The name "Laurain" is consistently misspelled as "Laurin" in the abstract, the introduction, and parts of the text; please correct this in all occurrences.
  5. [Appendix A] The equation labels (1.42), (1.43), (1.45), and (1.50) should be numbered within the appendix (e.g., (A.2), (A.3), etc.) to avoid confusion with equations in the introduction.
  6. [Section 3.1, definition of F3] The symbol d is overloaded: in (3.25) and the definition of F3, d^T appears to denote the k-vector (1,d_2,...,d_k), whereas elsewhere d denotes (d_2,...,d_k). Please clarify the notation.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the rate formula (1.5) is derived from solving the reduced system for the unknowns (λ,d,ξ), not imposed or fitted; self-citations are standard, non-load-bearing technical inputs.

full rationale

The paper's central derivation is self-contained with respect to the circularity criteria. The unknowns λ > 0, d_i > 0, and ξ are solved from the reduced system (2.24): λ is not a pre-fitted parameter, but is determined by the equations F1 = 0 and F2 = 0, which force λ(ξ) = Λ1(ξ) and d(ξ) via (3.27)-(3.28). The claimed concentration rate (1.5), ε log δ_{i,ε}^{-1} → Λ1(ξ0), is a consequence of this solution of the reduced problem, not an input to it. The ansatz (2.13) uses the exponential form known from the necessary conditions of König-Laurain, but using a necessary-rate ansatz to construct solutions is not circular: the existence proof still has to solve the reduced system, which it does. Prior expansions from Rey and from Musso-Pistoia are used as computational inputs for the error estimates and reduced equations (Propositions 2.1 and 2.2); these are external, independently established asymptotic facts, not the conclusion being proved, and they are not the load-bearing source of the existence claim. The stable-critical-set hypothesis and the assumption Λ1 > 0 are explicit hypotheses of Theorem 1.2, and the paper honestly limits its verification to the dumbbell and to k = 2, 4 in the thin annulus (Section 4.2), even conjecturing the general case (4.41). These are limitations on scope, not circular steps. The skeptic's objection that (3.31) is false as written is a technical gap in the Brouwer-degree computation: the displayed computation above (3.31) gives positive diagonal multipliers relating F3 to ∇Λ1, so the zero set and the degree conclusion are repairable. That is a correctness risk, not a circularity. No step reduces a prediction to a fitted parameter, a self-citation is not used to forbid alternatives, and no known result is merely renamed. Therefore the circularity score is 0.

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

The proof introduces no fitted constants and no new physical or mathematical entities. The unknowns λ and d are solved from the reduced system. The main external inputs are standard PDE expansions and known Green's-function formulas for the annulus.

assumptions (6)
  • domain assumption The projected bubble expansions (2.10) and (2.11) hold uniformly on compact sets.
    Taken from Rey [30] and used throughout the error estimates and Proposition 2.2.
  • domain assumption The domain Ω is bounded and regular, and blow-up points are separated from the boundary and from each other through the set Dρ.
    This separation controls interactions between bubbles and is assumed in the reduction and in the examples.
  • standard math The Perron-Frobenius theorem guarantees that the smallest eigenvalue Λ1(ξ) is simple with a positive eigenvector.
    Used to define e(ξ) and to prove invertibility of M(ξ)-Λ1(ξ)I in Section 3.1.
  • standard math Stable critical sets with nonzero Brouwer degree persist under small uniform perturbations.
    This is the mechanism at the end of Section 3 that turns the unperturbed degree condition into a solution of the reduced system with o(1) terms.
  • domain assumption The explicit Green and Robin function formulas for the annulus from Grossi and Vujadinović [14] are correct.
    These formulas are used in Proposition 4.2 to verify positivity of Λ1 for k=2 and k=4.
  • ad hoc to paper Λ1(ξ)>0 on the critical set K.
    This is assumed in Theorem 1.2 and only verified in the two examples; no general criterion is proved.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Multi-bubble solutions for the Brezis-Nirenberg problem in four dimensions." pith.science (2026). https://pith.science/paper/TOLOISFB

@misc{pith2026250524387,
  author       = {Pith},
  title        = {Pith review of: Multi-bubble solutions for the Brezis-Nirenberg problem in four dimensions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TOLOISFB}},
  note         = {Machine review of arXiv:2505.24387}
}
read the original abstract

The paper addresses the existence of multi-bubble solutions for the well-known Brezis-Nirenberg problem. Although there is extensive literature on the subject, the existence of solutions that blow up at multiple points in a 4D bounded domain remains an open problem. The goal of the present paper is to resolve this longstanding issue. In particular, we exhibit examples of domains where a large number of multi-bubble solutions exist. Our result can also be seen as the counterpart of the asymptotic analysis carried out by Konig and Laurin in Ann. Inst. H. Poincar\`e C Anal. Non Lin\`eaire, 2024.

Discussion (0). Sign in to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Sharp quantitative stability estimates for the Brezis-Nirenberg problem

    math.AP 2025-06 conditional novelty 7.0 of 10

    Nearly stationary functions for the Brezis-Nirenberg problem on bounded domains lie within a sharp, dimension-dependent distance of a solution plus bubbles, and the optimal exponents are identified.

Reference graph

Works this paper leans on

32 extracted references · 30 canonical work pages · cited by 1 Pith paper

  1. [21]

    F. Liu, G. Vaira, J. Wei, Y. Wu, Construction of bubbling solutions of the Brezis-Nirenberg problem in general bounded domains (I): the dimensions 4 and 5, preprint arXiv:2503.09250 (2025)

    arXiv 2025

  2. [1]

    A. L. Amadori, F. Gladiali, M. Grossi, A. Pistoia, G. Vaira, A complete scenario on nodal radial solutions to the Brezis - Nirenberg problem in low dimensions,Nonlinearity, 34 (11), 8055–8093 (2021)

    work page 2021

  3. [2]

    Aubin, Espaces de Sobolev sur les vari ´et´es riemanniennes, Bull

    T. Aubin, Espaces de Sobolev sur les vari ´et´es riemanniennes, Bull. Sci. Math.,2, 149–173 (1976)

    work page 1976

  4. [3]

    F. V . Atkinson, H. Brezis, L. A. Peletier, Nodal solutions of elliptic equations with critical Sobolev exponents,J. Diff. Eqs., 85, 151–170 (1990)

    work page 1990

  5. [4]

    F. V . Atkinson, L. A. Peletier, Large solutions of elliptic equations involving critical exponents, Asymptotic Analysis, 1, 139–160 (1988)

    work page 1988

  6. [5]

    Axler, P

    S. Axler, P . Bourdon, W. Ramey, Harmonic function theorySpringer-Verlag New York, Inc.Second Edition, (2000)

    work page 2000

  7. [6]

    Bahri and J.-M

    A. Bahri and J.-M. Coron. On a nonlinear elliptic equation involving the critical Sobolev exponent: the effect of the topology of the domain. Comm. Pure Appl. Math., 41 (3), 253–294 (1988)

    work page 1988

  8. [7]

    Bahri, Y

    A. Bahri, Y. Li, and O. Rey, On a variational problem with lack of compactness: The topolog- ical effect of the critical points at infinity. Calc. Var. Partial Differential Equations3, no. 1, 67–93 (1995)

    work page 1995

Show all 32 references
  1. [8]

    Brezis, L

    H. Brezis, L. Nirenberg, Positive solutions of nonlinear elliptic equations involving critical Sobolev exponents, Comm. Pure Appl. Math., 36, no. 4, 437–477 (1983)

    1983

  2. [9]

    Brezis, L

    H. Brezis, L. Peletier, Asymptotics for elliptic equations involving critical growth, pp. 149–192 in Partial Dif- ferential Equations and the Calculus of Variations, Progress in Nonlinear Differential Equations Applications 1, Birkh¨auser, Boston, 1989

    1989

  3. [10]

    L. A. Caffarelli, B. Gidas, J. Spruck, Asymptotic symmetry and local behavior of semilinear elliptic equations with critical Sobolev growth,Comm. Pure Appl. Math., 42 no.3, 271–297 (1989)

    1989

  4. [11]

    del Pino, P

    M. del Pino, P . Felmer, M. Musso, Multi-bubble solutions for slightly super-critical elliptic problems in domains with symmetries. Bull. London Math. Soc.35 no. 4, 513–521 (2003)

    2003

  5. [12]

    Druet, Elliptic equations with critical Sobolev exponents in dimension 3, Ann

    O. Druet, Elliptic equations with critical Sobolev exponents in dimension 3, Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire,19 (2) 125–142 (2002)

    2002

  6. [13]

    Esposito, N

    P . Esposito, N. Ghoussoub, A. Pistoia, and G. Vaira. Sign-changing solutions for critical equations with Hardy potential, Anal. PDE, 14 2, 533–566 (2021)

    2021

  7. [14]

    Grossi, D

    M. Grossi, D. Vujadinovi ´c. On the Green function of the annulus, Anal. Theory Appl., 32,1, 52–64. (2016)

    2016

  8. [15]

    Z.-C. Han. Asymptotic approach to singular solutions for nonlinear elliptic equations involving critical Sobolev exponent,Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire, 8, 2 159–174 (1991). MULTI-BUBBLE SOLUTIONS FOR THE BREZIS-NIRENBERG PROBLEM IN FOUR DIMENSIONS 23

    1991

  9. [16]

    K ¨onig, P

    T. K ¨onig, P . Laurain. Fine multibubble analysis in the higher-dimensional Brezis- Nirenberg problem,Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire,41 , no. 5, 1239–1287 (2024)

    2024

  10. [17]

    K ¨onig, P

    T. K ¨onig, P . Laurain, Multibubble blow-up analysis for the Brezis-Nirenberg problem in three dimensions, ArXiv 2208.12337 (2022)

    2022

  11. [18]

    Iacopetti, G

    A. Iacopetti, G. Vaira, Sign-changing tower of bubbles for the Brezis-Nirenberg problem, Commun. Contemp. Math., 18 , no. 1, 1550036, 53 pp. (2016)

    2016

  12. [19]

    Iacopetti, G

    A. Iacopetti, G. Vaira, Sign-changing blowing-up solutions for the Brezis-Nirenberg problem in dimension four and five, Ann. Sc. Norm. Pisa Cl. Sci., 18, 1–38 (2018)

    2018

  13. [20]

    Li, Yan Yan On a singularly perturbed equation with Neumann boundary condition. Comm. Partial Differential Equations 23 (1998), no. 3-4, 487–545

    1998

  14. [22]

    Morabito, A

    F. Morabito, A. Pistoia, and G. Vaira, Towering phenomena for the Yamabe equation on symmetric manifolds. Potential Anal., 47 (1) 53–102,( 2017)

    2017

  15. [23]

    Musso, A

    M. Musso, A. Pistoia, Multispike solutions for a nonlinear elliptic problem involving the critical Sobolev expo- nent, Indiana Univ. Math. J.51 , no. 3, 541–579 (2002)

    2002

  16. [24]

    Musso, D

    M. Musso, D. Salazar, Multispike solutions for the Brezis-Nirenberg problem in dimension three J. Differential Equations 264 no. 11, 6663–6709 (2018)

    2018

  17. [25]

    Pistoia, O

    A. Pistoia, O. Rey, Multiplicity of solutions to the supercritical Bahri-Coron’s problem in pierced domains.Adv. Differential Equations 11 no. 6, 647–666 (2006)

    2006

  18. [26]

    Pistoia, S

    A. Pistoia, S. Rocci, The Brezis-Nirenberg problem in 4D, Discrete Contin. Dyn. Syst. Ser. S 17 no. 4, 1562–1572 (2024)

    2024

  19. [27]

    Pistoia, G

    A. Pistoia, G. Vaira, Nodal solutions of the Brezis-Nirenberg problem in dimension 6, Anal. Theory Appl., 38 (2022), no. 1, 1–25

    2022

  20. [28]

    Pistoia and G

    A. Pistoia and G. Vaira, Clustering phenomena for linear perturbation of the Yamabe equation. In Partial dif- ferential equations arising from physics and geometry, volume 450 of London Math. Soc. Lecture Note Ser., pages 311–331. Cambridge Univ. Press, Cambridge, 2019

    2019

  21. [29]

    Premoselli

    B. Premoselli. Towers of bubbles for Yamabe-type equations and for the Brezis-Nirenberg problem in dimen- sions n ≥ 7,J. Geom. Anal., 32 (3):Paper No. 73, 65, (2022)

    2022

  22. [30]

    Rey, The role of the Green’s function in a nonlinear elliptic equation involving the critical Sobolev exponent, J

    O. Rey, The role of the Green’s function in a nonlinear elliptic equation involving the critical Sobolev exponent, J. Funct. Anal., 89, no. 1, 1–52 (1990)

    1990

  23. [31]

    Stein, G

    E. Stein, G. Weiss, Introduction to Fourier Analysis on Euclidean Spaces Princeton University Press, (1971)

    1971

  24. [32]

    Talenti, Best constant in Sobolev inequality,Ann

    G. Talenti, Best constant in Sobolev inequality,Ann. Mat. Pura Appl., (4) 353–372 (1976). DIPARTIMENTO DI SCIENZE DI BASE E APPLICATE PER L ’I NGEGNERIA , U NIVERSIT `A DEGLI STUDI DI ROMA SAPIENZA Email address: angela.pistoia@uniroma1.it DIPARTIMENTO DI MATEMATICA , U NIVERS...

    1976

Pith tools

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