Extremal Eigenvalues of Weighted Steklov Problems

We study the optimization of Steklov eigenvalues with respect to a boundary density function $ρ$ on a bounded Lipschitz domain $Ω\subset \mathbb{R}^N$. We investigate the minimization and maximization of $λ_k(ρ)$, 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 $ρ\mapsto λ_k(ρ)$ and $ρ\mapsto 1/λ_k(ρ)$ 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échet 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.