Extremal Eigenvalues of Weighted Steklov Problems
Pith reviewed 2026-05-18 12:05 UTC · model grok-4.3
The pith
Optimal boundary densities for Steklov eigenvalue extrema exist and satisfy specific structural properties under integral and pointwise constraints.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Existence of optimal densities is established for both min and max of λ_k(ρ). Minimizers are bang-bang and may have disconnected support; maximizers are not necessarily bang-bang. On circular domains infinitely many distinct minimizers arise from rotational symmetry while infinitely many distinct maximizers arise without symmetry. The objective functionals are neither convex nor concave, so a Fréchet-differentiable surrogate is introduced to obtain first-order optimality conditions and to design an efficient numerical algorithm.
What carries the argument
The weighted Steklov eigenvalue λ_k(ρ) defined by the boundary integral constraint ∫_∂Ω ρ u v dσ together with the harmonic extension of u inside Ω, optimized over densities ρ that satisfy m ≤ ρ ≤ M and ∫_∂Ω ρ dσ = constant.
If this is right
- Classical convex-optimization methods cannot be applied directly because the maps ρ ↦ λ_k(ρ) and ρ ↦ 1/λ_k(ρ) are neither convex nor concave.
- Minimizers on the circle are generated by arbitrary rotations of any single bang-bang density.
- Maximization on the circle produces infinitely many geometrically distinct optimal densities unrelated by symmetry.
- The surrogate functional permits derivation of optimality conditions even when the original objective is non-differentiable.
Where Pith is reading between the lines
- The same surrogate approach could be used to obtain optimality conditions for other boundary-weighted spectral problems on non-smooth domains.
- Disconnected support of minimizers suggests that optimal densities may concentrate on subsets of the boundary whose geometry is dictated by nodal sets of the eigenfunctions.
- The observed oscillations in numerical recoveries indicate that mesh refinement or regularization may be needed when the true optimum is bang-bang.
Load-bearing premise
The domain must be bounded and Lipschitz and the admissible densities must obey fixed pointwise bounds together with a fixed integral constraint.
What would settle it
A concrete numerical or analytic example on a Lipschitz domain in which every minimizer is smooth and connected or in which a maximizer is strictly bang-bang would refute the structural claims.
Figures
read the original abstract
We study the optimization of Steklov eigenvalues with respect to a boundary density function $\rho$ on a bounded Lipschitz domain $\Omega \subset \mathbb{R}^N$. We investigate the minimization and maximization of $\lambda_k(\rho)$, the $k$th Steklov eigenvalue, over admissible densities satisfying pointwise bounds and a fixed integral constraint. Our analysis covers both first and higher-order eigenvalues and applies to general, not necessarily convex or simply connected, domains. We establish the existence of optimal solutions and provide structural characterizations: minimizers are bang--bang functions and may have disconnected support, while maximizers are not necessarily bang--bang. On circular domains, the minimization problem admits infinitely many minimizers generated by rotational symmetry, while the maximization problem has infinitely many distinct maximizers that are not symmetry-induced. We also show that the maps $\rho \mapsto \lambda_k(\rho)$ and $\rho \mapsto 1/\lambda_k(\rho)$ are generally neither convex nor concave, limiting the use of classical convex optimization tools. To address these challenges, we analyze the objective functional and introduce a Fr\'echet differentiable surrogate that enables the derivation of optimality conditions. We further design an efficient numerical algorithm, with experiments illustrating the difficulty of recovering optimal densities when they lack smoothness or exhibit oscillations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies extremal problems for the k-th weighted Steklov eigenvalue λ_k(ρ) on a bounded Lipschitz domain Ω ⊂ R^N, where the weight ρ belongs to an admissible set defined by pointwise bounds a ≤ ρ ≤ b and a fixed integral constraint ∫_∂Ω ρ = m. The authors prove existence of minimizers and maximizers, establish that minimizers are bang-bang (and may have disconnected support) while maximizers need not be, show that both ρ ↦ λ_k(ρ) and ρ ↦ 1/λ_k(ρ) are in general neither convex nor concave, and introduce a Fréchet-differentiable surrogate functional to derive first-order optimality conditions when λ_k has multiplicity greater than one. They also treat the circular case, where rotational symmetry produces infinitely many distinct minimizers and maximizers, and present a numerical algorithm together with experiments that illustrate recovery difficulties for non-smooth or oscillatory optima.
Significance. If the central claims hold, the work supplies concrete structural information (bang-bang minimizers, possible disconnection, non-bang-bang maximizers) for a class of Steklov optimization problems that had previously been treated mainly for k=1 or convex domains. The surrogate approach and the explicit treatment of multiplicity are technically useful for extending optimality conditions beyond the simple eigenvalue case. The observation that the maps are neither convex nor concave correctly limits the applicability of convex-relaxation techniques and motivates the surrogate construction. The numerical examples on the disk usefully demonstrate both symmetry-induced multiplicity of optima and the practical challenge of recovering bang-bang densities.
major comments (2)
- [§4] §4 (surrogate construction and optimality conditions): the stationarity condition obtained from the Fréchet-differentiable surrogate must be shown to imply the corresponding first-order condition for the original min-max functional λ_k when the k-th eigenvalue has multiplicity >1. A limiting argument or exact equivalence between critical points of the surrogate and those of λ_k is required; without it the bang-bang characterization for minimizers rests on an unverified transfer.
- [Existence theorem] Theorem on existence (likely §3): the compactness argument uses boundedness in L^∞ together with the L^1 constraint and weak-* convergence; the proof should explicitly verify that the Steklov eigenvalues are continuous with respect to this topology on the admissible set, especially when the boundary is only Lipschitz.
minor comments (2)
- [Numerical experiments] The numerical section should specify the finite-element or boundary-element discretization used for the Steklov eigenvalue problem and the quadrature rule for the integral constraint.
- [Notation and preliminaries] Notation for the admissible set (a, b, m) and the surrogate parameter should be introduced once and used consistently; a short table summarizing the constants would improve readability.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. We respond to each major comment below and indicate the planned revisions.
read point-by-point responses
-
Referee: [§4] §4 (surrogate construction and optimality conditions): the stationarity condition obtained from the Fréchet-differentiable surrogate must be shown to imply the corresponding first-order condition for the original min-max functional λ_k when the k-th eigenvalue has multiplicity >1. A limiting argument or exact equivalence between critical points of the surrogate and those of λ_k is required; without it the bang-bang characterization for minimizers rests on an unverified transfer.
Authors: We thank the referee for highlighting this point. The surrogate is designed as a differentiable approximation to handle multiplicity. In the revision we will add a rigorous limiting argument showing that stationarity conditions for the surrogate converge to the first-order optimality conditions of the original functional λ_k as the regularization parameter tends to zero, thereby justifying the transfer to the bang-bang characterization. revision: yes
-
Referee: [Existence theorem] Theorem on existence (likely §3): the compactness argument uses boundedness in L^∞ together with the L^1 constraint and weak-* convergence; the proof should explicitly verify that the Steklov eigenvalues are continuous with respect to this topology on the admissible set, especially when the boundary is only Lipschitz.
Authors: We agree that an explicit verification strengthens the argument. While continuity of Steklov eigenvalues under weak-* convergence of densities on Lipschitz domains follows from standard spectral perturbation results for the boundary integral operator, we will insert a dedicated lemma or remark in the revised manuscript that directly establishes this continuity with respect to the topology used in the compactness argument. revision: yes
Circularity Check
No circularity: variational existence and surrogate-based optimality conditions are self-contained
full rationale
The paper derives existence of optimal densities via compactness arguments in the admissible set of bounded functions with fixed integral constraint, then introduces a Fréchet-differentiable surrogate solely to obtain first-order conditions when the original eigenvalue map is non-differentiable due to multiplicity. No claimed result (existence, bang-bang structure of minimizers, or non-bang-bang maximizers) is shown to reduce by construction to a fitted parameter, a self-definition, or a load-bearing self-citation; the surrogate serves as an auxiliary device whose critical points are transferred to the original problem via explicit limiting arguments within the given Lipschitz-domain setting. The analysis therefore remains independent of its inputs and does not exhibit any of the enumerated circular patterns.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Ω is a bounded Lipschitz domain in R^N
- domain assumption Admissible densities satisfy pointwise bounds and a fixed integral constraint
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
To address these challenges, we analyze the objective functional and introduce a Fréchet differentiable surrogate that enables the derivation of optimality conditions.
-
IndisputableMonolith/Foundation/BranchSelection.leanbranch_selection unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
minimizers are bang–bang functions and may have disconnected support
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Eldar Akhmetgaliyev, Chiu-Yen Kao, and Braxton Osting. Computational methods for extremal Steklov prob- lems.SIAM Journal on Control and Optimization, 55(2):1226–1240, 2017
work page 2017
-
[2]
Pedro Ricardo Sim˜ ao Antunes and Filippo Gazzola. Convex shape optimization for the least biharmonic Stekloveigenvalue.ESAIM: Control, Optimisation and Calculus of Variations, 19(2):385–403, 2013
work page 2013
-
[3]
Pedro RS Antunes and Beniamin Bogosel. Optimization of the Steklov-Lam´ e eigenvalues with respect to the domain.Journal of Differential Equations, 426:1–35, 2025
work page 2025
-
[4]
Pedro RS Antunes, Seyyed Abbas Mohammadi, and Heinrich Voss. A nonlinear eigenvalue optimization problem: Optimal potential functions.Nonlinear Analysis: Real World Applications, 40:307–327, 2018
work page 2018
-
[5]
Springer Science & Business Media, 2009
Kendall Atkinson and Weimin Han.Theoretical Numerical Analysis: A Functional Analysis Framework, vol- ume 39. Springer Science & Business Media, 2009
work page 2009
-
[6]
Steklov eigenproblems and the representation of solutions of elliptic boundary value problems
Giles Auchmuty. Steklov eigenproblems and the representation of solutions of elliptic boundary value problems. Numerical Functional Analysis and Optimization, 25(3-4):321–348, 2005
work page 2005
-
[7]
Catherine Bandle.Isoperimetric inequalities and applications, volume 7. Pitman Publishing, 1980. 22
work page 1980
-
[8]
Farid Bozorgnia. Optimal partitions for first eigenvalues of the Laplace operator.Numerical Methods for Partial Differential Equations, 31(3):923–949, 2015
work page 2015
-
[9]
Farid Bozorgnia, Leon Bungert, and Daniel Tenbrinck. The infinity Laplacian eigenvalue problem: reformulation and a numerical scheme.Journal of Scientific Computing, 98(2):40, 2024
work page 2024
-
[10]
Dorin Bucur and Giuseppe Buttazzo.Variational methods in shape optimization problems. Springer, 2005
work page 2005
-
[11]
Sagun Chanillo, Daniel Grieser, Masaki Imai, Kazuhiro Kurata, and Isamu Ohnishi. Symmetry breaking and other phenomena in the optimization of eigenvalues for composite membranes.Communications in Mathematical Physics, 214(2):315–337, 2000
work page 2000
-
[12]
Bruno Colbois, Alexandre Girouard, Carolyn Gordon, and David Sher. Some recent developments on the Steklov eigenvalue problem.Revista Matem´ atica Complutense, 37(1):1–161, 2024
work page 2024
-
[13]
Carlos Conca, Antoine Laurain, and Rajesh Mahadevan. Minimization of the ground state for two phase con- ductors in low contrast regime.SIAM Journal on Applied Mathematics, 72(4):1238–1259, 2012
work page 2012
-
[14]
Steven Cox and Robert Lipton. Extremal eigenvalue problems for two-phase conductors.Archive for rational mechanics and analysis, 136(2):101–118, 1996
work page 1996
-
[15]
Steven J Cox. The two phase drum with the deepest bass note.Japan journal of industrial and applied mathe- matics, 8(3):345–355, 1991
work page 1991
-
[16]
Elisa Davoli, Idriss Mazari, and Ulisse Stefanelli. Spectral optimization of inhomogeneous plates.SIAM Journal on Control and Optimization, 61(2):852–871, 2023
work page 2023
-
[17]
Behrouz Emamizadeh and Mohsen Zivari-Rezapour. Rearrangements and minimization of the principal eigen- value of a nonlinear Steklov problem.Nonlinear Analysis: Theory, Methods & Applications, 74(16):5697–5704, 2011
work page 2011
-
[18]
Ailana Fraser and Richard Schoen. Sharp eigenvalue bounds and minimal surfaces in the ball.Inventiones mathematicae, 203(3):823–890, 2016
work page 2016
-
[19]
Shmuel Friedland. Extremal eigenvalue problems defined for certain classes of functions.Archive for Rational Mechanics and Analysis, 67(1):73–81, 1977
work page 1977
-
[20]
David Gilbarg and Neil S Trudinger.Elliptic partial differential equations of second order. springer, 2015
work page 2015
-
[21]
On the Hersch-Payne-Schiffer inequalities for Steklov eigenvalues
Alexandre Girouard and Iosif Polterovich. On the Hersch-Payne-Schiffer inequalities for Steklov eigenvalues. Functional Analysis and its Applications, 44(2):106–117, 2010
work page 2010
-
[22]
Spectral geometry of the Steklov problem.Journal of Spectral Theory, 7(2):321–359, 2017
Alexandre Girouard and Iosif Polterovich. Spectral geometry of the Steklov problem.Journal of Spectral Theory, 7(2):321–359, 2017
work page 2017
-
[23]
Jens A Griepentrog and Lutz Recke. Linear elliptic boundary value problems with non–smooth data: normal solvability on sobolev–campanato spaces.Mathematische Nachrichten, 225(1):39–74, 2001
work page 2001
-
[24]
Springer Science & Business Media, 2006
Antoine Henrot.Extremum problems for eigenvalues of elliptic operators. Springer Science & Business Media, 2006
work page 2006
-
[25]
Antoine Henrot.Shape optimization and spectral theory. De Gruyter Open, 2017
work page 2017
-
[26]
European Mathematical Society (EMS) Publishing House, Z¨ urich, Switzerland, 2018
Antoine Henrot and Michel Pierre.Shape Variation and Optimization: A Geometrical Analysis, volume 28 of EMS Tracts in Mathematics. European Mathematical Society (EMS) Publishing House, Z¨ urich, Switzerland, 2018
work page 2018
-
[27]
Di Kang, Patrick Choi, and Chiu-Yen Kao. Minimization of the first nonzero eigenvalue problem for two-phase conductors with Neumann boundary conditions.SIAM Journal on Applied Mathematics, 80(4):1607–1628, 2020
work page 2020
-
[28]
Chiu Y Kao, Stanley Osher, and Eli Yablonovitch. Maximizing band gaps in two-dimensional photonic crystals by using level set methods.Applied Physics B, 81(2):235–244, 2005
work page 2005
-
[29]
Chiu-Yen Kao, Junshan Lin, and Braxton Osting. A semi-definite optimization method for maximizing the shared band gap of topological photonic crystals.Journal of Computational Physics, 521:113538, 2025
work page 2025
-
[30]
Principal eigenvalue for an elliptic problem with indefinite weight on cylindrical domains.Math
Chiu-Yen Kao, Yuan Lou, and Eiji Yanagida. Principal eigenvalue for an elliptic problem with indefinite weight on cylindrical domains.Math. Biosci. Eng, 5(2):315–335, 2008
work page 2008
-
[31]
Chiu-Yen Kao and Shu Su. Efficient rearrangement algorithms for shape optimization on elliptic eigenvalue problems.Journal of Scientific Computing, 54(2):492–512, 2013
work page 2013
-
[32]
Mark G. Krein.On Certain Problems on the Maximum and Minimum of Characteristic Values and on the Lyapunov Zones of Stability. American Mathematical Society translations. American Mathematical Society, 1955
work page 1955
-
[33]
The legacy of Vladimir Andreevich Steklov.Notices of the AMS, 61(1):190, 2014
Nikolay Kuznetsov, Tadeusz Kulczycki, M Kwa´ snicki, Alexander Nazarov, Sergey Poborchi, Iosif Polterovich, and Bart lomiej Siudeja. The legacy of Vladimir Andreevich Steklov.Notices of the AMS, 61(1):190, 2014
work page 2014
-
[34]
Viewing the Steklov eigenvalues of the Laplace operator as critical Neumann eigenvalues
Pier Domenico Lamberti and Luigi Provenzano. Viewing the Steklov eigenvalues of the Laplace operator as critical Neumann eigenvalues. InCurrent trends in analysis and its applications, pages 171–178. Springer, 2015
work page 2015
-
[35]
American Mathematical Soc., 2001
Elliott H Lieb and Michael Loss.Analysis, volume 14. American Mathematical Soc., 2001
work page 2001
-
[36]
Andrea Manzoni, Alfio Quarteroni, and Sandro Salsa.Optimal control of partial differential equations. Springer, 2021. 23
work page 2021
-
[37]
Idriss Mazari, Gr´ egoire Nadin, and Yannick Privat. Shape optimization of a weighted two-phase Dirichlet eigen- value.Archive for Rational Mechanics and Analysis, 243(1):95–137, 2022
work page 2022
-
[38]
Idriss Mazari and Yannick Privat. Qualitative analysis of optimisation problems with respect to non-constant robin coefficients.Annali della Scuola Normale Superiore di Pisa, Classe di Scienze, 2022
work page 2022
-
[39]
Idriss Mazari and Yannick Privat. Qualitative analysis of optimisation problems with respect to non-constant robin coefficients.ANNALI SCUOLA NORMALE SUPERIORE-CLASSE DI SCIENZE, pages 1309–1356, 2024
work page 2024
-
[40]
Abbasali Mohammadi, Fariba Bahrami, and Hakimeh Mohammadpour. Shape dependent energy optimization in quantum dots.Applied Mathematics Letters, 25(9):1240–1244, 2012
work page 2012
-
[41]
Abbasali Mohammadi and Mohsen Yousefnezhad. Optimal ground state energy of two-phase conductors.Elec- tronic Journal of Differential Equations, 2014(171):1–8, 2014
work page 2014
-
[42]
Extremal principal eigenvalue of the bi-Laplacian operator
Seyyed Abbas Mohammadi and Fariba Bahrami. Extremal principal eigenvalue of the bi-Laplacian operator. Applied Mathematical Modelling, 40(3):2291–2300, 2016
work page 2016
-
[43]
Seyyed Abbas Mohammadi, Farid Bozorgnia, and Heinrich Voss. Optimal shape design for the p-Laplacian eigenvalue problem.Journal of Scientific Computing, 78(2):1231–1249, 2019
work page 2019
-
[44]
Seyyed Abbas Mohammadi and Heinrich Voss. A minimization problem for an elliptic eigenvalue problem with nonlinear dependence on the eigenparameter.Nonlinear Analysis: Real World Applications, 31:119–131, 2016
work page 2016
-
[45]
Robin Nittka. Regularity of solutions of linear second order elliptic and parabolic boundary value problems on Lipschitz domains.Journal of Differential Equations, 251(4-5):860–880, 2011
work page 2011
-
[46]
Edouard Oudet, Chiu-Yen Kao, and Braxton Osting. Computation of free boundary minimal surfaces via ex- tremal Steklov eigenvalue problems.ESAIM: Control, Optimisation and Calculus of Variations, 27:34, 2021
work page 2021
-
[47]
Abdelkader Al Sayed, Beniamin Bogosel, Antoine Henrot, and Florent Nacry. Maximization of the Steklov eigenvalues with a diameter constraint.SIAM Journal on Mathematical Analysis, 53(1):710–729, 2021
work page 2021
-
[48]
Nathan Schroeder, Weaam Alhejaili, and Chiu-Yen Kao. Steklov eigenvalue problems on nearly spherical and nearly annular domains.arXiv preprint arXiv:2310.19585, 2023
-
[49]
American Mathematical Soc., 2010
Fredi Tr¨ oltzsch.Optimal control of partial differential equations: theory, methods, and applications, volume 112. American Mathematical Soc., 2010
work page 2010
-
[50]
Robert Weinstock. Inequalities for a classical eigenvalue problem.Indiana University Mathematics Journal, 3(6):745–753, 1954. Department of Mathematical Sciences, Claremont McKenna College, Claremont, CA 91711 Email address:ckao@cmc.edu Division of Mathematics, University of Dundee, Dundee DD1 4HN, United Kingdom; School of Computer Science and Applied Ma...
work page 1954
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.