---
abstract: 'The generalized Lax conjecture asserts that each hyperbolicity cone is a linear slice of the cone of positive semidefinite matrices.  We prove the conjecture for a multivariate generalization of the matching polynomial.  This is further extended (albeit in a weaker sense) to a multivariate version of the independence polynomial for simplicial graphs.  As an application we give a new proof of the conjecture for elementary symmetric polynomials (originally due to Brändén).  Finally we consider a hyperbolic convolution of determinant polynomials generalizing an identity of Godsil and Gutman.' 
address: 'Department of Mathematics, Royal Institute of Technology, SE-100 44 Stockholm, Sweden'
author:
- Nima Amini
title: Spectrahedrality of hyperbolicity cones of multivariate matching polynomials
---

Introduction {#sec::intro}
============

A homogeneous polynomial $h({\mathbf{x}}) \in {\mathbb{R}}[x_1, \ldots, x_n]$ is *hyperbolic* with respect to a vector ${\mathbf{e}}\in {\mathbb{R}}^n$ if $h({\mathbf{e}}) \neq 0$, and if for all ${\mathbf{x}}\in {\mathbb{R}}^n$ the univariate polynomial $t \mapsto h(t{\mathbf{e}}-{\mathbf{x}})$ has only real zeros.  Note that if $h$ is a hyperbolic polynomial of degree $d$, then we may write $$\begin{aligned}
\displaystyle h(t\mathbf{e} - \mathbf{x}) = h(\mathbf{e}) \prod_{j=1}^d (t - \lambda_j(\mathbf{x})), \end{aligned}$$ where $${\lambda_{\rm max}}({\mathbf{x}})=\lambda_1(\mathbf{x}) \geq \cdots \geq \lambda_d(\mathbf{x})={\lambda_{\rm min}}({\mathbf{x}})$$ are called the *eigenvalues* of $\mathbf{x}$ (with respect to $\mathbf{e}$).  The *hyperbolicity cone* of $h$ with respect to $\mathbf{e}$ is the set $\Lambda_{+}(h, \mathbf{e}) = \{ \mathbf{x} \in \mathbb{R}^n : {\lambda_{\rm min}}(\mathbf{x}) \geq 0 \}$. If ${\mathbf{v}}\in \Lambda_{+}(h, {\mathbf{e}})$, then $h$ is hyperbolic with respect to ${\mathbf{v}}$ and $\Lambda_{+}(h,{\mathbf{v}}) = \Lambda_{+}(h, {\mathbf{e}})$. For this reason we usually abbreviate and write $\Lambda_{+}(h)$ if there is no risk for confusion.  We denote by $\Lambda_{++}(h)$ the interior of $\Lambda_{+}(h)$. The cone $\Lambda_{++}(h)$ is convex and can be characterized as the connected component of the set $\{{\mathbf{x}}\in {\mathbb{R}}^n :h({\mathbf{x}}) \neq 0 \}$ containing ${\mathbf{e}}$. These are all facts due to Grding [@Gar]. 

\[detta\] An important example of a hyperbolic polynomial is $\det(X)$, where $X = (x_{ij})_{i,j =1}^n$ is a matrix of variables where we impose $x_{ij} = x_{ji}$. Note that $t \mapsto \det(tI - X)$ where $I = \text{diag}(1,\dots, 1)$, is the characteristic polynomial of a symmetric matrix so it has only real zeros.  Hence $\det(X)$ is a hyperbolic polynomial with respect to $I$, and its hyperbolicity cone is the cone of positive semidefinite matrices. 

Denote the directional derivative of $h({\mathbf{x}}) \in {\mathbb{R}}[x_1, \dots, x_n]$ with respect to ${\mathbf{v}}= (v_1, \dots, v_n)^T \in {\mathbb{R}}^n$ by $$\displaystyle D_{{\mathbf{v}}}h({\mathbf{x}}) = \sum_{k=1}^n v_k \frac{\partial h}{\partial x_k} ({\mathbf{x}}).$$ The following lemma is well-known and essentially follows from the identity $D_{{\mathbf{v}}}h(t) = \frac{d}{dt}h(t {\mathbf{v}}+ {\mathbf{x}}) |_{t=0}$ together with Rolle’s theorem (see [@Gar] [@R]). 

\[derivhyp\] Let $h$ be a hyperbolic polynomial and let ${\mathbf{v}}\in \Lambda_+$ be such that $D_{{\mathbf{v}}}h \not \equiv 0$. Then $D_{{\mathbf{v}}} h$ is hyperbolic with $\Lambda_{+}(h, {\mathbf{v}}) \subseteq \Lambda_{+}(D_{{\mathbf{v}}} h, {\mathbf{v}})$.

A class of polynomials which is intimately connected to hyperbolic polynomials is the class of stable polynomials.  A polynomial $P({\mathbf{x}}) \in \mathbb{C}[x_1, \dots, x_n]$ is *stable* if $P(z_1,\ldots,z_n) \neq 0$ whenever ${{\rm Im}}(z_j)>0$ for all $1\leq j \leq n$. Hyperbolic and stable polynomials are related as follows, see [@PLMS Prop.  1.1]. 

\[hypbas\] Let $P \in \mathbb{R}[x_1, \dots, x_n]$ be a homogenous polynomial.  Then $P$ is stable if and only if $P$ is hyperbolic with $\mathbb{R}_+^n \subseteq \Lambda_+(P)$.

The next theorem which follows (see [@LPR]) from a theorem of Helton and Vinnikov [@HV] proved the Lax conjecture (after Peter Lax $1958$ [@Lax]). 

\[heltonvinnikov\] Suppose that $h(x,y,z)$ is of degree $d$ and hyperbolic with respect to $e = (e_1,e_2,e_3)^T$. Suppose further that $h$ is normalized such that $h(e) = 1$. Then there are symmetric $d \times d$ matrices $A, B, C$ such that $e_1A+e_2B+ e_3C = I$ and $$\begin{aligned}
\displaystyle h(x,y,z) = \det(xA + yB + zC).\end{aligned}$$

The exact analogue of Theorem \[heltonvinnikov\] fails for $n>3$ variables.  This may be seen by comparing dimensions.  The set of polynomials on $\mathbb{R}^n$ of the form $\det(x_1A_1 + \cdots x_nA_n)$ with $A_i$ symmetric for $1 \leq i \leq n$, has dimension at most $n \binom{d}{2}$ whereas the set of hyperbolic polynomials on $\mathbb{R}^n$ has dimension $\binom{n+d-1}{d}$.

A convex cone in ${\mathbb{R}}^n$ is *spectrahedral* if it is of the form $$\begin{aligned}
\displaystyle \left \{ \mathbf{x} \in \mathbb{R}^n : \sum_{i = 1}^n x_i A_i \thickspace \text{ is positive semidefinite} \right \}\end{aligned}$$ where $A_i$, $i = 1, \dots, n$ are symmetric matrices such that there exists a vector $(y_1, \dots, y_n) \in \mathbb{R}^n$ with $\sum_{i=1}^n y_i A_i$ positive definite.  It is easy to see that spectrahedral cones are hyperbolicity cones.  A major open question asks if the converse is true. 

\[glc\] All hyperbolicity cones are spectrahedral. 

An important consequence of Conjecture \[glc\] in the field of optimization is that hyperbolic programming [@R] is the same as semidefinite programming. 

We may reformulate Conjecture \[glc\] as follows, see [@HV; @Vin].  The hyperbolicity cone of $h(\mathbf{x})$ with respect to $\mathbf{e} = (e_1, \dots ,e_n)$ is spectrahedral if there is a homogeneous polynomial $q(\mathbf{x})$ and real symmetric matrices $A_1, \ldots, A_n$ of the same size such that $$\begin{aligned}
\label{qhd}
\displaystyle q(\mathbf{x})h(\mathbf{x}) = \det \left ( \sum_{i = 1}^n x_i A_i \right )\end{aligned}$$ where $\Lambda_{++}(h, \mathbf{e}) \subseteq \Lambda_{++}(q, \mathbf{e})$ and $\sum_{i=1}^n e_iA_i$ is positive definite.  If we can choose $q({\mathbf{x}}) \equiv 1$, then we say that $h({\mathbf{x}})$ admits a *definite determinantal representation*. 

-   Conjecture \[glc\] is true for $n=3$ by Theorem \[heltonvinnikov\],

-   Conjecture \[glc\] is true for homogeneous cones [@Chua], i.e., cones for which the automorphism group acts transitively on its interior,

-   Conjecture \[glc\] is true for quadratic polynomials, see e.g. [@NT],

-   Conjecture \[glc\] is true for elementary symmetric polynomials, see [@B2],

-   Weaker versions of Conjecture \[glc\] are true for smooth hyperbolic polynomials, see [@Kum2; @NS]. 

-   Stronger algebraic versions of Conjecture \[glc\] are false, see [@AB; @B]. 

The paper is organized as follows.  In Section \[sec::matching\] we prove Conjecture \[glc\] for a multivariate generalization of the matching polynomial (Theorem \[matchconespec\]).  We also show that this implies Conjecture \[glc\] for elementary symmetric polynomials (Theorem \[elemspec\]).  Our result may therefore be viewed as a generalization of [@B2].  In Section \[sec::indep\] we generalize further to a multivariate version of the independence polynomial using a recent divisibility relation of Leake and Ryder [@LR] (Theorem \[indepspec\]).  The variables of the homogenized independence polynomial do not fully correspond combinatorially (under the line graph operation) to the more refined homogeneous matching polynomial.  The restriction of Theorem \[indepspec\] to line graphs is therefore weaker than Theorem \[matchconespec\].  Finally, in Section \[sec::conv\] we consider a hyperbolic convolution of determinant polynomials generalizing an identity of Godsil and Gutman [@GG] which asserts that the expected characteristic polynomial of a random signing of the adjacency matrix of a graph is equal to its matching polynomial. 

Unless stated otherwise, $G = (V(G),E(G))$ denotes a simple undirected graph.  We shall adopt the following notational conventions. 

-   $\text{Sym}(S)$ denotes the symmetric group on the set $S$. Write $\mathfrak{S}_n = \text{Sym}([n])$.

-   $N_G[v]$ (resp.  $N_G(v)$) denotes the closed (resp.  open) neighbourhood of $v$.

-   If $S \subseteq V(G)$, then $G[S]$ denotes the subgraph of $G$ induced by $S$.

-   $G \sqcup H$ denotes the disjoint union of the graphs $G$ and $H$.

-   ${\mathbb{R}}^S = \{ (a_s)_{s \in S} : a_s \in {\mathbb{R}}\} \cong {\mathbb{R}}^{|S|}$.

-   ${\mathbb{R}}^G = {\mathbb{R}}^{V(G)} \times {\mathbb{R}}^{E(G)}$.

Hyperbolicity cones of multivariate matching polynomials {#sec::matching}
========================================================

A *$k$-matching* in $G$ is a subset $M \subseteq E(G)$ of $k$ edges, no two of which have a vertex in common.  Let $\mathcal{M}(G)$ denote the set of all matchings in $G$ and let $m(G,k)$ denote the number of $k$-matchings in $G$. By convention $m(G,0) = 1$. We denote by $V(M)$ the set of vertices contained in the matching $M$. If $|V(M)| = |V(G)|$, then we call $M$ a *perfect matching*.  The (univariate) *matching polynomial* is defined by $$\begin{aligned}
\displaystyle \mu(G, t)= \sum_{k \geq 0} (-1)^k m(G,k) t^{|V(G)|-2k}.\end{aligned}$$ Note that this is indeed a polynomial since $m(G,k) = 0$ for $k > \frac{|V(G)|}{2}$. Heilmann and Lieb [@HL] studied the following multivariate version of the matching polynomial with variables $\mathbf{x} = (x_i)_{i \in V}$ and non-negative weights $\boldsymbol{\lambda} = (\lambda_e)_{e \in E}$, $$\begin{aligned}
\displaystyle \mu_{\boldsymbol{\lambda}}(G, \mathbf{x}) = \sum_{M \in \mathcal{M}(G)} (-1)^{|M|} \prod_{ij \in M} \lambda_{ij} x_i x_j.\end{aligned}$$

\[specrem\] Note that $\displaystyle t^{|V(G)|} \mu_{\mathbf{1}}(G, t^{-1} \mathbf{1}) = \mu(G,t)$, where $\mathbf{1} = (1, \dots, 1)$.

\[heilmannlieb\] If $\boldsymbol{\lambda} = (\lambda_e)_{e \in E}$ is a sequence of non-negative edge weights, then $\mu_{\boldsymbol{\lambda}}(G, \mathbf{x})$ is stable. 

A quick way to see Theorem \[heilmannlieb\] is to observe that $$\begin{aligned}
\displaystyle \text{MAP} \left ( \prod_{e = (i,j) \in E(G)} (1- \lambda_{e} x_i x_j) \right ) = \mu_{\boldsymbol{\lambda}}(G, \mathbf{x})\end{aligned}$$ where $\text{MAP}: \mathbb{C}[z_1, \dots, z_n] \to \mathbb{C}[z_1, \dots, z_n]$ is the stability preserving linear map taking a multivariate polynomial to its multiaffine part (see [@BB]).  Since real stable univariate polynomials are real-rooted the Heilmann-Lieb theorem (together with Remark \[specrem\]) implies the real-rootedness of $\mu(G, t)$.

We will consider the following homogeneous multivariate version of the matching polynomial

Let $\mathbf{x} = (x_v)_{v \in V}$ and $\mathbf{w} = (w_e)_{e \in  E}$ be indeterminates.  Define the *homogeneous multivariate matching polynomial* $\mu(G, \mathbf{x} \oplus \mathbf{w} ) \in \mathbb{R}[\mathbf{x}, \mathbf{w}]$ by $$\begin{aligned}
\displaystyle \mu(G, \mathbf{x} \oplus \mathbf{w} ) = \sum_{M \in \mathcal{M}(G)} (-1)^{|M|} \prod_{v \not \in V(M)} x_v \prod_{e \in M} w_e^2.\end{aligned}$$

(x1) at (2,0) ;

(x4) at (-1,0) ;

(x3) at (-2,-2) ;

(x2) at (1,-2) ;

\(e) at (0.3, 0.2) [$e$]{}; (d) at (-1.8, -1) [$d$]{}; (b) at (-0.5, -2.3) [$b$]{}; (a) at (1.7, -1) [$a$]{}; (c) at (0.2, -1.2) [$c$]{};

(x4) – (x1); (x4) – (x3); (x3) – (x2); (x2) – (x1); (x3) – (x1);

\[MatchignEx\] The homogeneous multivariate matching polynomial of the graph $G$ in Figure \[MatchingExFig\] is given by $$\mu(G, {\mathbf{x}}\oplus {\mathbf{w}}) = x_1x_2x_3x_4 - x_3x_4w_a^2 - x_1x_4w_b^2 - x_2x_4w_c^2 - x_1x_2w_d^2 - x_2x_3w_e^2 + w_a^2w_d^2 + w_b^2w_e^2.$$

Note that $\displaystyle \mu(G, t \mathbf{1} \oplus \mathbf{1}) = \mu(G, t)$ and that the perfect matching polynomial is given by $\mu(G, \mathbf{0} \oplus \mathbf{w})$.

In this section we prove Conjecture \[glc\] in the affirmative for the polynomials $\mu(G, \mathbf{x} \oplus \mathbf{w} )$.  We first assert that $\mu(G, \mathbf{x} \oplus \mathbf{w} )$ is indeed a hyperbolic polynomial. 

The polynomial $\mu(G, \mathbf{x} \oplus \mathbf{w} )$ is hyperbolic with respect to $\mathbf{e} = \mathbf{1} \oplus \mathbf{0}$.

Clearly $\mu(G, \mathbf{1} \oplus \mathbf{0} ) = 1 \neq 0$. Let $\mathbf{x} \oplus \mathbf{w} \in {\mathbb{R}}^G$ and $\lambda_e = w_{e}^2$ for all $e \in  E(G)$. Then $$\begin{aligned}
\displaystyle \mu(G, t \mathbf{e} - \mathbf{x} \oplus \mathbf{w} ) = \left ( \prod_{v \in V} (t - x_v) \right ) \mu_{\boldsymbol{\lambda}}(G, (t\mathbf{1} - \mathbf{x})^{-1}).\end{aligned}$$ Since $\mu_{\boldsymbol{\lambda}}(G, \mathbf{x})$ is real stable by Heilmann-Lieb theorem it follows that the right hand side is real-rooted.  Hence $\mu(G, \mathbf{x} \oplus \mathbf{w} )$ is hyperbolic with respect to $\mathbf{e} = \mathbf{1} \oplus \mathbf{0}$.

Analogues of the standard recursions for the univariate matching polynomial (see [@God Thm 1.1]) also hold for $\mu(G, \mathbf{x} \oplus \mathbf{w} )$.  In particular the following recursion is used frequently so we give details. 

\[matchrec\] Let $u \in V(G)$. Then the homogeneous multivariate matching polynomial satisfies the recursion $$\begin{aligned}
\displaystyle \mu(G, \mathbf{x} \oplus \mathbf{w} ) = x_u \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w} ) - \sum_{v \in N(u)} w_{uv}^2 \mu((G \setminus u) \setminus v, \mathbf{x} \oplus \mathbf{w} ).\end{aligned}$$

The identity follows by partitioning the matchings $M \in \mathcal{M}(G)$ into two parts depending on whether $u \in V(M)$ or $u \not \in V(M)$. Let $f_G(M) = \prod_{v \not \in V(M)} x_v \prod_{e \in M} w_e^2$. Then $$\begin{aligned}
\displaystyle \mu(G, \mathbf{x} \oplus \mathbf{w}) &= \sum_{M \in \mathcal{M}(G)} (-1)^{|M|} f_G(M) \\ &= \sum_{\substack{M \in \mathcal{M}(G) \\ u \not \in V(M) }} (-1)^{|M|} f_{G}(M) + \sum_{\substack{ M \in \mathcal{M}(G) \\ u \in V(M)} } (-1)^{|M|} f_G(M) \\ &= x_u \sum_{M \in M(G \setminus u)} (-1)^{|M|} f_{G \setminus u}(M) + \sum_{v \in N(u)} \sum_{\substack{M \in \mathcal{M}(G) \\ uv \in M }} (-1)^{|M|} f_G(M) \\ &= x_u \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w} ) - \sum_{v \in N(u)} w_{uv}^2 \sum_{M \in M( (G \setminus u) \setminus v)} (-1)^{|M|} f_{(G \setminus u)\setminus v}(M) \\ &= x_u \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w} ) - \sum_{v \in N(u)} w_{uv}^2 \thinspace \mu((G \setminus u) \setminus v, \mathbf{x} \oplus \mathbf{w}).\end{aligned}$$

Let $G$ be a graph and $u \in V(G)$. The *path tree* $T(G,u)$ is the tree with vertices labelled by paths in $G$ starting at $u$ and where two vertices are joined by an edge if one vertex is labelled by a maximal subpath of the other. 

[ cc ]{} $G$ & $T(G,1)$\
\

(x1) at (-2.5,0) ;

(x3) at (0,0) ;

(x5) at (2.5,0) ;

(x2) at (-2.5,-2.5) ;

(x4) at (0,-2.5) ;

(x6) at (2.5,-2.5) ;

(x1) – (x3); (x3) – (x5); (x5) – (x6); (x6) – (x4); (x4) – (x2); (x2) – (x1); (x3) – (x4);

& ![ image](PathTreeExSt.pdf){width="65mm"}

\[MatchignExFig\]

Let $G$ be a graph and $u \in V(G)$. Let $\phi: \mathbb{R}^{T(G,u)} \to \mathbb{R}^{G}$ denote the linear change of variables defined by $$\begin{aligned}
\displaystyle x_{p} &
\mapsto x_{i_k}, \\ w_{pp'} &\mapsto w_{i_k i_{k+1}},\end{aligned}$$ where $p = i_1 \cdots i_k$ and $p' = i_1 \cdots i_k i_{k+1}$ are adjacent vertices in $T(G,u)$. For every subforest $T \subseteq T(G,u)$, define the polynomial $$\eta(T, {\mathbf{x}}\oplus {\mathbf{w}}) = \mu(T, \phi({\mathbf{x}}' \oplus {\mathbf{w}}'))$$ where ${\mathbf{x}}' = (x_p)_{p \in V(T)}$ and ${\mathbf{w}}' = (w_e)_{e \in E(T)}$.

Note that $\eta(T, {\mathbf{x}}\oplus {\mathbf{w}})$ is a polynomial in variables ${\mathbf{x}}= (x_v)_{v \in V(G)}$ and ${\mathbf{w}}= (w_e)_{e \in E(G)}$.

For the univariate matching polynomial we have the following rather unexpected divisibility relation due to Godsil [@God2], $$\begin{aligned}
\frac{\mu(G \setminus u, t)}{\mu(G, t)} = \frac{\mu(T(G,u) \setminus u, t)}{\mu(T(G,u), t)}.\end{aligned}$$ Below we prove a multivariate analogue of this fact.  A similar multivariate analogue was also noted independently by Leake and Ryder [@LR].  In fact they were able to find a further generalization to independence polynomials of simplicial graphs.  We will revisit their results in Section \[sec::indep\].  The arguments all closely resemble Godsil’s proof for the univariate matching polynomial.  For the convenience of the reader we provide the details in our setting. 

\[matchquot\] Let $u \in V(G)$. Then $$\begin{aligned}
\displaystyle \frac{\mu(G \setminus u, \mathbf{x} \oplus \mathbf{w} )}{ \mu(G, \mathbf{x} \oplus \mathbf{w})} = \frac{\eta(T(G,u) \setminus u, \mathbf{x} \oplus \mathbf{w})}{ \eta(T(G,u), \mathbf{x} \oplus \mathbf{w})}.\end{aligned}$$

If $G$ is a tree, then $\mu(G, {\mathbf{x}}\oplus {\mathbf{w}}) = \eta(T(G,u), {\mathbf{x}}\oplus {\mathbf{w}})$ and $\mu(G \setminus u, {\mathbf{x}}\oplus {\mathbf{w}}) = \eta(T(G,u) \setminus u, {\mathbf{x}}\oplus {\mathbf{w}})$ so the lemma holds.  In particular the lemma holds for all graphs with at most two vertices.  We now argue by induction on the number of vertices of $G$. We first claim that $$\frac{\eta(T(G,u)\setminus \{ u, uv \}, {\mathbf{x}}\oplus {\mathbf{w}})}{\eta(T(G,u) \setminus u, {\mathbf{x}}\oplus {\mathbf{w}})} = \frac{\eta(T(G \setminus u, v)\setminus v, {\mathbf{x}}\oplus {\mathbf{w}}) }{\eta(T(G \setminus u, v), {\mathbf{x}}\oplus {\mathbf{w}})}.$$ Let $v \in N(u)$. By examining the path tree $T(G,u)$ we note the following isomorphisms $$\begin{aligned}
T(G, u) \setminus u &\cong \bigsqcup_{n \in N(u)} T(G \setminus u, n), \\ T(G, u) \setminus \{ u,uv  \} &\cong \left ( \bigsqcup_{\substack{n \in N(u) \\ n \neq v}} T(G \setminus u, n) \right ) \sqcup T(G \setminus u,v) \setminus v, \end{aligned}$$ following from the fact that $T(G\setminus u, n)$ is isomorphic to the connected component of $T(G,u)\setminus u$ which contains the path $un$ in $G$. By the definition of $\phi$ and the general multiplicative identity $$\begin{aligned}
\displaystyle \mu(G \sqcup H, \mathbf{x} \oplus \mathbf{w}) = \mu(G, \mathbf{x} \oplus \mathbf{w})\mu(H, \mathbf{x} \oplus \mathbf{w}),\end{aligned}$$ the above isomorphisms translate to the following identities $$\begin{aligned}
\displaystyle \eta(T(G, u) \setminus u, {\mathbf{x}}\oplus {\mathbf{w}}) &= \prod_{n \in N(u)} \eta(T(G \setminus u, n), {\mathbf{x}}\oplus {\mathbf{w}}), \\  \eta(T(G,u) \setminus \{u,uv\}, {\mathbf{x}}\oplus {\mathbf{w}}) &= \eta(T(G\setminus u,v) \setminus v, {\mathbf{x}}\oplus {\mathbf{w}}) \prod_{\substack{n \in N(u) \\ n \neq v}} \eta(T(G \setminus u, n), {\mathbf{x}}\oplus {\mathbf{w}}),\end{aligned}$$ from which the claim follows.  By Lemma \[matchrec\], induction, above claim and the definition of $\phi$ we finally get $$\begin{aligned}
\displaystyle \frac{\mu(G, {\mathbf{x}}\oplus {\mathbf{w}})}{\mu(G \setminus u, {\mathbf{x}}\oplus {\mathbf{w}})} &= \frac{x_u \mu(G\setminus u, {\mathbf{x}}\oplus {\mathbf{w}}) - \sum_{v \in N(u)}w_{uv}^2 \mu(G \setminus \{u,v \}, {\mathbf{x}}\oplus {\mathbf{w}})  }{\mu(G \setminus u, {\mathbf{x}}\oplus {\mathbf{w}})} \\ &= x_u - \sum_{v \in N(u)}w_{uv}^2 \frac{\mu((G \setminus u) \setminus  v, {\mathbf{x}}\oplus {\mathbf{w}})}{\mu(G \setminus u, {\mathbf{x}}\oplus {\mathbf{w}})} \\ &= x_u - \sum_{v \in N(u)} w_{uv}^2 \frac{\eta(T(G \setminus u, v) \setminus v, {\mathbf{x}}\oplus {\mathbf{w}})}{\eta(T(G \setminus u, v), {\mathbf{x}}\oplus {\mathbf{w}})} \\ &= x_u - \sum_{v \in N(u)} w_{uv}^2 \frac{\eta(T(G,u)\setminus \{ u, uv \}, {\mathbf{x}}\oplus {\mathbf{w}})}{\eta(T(G,u) \setminus u, {\mathbf{x}}\oplus {\mathbf{w}})} \\ &= \frac{\eta(T(G,u), {\mathbf{x}}\oplus {\mathbf{w}})}{\eta(T(G,u) \setminus u, {\mathbf{x}}\oplus {\mathbf{w}})}\end{aligned}$$ which is the reciprocal of the desired identity. 

\[matchdiv\] Let $u \in V(G)$. Then $\mu(G, {\mathbf{x}}\oplus {\mathbf{w}})$ divides $\eta(T(G,u), {\mathbf{x}}\oplus {\mathbf{w}})$.

The argument is by induction on the number of vertices of $G$. Deleting the root $u$ of $T(G,u)$ we get a forest with $|N(u)|$ disjoint components isomorphic to $T(G \setminus u, v)$ respectively for $v \in N(u)$. This gives $$\begin{aligned}
 \label{proddecomp}
\displaystyle \eta(T(G,u) \setminus u, \mathbf{x} \oplus \mathbf{w}) = \prod_{v \in N(u)} \eta(T(G \setminus u,v), \mathbf{x} \oplus \mathbf{w}).\end{aligned}$$ Therefore $\eta(T(G \setminus u,v), \mathbf{x} \oplus \mathbf{w})$ divides $\eta(T(G,u) \setminus u, \mathbf{x} \oplus \mathbf{w})$ for all $v \in N(u)$. By induction $\mu(G \setminus u, \mathbf{x} \oplus \mathbf{w})$ divides $\eta(T(G \setminus u,v), \mathbf{x} \oplus \mathbf{w})$ for all $v \in N(u)$. Hence $\mu(G \setminus u, \mathbf{x} \oplus \mathbf{w})$ divides $\eta(T(G,u) \setminus u, \mathbf{x} \oplus \mathbf{w})$, so by Lemma \[matchquot\], $\mu(G, \mathbf{x} \oplus \mathbf{w})$ divides $\eta(T(G,u), \mathbf{x} \oplus \mathbf{w})$.

In [@GG] Godsil and Gutman proved the following relationship between the univariate matching polynomial $\mu(G,t)$ of a graph $G$ and the characteristic polynomial $\chi(A,t)$ of its adjacency matrix $A$ $$\begin{aligned}
\displaystyle \chi(A,t) = \sum_{C} (-2)^{\text{comp}(C)} \mu(G \setminus C, t),\end{aligned}$$ where the sum ranges over all subgraphs $C$ (including $C = \emptyset$) in which each component is a cycle of degree $2$ and $\text{comp}(C)$ is the number of connected components of $C$.  In particular if $T$ is a tree, then the only such subgraph is $C = \emptyset$ and therefore $$\begin{aligned}
\displaystyle \chi(A,t) = \mu(T,t).\end{aligned}$$ Next we will derive a multivariate analogue of this relationship for trees. 

\[treedetrep\] Let $T = (V,E)$ be a tree.  Then $\mu(T, \mathbf{x} \oplus \mathbf{w})$ has a definite determinantal representation. 

Let $X = \text{diag}(\mathbf{x})$ and $A = (A_{ij})$ be the matrix $$A_{ij} = \begin{cases} w_{ij} & \text{ if } ij \in E(T) \\ 0 & \text{ otherwise} \end{cases}$$ for all $i,j \in V(T)$. If $\sigma \in \text{Sym}(V(T))$ is an involution (i.e $\sigma^2 = id$), then clearly $A_{j \sigma(j)} = w_{j\sigma(j)} = A_{\sigma(j) \sigma^2(j)}$ since $A$ is symmetric.  Hence by acyclicity of trees we have that $$\begin{aligned}
\displaystyle  \det(X + A) &= \sum_{\sigma \in \text{Sym}(V(T))} \text{sgn}(\sigma) \prod_{i \in V(T)} (X_{i \sigma(i)} + A_{i \sigma(i)}) \\ &=  \sum_{S \subseteq V(T)} \prod_{i \in V(T) \setminus S} x_i \sum_{\substack{\sigma \in \text{Sym}(S) \\ \sigma(j) \neq j \thinspace \forall j \in S \\ \sigma^2 = \text{id}}} \text{sgn}(\sigma) \prod_{j \in S} A_{j \sigma(j)} \\ &= \sum_{S \subseteq V(T)} \prod_{i \in V(T) \setminus S} x_i \sum_{\substack{M \in \mathcal{M}(T[S])\\ M \text{ perfect}}}(-1)^{|M|} \prod_{jk \in M} w_{jk}^2  \\ &= \sum_{M \in \mathcal{M}(T)} (-1)^{|M|} \prod_{i \not \in M} x_i \prod_{jk \in M} w_{jk}^2 && \\ &= \mu(T, \mathbf{x} \oplus \mathbf{w}).\end{aligned}$$ Write $$\begin{aligned}
\displaystyle X + A =  \sum_{i \in V(T)} x_i E_{ii} + \sum_{ij \in E(T)} w_{ij} (E_{ij} + E_{ji}),\end{aligned}$$ where $\{E_{ij} : i,j \in V(T)\}$ denotes the standard basis for the vector space of all real $|V(T)| \times |V(T)|$ matrices.  Evaluated at ${\mathbf{e}}= \mathbf{1} \oplus \mathbf{0}$ we obtain the identity matrix $I$ which is positive definite. 

The proof of Lemma \[treedetrep\] is not dependent on $T$ being connected so the statement remains valid for arbitrary undirected acyclic graphs (i.e. forests). 

We now have all the ingredients to prove our main theorem. 

\[matchconespec\] The hyperbolicity cone of $\mu(G, \mathbf{x} \oplus \mathbf{w})$ is spectrahedral. 

The proof is by induction on the number of vertices of $G$. For the base case we have $\mu(G, {\mathbf{x}}\oplus {\mathbf{w}}) = x_v$, so $\Lambda_{+} = \{ x \in {\mathbb{R}}: x \geq 0 \}$ which is clearly spectrahedral.  Assume $G$ contains more than one vertex.  If $G = G_1 \sqcup G_2$ for some non-empty graphs $G_1,G_2$, then $\Lambda_{++}(\mu(G_i, {\mathbf{x}}, \oplus {\mathbf{w}}))$ is spectrahedral by induction for $i = 1,2$. Therefore $$\begin{aligned}
\displaystyle \Lambda_{++}(\mu(G, {\mathbf{x}}, \oplus {\mathbf{w}})) &= \Lambda_{++}(\mu(G_1 \sqcup G_2, {\mathbf{x}}, \oplus {\mathbf{w}})) \\&= \Lambda_{++} \left (\mu(G_1, {\mathbf{x}}, \oplus {\mathbf{w}})\mu(G_1, {\mathbf{x}}, \oplus {\mathbf{w}}) \right )  \\ &= \Lambda_{++}(\mu(G_1, {\mathbf{x}}, \oplus {\mathbf{w}}))  \cap \Lambda_{++}(\mu(G_2, {\mathbf{x}}, \oplus {\mathbf{w}}))\end{aligned}$$ showing that $\Lambda_{++}(\mu(G, {\mathbf{x}}, \oplus {\mathbf{w}}))$ is spectrahedral.  We may therefore assume $G$ is connected.  Let $u \in V(G)$. Since $G$ is connected and has size greater than one, $N(u) \neq \emptyset$. By Lemma \[matchdiv\] we may define the polynomial $$q_{G,u}({\mathbf{x}}\oplus {\mathbf{w}}) = \frac{\eta(T(G,u), {\mathbf{x}}\oplus {\mathbf{w}})}{\mu(G, {\mathbf{x}}\oplus {\mathbf{w}})}$$ for each graph $G$ and $u \in V(G)$. We want to show that $$\Lambda_{++}(\mu(G, {\mathbf{x}}\oplus {\mathbf{w}})) \subseteq \Lambda_{++}(q_{G,u}({\mathbf{x}}\oplus {\mathbf{w}})).$$ By Lemma \[matchquot\] we have that $$\begin{aligned}
\displaystyle q_{G,u}(\mathbf{x} \oplus \mathbf{w}) \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w} ) = \eta(T(G,u) \setminus u, \mathbf{x} \oplus \mathbf{w}).\end{aligned}$$ Fixing $v \in N(u)$ it follows using (\[proddecomp\]) that $$\begin{aligned}
\displaystyle \frac{q_{G,u}(\mathbf{x} \oplus \mathbf{w}) }{q_{G \setminus u,v}(\mathbf{x} \oplus \mathbf{w})} &=  \frac{q_{G,u}(\mathbf{x} \oplus \mathbf{w}) \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w}) }{q_{G \setminus u,v}(\mathbf{x} \oplus \mathbf{w}) \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w})} \\ &=  \frac{\eta(T(G,u) \setminus u, \mathbf{x} \oplus \mathbf{w} )}{ \eta(T(G \setminus u, v), \mathbf{x} \oplus \mathbf{w})} \\ &= \prod_{w \in N(u) \setminus v} \eta(T(G \setminus u, w), \mathbf{x} \oplus \mathbf{w}) \\ &= \prod_{w \in N(u) \setminus v} q_{G \setminus u,w}(\mathbf{x} \oplus \mathbf{w}) \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w}).\end{aligned}$$ Note that $$\frac{\partial}{\partial x_u} \mu(G, \mathbf{x} \oplus \mathbf{w}) = \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w}).$$ Therefore by Lemma \[derivhyp\], $$\begin{aligned}
\displaystyle \Lambda_{++}(\mu(G, \mathbf{x} \oplus \mathbf{w})) \subseteq \Lambda_{++}(\mu(G \setminus u, \mathbf{x} \oplus \mathbf{w})) \subseteq \Lambda_{++}(q_{G \setminus u,w}(\mathbf{x} \oplus \mathbf{w}))\end{aligned}$$ for all $w \in N(u)$ where the last inclusion follows by inductive hypothesis.  Hence $$\begin{aligned}
\displaystyle \Lambda_{++}(\mu(G, \mathbf{x} \oplus \mathbf{w})) &\subseteq \bigcap_{w \in N(u)} \Lambda_{++} (q_{G \setminus u,w}({\mathbf{x}}\oplus {\mathbf{w}}) \cap \Lambda_{++} \left (\mu(G \setminus u, {\mathbf{x}}\oplus {\mathbf{w}}\right ) \\ &=  \Lambda_{++} \left ( q_{v, G \setminus u}(\mathbf{x} \oplus \mathbf{w}) \prod_{w \in N(u) \setminus v} q_{w, G \setminus u}(\mathbf{x} \oplus \mathbf{w}) \mu(G \setminus u, \mathbf{x} \oplus \mathbf{w})  \right ) \\ &= \Lambda_{++}(q_{G,u}({\mathbf{x}}\oplus {\mathbf{w}})).\end{aligned}$$ Finally by Lemma \[treedetrep\], $\eta(T(G,u), \mathbf{x} \oplus \mathbf{w})$ has a definite determinantal representation.  Hence the theorem follows by induction. 

To show that a hyperbolic polynomial $h$ has a spectrahedral hyperbolicity cone it is by Theorem \[matchconespec\] sufficient to show that $h$ can be realized as a factor of a matching polynomial $\mu(G, \mathbf{x} \oplus \mathbf{w})$ with $\Lambda_{++}(h, \mathbf{e}) \subseteq \Lambda_{++} \left ( \frac{\mu(G, \mathbf{x} \oplus \mathbf{w})}{h}, \mathbf{e} \right )$ (possibly after a linear change of variables). 

\(p) ;

(x1) ;

(x2) at (p.corner 1) ;

(xn) at (p.corner 2) ;

(x6) at (p.corner 3) ;

(x5) at (p.corner 4) ;

(x4) at (p.corner 5) ;

(x3) at (p.corner 6) ;

(anchor1) at (-0.4,0.7); (anchor2) at (-0.8,0);

(w1) at (0.9, 1) [$w_1$]{}; (w2) at (1.5, -0.2) [$w_2$]{}; (w3) at (0.5, -1.4) [$w_3$]{}; (w4) at (-1, -1.2) [$w_4$]{}; (w5) at (-1.5, 0.15) [$w_5$]{}; (wn) at (-0.4, 1.2) [$w_n$]{};

(x1) – (x2); (x1) – (x3); (x1) – (x4); (x1) – (x5); (x1) – (x6); (x1) – (xn); (anchor1) to\[out=140,in=140\] (anchor2);

The *elementary symmetric polynomial* $e_d({\mathbf{x}}) \in {\mathbb{R}}[x_1, \dots, x_n]$ of degree $d$ in $n$ variables is defined by $$e_d({\mathbf{x}}) = \sum_{\substack{S \subseteq [n] \\ |S| = d}} \prod_{i \in S} x_i.$$ The polynomials $e_d({\mathbf{x}})$ are hyperbolic (in fact stable) as a consequence of e.g Grace-Walsh-Szegő theorem (see [@Mar Thm 15.4]). 

The star graph, denoted $S_n$, is given by the complete bipartite graph $K_{1,n}$ with $n+1$ vertices.  As an application of Theorem \[matchconespec\] we show that several well-known instances of hyperbolic polynomials have spectrahedral hyperbolicity cones by realizing them as factors of the multivariate matching polynomial of $S_n$ under some linear change of variables.  With notation as in Figure \[stargraph\], using the recursion in Lemma \[matchrec\], the multivariate matching polynomial of $S_n$ is given by $$\begin{aligned}
\displaystyle \mu(S_n, \mathbf{x} \oplus \mathbf{w}) = \prod_{i = 1}^{n+1} x_i - \sum_{i=1}^n w_i^2 \prod_{\substack{j = 1 \\ j \neq i}}^{n} x_j.  \end{aligned}$$

1.   For $h(\mathbf{x}) = e_{n-1}(\mathbf{x})$ consider the linear change of variables $x_n \mapsto -x_n$ and $w_i \mapsto x_n$ for $i = 1, \dots, n-1$. Then $\mu(S_{n-1}, \mathbf{x} \oplus \mathbf{w}) \mapsto -x_n e_{n-1}(\mathbf{x})$. Clearly $\Lambda_{++}(e_{n-1}(\mathbf{x}), \mathbf{1}) \subseteq \Lambda_{++}(x_n, \mathbf{1})$. The spectrahedrality of $\Lambda_{++}(e_{n-1}({\mathbf{x}}), \mathbf{1})$ was first proved by Sanyal in [@S]. 

2.   For $h(\mathbf{x}) = e_{2}(\mathbf{x})$ consider the linear change of variables $x_i \mapsto e_1(\mathbf{x})$ and $w_i \mapsto x_i$ for $i = 1, \dots, n+1$. Then $\mu(S_{n}, \mathbf{x} \oplus \mathbf{w}) \mapsto 2e_1(\mathbf{x})^{n-1}e_{2}(\mathbf{x})$. Since $D_{\mathbf{1}} e_2({\mathbf{x}}) = (n-1)e_1({\mathbf{x}})$, Lemma \[derivhyp\] implies that $\Lambda_{++}(e_2({\mathbf{x}}), \mathbf{1}) \subseteq \Lambda_{++}(e_1({\mathbf{x}}), \mathbf{1})$. Hence $\Lambda_{++}(e_2({\mathbf{x}}), \mathbf{1})$ is spectrahedral. 

3.   Let $h(\mathbf{x}) = x_n^2 - x_{n-1}^2 - \cdots - x_{1}^2$. Recall that $\Lambda_{++}(h, \mathbf{e})$ is the *Lorentz cone* where $\mathbf{e} = (0, \dots, 0,1)$. Consider the linear change of variables $x_i \mapsto x_n$ and $w_i \mapsto x_{i}$ for $i = 1, \dots, n$. Then $\mu(S_{n-1}, \mathbf{x} \oplus \mathbf{w}) \mapsto x_n^n - \sum_{i = 1}^{n-1} x_{i}^2 x_n^{n-2} = x_n^{n-2}h(\mathbf{x})$. Clearly $\Lambda_{++}(h, \mathbf{e}) \subseteq \Lambda_{++}(x_n^{n-2}, \mathbf{e})$. Hence the Lorentz cone is spectrahedral.  Of course this (and the preceding example) also follow from the fact that all quadratic hyperbolic polynomials have spectrahedral hyperbolicity cone [@NT]. 

Hyperbolicity cones of elementary symmetric polynomials have been studied by Zinchenko [@Z], Sanyal [@S] and Brändén [@B2].  Brändén proved that all hyperbolicity cones of elementary symmetric polynomials are spectrahedral.  As an application of Theorem \[matchconespec\] we give a new proof of this fact using matching polynomials. 

(v1) at (0,0) ;

(v12) at (2.5,0) ;

(v13) at (0,-2.5) ;

(v11) at (0,2.5) ;

(v1n) at (-2.5,0) ;

(p1) at (v12) ;

(v21) at (p1.corner 1) ;

(t1) at (v21) ;

(T21) at (\[xshift=3pt\]t1) [$T_{n-2,k-2}$]{};

(v22) at (p1.corner 8) ;

(t2) at (v22) ;

(T22) at (\[yshift=-3.5pt\]t2) [$T_{n-2,k-2}$]{};

(v23) at (p1.corner 7) ;

(t3) at (v23) ;

(T23) at (\[xshift=-1.5pt\]t3) [$T_{n-2,k-2}$]{};

(e11) at (\[xshift= 8pt, yshift=-35\]v11) [$x_1$]{};

(e12) at (\[xshift=-35pt, yshift=-5\]v12) [$x_2$]{};

(e13) at (\[xshift=8pt, yshift=35\]v13) [$x_{3}$]{};

(e1n) at (\[xshift=35pt, yshift=-5\]v1n) [$x_{n-1}$]{};

(e21) at (\[xshift=-17pt, yshift=-25\]v21) [$x_1$]{};

(e22) at (\[xshift=-20pt, yshift=-14\]v22) [$x_3$]{};

(e23) at (\[xshift=-15pt, yshift=12\]v23) [$x_{n-1}$]{};

(p2) at (v1n) ;

(vn1) at (p2.corner 5) ;

(t4) at (vn1) ;

(vn2) at (\[xshift=-10pt\] p2.corner 4) ;

(t5) at (vn2) ;

(vn3) at (p2.corner 3) ;

(t6) at (vn3) ;

(en1) at (\[xshift= 18pt, yshift=25\]vn1) [$x_1$]{};

(en2) at (\[xshift= 30pt, yshift=15\]vn2) [$x_2$]{};

(en3) at (\[xshift= 20pt, yshift=-15\]vn3) [$x_{n-2}$]{};

(Tn1) at (\[yshift=-3.5pt, xshift=-3pt\]t4) [$T_{n-2,k-2}$]{};

(Tn2) at (\[yshift=1pt, xshift=1pt\]t5) [$T_{n-2,k-2}$]{};

(Tn3) at (\[yshift=-5pt, xshift=0pt\]t6) [$T_{n-2,k-2}$]{};

(p3) at (v11) ;

(v111) at (p3.corner 3) ;

(e111) at (\[xshift= 20pt, yshift=-15\]v111) [$x_2$]{};

(t7) at (v111) ;

(T11) at (\[yshift=-3.5pt, xshift=0pt\]t7) [$T_{n-2,k-2}$]{};

(v112) at (p3.corner 2) ;

(e112) at (\[xshift= 17pt, yshift=-25\]v112) [$x_3$]{};

(t8) at (v112) ;

(T12) at (\[yshift=-2pt, xshift=-2pt\]t8) [$T_{n-2,k-2}$]{};

(v113) at (p3.corner 1) ;

(e113) at (\[xshift= 0pt, yshift=-30\]v113) [$x_{n-1}$]{};

(t9) at (v113) ;

(T13) at (\[yshift=0pt, xshift=2pt\]t9) [$T_{n-2,k-2}$]{};

(p4) at (v13) ;

(v131) at (p4.corner 7) ;

(e131) at (\[xshift= -18pt, yshift=12\]v131) [$x_{1}$]{};

(t10) at (v131) ;

(T31) at (\[yshift=3pt, xshift=-2pt\]t10) [$T_{n-2,k-2}$]{};

(v132) at (p4.corner 6) ;

(e132) at (\[xshift= -16pt, yshift=25\]v132) [$x_{2}$]{};

(t11) at (v132) ;

(T32) at (\[yshift=3pt, xshift=3pt\]t11) [$T_{n-2,k-2}$]{};

(v133) at (p4.corner 5) ;

(e133) at (\[xshift= -3pt, yshift=22\]v133) [$x_{n-1}$]{};

(t12) at (v133) ;

(T33) at (\[yshift=0pt, xshift=-3pt\]t12) [$T_{n-2,k-2}$]{};

(v1) – (v12); (v1) – (v13); (v1) – (v11); (v1) – (v1n);

(v12) – (v21); (v12) – (v22); (v12) – (v23);

(v1n) – (vn1); (v1n) – (vn2); (v1n) – (vn3);

(v11) – (v111); (v11) – (v112); (v11) – (v113);

(v13) – (v131); (v13) – (v132); (v13) – (v133);

plot ([cos()]{}, [-3.3 + sin()]{});

plot ([cos()]{}, [3.3 + sin()]{});

plot ([3.3 + cos()]{}, [sin()]{});

plot ([2\*cos()]{}, [2\*sin()]{});

plot ([-3.3 + cos()]{}, [sin()]{});

\[elemspec\] Hyperbolicity cones of elementary symmetric polynomials are spectrahedral. 

For a subset $S \subseteq [n]$ we shall use the notation $$\begin{aligned}
\displaystyle e_k(S) = \sum_{\substack{T \subseteq S \\ |T| = k}} \prod_{j \in T} x_j.\end{aligned}$$ We show that $e_k(\mathbf{x}) = e_k([n])$ divides the multivariate matching polynomial of the length $k$-truncated path tree $T_{n,k}$ of the complete graph $K_n$ rooted at a vertex $v$ after a linear change of variables.  Let $(C_k)_{k \geq 0}$ denote the real sequence defined by $$\begin{aligned}
\displaystyle C_{0} = 1, \thickspace C_{1} = 1, \thickspace C_k = \prod_{j = 0}^{\lfloor k/2 \rfloor - 1} \frac{k - 2j}{k - 2j-1} \hspace{0.2cm} \text{ for } k \geq 2,\end{aligned}$$ so that $$\begin{aligned}
\displaystyle C_{k}C_{k-1} = k \text{ for all } k \geq 1.\end{aligned}$$ Consider the family $$\begin{aligned}
\displaystyle M_{S,k,i} = \mu(T_{|S|, k}, \phi_{S,k,i}(\mathbf{x} \oplus \mathbf{w})) \end{aligned}$$ of multivariate matching polynomials where $i \in S, \thinspace k \in \mathbb{N}$ and $\phi_{S,k,i}$ is the linear change of variables defined recursively (see Fig \[truncpathtree\]) via

1.   $\phi_{S,0,i}$ is the map $x_v \mapsto e_1(S)$ for all $S \subseteq [n]$ and $i \in S$.

2.   $x_{v} \mapsto L_{S,k,i}$ if $k \geq 1$ where $$L_{S,k,i} = \frac{1}{C_{k-1}}e_1(S \setminus i) + C_k x_i$$ and $x_v$ is the variable corresponding to the root of $T_{n,k}$.

3.   $w_{e_j} \mapsto x_j$ for $j \in S \setminus i$ where $w_{e_j}$ are the variables corresponding to the edges $e_j$ incident to the root of $T_{n,k}$.

4.   For each $j \in S \setminus i$ make recursively the linear substitutions $\phi_{S \setminus i, k-1, j}$ respectively to the variables corresponding to the $j$-indexed copies of the subtrees of $T_{n,k}$ isomorphic to $T_{n-1,k-1}$.

We claim $$\begin{aligned}
\displaystyle
M_{S,0,i} &= e_1(S), \\   
M_{S,k,i} &= \frac{C_k e_k(S)}{e_{k-1}(S \setminus \{i\})} \prod_{j \in S \setminus \{i\}} M_{S \setminus i, k-1, j}\end{aligned}$$ for all $S \subseteq [n], \thinspace i \in S$ and $k \in \mathbb{N}$ by induction on $k$. Clearly $M_{S,0,i} = e_1(S)$ since $\mu(T_{n,0}, \mathbf{x} \oplus \mathbf{w}) = x_v$. By Lemma \[matchrec\] and induction we have $$\begin{aligned}
\displaystyle &M_{S,k,i} \\ &= L_{S,k,i} \prod_{s \in S \setminus \{i \}} M_{S \setminus i, k - 1, s} - \sum_{j \in S \setminus i} x_j^2 \prod_{s \in S \setminus \{i,j\}} M_{S \setminus i, k - 1, s} \thinspace M_{S \setminus \{i,j \}, k - 2, s}
\\ &= \left ( \frac{1}{C_{k-1}}e_1(S \setminus i) + C_kx_i - \sum_{j \in S \setminus \{i\}} x_j^2 \thinspace \frac{e_{k-2}(S \setminus \{i,j \})}{C_{k-1}e_{k-1}(S \setminus i)} \right )  \prod_{s \in S \setminus \{i\}} M_{S \setminus i, k-1, s} \\ &=
 \frac{1}{e_{k-1}(S \setminus i)} \left ( \left (\frac{1}{C_{k-1}}e_1(S \setminus i) + C_k x_i \right ) e_{k-1}(S \setminus i) - \frac{1}{C_{k-1}}\sum_{j \in S \setminus i} x_j^2 e_{k-2}(S \setminus \{i,j \})  \right ) \\ & \times  \prod_{s \in S \setminus \{i\}} M_{S \setminus i, k-1, s} \\ &= \frac{1}{e_{k-1}(S \setminus i)} \left ( \frac{k}{C_{k-1}} e_{k}(S \setminus i) + C_kx_i e_{k-1}(S \setminus i) \right ) \prod_{s \in S \setminus \{i\}} M_{S \setminus i, k-1, s} \\ &=  \frac{C_k e_k(S)}{e_{k-1}(S \setminus i)} \prod_{s \in S \setminus \{i\}} M_{S \setminus i, k-1, s}.\end{aligned}$$ Unwinding the above recursion it follows that $M_{S,k,i}$ is of the form $$\begin{aligned}
\displaystyle M_{S,k,i} = C e_k(S) \prod_{\substack{T \subseteq S \setminus i \\ |T| > |S| - k}} e_{k + |T| - |S|}(T)^{\alpha_T}\end{aligned}$$ for some constant $C$ and exponents $\alpha_T \in \mathbb{N}$ .  Taking $S = [n]$ we thus see that $e_k(\mathbf{x})$ is a factor of the multivariate matching polynomial $M_{[n], k, n}$. It remains to show that $$\begin{aligned}
\displaystyle \Lambda_{++}(e_k(\mathbf{x}), \mathbf{1}) \subseteq \Lambda_{++} \left (\frac{M_{[n],k,n}}{e_k(\mathbf{x})}, \mathbf{1} \right ).\end{aligned}$$ for all $k \leq n$. By Lemma \[derivhyp\] above inclusion follows from the fact that $$\begin{aligned}
\displaystyle \Lambda_{++}(e_k(S), \mathbf{1}) \subseteq \Lambda_{++}(e_{k-1}(S), \mathbf{1})\end{aligned}$$ for all $k \geq 1$ since $D_{\mathbf{1}} e_{k}(S) = (|S|-k)e_{k-1}(S)$, and from the fact that $$\begin{aligned}
\displaystyle \Lambda_{++}(e_k(S), \mathbf{1}) \subseteq \Lambda_{++} \left (e_k(T), \mathbf{1} \right )\end{aligned}$$ for all $T \subseteq S$ since $e_k(T) = \left ( \prod_{i \in S \setminus T} \frac{\partial}{\partial x_i} \right ) e_k(S)$. Hence $\Lambda_{++}(e_k(\mathbf{x}), \mathbf{1})$ is spectrahedral by Theorem \[matchconespec\]. 

Hyperbolicity cones of multivariate independence polynomials {#sec::indep}
============================================================

A subset $I \subseteq V(G)$ is *independent* if no two vertices of $I$ are adjacent in $G$. Let $\mathcal{I}(G)$ denote the set of all independent sets in $G$ and $i(G,k)$ denote the number of independent sets in $G$ of size $k$. By convention $i(G,0) = 1$. The (univariate) *independence polynomial* is defined by $$I(G,t) = \sum_{k \geq 0} i(G,k) t^k.$$ The *line graph* $L(G)$ of $G$ is the graph having vertex set $E(G)$ and where two vertices in $L(G)$ are adjacent if and only if the corresponding edges in $G$ are incident.  It follows that $\mu(G,t) = t^{|V(G)|}I(L(G),-t^{-2})$. Therefore the independence polynomial can be viewed as a generalization of the matching polyomial.  In contrast to the matching polynomial, the independence polynomial of a graph is not real-rooted in general.  However Chudnovsky and Seymour [@CS] proved that $I(G,t)$ is real-rooted if $G$ is claw-free, that is, if $G$ has no induced subgraph isomorphic to the complete bipartite graph $K_{1,3}$. The theorem was later generalized by Engström to graphs with weighted vertices. 

\[engstrom\] Let $G$ be a claw-free graph and $\boldsymbol{\lambda} = (\lambda_v)_{v \in  V(G)}$ a sequence of non-negative vertex weights.  Then the polynomial $$I_{\boldsymbol{\lambda}}(G,t) =  \sum_{I \in \mathcal{I}(G)} \left ( \prod_{v \in I} \lambda_v \right ) t^{|I|}$$ is real-rooted. 

A full characterization of the graphs for which $I(G,t)$ is real-rooted remains an open problem. 

A natural multivariate analogue of the independence polynomial is given by $$I(G, {\mathbf{x}}) = \sum_{I \in \mathcal{I}(G)}  \prod_{v \in I} x_v.$$ Leake and Ryder [@LR] define a strictly weaker notion of stability which they call same-phase stability.  A polynomial $p({\mathbf{z}}) \in {\mathbb{R}}[z_1, \dots, z_n]$ is *(real) same-phase stable* if for every ${\mathbf{x}}\in {\mathbb{R}}_+^n$, the univariate polynomial $p(t {\mathbf{x}})$ is real-rooted.  The authors prove that $I(G, {\mathbf{x}})$ is same-phase stable if and only if $G$ is claw-free.  In fact the same-phase stability of $I(G, {\mathbf{x}})$ is an immediate consequence of Theorem \[engstrom\]. 

The added variables in a homogeneous multivariate independence polynomial should preferably have labels carrying combinatorial meaning in the graph.  For line graphs it is additionally desirable to maintain a natural correspondence with the homogeneous multivariate matching polynomial $\mu(G, {\mathbf{x}}\oplus {\mathbf{w}})$. Unfortunately we have not found a hyperbolic definition that satisfies both of the above properties.  We have thus settled for the following definition. 

Let ${\mathbf{x}}= (x_v)_{v \in V}$ and $t$ be indeterminates.  Define the *homogeneous multivariate independence polynomial* $I(G,{\mathbf{x}}\oplus t) \in {\mathbb{R}}[{\mathbf{x}},t]$ by $$I(G,{\mathbf{x}}\oplus t) = \sum_{I \in \mathcal{I}(G)} (-1)^{|I|} \left ( \prod_{v \in I} x_v^2 \right ) t^{2|V(G)| -2 |I|}.$$

\[indephyp\] If $G$ is a claw-free graph, then $I(G, {\mathbf{x}},t)$ is a hyperbolic polynomial with respect to ${\mathbf{e}}= (0, \dots, 0,1) \in {\mathbb{R}}^{V(G)} \times {\mathbb{R}}$.

First note that $I(G, {\mathbf{e}}) = 1 \neq 0$. Let ${\mathbf{x}}\oplus t \in {\mathbb{R}}^{V(G)} \times {\mathbb{R}}$ and $\lambda_v = x_v^2$ for all $v \in V(G)$. Then $$I(G, s{\mathbf{e}}- {\mathbf{x}}\oplus t) = (s-t)^{2|V(G)|} I_{\boldsymbol{\lambda}}(G, -(s-t)^{-2}).$$ By Theorem \[engstrom\] the polynomial $I_{\boldsymbol{\lambda}}(G, s)$ is real-rooted.  Clearly all roots are negative which implies $I_{\boldsymbol{\lambda}}(G, -s^{-2})$ is real-rooted.  Hence the univariate polynomial $s \mapsto I(G, s{\mathbf{e}}- {\mathbf{x}}\oplus t)$ is real-rooted which shows that $I(G, {\mathbf{x}}\oplus t)$ is hyperbolic with respect to ${\mathbf{e}}$.

An induced clique $K$ in $G$ is called a *simplicial clique* if for all $u \in K$ the induced subgraph $N[u] \cap (G \setminus K)$ of $G \setminus K$ is a clique.  In other words the neighbourhood of each $u \in K$ is a disjoint union of two induced cliques in $G$. Furthermore, a graph $G$ is said to be *simplicial* if $G$ is claw-free and contains a simplicial clique. 

In this section we prove Conjecture \[glc\] for the polynomial $I(G, {\mathbf{x}},t)$ when $G$ is simplicial.  The proof unfolds in a parallel manner to Theorem \[matchconespec\] by considering a different kind of path tree.  Before the results can be stated we must outline the necessary definitions from [@LR]. 

A connected graph $G$ is a *block graph* if each $2$-connected component is a clique.  Given a simplicial graph $G$ with a simplicial clique $K$ we recursively define a block graph $T^{\boxtimes}(G,K)$ called the *clique tree* associated to $G$ and rooted at $K$ (see Figure \[clfig\]). 

We begin by adding $K$ to $T^{\boxtimes}(G,K)$. Let $K_u = N[u] \setminus K$ for each $u \in K$. Attach the disjoint union $\bigsqcup_{u \in K} K_u$ of cliques to $T^{\boxtimes}(G,K)$ by connecting $u \in K$ to every $v \in K_u$. Finally recursively attach $T^{\boxtimes}(G\setminus K,K_u)$ to the clique $K_u$ in $T^{\boxtimes}(G,K)$ for every $u \in K$. Note that the recursion is made well-defined by the following lemma. 

Let $G$ be a clawfree graph and let $K$ be a simplicial clique in $G$. Then $N[u] \setminus K$ is a simplicial clique in $G \setminus K$ for all $u \in K$.

[ cc ]{}

\(p) ;

\(a) at (p.corner 1) ;

\(b) at (p.corner 2) ;

\(c) at (p.corner 3) ;

\(d) at (p.corner 4) ;

\(e) at (p.corner 5) ;

\(f) at (p.corner 6) ;

\(a) – (c); (a) – (b); (b) – (c); (c) – (d); (d) – (e); (e) – (f); (a) – (f);

\(a) – (e); (b) – (d); (b) – (f); (c) – (e); (d) – (f);

\
$G$\
\

(t1) ; (s1) at (t1.corner 1) ; (b1) at (t1.corner 2) ; (a1) at (t1.corner 3) ;

(t2) at (0,2) ; (d3) at (t2.corner 2) ; (f3) at (t2.corner 3) ;

(t3) at (1.75,-1) ; (e2) at (t3.corner 1) ; (f2) at (t3.corner 2) ;

(t4) at (-1.75,-1) ; (e1) at (t4.corner 1) ; (d1) at (t4.corner 3) ;

(fd) at (\[yshift = 0.05pt, xshift=-45pt\]d1) ; (d1) – (fd);

(fe) at (-2.4, -3.2) ; (e1) – (fe);

(de) at (2.4, -3.2) ; (e2) – (de);

(df) at (\[yshift = 0.05pt, xshift=45pt\]f2) ; (f2) – (df);

(ed) at (1.7, 4) ; (d3) – (ed);

(ef) at (-1.7, 4) ; (f3) – (ef);

\
$T^\boxtimes(G, \{a,b,c\})$

It is well-known that a graph is the line graph of a tree if and only if it is a claw-free block graph [@Har Thm 8.5].  In [@LR] it was demonstrated that the block graph $T^{\boxtimes}(G,K)$ is the line graph of a certain *induced path tree* $T^{\angle}(G,K)$. Its precise definition is not important to us, but we remark that it is a subtree of the usual path tree defined in Section \[sec::matching\] that avoids traversed neighbours.  This enables us to find a definite determinantal representation of $I(T^{\boxtimes}(G,K), {\mathbf{x}}\oplus t)$ via Lemma \[treedetrep\].  The second important fact is that $I(G, {\mathbf{x}})$ divides $I(T^{\boxtimes}(G,K), {\mathbf{x}})$ where $T^{\boxtimes}(G,K)$ is relabelled according to the natural graph homomorphism $\phi_K: T^{\boxtimes}(G,K) \to G$. Hence using the recursion provided by the simplicial structure of $G$ we have almost all the ingredients to finish the proof of Conjecture \[glc\] for $I(G, {\mathbf{x}}\oplus t)$.

\[linecltree\] For any simplicial graph $G$, and any simplicial clique $K \leq G$, we have $$L(T^{\angle}(G,K)) \cong T^{\boxtimes}(G,K).$$

The following theorem is a generalization of Godsil’s divisibility theorem for matching polynomials.  It can be proved in a similar manner by induction using the recursive structure of simplicial graphs and removing cliques instead of vertices.  For the proof to go through in the homogeneous setting we must replace the usual recursion by $$I(G, {\mathbf{x}}\oplus t) = t^{2|K|} I(G \setminus K, {\mathbf{x}}\oplus t) - \sum_{v \in K} t^{2|N(v)|} x_v^2 I(G \setminus N[v], {\mathbf{x}}\oplus t).$$

\[Leake-Ryder [@LR])\] \[indepdiv\] Let $K$ be a simplicial clique of the simplicial graph $G$. Then $$\frac{I(G, {\mathbf{x}}\oplus t)}{I(G \setminus K, {\mathbf{x}}\oplus t)} = \frac{I(T^{\boxtimes}(G,K), {\mathbf{x}}\oplus t)}{I(T^{\boxtimes}(G,K) \setminus K, {\mathbf{x}}\oplus t)},$$ where $T^{\boxtimes}(G,K)$ is relabelled according to the natural graph homomorphism $\phi_K: T^{\boxtimes}(G,K) \to G$. Moreover $I(G, {\mathbf{x}}\oplus t)$ divides $I(T^{\boxtimes}(G,K), {\mathbf{x}}\oplus t)$.

The following lemma ensures the hyperbolicity cones behave well under vertex deletion. 

\[hypconedel\] Let $v \in V(G)$. Then $\Lambda_{++}(I(G, {\mathbf{x}}\oplus t)) \subseteq \Lambda_{++}(I(G \setminus v, {\mathbf{x}}\oplus t))$.

Let ${\mathbf{x}}\oplus t \in {\mathbb{R}}^{V(G)} \times {\mathbb{R}}$ and ${\mathbf{e}}= (0, \dots, 0,1)$. By Lemma \[indephyp\] the polynomials $s\mapsto I(G, s{\mathbf{e}}- {\mathbf{x}}\oplus t)$ and $s \mapsto I(G\setminus v, s{\mathbf{e}}- {\mathbf{x}}\oplus t)$ are both real-rooted.  Denote their roots by $\alpha_1, \dots, \alpha_{2n}$ and $\beta_1, \dots, \beta_{2n-2}$ respectively where $n = |V(G)|$. We claim that $$\min_i \alpha_i \leq \min_{i} \beta_i \leq \max_{i} \beta_i \leq \max_i \alpha_i$$ by induction on the number of vertices of $G$. Indeed the claim is vacuously true if $|V(G)| = 1$. Suppose therefore $|V(G)| > 1$. If $G$ is not connected, then $G = G_1 \sqcup G_2$ for some non-empty graphs $G_1,G_2$. Without loss assume $v \in G_1$. Then $G \setminus v = (G_1\setminus v) \sqcup G_2$. By induction the claim holds for the pair $G_1$ and $G_1 \setminus v$. This implies the claim for $G$ and $G \setminus v$ since $I(G, {\mathbf{x}}\oplus t)$ is multiplicative with respect to disjoint union.  We may therefore assume $G$ is connected.  Thus $G \setminus N[v]$ is of strictly smaller size than $G \setminus v$. We have $$\begin{aligned}
 \label{indeprec}
\displaystyle I(G, {\mathbf{x}}\oplus t) = t^2 I(G \setminus v, {\mathbf{x}}\oplus t) - x_v^2 t^{2|N(v)|} I(G \setminus N[v], {\mathbf{x}}\oplus t).\end{aligned}$$ By induction, the maximal root $\gamma$ of $I(G \setminus N[v], s{\mathbf{e}}- {\mathbf{x}}\oplus t)$ is less than the maximal root $\beta$ of $I(G \setminus v, s{\mathbf{e}}- {\mathbf{x}}\oplus t)$. Since $I(G \setminus N[v], s{\mathbf{e}}- {\mathbf{x}}\oplus t)$ is an even degree polynomial with positive leading coefficient we have that $I(G \setminus N[v], s {\mathbf{e}}- {\mathbf{x}}\oplus t) \geq 0$ for all $s \geq \gamma$. By (\[indeprec\]) this implies that $I(G, \beta {\mathbf{e}}- {\mathbf{x}}\oplus t) \leq 0$. Hence $\max_{i} \beta_i \leq \max_{i} \alpha_i$ since $I(G, s {\mathbf{e}}- {\mathbf{x}}\oplus t) \to \infty$ as $s \to \infty$. Since each of the terms involved in the polynomials $I(G, s {\mathbf{e}}- {\mathbf{x}}\oplus t)$ and $I(G \setminus v, s {\mathbf{e}}- {\mathbf{x}}\oplus t)$ have even degree in $s-t$, their respective roots are symmetric about $s = t$. Hence $\min_i \alpha_i \leq \min_i \beta_i$ proving the claim.  Finally if ${\mathbf{x}}_0 \oplus t_0 \in \Lambda_{++}(I(G, {\mathbf{x}}\oplus t))$, then $\min_i \alpha_i > 0$ so by the claim $\min_i \beta_i > 0$ showing that ${\mathbf{x}}_0 \oplus t_0 \in \Lambda_{++}(I(G \setminus v, {\mathbf{x}}\oplus t))$. This proves the lemma. 

Since $$I(G, {\mathbf{x}}\oplus t) \biggr |_{x_v = 0} = t^{2} I(G\setminus v, {\mathbf{x}}\oplus t),$$ we see by Lemma \[hypconedel\] that setting vertex variables equal to zero relaxes the hyperbolicity cone. 

\[indepspec\] If $G$ is a simplicial graph, then the hyperbolicity cone of $I(G, {\mathbf{x}}\oplus t)$ is spectrahedral. 

Let $K$ be a simplicial clique of $G$. Arguing by induction as in Theorem \[matchconespec\], using the clique tree $T^{\boxtimes}(G,K)$ instead of the path tree $T(G,u)$, and invoking Theorem \[indepdiv\] we get a factorization $$\begin{aligned}
 \label{clfac}
\displaystyle q_{G,K}({\mathbf{x}}\oplus t) = q_{G \setminus K, K_v}({\mathbf{x}}\oplus t) \prod_{w \in K \setminus v} q_{G\setminus K, K_w}({\mathbf{x}}\oplus t) I(G \setminus K, {\mathbf{x}}\oplus t),\end{aligned}$$ where $v \in K$ is fixed, $K_w = N[w]\setminus K$ and $$\begin{aligned}
q_{G, K}({\mathbf{x}}\oplus t)I(G, {\mathbf{x}}\oplus t) &= I(T^{\boxtimes}(G,K), {\mathbf{x}}\oplus t), \\
q_{G \setminus K, K_w}({\mathbf{x}}\oplus t) I(G \setminus K, {\mathbf{x}}\oplus t) &= I(T^{\boxtimes}(G \setminus K,K_w), {\mathbf{x}}\oplus t)\end{aligned}$$ for $w \in K$. Repeated application of Lemma \[hypconedel\] gives $$\Lambda_{++}(I(G, {\mathbf{x}}\oplus t)) \subseteq \Lambda_{++}(I(G \setminus K, {\mathbf{x}}\oplus t)).$$ By the factorization (\[clfac\]) and induction we hence get the desired cone inclusion $$\Lambda_{++}(I(G, {\mathbf{x}}\oplus t)) \subseteq \Lambda_{++}(q_{G,K}({\mathbf{x}}\oplus t)).$$ Since $L(T^{\angle}(G,K)) \cong T^{\boxtimes}(G,K)$ by Lemma \[linecltree\] we see that $$I(T^{\boxtimes}(G,K), {\mathbf{x}}\oplus t) = \mu(T^{\angle}(G,K), t\mathbf{1} \oplus {\mathbf{x}}).$$ Hence $I(T^{\boxtimes}(G,K), {\mathbf{x}}\oplus t)$ has a definite determinantal representation by Lemma \[treedetrep\] proving the theorem. 

Convolutions {#sec::conv}
============

If $G$ is a simple undirected graph with adjacency matrix $A = (a_{ij})$, then we may associate a signing ${\mathbf{s}}= (s_{ij}) \in \{\pm 1 \}^{E(G)}$ to its edges.  The symmetric adjacency matrix $A^{{\mathbf{s}}} = (a_{ij}^s)$ of the resulting graph is given by $a_{ij}^{{\mathbf{s}}} = s_{ij} a_{ij}$ for $ij \in E(G)$ and $a_{ij}^{{\mathbf{s}}} = 0$ otherwise.  Godsil and Gutman [@GG2] proved that $$\begin{aligned}
 \label{godsilgutman}
\mathop{\mathbb{E}}_{{\mathbf{s}}\in \{ \pm 1 \}^{E(G)}} \det \left ( tI - A^{{\mathbf{s}}} \right ) = \mu(G,t).\end{aligned}$$ In other words, the expected characteristic polynomial of an independent random signing of the adjacency matrix of a graph is equal to its matching polynomial.  Therefore the expected characteristic polynomial is real-rooted.  This was one of the facts used by Marcus, Spielman and Srivastava [@MSS1] in proving that there exist infinite families of regular bipartite Ramanujan graphs.  Since then, several other families of characteristic polynomials have been identified with real-rooted expectation (see e.g [@MSS3],[@HDS]).  Such families go under the name *interlacing families*, based on the fact that there exists a common root interlacing polynomial if and only if every convex combination of the family is real-rooted.  The method of interlacing families have been successfully applied to other contexts, in particular to the affirmative resolution of the Kadison-Singer problem [@MSS2]. 

In this section we define a convolution of multivariate determinant polynomials and show that it is hyperbolic as a direct consequence of a more general theorem by Brändén [@B3].  In particular this convolution can be viewed as a generalization of the fact that the expectation in (\[godsilgutman\]) is real-rooted.  Namely, we show that the expected characteristic polynomial over any finite set of independent random edge weightings is real-rooted barring certain adjustments to the weights of the loop edges. 

Recall that every symmetric matrix may be identified with the adjacency matrix of an undirected weighted graph (with loops). 

Let $W \subseteq {\mathbb{R}}$ be a finite set and let $A = (a_{ij})_{i,j = 1}^n$ be a real symmetric matrix.  Define a *weighting* ${\mathbf{w}}= (w_{ij})_{i < j} \in W^{\binom{n}{2}}$ of $A$ to be a symmetric matrix $A^{{\mathbf{w}}} = (a_{ij}^{{\mathbf{w}}})$ given by $$a_{ij}^{{\mathbf{w}}} = \begin{cases} w_{ij} a_{ij} & \text{ if } i < j \\ \sum_{k = 1}^i a_{ik} + \sum_{k = i+1}^n w_{ik}^2 a_{ik} & \text{ if } i = j \end{cases}.$$

Let $X = (x_{ij})_{i,j=1}^n$ and $Y = (y_{ij})_{i,j=1}^n$ be symmetric matrices in variables ${\mathbf{x}}= (x_{ij})_{i \leq j}$ and ${\mathbf{y}}= (y_{ij})_{i \leq j}$ respectively.  Let $W \subseteq {\mathbb{R}}$ be a finite set.  We define the convolution $$\begin{aligned}
\displaystyle \det(X) \operatorname{\ast}_{W} \det(Y) =  \mathop{\mathbb{E}}_{{\mathbf{w}}_1, {\mathbf{w}}_2 \in W^{\binom{n}{2}}} \det( X^{{\mathbf{w}}_1} + Y^{{\mathbf{w}}_2}) \in \mathbb{R}[{\mathbf{x}}, {\mathbf{y}}].\end{aligned}$$

We have the following general fact about hyperbolic polynomials. 

\[compatiblefamilies\] Let $h({\mathbf{x}})$ be a hyperbolic polynomial with respect to ${\mathbf{e}}\in \mathbb{R}^n$, let $V_1, \dots, V_m$ be finite sets of vectors of rank at most one in $\Lambda_+$. For $\mathbf{V} = ({\mathbf{v}}_1, \dots, {\mathbf{v}}_m) \in V_1 \times \cdots \times V_m$, let $$\begin{aligned}
\displaystyle g(\mathbf{V};t) = h(t{\mathbf{e}}+ {\mathbf{u}}-\alpha_1 {\mathbf{v}}_1 - \cdots - \alpha_m {\mathbf{v}}_m)\end{aligned}$$ where ${\mathbf{u}}\in \mathbb{R}^n$ and $(\alpha_1, \dots, \alpha_m) \in \mathbb{R}^{m}$. Then $\displaystyle \mathop{\mathbb{E}}_{\mathbf{V} \in V_1 \times \cdots \times V_m} \thinspace g(\mathbf{V};t)$ is real-rooted. 

\[convhyp\] Let $W \subseteq \mathbb{R}$ be a finite subset.  Then $\det(X) \operatorname{\ast}_W \det(Y)$ is hyperbolic with respect to ${\mathbf{e}}= I \oplus \mathbf{0}$ where $I$ denotes the identity matrix. 

Let $h(X \oplus Y) = \det(X) \operatorname{\ast}_{W} \det(Y)$. We note that $h({\mathbf{e}}) = 1 \neq 0$. Let $\delta_1, \dots, \delta_n$ denote the standard basis of $\mathbb{R}^n$. Put $$V_{ij} = \{ {\mathbf{v}}_{ijw} : w \in W \}$$ where ${\mathbf{v}}_{ijw} =  (\delta_{i} + w \delta_j) (\delta_{i} + w \delta_j)^{T}$ for $i < j$ and $w \in W$. Note that ${\mathbf{v}}_{ijw}$ is a rank one matrix belonging to the hyperbolicity cone of positive semidefinite matrices (with non-zero eigenvalue $w^2+1$).  Letting ${\mathbf{u}}= \mathbf{0}$ and $\alpha_{ij}^X = x_{ij}, \alpha_{ij}^Y = y_{ij}$ for $i < j$ we see that $$\begin{aligned}
\displaystyle h(t {\mathbf{e}}- X \oplus Y) &= \mathop{\mathbb{E}}_{{\mathbf{w}}_1, {\mathbf{w}}_2} \det(tI - X^{{\mathbf{w}}_1} - Y^{{\mathbf{w}}_2}) \\ &= \mathop{\mathbb{E}}_{\substack{{\mathbf{v}}_{ijw_1}, {\mathbf{v}}_{ijw_2} \in V_{ij} \\ i < j}} \det \left (tI + {\mathbf{u}}- \sum_{i < j} (\alpha_{ij}^X {\mathbf{v}}_{ijw_1} +  \alpha_{ij}^Y {\mathbf{v}}_{ijw_2}) \right ),\end{aligned}$$ where the right hand side is a real-rooted polynomial in $t$ by Theorem \[compatiblefamilies\].  Hence $\det(X) \operatorname{\ast}_W \det(Y)$ is hyperbolic with respect to ${\mathbf{e}}$.

\[diagrmk\] Taking $W = \{ \pm 1 \}$ we see that the diagonal adjustment in the weighting is constant.  Therefore setting ${\mathbf{u}}= \text{diag}(d_1, \dots, d_n)$ where $d_i = \sum_{j \neq i} (x_{ij} + y_{ij})$ in the proof of Corollary \[convhyp\], we get that $$\begin{aligned}
 \label{sgnconv}
\displaystyle \mathop{\mathbb{E}}_{{\mathbf{s}}_1, {\mathbf{s}}_2} \det(X^{{\mathbf{s}}_1} + Y^{{\mathbf{s}}_2})\end{aligned}$$ is hyperbolic, where the expectation is taken over independent random signings of the matrices $X$ and $Y$ as in (\[godsilgutman\]) without diagonal adjustment.  This shows in particular that the expectation in (\[godsilgutman\]) is real-rooted. 

Let $W \subseteq {\mathbb{R}}$ be a finite subset and $A$ a real symmetric $n \times n$ matrix.  Then $$\begin{aligned}
\displaystyle \mathop{\mathbb{E}}_{{\mathbf{w}}\in W^{\binom{n}{2}}} \det(tI - A^{{\mathbf{w}}})\end{aligned}$$ is real-rooted. 

By Corollary \[convhyp\] the polynomial $\det(Y) \operatorname{\ast}_W \det(X)$ is hyperbolic, so in particular $\displaystyle t \mapsto \mathop{\mathbb{E}}_{{\mathbf{w}}} \det(tI - A^{{\mathbf{w}}})$ is real-rooted with $X = \mathbf{0}$ and $Y = A$. 

Next we see that the convolution (\[sgnconv\]) over independent random signings can be realized as a convolution of multivariate matching polynomials.  The proof is similar to that of the univariate identity (\[godsilgutman\]) (cf [@GG2]).  Let $G_X$ and $G_Y$ denote the weighted graphs corresponding to the symmetric matrices $X$ and $Y$.

\[boxplus\] Let $X = (x_{ij})_{i,j = 1}^n$ and $Y = (y_{ij})_{i,j = 1}^n$ be symmetric matrices in variables ${\mathbf{x}}= (x_{ij})_{i \leq j}$ and ${\mathbf{y}}= (y_{ij})_{i \leq j}$. Then $$\begin{aligned}
\displaystyle &\mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}, {\mathbf{s}}^{(2)}} \det(X^{{\mathbf{s}}^{(1)}} + Y^{{\mathbf{s}}^{(2)}}) \\ &= \sum_{S \subseteq [n]} (-1)^{|S|/2} \prod_{i \not \in S} (x_{ii} + y_{ii}) \sum_{S_1 \sqcup S_2 = S} \mu(G_X[S_1], \mathbf{0} \oplus \mathbf{x}) \mu(G_Y[S_2], \mathbf{0} \oplus \mathbf{y})\end{aligned}$$ where the expectation is taken over independent random signings as in (\[godsilgutman\]). 

Expanding the convolution from the definition of the determinant we have $$\begin{aligned}
\displaystyle  &\mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}, {\mathbf{s}}^{(2)}} \det( X^{{\mathbf{s}}^{(1)}} + Y^{{\mathbf{s}}^{(2)}}) \\  &= \mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}, {\mathbf{s}}^{(2)}} \sum_{\sigma \in \mathfrak{S}_{n}} \text{sgn}(\sigma) \prod_{i = 1}^{n} \left (X^{{\mathbf{s}}^{(1)}} + Y^{{\mathbf{s}}^{(2)}} \right )_{i \sigma(i)} \\ & = \mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}, {\mathbf{s}}^{(2)}} \sum_{S \subseteq [n]} \prod_{i \not \in S} (x_{ii} + y_{ii}) \sum_{\substack{\sigma \in \text{Sym}(S) \\ \sigma(j) \neq j \thickspace \forall j \in S}} \text{sgn}(\sigma) \prod_{j \in  S} \left (  s_{j \sigma(j)}^{(1)}  x_{j \sigma(j)} + s_{j \sigma(j)}^{(2)} y_{j \sigma(j)} \right ) \\ &= 
\sum_{S \subseteq [n]} \prod_{i \not \in S} (x_{ii} + y_{ii}) \sum_{\substack{\sigma \in \text{Sym}(S) \\ \sigma(j) \neq j \thickspace \forall j \in S}} \text{sgn}(\sigma) \sum_{S_1 \sqcup S_2 = S} \mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}} \prod_{j \in  S_1} s_{j \sigma(j)}^{(1)} x_{j \sigma(j)} \mathop{\mathbb{E}}_{{\mathbf{s}}^{(2)}} \prod_{j \in  S_2} s_{j \sigma(j) }^{(2)} y_{j \sigma(j)}\end{aligned}$$ Note the following regarding the random variables $s_{ij}^{(k)}$, $k = 1,2$:

1.   $s_{ij}^{(k)}$ appears with power at most two in each of the products. 

2.   The random variables $s_{ij}^{(k)}$ are independent. 

3.   $\mathop{\mathbb{E}} s_{ij}^{(k)}  = 0$.

4.   $\mathop{\mathbb{E}} (s_{ij}^{(k)})^2 = 1$.

As a consequence, permutations with the following characteristics may be eliminated since they produce factors $s_{ij}^{(k)}$ of power one making the term vanish:

1.   $\sigma \in \mathfrak{S}_n$ having no factorization $\sigma = \sigma_1 \sigma_2$ for $\sigma_i \in \text{Sym}(S_i)$, $i = 1,2$.

2.   $\sigma \in \mathfrak{S}_n$ such that $\sigma$ is not a complete product of disjoint transpositions. 

This leaves us with products of fixed-point-free involutions in $\text{Sym}(S_1)$ and $\text{Sym}(S_2)$. Thus the non-vanishing terms are those corresponding to perfect matchings on $G_X[S_1]$ and $G_Y[S_2]$. Hence $$\begin{aligned}
\displaystyle \mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}, {\mathbf{s}}^{(2)}} \det( X^{{\mathbf{s}}^{(1)}} + Y^{{\mathbf{s}}^{(2)}}) = \sum_{S \subseteq [n]} \prod_{i \not \in S} (x_{ii} + y_{ii}) \sum_{S_1 \sqcup S_2 = S} P_1({\mathbf{x}}) P_2({\mathbf{y}})\end{aligned}$$ where $$\begin{aligned}
\displaystyle P_1({\mathbf{x}}) &= \sum_{\substack{ \sigma_1 \in \text{Sym}(S_1) \\ \sigma_1(j) \neq j \thickspace \forall j \in S_1}} \text{sgn}(\sigma_1) \mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}} \prod_{i \in S_1} s_{i \sigma_1(i) }^{(1)} x_{i \sigma_1(i)} \\ &= \sum_{ \substack{ M \in \mathcal{M}(G_X[S_1]) \\ M \text{ perfect }}} (-1)^{|S_1|/2} \prod_{ij \in M} \mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}} (s_{ i j }^{(1)})^2 x_{i j}^2 \\ &= (-1)^{|S_1|/2} \sum_{\substack{ M \in \mathcal{M}(G_X[S_1]) \\ M \text{ perfect }}} \prod_{ij \in M} x_{ij}^2 \\ &= (-1)^{|S_1|/2} \mu(G_X[S_1], \mathbf{0} \oplus \mathbf{x})\end{aligned}$$ and similarly for $P_2({\mathbf{y}})$.

The expression in Proposition \[boxplus\] may also be written $$\begin{aligned}
\displaystyle \displaystyle \mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}, {\mathbf{s}}^{(2)}} \det( X^{{\mathbf{s}}^{(1)}} + Y^{{\mathbf{s}}^{(2)}}) = \sum_{M \in \mathcal{M}(K_n)} (-1)^{|M|} \prod_{i \not \in V(M)} (x_{ii} + y_{ii}) \prod_{jk \in M} (x_{jk}^2 + y_{jk}^2).\end{aligned}$$

1.   Let $A$ be the adjacency matrix of a simple undirected graph $G$. Under the specialization $X = tI$ and $Y = -A$ in Proposition \[boxplus\] we recover the identity (\[godsilgutman\]) of Godsil and Gutman. 

2.   Let $A$ and $B$ both be adjacency matrices of the complete graph $K_n$. It is well-known (see e.g. [@God]) that the number of perfect matchings in $K_n$ is given by $(n-1)!!$ if $n$ is even and $0$ otherwise, where $(n)!!  = n(n-2)(n-4)\cdots$. By Proposition \[boxplus\] and a simple calculation it follows that $$\begin{aligned}
    \displaystyle \mathop{\mathbb{E}}_{{\mathbf{s}}^{(1)}, {\mathbf{s}}^{(2)}} \det( tI + A^{{\mathbf{s}}^{(1)}} + B^{{\mathbf{s}}^{(2)}}) &= \sum_{k = 0}^{\lfloor n/2 \rfloor} t^{n-2k} (-1)^k \binom{n}{2k} \sum_{i+j = k} \binom{2k}{2i} (2i-1)!!(2j-1)!! \\ & = \sum_{k = 0}^{\lfloor n/2 \rfloor} t^{n-2k} (-1)^k \binom{n}{2k} (2k-1)!!  \left ( \frac{3}{2} \right )^k \\ &= t^n \mu_{\frac{3}{2} \mathbf{1}}(K_n, t^{-1}\mathbf{1}).\end{aligned}$$

Final remarks {#final-remarks .unnumbered}
=============

In Theorem \[indepspec\] we proved Conjecture \[glc\] for $I(G, {\mathbf{x}}\oplus t)$ whenever $G$ is a simplicial graph.  An extension of the divisibility relation in Theorem \[indepdiv\] to all claw-free graphs would immediately extend Theorem \[indepspec\] to all claw-free graphs. 

An interesting extension of this work would be to study a family of stable graph polynomials introduced by Wagner [@Wag] in a general effort to prove Heilmann-Lieb type theorems.  Let $G = (V,E)$ be a graph.  For $H \subseteq E$, let $\deg_H: V \to \mathbb{N}$ denote the degree function of the subgraph $(V,H)$. Furthermore let $$\begin{aligned}
\displaystyle {\mathbf{u}}^{(v)} = (u_0^{(v)}, u_1^{(v)}, \dots, u_d^{(v)})\end{aligned}$$ denote a sequence of *activities* at each vertex $v \in V$ where $d = \deg_G(v)$. Define the polynomial $$\begin{aligned}
\displaystyle Z(G, {\boldsymbol{\lambda}}, {\mathbf{u}}; {\mathbf{x}}) = \sum_{H \subseteq E} (-1)^{|H|} {\boldsymbol{\lambda}}^H {\mathbf{u}}_{\deg_H} {\mathbf{x}}^{\deg_H}\end{aligned}$$ where ${\boldsymbol{\lambda}}= \{\lambda_e\}_{e \in  E}$ are edge weights and $$\begin{aligned}
\displaystyle {\boldsymbol{\lambda}}^H = \prod_{e \in H} \lambda_e, \hspace{0.3cm} {\mathbf{u}}_{\deg_H} = \prod_{v \in V} u_{\deg_H(v)}^{(v)}, \hspace{0.3cm} {\mathbf{x}}^{\deg_H} = \prod_{v \in .V} x_v^{\deg_H(v)}.\end{aligned}$$ Wagner proves that $Z(G, \boldsymbol{\lambda}, {\mathbf{u}}, {\mathbf{x}})$ is stable whenever $\lambda_e \geq 0$ for all $e \in E$ and the univariate *key-polynomial* $K_v(z) = \sum_{j=0}^d \binom{d}{j} u_j^{(v)} z^j$ is real-rooted for all $v \in V$ (cf [@Wag Thm 3.2]).  We note in particular that if $u_0^{(v)} = u_1^{(v)} = 1$, $u_k^{(v)} = 0$ for all $k > 1$ and $\thickspace v \in V$, then $Z(G, {\boldsymbol{\lambda}}, {\mathbf{u}}; {\mathbf{x}}) = \mu_{{\boldsymbol{\lambda}}}(G,{\mathbf{x}})$ where $\mu_{{\boldsymbol{\lambda}}}(G,{\mathbf{x}})$ is the weighted multivariate matching polynomial studied by Heilmann and Lieb [@HL].  An appropriate homogenization of $Z(G, {\boldsymbol{\lambda}}, {\mathbf{u}}; {\mathbf{x}})$ could be defined as $$W(G, \mathbf{u}; {\mathbf{x}}\oplus {\mathbf{w}}) = \sum_{H \subseteq E} (-1)^{|H|} {\mathbf{u}}_{\deg_H} {\mathbf{w}}^{2H} {\mathbf{x}}^{\deg_G - \deg_H}.$$ Since $W(G, \mathbf{u}; {\mathbf{x}}\oplus {\mathbf{w}}) = {\mathbf{x}}^{\deg_G} Z(G, {\mathbf{w}}^2, {\mathbf{u}}; {\mathbf{x}}^{-1} )$ we see that $W(G, \mathbf{u}; {\mathbf{x}}\oplus {\mathbf{w}})$ is hyperbolic with respect to ${\mathbf{e}}= \mathbf{1} \oplus \mathbf{0}$ whenever $K_v(z)$ is real-rooted for all $v \in  V$. We also note the following edge and node recurrences for $e \in E$ and $v \in V$, $$\begin{aligned}
&W(G, {\mathbf{u}}; {\mathbf{x}}\oplus {\mathbf{w}}) \\ &= {\mathbf{x}}^e W(G \setminus e, {\mathbf{u}}; {\mathbf{x}}\oplus {\mathbf{w}}) - {\mathbf{w}}^{2e} W(G \setminus e, {\mathbf{u}}\ll e; {\mathbf{x}}\oplus {\mathbf{w}}) \\
&= \sum_{S \subseteq N(v)} (-1)^{|S|} u_{|S|}^{(v)} {\mathbf{w}}^{2E(S,v)} x_v^{\deg_G(v) - |S|}{\mathbf{x}}^{N(v) \setminus S}W(G \setminus v,{\mathbf{u}}\ll S; {\mathbf{x}}\oplus {\mathbf{w}})\end{aligned}$$ where $E(S,v) = \{ sv \in E : s \in S \}$ and $\displaystyle ({\mathbf{u}}\ll S)^{(v)} = \begin{cases} (u_1^{(v)}, \dots, u_d^{(v)}), & v \in S \\ {\mathbf{u}}^{(v)}, & v \not \in S \end{cases}$

Although it is not clear in general how to find a definite determinantal representation of $W(G, \mathbf{u}; {\mathbf{x}}\oplus {\mathbf{w}})$, it may be possible to consider special form activity vectors and obtain a reduction by constructing divisibility relations in the spirit of Lemma \[matchquot\] and Theorem \[indepdiv\].  This may also be of independent interest for studying root bounds of their univariate specializations. ** Acknowledgements.**  The author is grateful to Petter Brändén. 

[99]{} N. Amini, P. Brändén, Non-representable hyperbolic matroids, arXiv:1512.05878 (2015).  J. Borcea, P. Brändén, [The Lee-Yang and Pólya-Schur programs.  II.  Theory of stable polynomials and applications,]{} Comm.  Pure Appl.  Math. [** 62**]{} (2009), no.  12, 1595–1631.  J. Borcea, P. Brändén, Multivariate Pólya–Schur classification problems in the Weyl algebra, Proc.  Lond.  Math.  Soc. [** 101**]{} (2010), 73–104.  P. Brändén, [Obstructions to determinantal representability,]{} Adv.  Math. [** 226**]{} (2011), 1202-1212.  P. Brändén, [Hyperbolicity cones of elementary symmetric polynomials are spectrahedral,]{} Optim.  Lett. [** 8**]{} (2014), 1773-1782.  P. Brändén, [Hyperbolic polynomials and the Marcus-Spielman-Srivastava theorem,]{} arXiv:1412.0245 (2014) P. Brändén, R. S. Gonzáles D’Léon, [On the half-plane property and the Tutte group of a matroid,]{} J. Combin.  Theory Ser.  B 100. [** 5**]{} (2010), 485-492.  Y. Choe, J. Oxley, A. Sokal, D. Wagner, [Homogeneous multivariate polynomials with the half-plane property,]{} Adv.  in Appl.  Math. [** 32**]{} (2004), no 1-2, 88-187.  C. B. Chua, Relating homogeneous cones and positive definite cones via T -algebras, SIAM J. Optim. [** 14**]{} (2003), 500–506.  M. Chudnovsky, P. Seymour, [The roots of the independence polynomial of a clawfree graph]{}, J. Combin.  Theory Ser.  B [**97**]{} (2007), 350–357.  A. Engström, [Inequalities on well-distributed point sets on circles,]{} JIPAM J. Inequal.  Pure Appl.  Math. [** 8**]{} (2007), no.2, Article 34, 5pp.  C. Godsil, [Matchings and walks in graphs,]{} Journal of Graph Theory [**5**]{} (1981), 285–-297.  C. Godsil, [Algebraic Combinatorics,]{} Chapman and Hall, (1993).  C. Godsil, I. Gutman [On the theory of the matching polynomial,]{} J. Graph Theory [**5**]{} (1981), 137–144.  C. Godsil, I. Gutman [On the matching polynomial of a graph,]{} Algebraic Methods in graph theory, volume I of Colloquia Mathematica Societatis János Bolyai [**25**]{} (1981), 241–249.  L. Gurvits, [Combinatorial and algorithmic aspects of hyperbolic polynomials,]{} arXiv:math/0404474 (2005).  L. Grding, [An inequality for hyperbolic polynomials,]{} J. Math.  Mech [**8**]{} (1959), 957-965.  C. Hall, D. Puder, W. F. Sawin [Ramanujan coverings of graphs,]{} arXiv:1506.02335 (2015).  F. Harary, [Graph Theory,]{} Addison-Wesley, Massachusetts, (1972).  Commun.  Pure Appl.  Math. [** 11**]{} (1958), 175–194.  O. J. Heilmann, E. H. Lieb [Theory of monomer-dimer systems,]{} Comm.  Math.  Phys. [** 25**]{} (1972), 190-232.  J. W. Helton, V. Vinnikov, [Linear matrix inequality representation of sets,]{} Comm.  Pure Appl.  Math. [** 60**]{} (2007), 654-674. 

M. Kummer, [A note on the hyperbolicity cone of the specialized Vámos polynomial,]{} arXiv:1306.4483 (2013).  M. Kummer, Determinantal Representations and Bézoutians, arXiv:1308.5560 (2013).  M. Kummer, D. Plaumann, C. Vinzant, [Hyperbolic polynomials, interlacers, and sums of squares,]{} arXiv:1212.6696 (2013).  P. Lax, Differential equations, difference equations and matrix theory J. Leake, N. Ryder, [Generalizations of the Matching Polynomial to the Multivariate Independence Polynomial,]{} arXiv:1610.00805 (2016).  A. Lewis, P. Parrilo, M. Ramana, The Lax conjecture is true, Proc.  Amer.  Math.  Soc. [** 133**]{} (2005), 2495–2499.  A. W.  Marcus, D. A. Spielman, N. Srivastava [Interlacing families I: Bipartite Ramanujan graphs of all degrees,]{} Annals of Mathematics, [**182**]{} (2015).  A. W.  Marcus, D. A. Spielman, N. Srivastava [Interlacing families II: Mixed characteristic polynomials and the Kadison–Singer problem,]{} Annals of Mathematics, [**182**]{} (2015).  A. W.  Marcus, D. A. Spielman, N. Srivastava [Finite free convolutions of polynomials,]{} arXiv:1504.00350 (2015).  M. Marden, Geometry of Polynomials, Amer.  Math.  Soc.,  Providence, RI (1966).  T. Netzer, R. Sanyal, [Smooth hyperbolicity cones are spectrahedral shadows,]{} arXiv:1208.0441 (2012).  T. Netzer, A. Thom, Polynomials with and without determinantal representations, Linear Algebra Appl. [** 437**]{} (2012), 1579–1595.  J. Renegar, [Hyperbolic programs, and their derivative relaxations,]{} Found.  Comput.  Math., [** 6**]{} (2006), 59-79.  R.Sanyal, [On the derivative cones of polyhedral cones,]{} Adv.  Geom. [** 13**]{} (2013), 315-321.  V. Vinnikov, LMI representations of convex semialgebraic sets and determinantal representations of algebraic hypersurfaces: past, present, and future, Mathematical methods in systems, optimization, and control, 325–349, Oper.  Theory Adv.  Appl., [** 222**]{}, Birkhauser/Springer Basel AG, Basel, (2012).  D. G. Wagner, [Weighted enumeration of spanning subgraphs with degree constraints,]{} J. Combin.  Theory Ser.  B [**99**]{} (2009), 347-357.  Y.Zinchenko, [On hyperbolicity cones associated with elementary symmetric polynomials,]{} Optim.  Lett. [** 2**]{} (2008), 389-402.