Nonlinear Rayleigh quotient optimization
Flavio Salizzoni, Luca Sodomaco, Julian Weigert
TL;DR
This work extends Rayleigh quotient optimization from quadratic forms to homogeneous polynomials of degree $\omega$ constrained to the affine cone over a projective variety $X$, introducing the X-eigenpoints and the Rayleigh-Ritz degree $\mathrm{RRdeg}_{\omega}(X)$. It establishes a deep link with Euclidean distance geometry by proving $\mathrm{RRdeg}_{\omega}(X)=\mathrm{DD}(\nu_{\omega}(X),Q_{BW})$ via the Bombieri-Weyl inner product and the Veronese embedding, enabling algebraic-geometric computation of critical points. The authors develop Porteous-type degeneracy formulas to compute or bound $\mathrm{RRdeg}_{\omega}(X)$ for implicit complete intersections and polynomial images, with explicit results for hypersurfaces and Segre-Veronese embeddings, complemented by Macaulay2 implementations. They also present a general-position theory that yields closed-form RR-degree expressions in terms of Chern classes and the generic-distance-degree framework, and a Segre-Veronese singular-tuple perspective that connects $X$-eigenvectors to singular $k$-tuples of multihomogeneous forms. The paper concludes with real-case considerations showing that RRdeg bounds may not tightly bound the number of real eigenpoints and provides concrete examples in rank-one tensor settings that illustrate the computational utility of the formulas.
Abstract
Rayleigh quotient minimization deals with optimizing a quadratic homogeneous function over a sphere. Its critical points correspond to the normalized eigenvectors of the symmetric matrix associated with the quadratic form. In this paper, we consider a homogeneous polynomial objective function $f$ over a sphere, a projective algebraic variety $X$, and we study the $X$-eigenpoints of $f$, which are classes of critical points of $f$ constrained to the sphere and the affine cone over $X$. The number of $X$-eigenpoints of a generic polynomial $f$ is the Rayleigh-Ritz degree of $X$. This invariant is a version of the Euclidean distance degree of a Veronese embedding of $X$. We provide concrete formulas in various scenarios, including those involving varieties of rank-one tensors.
