Source author record

Grigoriy Blekherman

Grigoriy Blekherman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

A Path Forward: Tropicalization in Extremal Combinatorics

Many important problems in extremal combinatorics can be be stated as proving a pure binomial inequality in graph homomorphism numbers, i.e., proving that hom$(H_1,G)^{a_1}\cdots$hom$(H_k,G)^{a_k}\geq$hom$(H_{k+1},G)^{a_{k+1}}\cdots$hom$(H_m,G)^{a_m}$ holds for some fixed graphs $H_1,\dots,H_m$ and all graphs $G$. One prominent example is Sidorenko's conjecture. For a fixed collection of graphs $\mathcal{U}=\{H_1,\dots,H_m\}$, the exponent vectors of valid pure binomial inequalities in graphs of $\mathcal{U}$ form a convex cone. We compute this cone for several families of graphs including complete graphs, even cycles, stars and paths; the latter is the most interesting and intricate case that we compute. In all of these cases, we observe a tantalizing polyhedrality phenomenon: the cone of valid pure binomial inequalities is actually rational polyhedral, and therefore all valid pure binomial inequalities can be generated from the finite collection of exponent vectors of the extreme rays. Using the work of Kopparty and Rossman, we show that the cone of valid inequalities is indeed rational polyhedral when all graphs $H_i$ are series-parallel and chordal, and we conjecture that polyhedrality holds for any finite collection $\mathcal{U}$. We demonstrate that the polyhedrality phenomenon also occurs in matroids and simplicial complexes. Our description of the inequalities for paths involves a generalization of the Erdős-Simonovits conjecture recently proved in its original form by Sağlam and a new family of inequalities not observed previously. We also solve an open problem of Kopparty and Rossman on the homomorphism domination exponent of paths. One of our main tools is tropicalization, a well-known technique in complex algebraic geometry. We prove several results about tropicalizations which may be of independent interest.

preprint2022arXiv

Tropicalization of Graph Profiles

A graph profile records all possible densities of a fixed finite set of graphs. Profiles can be extremely complicated; for instance the full profile of any triple of connected graphs is not known, and little is known about hypergraph profiles. We introduce the tropicalization of graph and hypergraph profiles. Tropicalization is a well-studied operation in algebraic geometry, which replaces a variety (the set of real or complex solutions to a finite set of algebraic equations) with its "combinatorial shadow". We prove that the tropicalization of a graph profile is a closed convex cone, which still captures interesting combinatorial information. We explicitly compute these tropicalizations for arbitrary sets of complete and star hypergraphs. We show they are rational polyhedral cones even though the corresponding profiles are not even known to be semialgebraic in some of these cases. We then use tropicalization to prove strong restrictions on the power of the sums of squares method, equivalently Cauchy-Schwarz calculus, to test (which is weaker than certification) the validity of graph density inequalities. In particular, we show that sums of squares cannot test simple binomial graph density inequalities, or even their approximations. Small concrete examples of such inequalities are presented, and include the famous Blakley-Roy inequalities for paths of odd length. As a consequence, these simple inequalities cannot be written as a rational sum of squares of graph densities.

preprint2022arXiv

Undecidability of polynomial inequalities in weighted graph homomorphism densities

Many problems and conjectures in extremal combinatorics concern polynomial inequalities between homomorphism densities of graphs where we allow edges to have real weights. Using the theory of graph limits, we can equivalently evaluate polynomial expressions in homomorphism densities on kernels $W$, i.e., symmetric, bounded, and measurable functions $W$ from $[0,1]^2 \to \mathbb{R}$. In 2011, Hatami and Norin proved a fundamental result that it is undecidable to determine the validity of polynomial inequalities in homomorphism densities for graphons (i.e., the case where the range of $W$ is $[0,1]$, which corresponds to unweighted graphs, or equivalently, to graphs with edge weights between $0$ and $1$). The corresponding problem for more general sets of kernels, e.g., for all kernels or for kernels with range $[-1,1]$, remains open. For any $a > 0$, we show undecidability of polynomial inequalities for any set of kernels which contains all kernels with range $\{0,a\}$. This result also answers a question raised by Lovász about finding computationally effective certificates for the validity of homomorphism density inequalities in kernels.

preprint2022arXiv

Weak and Strong Extremal Biquadratics

We study quasiconvex quadratic forms on $n \times m$ matrices which correspond to nonnegative biquadratic forms in $(n,m)$ variables. We disprove a conjecture stated by Harutyunyan--Milton (Comm. Pure Appl. Math. 70(11), 2017) as well as Harutyunyan--Hovsepyan (Arch. Ration. Mech. Anal. 244, 2022) that extremality in the cone of quasiconvex quadratic forms on $3\times 3$ matrices can follow only from the extremality of the determinant of its acoustic tensor, using previous work by Buckley--Šivic (Linear Algebra Appl. 598, 2020). Our main result is to establish a conjecture of Harutyunyan--Milton (Comm. Pure Appl. Math. 70(11), 2017) that weak extremal quasiconvex quadratics on $3 \times 3$ matrices are strong extremal. Our main technical ingredient is a generalization of the work of Kunert--Scheiderer on extreme nonnegative ternary sextics (Trans. Amer. Math. Soc. 370(6), 2018). Specifically, we show that a nonnegative ternary sextic, which is not a square, is extremal if and only if its variety (over the complex numbers) is a rational curve and all its singularities are real.

preprint2020arXiv

Sparse PSD approximation of the PSD cone

While semidefinite programming (SDP) problems are polynomially solvable in theory, it is often difficult to solve large SDP instances in practice. One technique to address this issue is to relax the global positive-semidefiniteness (PSD) constraint and only enforce PSD-ness on smaller $k\times k$ principal submatrices --- we call this the sparse SDP relaxation. Surprisingly, it has been observed empirically that in some cases this approach appears to produce bounds that are close to the optimal objective function value of the original SDP. In this paper, we formally attempt to compare the strength of the sparse SDP relaxation vis-à-vis the original SDP from a theoretical perspective. In order to simplify the question, we arrive at a data independent version of it, where we compare the sizes of SDP cone and the $k$-PSD closure, which is the cone of matrices where PSD-ness is enforced on all $k\times k$ principal submatrices. In particular, we investigate the question of how far a matrix of unit Frobenius norm in the $k$-PSD closure can be from the SDP cone. We provide two incomparable upper bounds on this farthest distance as a function of $k$ and $n$. We also provide matching lower bounds, which show that the upper bounds are tight within a constant in different regimes of $k$ and $n$. Other than linear algebra techniques, we extensively use probabilistic methods to arrive at these bounds. One of the lower bounds is obtained by observing a connection between matrices in the $k$-PSD closure and matrices satisfying the restricted isometry property (RIP).

preprint2020arXiv

Sums of Squares and Quadratic Persistence on Real Projective Varieties

We bound the Pythagoras number of a real projective subvariety: the smallest positive integer $r$ such that every sum of squares of linear forms in its homogeneous coordinate ring is a sum of at most $r$ squares. Enhancing existing methods, we exhibit three distinct upper bounds involving known invariants. In contrast, our lower bound depends on a new invariant of a projective subvariety called the quadratic persistence. Defined by projecting away from points, this numerical invariant is closely related to the linear syzygies of the variety. In addition, we classify the projective subvarieties of maximal and almost-maximal quadratic persistence, and determine their Pythagoras numbers.

preprint2020arXiv

Symmetric nonnegative forms and sums of squares

We study symmetric nonnegative forms and their relationship with symmetric sums of squares. For a fixed number of variables $n$ and degree $2d$, symmetric nonnegative forms and symmetric sums of squares form closed, convex cones in the vector space of $n$-variate symmetric forms of degree $2d$. Using representation theory of the symmetric group we characterize both cones in a uniform way. Further, we investigate the asymptotic behavior when the degree $2d$ is fixed and the number of variables $n$ grows. Here, we show that, in sharp contrast to the general case, the difference between symmetric nonnegative forms and sums of squares does not grow arbitrarily large for any fixed degree $2d$. We consider the case of symmetric quartic forms in more detail and give a complete characterization of quartic symmetric sums of squares. Furthermore, we show that in degree $4$ the cones of nonnegative symmetric forms and symmetric sums of squares approach the same limit, thus these two cones asymptotically become closer as the number of variables grows. We conjecture that this is true in arbitrary degree $2d$.

preprint2018arXiv

Generalized eigenvalue methods for Gaussian quadrature rules

A quadrature rule of a measure $μ$ on the real line represents a convex combination of finitely many evaluations at points, called nodes, that agrees with integration against $μ$ for all polynomials up to some fixed degree. In this paper, we present a bivariate polynomial whose roots parametrize the nodes of minimal quadrature rules for measures on the real line. We give two symmetric determinantal formulas for this polynomial, which translate the problem of finding the nodes to solving a generalized eigenvalue problem.

preprint2016arXiv

Do Sums of Squares Dream of Free Resolutions?

We associate to a real projective variety $X$ two convex cones which are fundamental in real algebraic geometry: the cone $P_X$ of quadratic forms nonnegative on $X$, and the cone $Σ_X$ of sums of squares of linear forms. The dual cone $Σ_X^\ast$ is a spectrahedron and we show that its convexity properties are closely related to homological properties of $X$. For instance, we show that all extreme rays of $Σ_X^\ast$ have rank one if and only if X has Castelnuovo-Mumford regularity two. More generally, if $Σ_X^\ast$ has an extreme ray of rank $p > 1$, then $X$ does not satisfy the property $N_{2,p}$. We show that the converse also holds in a wide variety of situations: the smallest $p$ for which property $N_{2,p}$ does not hold is equal to the smallest rank of an extreme ray of $Σ_X^\ast$ greater than one. These results allow us to generalize the work of Blekherman-Smith-Velasco on equality of nonnegative polynomials and sums of squares from irreducible varieties to reduced schemes and to classify all spectrahedral cones with only rank one extreme rays. Our results have applications to the positive semidefinite matrix completion problem and to the truncated moment problem on projective varieties.

preprint2015arXiv

Extreme Rays of the Hankel Spectrahedra for Ternary Forms

Hankel spectrahedra are the dual convex cones to the cone of sums of squares of real polynomials, and we study them from the point of view of convex algebraic geometry. We show that the Zariski closure of the union of all extreme rays of Hankel spectrahedra for ternary forms is an irreducible variety of codimension 10. It is the variety of all Hankel (or middle Catalecticant) matrices of corank at least 4. We explicitly construct a rational extreme ray of maximal rank using the Cayley-Bacharach Theorem for plane curves. We work out the rank stratification of the semi-algebraic set of extreme rays of Hankel spectrahedra in the first three nontrivial cases d = 3, 4, 5. Dually, we get a characterisation of the algebraic boundary of the cone of sums of squares via projective duality theory, extending previous work of Blekherman, Hauenstein, Ottem, Ranestad, and Sturmfels.

preprint2015arXiv

On real typical ranks

We study typical ranks with respect to a real variety $X$. Examples of such are tensor rank ($X$ is the Segre variety) and symmetric tensor rank ($X$ is the Veronese variety). We show that any rank between the minimal typical rank and the maximal typical rank is also typical. We investigate typical ranks of $n$-variate symmetric tensors of order $d$, or equivalently homogeneous polynomials of degree $d$ in $n$ variables, for small values of $n$ and $d$. We show that $4$ is the unique typical rank of real ternary cubics, and quaternary cubics have typical ranks $5$ and $6$ only. For ternary quartics we show that $6$ and $7$ are typical ranks and that all typical ranks are between $6$ and $8$. For ternary quintics we show that the typical ranks are between $7$ and $13$.

preprint2015arXiv

Real Rank with Respect to Varieties

We study the real rank of points with respect to a real variety $X$. This is a generalization of various tensor ranks, where $X$ is in a specific family of real varieties like Veronese or Segre varieties. The maximal real rank can be bounded in terms of the codimension of $X$ only. We show constructively that there exist varieties $X$ for which this bound is tight. The same varieties provide examples where a previous bound of Blekherman-Teitler on the maximal $X$-rank is tight. We also give examples of varieties $X$ for which the gap between maximal complex and the maximal real rank is arbitrarily large. To facilitate our constructions we prove a conjecture of Reznick on the maximal real symmetric rank of symmetric bivariate tensors. Finally we study the geometry of the set of points of maximal real rank in the case of real plane curves.

preprint2014arXiv

Gap vectors of real projective varieties

Let $X\subseteq \mathbb{P}^m$ be a totally real, non-degenerate, projective variety and let $Γ\subseteq X(\mathbb{R})$ be a generic set of points. Let $P$ be the cone of nonnegative quadratic forms on $X$ and let $Σ$ be the cone of sums of squares of linear forms. We examine the dimensions of the faces $P(Γ)$ and $Σ(Γ)$ consisting of forms in $P$ and $Σ$, which vanish on $Γ$. As the cardinality of the set $Γ$ varies in $1,\dots,\rm{codim}(X)$, the difference between the dimensions of $P(Γ)$ and $Σ(Γ)$ defines a numerical invariant of $X$, which we call the gap vector of X. In this article we begin a systematic study of its fundamental properties. Our main result is a formula relating the components of the gap vector of $X$ and the quadratic deficiencies of $X$ and its generic projections. The quadratic deficiency is a fundamental numerical invariant of projective varieties introduced by F. L. Zak. The relationship between quadratic deficiency and gap vectors allows us to effectively compute the gap vectors of concrete varieties as well as to prove several general properties. We prove that gap vectors are weakly increasing, obtain upper bounds for their rate of growth and prove that these upper bounds are eventually achieved for all varieties. Moreover, we give a characterization of the varieties with the simplest gap vectors: We prove that the gap vector vanishes identically precisely for varieties of minimal degree, and characterize the varieties whose gap vector equals $(0,\dots, 0,1)$. In particular, our results give a new proof of the theorem of Blekherman, Smith and Velasco saying that there are nonnegative quadratic forms which are not sums of squares on every variety, which is not of minimal degree. Finally, we determine the gap vector of all Veronese embeddings of $\mathbb{P}^2$, generalizing work of the first three authors.

preprint2014arXiv

Sums of Squares on the Hypercube

Let X be a finite set of points in R^n. A polynomial p nonnegative on X can be written as a sum of squares of rational functions modulo the vanishing ideal I(X). From the point of view of applications, such as polynomial optimization, we are interested in rational function representations of small degree. We derive a general upper bound in terms of the Hilbert function of X, and we show that this upper bound is tight for the case of quadratic functions on the hypercube C={0,1}^n, a very well studied case in combinatorial optimization. Using the lower bounds for C we construct a family of globally nonnegative quartic polynomials, which are not sums of squares of rational functions of small degree. To our knowledge this is the first construction for Hilbert's 17th problem of a family of polynomials of bounded degree which need increasing degrees in rational function representations as the number of variables n goes to infinity. We note that representation theory of the symmetric group S_n play a crucial role in our proofs of the lower bounds.

preprint2013arXiv

Dimensional Differences Between Faces of the Cones of Nonnegative Polynomials and Sums of Squares

We study dimensions of the faces of the cone of nonnegative polynomials and the cone of sums of squares; we show that there are dimensional differences between corresponding faces of these cones. These dimensional gaps occur in all cases where there exist nonnegative polynomials that are not sums of squares. The gaps occur generically, they are not the product of selecting special faces of the cones. For ternary forms and quaternary quartics, we completely characterize when these differences are observed. Moreover, we provide an explicit description for these differences in the two smallest cases, in which the cone of nonnegative polynomials and the cone of sums of squares are different. Our results follow from more general results concerning the relationship between the second ordinary power and the second symbolic power of the vanishing ideal of points in projective space.

preprint2013arXiv

Sums of squares and varieties of minimal degree

Let X be a real nondegenerate projective subvariety such that its set of real points is Zariski dense. We prove that every real quadratic form that is nonnegative on X is a sum of squares of linear forms if and only if X is a variety of minimal degree. This substantially extends Hilbert's celebrated characterization of equality between nonnegative forms and sums of squares. We obtain a complete list for the cases of equality and also a classification of the lattice polytopes Q for which every nonnegative Laurent polynomial with support contained in 2Q is a sum of squares.

preprint2012arXiv

Nonnegative Polynomials and Sums of Squares

In the smallest cases where there exist nonnegative polynomials that are not sums of squares we present a complete explanation of this distinction. The fundamental reason that the cone of sums of squares is strictly contained in the cone of nonnegative polynomials is that polynomials of degree $d$ satisfy certain linear relations, known as the Cayley-Bacharach relations, which are not satisfied by polynomials of full degree 2d. For any nonnegative polynomial that is not a sum of squares we can write down a linear inequality coming from a Cayley-Bacharach relation that certifies this fact. We also characterize strictly positive sums of squares that lie on the boundary of the cone of sums of squares and extreme rays of the cone dual to the cone of sums of squares

preprint2012arXiv

Positive Gorenstein Ideals

We introduce positive Gorenstein ideals. These are Gorenstein ideals in the graded ring $\RR[x]$ with socle in degree 2d, which when viewed as a linear functional on $\RR[x]_{2d}$ is nonnegative on squares. Equivalently, positive Gorenstein ideals are apolar ideals of forms whose differential operator is nonnegative on squares. Positive Gorenstein ideals arise naturally in the context of nonnegative polynomials and sums of squares, and they provide a powerful framework for studying concrete aspects of sums of squares representations. We present applications of positive Gorenstein ideals in real algebraic geometry, analysis and optimization. In particular, we present a simple proof of Hilbert's nearly forgotten result on representations of ternary nonnegative forms as sums of squares of rational functions. Drawing on our previous work, our main tools are Cayley-Bacharach duality and elementary convex geometry.

preprint2012arXiv

Typical Real Ranks of Binary Forms

We prove a conjecture of Comon and Ottaviani that typical real Waring ranks of bivariate forms of degree $d$ take all integer values between $\lfloor \frac{d+2}{2}\rfloor$ and $d$. That is we show that for all $d$ and all $\lfloor \frac{d+2}{2}\rfloor \leq m \leq d$ there exists a bivariate form $f$ such that $f$ can be written as a linear combination of $m$ $d$-th powers of real linear forms and no fewer, and additionally all forms in an open neighborhood of $f$ also possess this property. Equivalently we show that for all $d$ and any $\lfloor \frac{d+2}{2}\rfloor \leq m \leq d$ there exists a symmetric real bivariate tensor $t$ of order $d$ such that $t$ can be written as a linear combination of $m$ symmetric real tensors of rank 1 and no fewer, and additionally all tensors in an open neighborhood of $t$ also possess this property.

preprint2011arXiv

ADAM: Analysis of Discrete Models of Biological Systems Using Computer Algebra

Background: Many biological systems are modeled qualitatively with discrete models, such as probabilistic Boolean networks, logical models, Petri nets, and agent-based models, with the goal to gain a better understanding of the system. The computational complexity to analyze the complete dynamics of these models grows exponentially in the number of variables, which impedes working with complex models. Although there exist sophisticated algorithms to determine the dynamics of discrete models, their implementations usually require labor-intensive formatting of the model formulation, and they are oftentimes not accessible to users without programming skills. Efficient analysis methods are needed that are accessible to modelers and easy to use. Method: By converting discrete models into algebraic models, tools from computational algebra can be used to analyze their dynamics. Specifically, we propose a method to identify attractors of a discrete model that is equivalent to solving a system of polynomial equations, a long-studied problem in computer algebra. Results: A method for efficiently identifying attractors, and the web-based tool Analysis of Dynamic Algebraic Models (ADAM), which provides this and other analysis methods for discrete models. ADAM converts several discrete model types automatically into polynomial dynamical systems and analyzes their dynamics using tools from computer algebra. Based on extensive experimentation with both discrete models arising in systems biology and randomly generated networks, we found that the algebraic algorithms presented in this manuscript are fast for systems with the structure maintained by most biological systems, namely sparseness, i.e., while the number of nodes in a biological network may be quite large, each node is affected only by a small number of other nodes, and robustness, i.e., small number of attractors.

preprint2008arXiv

Experimental study of energy-minimizing point configurations on spheres

In this paper we report on massive computer experiments aimed at finding spherical point configurations that minimize potential energy. We present experimental evidence for two new universal optima (consisting of 40 points in 10 dimensions and 64 points in 14 dimensions), as well as evidence that there are no others with at most 64 points. We also describe several other new polytopes, and we present new geometrical descriptions of some of the known universal optima.

preprint2003arXiv

There are Significantly More Nonnegative Polynomials than Sums of Squares

We investigate the quantitative relationship between nonnegative polynomials and sums of squares of polynomials. We show that if the degree is fixed and the number of variables grows then there are significantly more nonnegative polynomials than sums of squares. More specifically, we take compact bases of the cone of nonnegative polynomials and the cone of sums of squares and derive bounds for the volumes of the bases. If the degree is greater than 2 then we show that the ratio of the volumes of the bases, raised to the power reciprocal to the ambient dimension, tends to 0 as the number of variables tends to infinity.