---
abstract: 'A classical inequality of Szász bounds polynomials with no zeros in the upper half plane entirely in terms of their first few coefficients. Borcea-Brändén generalized this result to several variables as a piece of their characterization of linear maps on polynomials preserving stability. In this paper, we improve Szász’s original inequality, use determinantal representations to prove Szász type inequalities in two variables, and then prove that one can use the two variable inequality to prove an inequality for several variables.'
address: |
    Washington University in St. Louis\
    Department of Mathematics, One Brookings Drive, St. Louis, MO 63130
author:
- Greg Knese
title: Global bounds on stable polynomials
---

Introduction
============

We say $p \in {\mathbb{C}}[z]$ is *stable* if $p$ has no zeros in ${\mathbb{C}}_+:=\{z\in {\mathbb{C}}: {\text{Im}}z > 0\}$.

This note is about improvements and generalizations of the following classical inequality of O. Szász.

\[szasz\] If $p(z) = \sum_{j=0}^{d}c_j z^j \in {\mathbb{C}}[z]$ is stable and $p(0)=1$ then $$|p(z)| \leq \exp(|z||c_1| + 3|z|^2(|c_1|^2 + |c_2|)).$$

The purpose of the theorem is to prove $$\mathcal{F}_C = \{p\in {\mathbb{C}}[z]: p \text{ is stable}, p(0)=1, |p'(0)|, |p''(0)| \leq C\}$$ is a normal family whose local uniform limits are entire functions of order at most $2$. One can use this to give a complete characterization of the local uniform limits of stable polynomials. See Theorem 4 Chapter VIII (page 334) of [@Levin].

Only recently have multivariable Szász type inequalities been considered. We say $p \in {\mathbb{C}}[z_1,\dots, z_n]$ is *stable* if $p$ has no zeros in $({\mathbb{C}}_{+})^n$. In their groundbreaking characterization of linear operators on polynomials $T:{\mathbb{C}}[z_1,\dots, z_n] \to {\mathbb{C}}[z_1,\dots, z_n]$ that preserve stability, Borcea-Brändén [@BB] established a Szász type inequality. Its purpose was to prove that the *symbol* of $T$ $$\overline{G_T}(z,w) = \sum_{\alpha \in {\mathbb{N}}^n} (-1)^{\alpha} T(z^\alpha) \frac{w^\alpha}{\alpha!}$$ is actually an entire function. Formally, the symbol is given as $\overline{G_T}(z,w) = T[e^{-z\cdot w}]$.

We let $e_1,\dots, e_n$ be standard basis vectors of ${\mathbb{Z}}^n$.

\[BBthm\] Suppose that $p(z) = \sum_\beta a(\beta) z^{\beta} \in {\mathbb{C}}[z_1,\dots, z_n]$ is stable with $p(0)=1$. Let $$B = 2^{n-1} \frac{\sqrt{2e^2-e}}{e-1} = 2^{n-1}\cdot 2.0210\dots,$$ $$C=6e^2\left(\sum_{i=1}^{n}|a(e_i)|\right)^2 + 
4e^2\sum_{i,j=1}^{n} |a(e_i+e_j)|.$$ Then, $$|p(z)| \leq B\exp(C\|z\|_{\infty}^2).$$

Here $\|z\|_{\infty} = \max\{|z_j|: 1\leq j \leq n\}$.

The original proof of this theorem uses an inequality of Szász (actually inequality below) along with some linear operators that preserve stability to bound all of the coefficients $\{a(\beta)\}$ and then reassemble $p$ to obtain the bound above.

The goal of this paper is to improve both Theorem \[szasz\] and Theorem \[BBthm\] and to give a more geometric proof in the multivariable case.

Our strategy is to first make some minor improvements to Szász’s original argument, then to use determinantal representations of two variable stable polynomials to prove a version of Theorem \[BBthm\] in two variables, and finally to show that an inequality in the two variable case yields an inequality in several variables. We also show how to handle the case where $p(0)=0$.

First, Szász’s inequality can be slightly improved to the following.

\[szimprove\] Suppose $p(z) = \sum_{j=0}^{d}p_j z^j \in {\mathbb{C}}[z]$ is stable and $p(0)=1$. Then, $$|p(z)| \leq \exp( {\text{Re}}(p_1 z) + \frac{1}{2}(|p_1|^2 - 2{\text{Re}}(p_2))|z|^2).$$ The inequality is sharp on the imaginary axis for stable $p \in
{\mathbb{R}}[z]$. Specifically, given $c_1,c_2 \in {\mathbb{R}}$ with $\gamma := \frac{1}{2}(c_1^2-2c_2)> 0$ there exist stable polynomials $p_n \in {\mathbb{R}}[z]$ with $p_n=1+c_1z+c_2z^2
+ \dots$ such that $$\lim_{n\to \infty} |p_n(iy)| = \exp(\gamma y^2).$$

We are unsure about sharpness more generally. Note that necessarily $c_1^2-2c_2 \geq 0$ if $p=1+c_1z+c_2z^2+\dots
\in {\mathbb{R}}[z]$ is stable (e.g. examine $c_1,c_2$ in terms of roots) and $c_1^2-2c_2=0$ if and only if $p\equiv 1$.

The next theorem subsumes this theorem, however we feel it is instructive to go through the elementary one variable argument as a warm-up in Section \[sec:szimprove\].

Using determinantal representations for stable two variable polynomials we are able to offer the following improvement on Theorem \[BBthm\]. Below, $p_j = \frac{\partial p}{\partial z_j}$ and $p_{j,k} = \frac{\partial^2 p}{\partial z_j \partial z_k}$.

\[bisz\] Let $p\in {\mathbb{C}}[z_1,z_2]$ be stable. If $p(0)=1$, then $$\label{mainineq}
|p(z)|\leq\exp({\text{Re}}(\sum_{j=1}^{2} z_j p_j(0)) +
\frac{1}{2}\|z\|_{\infty}^2 (|\sum_{j=1}^{2} p_j(0)|^2-{\text{Re}}(\sum_{j,k=1}^{2} p_{j,k}(0)))).$$

For comparison, using $a \leq (1+a^2)/2$ for the first term in the exponential and writing $p = 1+\sum_{\beta\ne 0} a(\beta) z^{\beta}$ we get $$|p(z)| \leq \sqrt{e} \exp(C\|z\|_{\infty}^2)$$ $$C=(\sum_{j=1}^{2} |a(e_j)|)^2+\sum_{j,k=1}^{2} |{\text{Re}}[a(e_j+e_k)]|.$$ Here we used $p_{jj}(0)= 2 a(2e_j)$ and $p_{jk}(0) = a(e_j+e_k)$ for $j\ne k$.

It turns out that the two variable result can be used to prove an $n$-variable result.

\[msz\] Suppose $p\in {\mathbb{C}}[z_1,\dots,z_n]$ is stable. If $p(0)=1$ then $$|p(z)| \leq \exp(
\sqrt{2}|\nabla p(0)||z|
+(|\nabla p(0)|^2+\|{\text{Re}}Hp(0)\|)|z|^2)$$ where $Hp$ is the Hessian of $p$ and $\|\cdot \|$ denotes operator norm.

We also get an inequality more comparable to Theorem \[BBthm\]: $$|p(z)| \leq \sqrt{e}\exp(C \|z\|_{\infty}^2)$$ $$C = 2(\sum_{j=1}^{n} |a(e_j)|)^2
+2 \sum_{j,k=1}^{n} |{\text{Re}}[a(e_j+e_k)]|).$$

To close the introduction we discuss what happens when $p(0)=0$. In one variable there is only a minor issue because we can factor $p(z) = z^{k}q(z)= p_kz^k + p_{k+1}z^{k+1} +
p_{k+2}z^{k+2}+\dots$ and get a bound depending on $p_k\ne 0 ,p_{k+1},p_{k+2}$: $$|p(z)| \leq |p_k| \exp \left[k (|z|-1) + {\text{Re}}\left(\frac{p_{k+1}}{p_{k}}z\right) + 
\frac{1}{2}|z|^2\left( \left|\frac{p_{k+1}}{p_{k}}\right|^2 -
  2{\text{Re}}\left(\frac{p_{k+2}}{p_k} \right) \right) \right]$$ for stable $p \in {\mathbb{C}}[z]$ with a zero of order $k$ at $0$. We used the inequality $\log |z|^k \leq k(|z|-1)$.

In several variables the case $p(0)=0$ is a little more delicate. Borcea-Brändén covered this case as follows. Given $p(z) =\sum_{\alpha} a(\alpha)z^{\alpha} \in {\mathbb{C}}[z_1,\dots,z_n]$ let $\text{supp}(p) = \{\alpha \in {\mathbb{N}}^n: a(\alpha) \ne 0\}$ and let $\mathcal{M}(p)$ denote the set of minimal elements of $\text{supp}(p)$ with respect to the partial order $\leq$ on ${\mathbb{N}}^n$. Also, for fixed $\mathcal{M} \subset {\mathbb{N}}^n$ let $$\mathcal{M}_2 = \{\alpha+ \beta : \alpha\in \mathcal{M}, \beta \in
{\mathbb{N}}^n, |\beta| \leq 2\}$$

\[BBthm2\] Let $\mathcal{M} \subset {\mathbb{N}}^n$ be a finite nonempty set and $p(z) = \sum_{\alpha} a(\alpha) z^{\alpha} \in {\mathbb{C}}[z_1,\dots, z_n]$ be stable with $\mathcal{M}(p)= \mathcal{M}$. Then, there are constants $B$ and $C$ depending only on the coefficients $a(\alpha)$ with $\alpha \in
\mathcal{M}_2$ such that $$|p(z)| \leq B\exp(C\|z\|_{\infty}^2).$$ Moreover, $B$ and $C$ can be chosen so that they depend continuously on the aforementioned set of coefficients.

With our approach we are able to to get a more explicit estimate in two and several variables. Set $\vec{1} = (1,\dots, 1) \in {\mathbb{C}}^n$.

\[bisz2\] Let $p\in {\mathbb{C}}[z_1,z_2]$ be stable and assume $p$ vanishes to order $r$ at $0$. Write out the homogeneous expansion of $p$: $$p(z) = \sum_{j=r}^{d} P_j(z)$$ where $P_j$ is homogeneous of degree $j$. Then, $$|p(z)| \leq |P_r(\vec{1})|e^{-r/2}\exp\left[{\text{Re}}(\sum_{j=1}^{2} c_j z_j)
+B \|z\|_{\infty}^2 \right]$$ where $$c_j = \frac{1}{P_{r}(\vec{1})} \left[ \frac{\partial P_r}{\partial
      z_j}(\vec{1})\left(1-\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right)
 + \frac{\partial P_{r+1}}{\partial z_j}(\vec{1})\right]$$ $$B = \frac{1}{2}\left(\left|\frac{P_{r+1}(\vec{1})}{P_{r}(\vec{1})}\right|^2
  -2{\text{Re}}\left(\frac{P_{r+2}(\vec{1})}{P_r(\vec{1})}\right) + r\right).$$

This is proved in Section \[sec:bisz2\].

Finally, we present a multivariable Szász inequality for the case $p(0)=0$.

\[msz2\] Suppose $p\in {\mathbb{C}}[z_1,\dots,z_n]$ is stable and vanishes to order $r$ at $0$. If we write out the homogeneous expansion of $p$ $$p(z) = \sum_{j=r}^{d} P_j(z)$$ then $$|p(z)| \leq \|z\|_{\infty}^{r} |P_r(\vec{1})| \exp(C_0+ C_1
\|z\|_{\infty} + C_2 \|z\|_{\infty}^2)$$ where $C_0, C_1,C_2$ are constants depending on $r, P_r(\vec{1}),
P_{r+1}(\vec{1}), P_{r+2}(\vec{1}), \nabla P_{r}(\vec{1}), \nabla
P_{r+1}(\vec{1})$.

The constants $C_0,C_1,C_2$ along with the proof of this theorem are explicitly given in Section \[sec:msz2\].

Improved one variable inequality {#sec:szimprove}
================================

In this section we prove Theorem \[szimprove\].

\[squareslemma\] Suppose ${\alpha}_1,\dots, {\alpha}_d\in {\mathbb{C}}$ with ${\text{Im}}{\alpha}_j \leq 0$. Then, $$\sum_{j=1}^{d} |{\alpha}_j|^2 \leq |\sum_{j=1}^{d} {\alpha}_j|^2 - 2 {\text{Re}}\sum_{j<k} {\alpha}_j {\alpha}_k.$$ Equality holds if and only if either $d=1$ or ${\alpha}_j \in {\mathbb{R}}$ for all $j$.

Note $$|\sum_{j=1}^{d} {\alpha}_j|^2 = \sum_{j=1}^{d} |{\alpha}_j|^2 + 2 {\text{Re}}\sum_{j<k} {\alpha}_j\overline{{\alpha}_k}$$ So, our inequality reduces to showing the following is non-negative: $$2{\text{Re}}\sum_{j<k} ({\alpha}_j \overline{{\alpha}_k} - {\alpha}_j {\alpha}_k) = 2{\text{Re}}\sum_{j<k} {\alpha}_j(-2i){\text{Im}}{\alpha}_k
= 4\sum_{j<k} {\text{Im}}{\alpha}_j {\text{Im}}{\alpha}_k.$$ The last quantity is evidently non-negative and equals zero exactly when either $d=1$ (and the sum is empty) or ${\text{Im}}{\alpha}_j=0$ for all $j$.

Szász uses the inequality $|(1+z)e^{-z}|\leq e^{|z|^2}$ instead of the stronger inequality:

\[loglemma\] For $z \in {\mathbb{C}}$, $z\ne -1$ $$\log|1+z|\leq {\text{Re}}z + \frac{1}{2}|z|^2.$$

Since $\log (1+x) \leq x$ we have $$\begin{aligned}
\log|1+z|&=\frac{1}{2}\log|1+z|^2 \\
&= \frac{1}{2}(\log(1+2{\text{Re}}z + |z|^2)\\
&\leq \frac{1}{2}(2{\text{Re}}z + |z|^2).
\end{aligned}$$

Write $p(z) = \prod_{j=1}^{d}(1+{\alpha}_j z)$ where ${\text{Im}}{\alpha}_j \leq 0$. Note $\sum_{j} {\alpha}_j = p_1$ and $\sum_{j<k} {\alpha}_j {\alpha}_k = p_2$. By Lemma \[loglemma\] and Lemma \[squareslemma\] $$\begin{aligned}
\log |p(z)| & \leq \sum_{j} ({\text{Re}}({\alpha}_j z) + \frac{1}{2}|{\alpha}_j|^2|z|^2) \\
& = {\text{Re}}(p_1 z) + \frac{1}{2}(\sum_{j} |{\alpha}_j|^2) |z|^2 \\
& \leq {\text{Re}}(p_1 z) + \frac{1}{2}(|p_1|^2 - 2{\text{Re}}p_2)|z|^2.
\end{aligned}$$

Regarding sharpness define $\gamma = (c_1^2-2c_2)/2 >0$. Choose $n$ large enough that $d_n = \gamma - c_1^2/(2n)\geq 0.$ Then, the polynomial $$p_n(z) = \left(1+\frac{c_1z}{n}\right)^n\left(1-\frac{d_n z^2}{n}\right)^n$$ is stable, belongs to ${\mathbb{R}}[z]$, and has the correct normalizations. Since $p_n(z) \to \exp{(c_1z-\gamma z^2)}$ locally uniformly, we have $$\lim_{n\to \infty} |p_n(iy)| = \exp{\gamma y^2}$$ which is exactly what was claimed.

It is worth pointing out that Szász proves $$\sum_j |\alpha_j|^2 \leq 2|\sum_j \alpha_j|^2+|\sum_j \alpha_j^2|$$ and then converts $\sum_j \alpha_j^2 = p_1^2-2p_2$ to get the estimate $$\label{BBsz}
\sum_j |\alpha_j|^2 \leq 3|p_1|^2+2|p_2|.$$ By sidestepping the former inequality and estimating $\sum_j |{\alpha}_j|^2$ directly in terms of polynomial coefficients we get a better bound. The inequality is used in [@BB] to prove multivariable Szász inequalities. So, using Lemma \[squareslemma\] in their proof would already improve Theorem \[BBthm\].

Two variable Szász inequality
=============================

Using determinantal formulas it is possible to establish a Szász inequality for two variable polynomials.

\[detdef\] We shall say a stable polynomial $p\in {\mathbb{C}}[z_1,\dots,z_n]$ of total degree $d$ has a *determinantal representation* if there exist $d\times d$ matrices $A,B_1,\dots, B_n$ and a constant $c\in {\mathbb{C}}$ such that

1.  ${\text{Im}}A := \frac{1}{2i}(A-A^*) \geq 0$

2.  for all $j$, $B_j \geq 0$

3.  $\sum_{j=1}^{n} B_j=I$.

4.  $p(z) = c \det(A+\sum_{j=1}^{n} z_j B_j)$.

Theorem \[bisz\] will be broken into two theorems.

\[detrep\] If $p\in {\mathbb{C}}[z_1,z_2]$ is stable, then $p$ has a determinantal representation.

Several different determinantal representations are closely related to this one but not quite equivalent. There are determinantal representations for three variable hyperbolic polynomials, two variable real-zero polynomials, and two variable real-stable polynomials (see [@HV; @VV; @Grinshpan]). It turns out this formula can be derived from a determinantal representation for polynomials with no zeros on the bidisk ${\mathbb{D}}^2 = \{(z_1,z_2):
|z_1|,|z_2|<1\}$ from [@Grinshpan]. We show how to convert from the bidisk formula to Theorem \[detrep\] in Section \[sec:detrep\]. The method of conversion is a very slight modification of what is done in the paper [@GK]. We include the argument for the reader’s convenience; the essence of Section \[sec:detrep\] is not new.

In Section \[sec:ineqfordets\] we prove the following Szász inequality for stable polynomials with determinantal representations.

\[ineqfordets\] Suppose $p\in {\mathbb{C}}[z_1,\dots, z_n]$ has a determinantal representation as above. If $p(0)=1$, then $$|p(z)|\leq\exp({\text{Re}}(\sum_{j=1}^{n} z_j p_j(0)) +
\frac{1}{2}\|z\|_{\infty}^2 (|\sum_{j=1}^{n} p_j(0)|^2-{\text{Re}}(\sum_{j,k=1}^{n} p_{j,k}(0)))).$$

Theorems \[detrep\] and \[ineqfordets\] combine to give Theorem \[bisz\].

Determinantal representations {#sec:detrep}
=============================

In this section we prove Theorem \[detrep\]. We begin by recalling the following.

\[bidiskdetrep\] If $q\in {\mathbb{C}}[z_1,z_2]$ has no zeros in ${\mathbb{D}}^2$ and bidegree $(n,m)$, then there exists a constant $c$ and an $(n+m)\times(n+m)$ contractive matrix $D$ such that $$\label{qdet}
q(z) = c \det(I-D \Delta(z))$$ where $\Delta(z) = z_1P_1+z_2P_2$ and $$P_1 = \begin{pmatrix} I_n & 0 \\ 0 & O_m \end{pmatrix} \qquad 
P_2 = \begin{pmatrix} O_n & 0 \\ 0 & I_m \end{pmatrix}.$$

Let $p \in {\mathbb{C}}[z_1,z_2]$ be stable and have bidegree $(n,m)$. Define $\phi(\zeta) = i\frac{1+\zeta}{1-\zeta}$ and $$q(z_1,z_2) = p(\phi(z_1),\phi(z_2))\left(\frac{1-z_1}{2i}\right)^n \left(\frac{1-z_2}{2i}\right)^m.$$ One can calculate that $\phi^{-1}(\zeta) = \frac{\zeta - i}{\zeta +
  i}$ and $$p(z_1,z_2) = q(\phi^{-1}(z_1),\phi^{-1}(z_2)) (z_1+i)^n(z_2+i)^m.$$ Then, $q$ has no zeros in ${\mathbb{D}}^2$ and so the conclusion of Theorem \[bidiskdetrep\] holds. Then, converting to a formula for $p$ yields $$\begin{aligned}
p(z) &= c\det( (z_1+i)P_1 + (z_2+i)P_2 - D((z_1-i)P_1+(z_2-i)P_2))\\
&= c \det( (I-D)\Delta(z) +i (I+D))
\end{aligned}$$ Since $D$ is a contraction, the eigenspace corresponding to eigenvalue $1$ is reducing (if nontrivial). Thus, there exists a unitary $U$ such that $$D = U\begin{pmatrix} I & 0 \\ 0 & K \end{pmatrix} U^*$$ where $K$ is a contractive $k\times k$ matrix for which $1$ is not an eigenvalue. Here $k$ is the codimension of the eigenspace of $D$ corresponding to eigenvalue $1$.

Then, $$\begin{aligned}
p(z) &= c \det\left( \begin{pmatrix} 0 & 0 \\ 0 & I-K \end{pmatrix} U^* \Delta(z) U + 
i\begin{pmatrix} 2I & 0 \\ 0 & I+K \end{pmatrix}\right)\\
&= c \det(I-K)
\det\left( \begin{pmatrix} 0 & 0 \\ 0 & I \end{pmatrix} U^*\Delta(z) U + \begin{pmatrix} 2iI & 0 \\ 0 & A \end{pmatrix}\right)
\end{aligned}$$ where $A = i(I+K)(I-K)^{-1}$.

Let $B_j$ equal the bottom right $k\times k$ block of $U^*P_jU$. Then, $$p(z) = c\det(I-K) 
\det\begin{pmatrix} 2iI & 0 \\ * & A + \sum_j z_jB_j \end{pmatrix}
= c_0 \det(A+\sum_j z_j B_j)$$ where $c_0$ is a new constant (the $*$ denotes a block we are unconcerned with). Since $P_1+P_2 = I$, $B_1+B_2= I$. Also note that $p(t,t) = c_0 \det(A+tI)$ has degree $k$ so that $k \leq \deg p$. On the other hand, the determinantal formula for $p$ has total degree at most $k$, so that $\deg p \leq k$. Therefore the matrices in our formula have size matching the total degree of $p$. Finally, $${\text{Im}}A  = (I-K)^{-1}(I-KK^*)(I-K^*)^{-1} \geq 0.$$ This proves Theorem \[detrep\].

Szász inequality for determinantal polynomials {#sec:ineqfordets}
==============================================

In this section we prove Theorem \[ineqfordets\].

Suppose $p(z) = c \det(A+\sum_{j=1}^{n} z_j B_j)$ where $\sum_{j=1}^{n} B_j = I$, $B_j \geq 0$, ${\text{Im}}A \geq 0$, and $p(0)=1$. By the last normalization $A$ is invertible with $c\det A = 1$ so that $$p(z) = \det( I + \sum_{j=1}^{n} z_j X_j)$$ where $X_j = B_j A^{-1}$. As with complex numbers, ${\text{Im}}(A^{-1}) \leq 0$.

It helps to make note of a few formulas for the derivatives of $p$. Recall that if $A(t)$ is a differentiable matrix function then $$\frac{d}{dt} \det A(t) = {\text{tr}}( A'(t) A(t)^{-1})) \det A(t)$$ whenever $A(t)$ is invertible. Here ${\text{tr}}$ is the trace of a matrix.

Letting $X(z) = I+\sum_{j=1}^n z_j X_j$, whenever $p(z) \ne 0$ we have $$\begin{aligned}
p_j(z) &= \text{tr}(X_j(X(z))^{-1}) p(z)\\
p_{jk}(z) &= -\text{tr}(X_j (X(z))^{-1}X_k(X(z))^{-1})p(z) +
\text{tr}(X_j(X(z))^{-1})\text{tr}(X_k(X(z))^{-1})p(z)
\end{aligned}$$ so that $$p_j(0) = \text{tr} X_j \quad p_{jk}(0)=-\text{tr}(X_jX_k) + \text{tr}(X_j)\text{tr}(X_k).$$

For a positive definite matrix $P$ we have $\log P \leq P-I$ simply because the same inequality holds for the eigenvalues of $P$. Therefore,

$$\begin{aligned}
\log|p(z)| &= \frac{1}{2} \log \det(X(z)^*X(z)) \label{logp1} \\
&= (1/2) {\text{tr}}\log X(z)^*X(z) \label{logp2} \\
& \leq (1/2) {\text{tr}}(X(z)^*X(z) - I) \label{logp3} \\
& = (1/2) {\text{tr}}(2 {\text{Re}}(\sum_{j=1}^{n} z_j X_j)+ (\sum_{j=1}^{n} z_j X_j)^*(\sum_{k=1}^{n} z_k X_k)). \nonumber
$$

Now, $$\begin{aligned}
{\text{tr}}(\sum_{j}z_j X_j)^*(\sum_{k} z_k X_k) &= {\text{tr}}[(A^*)^{-1}
(\sum_{j}z_j B_j)^*(\sum_{k} z_k B_k) A^{-1}]\\
& = {\text{tr}}[ (\sum_{j}z_j
B_j)^*(\sum_{k} z_k B_k) A^{-1}(A^*)^{-1}].
\end{aligned}$$

By Lemma \[PM\], Lemma \[sumB\], and Lemma \[Im\] below we have $$\begin{aligned}
{\text{tr}}(\sum_{j}z_j X_j)^*(\sum_{k} z_k X_k) &\leq \|\sum_j z_j B_j\|^{2}
{\text{tr}}\left[(A^*)^{-1} A^{-1}\right] \\
&\leq \|z\|_{\infty}^2\left[ |{\text{tr}}(A^{-1})|^2 - {\text{Re}}(
({\text{tr}}A^{-1})^2 - {\text{tr}}A^{-2})\right].
\end{aligned}$$

Finally, since $\sum_j B_j = I$ we have $$\sum_j p_j(0) = {\text{tr}}A^{-1}$$ and $$\sum_{j,k} p_{jk}(0) = ({\text{tr}}A^{-1})^2 - {\text{tr}}A^{-2}.$$ Thus, $$\log|p(z)| \leq {\text{Re}}(\sum_{j=1}^{n} p_j(0)z_j) + \frac{1}{2} \|z\|_{\infty}^2 (|\sum_{j=1}^{n} p_j(0)|^2 - {\text{Re}}(\sum_{j,k=1}^{n} p_{jk}(0)))$$ which proves Theorem \[ineqfordets\] modulo the following three lemmas.

\[PM\] Let $P, M$ be $n\times n$ matrices. If $P \geq 0$, then $$|\text{tr}(MP)| \leq \|M\| \text{tr}P$$

Since $P\geq 0$, we can decompose $P = \sum_j v_jv_j^*$ where $v_j \in {\mathbb{C}}^n$. Then, $$|\text{tr}{MP}| \leq \sum_j |\text{tr}{Mv_jv_j^*}|
= \sum_j |\langle Mv_j, v_j \rangle|
\leq \|M\| \sum_j \|v_j\|^2 = \|M\| \text{tr}{P}.$$

The following is a standard result.

\[sumB\] Suppose $B_1,\dots, B_n$ are $N\times N$ matrices. Assume for all $j$, $B_j \geq 0$ and $\sum_j B_j = I$. Then, there exist pairwise orthogonal projection matrices $P_1,\dots, P_n$ of size $m\times m$ where $m= nN$ such that $$\label{bj}
 B_j = (I_N,0,\dots, 0) P_j (I_N,0,\dots, 0)^t.$$ In particular, for $z=(z_1,\dots, z_n) \in {\mathbb{C}}^n$ $$\|\sum_{j} z_j B_j\| \leq \|z\|_{\infty}.$$

We can factor $B_j = A_j^* A_j$ with $N\times N$ matrix $A_j$. The $nN\times N$ matrix $$T = \begin{pmatrix} A_1 \\ \vdots \\ A_n\end{pmatrix}$$ is an isometry from ${\mathbb{C}}^N$ to ${\mathbb{C}}^{nN}$ since $T^*T = \sum_j B_j = I$. We can extend $T$ to a $m\times m$ unitary $U$. Let $Q_j$ be the orthogonal projection onto the $j$-th block of ${\mathbb{C}}^{m} = {\mathbb{C}}^N\oplus \cdots \oplus {\mathbb{C}}^N$. Set $P_j = U^*Q_j U$. Then, holds and $$\|\sum_{j} z_j B_j\| \leq \|\sum_j z_j P_j \| \leq \|z\|_{\infty}.$$

The following lemma is an adaptation of our one variable argument.

\[Im\] If $M$ is a square matrix with ${\text{Im}}M \geq 0$, then $$\text{tr}M^*M \leq |\text{tr}M|^2 - {\text{Re}}( (\text{tr} M)^2 - \text{tr} M^2).$$

Write $A = {\text{Re}}M, B = {\text{Im}}M$. Then, $$\begin{aligned}
|{\text{tr}}M|^2 &= ({\text{tr}}A)^2+({\text{tr}}B)^2,\\
{\text{tr}}M^*M &= {\text{tr}}(A^2+B^2 + i(AB-BA)) = {\text{tr}}A^2 + {\text{tr}}B^2,\\
{\text{Re}}({\text{tr}}M)^2 &= ({\text{tr}}A)^2 - ({\text{tr}}B)^2,\\
{\text{Re}}[{\text{tr}}M^2] &= {\text{Re}}[{\text{tr}}( A^2-B^2 + i(AB+BA))] = {\text{tr}}A^2-{\text{tr}}B^2. 
\end{aligned}$$ Then, $$|\text{tr}M|^2-{\text{tr}}M^*M - {\text{Re}}( (\text{tr} M)^2 - \text{tr} M^2) = 2(({\text{tr}}B)^2 - {\text{tr}}B^2)$$ If $B$ has eigenvalues $\beta_j\geq 0$ then $$({\text{tr}}B)^2 - {\text{tr}}B^2 = (\sum \beta_j)^2 - \sum \beta_j^2 
= \sum_{j\ne k} \beta_j\beta_k \geq 0.$$ This proves the claimed inequality.

Szász inequality for determinants with $p(0)=0$ {#sec:bisz2}
===============================================

As with Theorem \[bisz\] we will prove a Szász inequality for polynomials with determinantal representations and Theorem \[bisz2\] will follow via Theorem \[detrep\].

\[szdetreps2\] Suppose $p\in {\mathbb{C}}[z_1,\dots,z_n]$ has a determinantal representation as in Definition \[detdef\]. Assume $p$ vanishes to order $r$ at $0$. Write out the homogeneous expansion of $p$: $$p(z) = \sum_{j=r}^{d} P_j(z)$$ where $P_j$ is homogeneous of degree $j$. Then, $$|p(z)| \leq |P_r(\vec{1})|e^{-r/2}\exp\left[{\text{Re}}(\sum_{j=1}^{n} c_j z_j)
+B \|z\|_{\infty}^2 \right]$$ where $$c_j = \frac{1}{P_{r}(\vec{1})} \left[ \frac{\partial P_r}{\partial
      z_j}(\vec{1})(1-\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})})
 + \frac{\partial P_{r+1}}{\partial z_j}(\vec{1})\right]$$ $$B = \frac{1}{2}\left(\left|\frac{P_{r+1}(\vec{1})}{P_{r}(\vec{1})}\right|^2
  -2{\text{Re}}\left(\frac{P_{r+2}(\vec{1})}{P_r(\vec{1})}\right) + r\right)$$ and $\vec{1} = (1,\dots,1) \in {\mathbb{C}}^n$.

Write $p(z) = c\det(A+\sum_{j=1}^{d} z_j B_j)$ as in Definition \[detdef\]. Since $p(0)=0$, $\det A = 0$. Since ${\text{Im}}A \geq 0$, the eigenspace corresponding to eigenvalue $0$ is reducing for $A$; see Lemma \[ImA\] below. Let $s$ equal the dimension of the kernel of $A$.

So, after conjugating by a unitary we can rewrite $p$ in the form $$p(z) = c\det\left( \begin{pmatrix} 0 & 0 \\ 0 & C \end{pmatrix} +
\sum_{j=1}^{d} z_j B_j \right)$$ where $C$ is an invertible $(d-s)\times(d-s)$ matrix with ${\text{Im}}C\geq 0$ and the $B_j$ are relabelled after conjugating (they satisfy all of the same properties as before). Define $X_j = B_j \begin{pmatrix} I & 0
  \\ 0 & C^{-1} \end{pmatrix}$ and $J = \begin{pmatrix} O_s & 0 \\ 0 &
  I_{d-s} \end{pmatrix}$. Then, $$p(z) = c_0 \det (J +\sum_{j=1}^{d} z_j X_j)$$ where $c_0 = c \det C$. Let $X(z) = J +\sum_{j=1}^{d} z_j X_j$. Let $X_s(z)$ be the top left $s\times s$ block of $X(z)$. Evaluating $\det
X$ starting with the top left $s\times s$ block gives $$\det X(z) = \det X_s(z) + \text{ higher order terms}.$$ Note that $\det X_s(z)$ is homogeneous of degree $s$ and $X_s(\vec{1}) = I_s$ since $\sum_j B_j=I$. This proves $s=r$.

We can follow some of the argument in Section \[sec:ineqfordets\]. Equations , , hold when $p(z)\ne 0$ but rearranges into $$\begin{aligned}
\log|p(z)/c_0| & \leq (1/2) \text{tr}\left( \begin{pmatrix} -I_r & 0 \\ 0 &
  0 \end{pmatrix} + 2{\text{Re}}(\sum_{j=1}^{n} z_j X_j) + (\sum_j z_j
X_j)^*(\sum_k z_k X_k)\right)\\
& \leq -r/2 + {\text{Re}}(\sum_{j=1}^{n} z_j {\text{tr}}(X_j)) + (1/2){\text{tr}}(\sum_{j}z_j
X_j)^*(\sum_k z_k X_k).
\end{aligned}$$

As before using Lemmas \[PM\], \[sumB\], \[Im\] we have $$\begin{aligned}
{\text{tr}}(\sum_{j}z_j
X_j)^*(\sum_k z_k X_k) &\leq  \|\sum_{j} z_j B_j\|^2 {\text{tr}}( \begin{pmatrix} I_r & 0 \\ 0 &
  (C^*)^{-1}C^{-1} \end{pmatrix}) \\
&\leq \|z\|_{\infty}^2( r + |{\text{tr}}(C^{-1})|^2 - {\text{Re}}(
({\text{tr}}C^{-1})^2 - {\text{tr}}C^{-2})).
\end{aligned}$$

Now we must relate these quantities to intrinsic quantities of $p$.

First, $p(t\vec{1}) = c_0 \det \begin{pmatrix} tI_r & 0 \\ 0 &
  I+tC^{-1} \end{pmatrix} = c_0 t^r \det(I+tC^{-1}).$ So, using this formula and the homogeneous expansion of $p$ we get $$\begin{aligned}
 t^{-r}p(t\vec{1}) \big|_{t=0} &=  c_0  = P_r(\vec{1}) \\ 
 \frac{d}{dt}\big|_{t=0} t^{-r}p(t\vec{1})  &= c_0 {\text{tr}}C^{-1} = P_{r+1}(\vec{1})\\
 \frac{d^2}{dt^2}\big|_{t=0} t^{-r} p(t\vec{1}) & = c_0( ({\text{tr}}C^{-1})^2 -
 {\text{tr}}C^{-2}) = 2P_{r+2}(\vec{1}).
\end{aligned}$$ It is more difficult to calculate ${\text{tr}}X_j$. Define $$q(s,t) = p(s e_j + t\vec{1}) = c_0 \det\left(J + s X_j + t \begin{pmatrix} I & 0 \\
  0 & C^{-1} \end{pmatrix}\right).$$ Note $\frac{\partial q}{\partial s}(0,t) = \frac{\partial p}{\partial z_j}(t\vec{1})$. Then, $$\begin{aligned}
\frac{\partial q}{\partial s}(0,t) &= {\text{tr}}\left( X_j \begin{pmatrix}
t^{-1}I_r & 0 \\ 0 & (I+tC^{-1})^{-1} \end{pmatrix}\right)  p(t\vec{1}) \\
&= {\text{tr}}\left(X_j \begin{pmatrix}
I_r & 0 \\ 0 & t(I+tC^{-1})^{-1} \end{pmatrix}\right) c_0 t^{r-1}
\det(I+tC^{-1}) 
\end{aligned}$$ Thus, we can do the following computation with matrices and also with the homogeneous expansion of $p$ $$t^{-r+1}\frac{\partial q}{\partial s}(0,t)\Big|_{t=0} = c_0 {\text{tr}}\left(X_j \begin{pmatrix} I_r & 0 \\ 0 & 0 \end{pmatrix}\right) = \frac{\partial
  P_r}{\partial z_j}(\vec{1})$$ Therefore, $$\frac{\partial q}{\partial s}(0,t) -
t^{-1}\frac{p(t\vec{1})}{P_r(\vec{1})} \frac{\partial P_r}{\partial
  z_j}(\vec{1}) =
{\text{tr}}\left( X_j \begin{pmatrix} 0 & 0 \\ 0 &
  t(I+tC^{-1})^{-1} \end{pmatrix}\right)c_0 t^{r-1} \det(I+tC^{-1})$$ which implies $$t^{-r}\left(\frac{\partial q}{\partial s}(0,t) -
t^{-1}\frac{p(t\vec{1})}{P_r(\vec{1})} \frac{\partial P_r}{\partial
  z_j}(\vec{1})\right)\Big|_{t=0} = {\text{tr}}\left(X_j \begin{pmatrix} 0 & 0 \\ 0 &
  I\end{pmatrix}\right) c_0 =\frac{\partial P_{r+1}}{\partial z_j}(\vec{1})
- \frac{P_{r+1}(\vec{1})}{P_r(\vec{1})} \frac{\partial P_r}{\partial z_j}(\vec{1})$$ Therefore, $$c_0 {\text{tr}}X_j = \left(1-\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right)
\frac{\partial P_r}{\partial z_j}(\vec{1}) + \frac{\partial
  P_{r+1}}{\partial z_j}(\vec{1}).$$ If we reassemble we get $$\log|p(z)/P_r(\vec{1})| \leq -r/2 + {\text{Re}}(\sum_{j=1}^{n} c_j z_j ) +
B \|z\|_\infty^2$$ where $$c_j = \frac{1}{P_r(\vec{1})} \left[\left(1-\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right)
\frac{\partial P_r}{\partial z_j}(\vec{1}) + \frac{\partial
  P_{r+1}}{\partial z_j}(\vec{1})\right]$$ $$B = \frac{1}{2}\left( r+ \left|
    \frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right|^2
  -2{\text{Re}}\left(\frac{P_{r+2}(\vec{1})}{P_r(\vec{1})}\right)\right)$$ and this concludes the proof.

\[ImA\] Suppose $A$ is a matrix with ${\text{Im}}A \geq 0$. If $0$ is an eigenvalue of $A$ with eigenspace of dimension $s$, then there exists a unitary $U$ such that $$U^* A U = \begin{pmatrix} O_{s} & 0 \\ 0 & C\end{pmatrix}$$ where ${\text{Im}}C \geq 0$ and $C$ is invertible.

If we write $A$ using an orthonormal basis for its kernel followed by an orthonormal basis for the orthogonal complement of its kernel, we can put $A$ into the form $$\begin{pmatrix} O_s & B \\ 0 & C\end{pmatrix}$$ by conjugating by a unitary. This matrix will still have positive semi-definite imaginary part: $$\begin{pmatrix} O_s & \frac{1}{2i} B \\ -\frac{1}{2i} B^* & {\text{Im}}C\end{pmatrix} \geq 0.$$ This implies $B=0$. Note $C$ is invertible because it cannot have $0$ as an eigenvalue.

Multivariable Szász inequalities
================================

Using the two variable Szász inequality we can establish the multivariable inequality Theorem \[msz\].

We will frequently use the component-wise partial order on ${\mathbb{R}}^n$: $x \geq y$ if and only if for all $j=1,\dots, n$, $x_j\geq y_j$.

For $z=x+iy \in {\mathbb{C}}_+^{n}$ we have $$|p(x_1+iy_1,\dots, x_n+iy_n)| \geq |p(x_1\pm i y_1, \dots, x_n\pm iy_n)|$$ for all independent choices of $\pm$ by Lemma \[Polya2\] below (more precisely, we can hold fixed any variables with a “$+$” and apply Lemma \[Polya2\] to the remaining variables). So, it is enough to prove Theorem \[msz\] for $z \in {\mathbb{C}}_+^n$.

By Lemma \[Polya3\] below, if $0\leq y\leq \tilde{y}$ then $$|p(x+iy)| \leq |p(x+i\tilde{y})|.$$ Define $\tilde{y}$ as the vector with $j$-th component $$\tilde{y}_j = \max(|x_j|, y_j).$$ Then, $\tilde{y} \geq \pm x$ and $\tilde{y}\geq y$.

Define $$q(w_1,w_2) = p(w_1(\tilde{y}+x)+w_2(\tilde{y}-x))$$ which has no zeros in ${\mathbb{C}}_+^2$ and $q(0)=1$. We will now apply Theorem \[bisz\] using all of the following computations. $$q\left(\frac{i+1}{2}, \frac{i-1}{2}\right) = p(x+i\tilde{y})$$ $$q_1(w) = \sum_j p_j(w_1(\tilde{y}+x)+w_2(\tilde{y}-x))(\tilde{y}_j+x_j)$$ $$q_2(w) = \sum_j p_j(w_1(\tilde{y}+x)+w_2(\tilde{y}-x))(\tilde{y}_j-x_j)$$ $$q_{11}(0) = \sum_{j,k} p_{jk}(0)(\tilde{y}_j+x_j)(\tilde{y}_k+x_k)$$ $$q_{12}(0) = \sum_{j,k} p_{jk}(0)(\tilde{y}_j+x_j)(\tilde{y}_k-x_k)$$ $$q_{22}(0) = \sum_{j,k} p_{jk}(0)(\tilde{y}_j-x_j)(\tilde{y}_k-x_k)$$ $$q_1(0)(i+1)/2 + q_2(0)(i-1)/2 = \nabla p(0)\cdot (x+i\tilde{y})$$ $$q_1(0)+q_2(0) = 2\nabla p(0)\cdot \tilde{y}$$ $$q_{11}(0)+2q_{12}(0)+q_{22}(0) = 4\sum_{j,k}p_{jk}(0)\tilde{y}_j\tilde{y}_k.$$

Thus, by Theorem \[bisz\]

$$\begin{gathered}
\label{worth}
\log\left|q\left(\frac{i+1}{2},\frac{i-1}{2}\right)\right|  \leq
{\text{Re}}(\nabla p(0)\cdot (x+i\tilde{y})) \\+\frac{1}{4}(|2\nabla p(0)\cdot \tilde{y}|^2 
-4Re(\sum_{j,k} p_{jk}(0)\tilde{y}_j\tilde{y}_k)) \\\leq 
\sqrt{2}|\nabla p(0)||z|
+(|\nabla p(0)|^2+\|{\text{Re}}Hp(0)\|)|z|^2 $$

where we have used $|x+i\tilde{y}| \leq \sqrt{2}|z|$ and $|\tilde{y}|\leq |z|$.

Since Theorem \[BBthm\] is an estimate on polydisks it is worth pointing out that yields $$\begin{aligned}
\log|p(z)| &\leq \|z\|_{\infty} \sqrt{2} \sum_j |p_j(0)| + \|z\|_{\infty}^2((\sum_j |p_j(0)|)^2
+ \sum |{\text{Re}}[p_{jk}(0)]|) \\
& \leq \|z\|_{\infty} \sqrt{2} \sum_j |a(e_j)| + \|z\|_{\infty}^2((\sum_j |a(e_j)|)^2
+2 \sum |{\text{Re}}[a(e_j+e_k)]|) \\
&\leq \frac{1}{2} + \|z\|_{\infty}^2(2(\sum_j |a(e_j)|)^2
+2 \sum |{\text{Re}}[a(e_j+e_k)]|)
\end{aligned}$$ where $p = \sum a(\beta) z^{\beta}$ and in the last line we used the inequality $a\leq (1+a^2)/2$. This gives $$|p(z)| \leq \sqrt{e}\cdot \exp(C\|z\|_{\infty}^2)$$ $$C = 2(\sum_j |a(e_j)|)^2
+2 \sum |{\text{Re}}[a(e_j+e_k)]|.$$

The following is a standard result. See Lemma 2.8 of [@BB] for instance.

\[Polya2\] If $p\in {\mathbb{C}}[z_1,\dots, z_n]$ has no zeros in ${\mathbb{C}}_+^n$ then for $z=x+iy\in {\mathbb{C}}_{+}^n$ $$|p(x+iy)| \geq |p(x-iy)|.$$

The one variable polynomial $q(\zeta) = p(x+\zeta y)$ has no zeros in ${\mathbb{C}}_+$. Then, $q$ can be factored as a product of terms of the form $(1+\alpha \zeta)$ where ${\text{Im}}\alpha \leq 0$. We can then check directly that $$|1+i \alpha| \geq |1-i\alpha|$$ which implies $|q(i)|\geq |q(-i)|$.

\[Polya3\] If $p\in {\mathbb{C}}[z_1,\dots, z_n]$ has no zeros in ${\mathbb{C}}_{+}^n$ and if $0\leq y \leq \tilde{y}$ then for any $x\in {\mathbb{R}}^n$ $$|p(x+iy)| \leq |p(x+i\tilde{y})|.$$

The one variable polynomial $q(\zeta) = p(x+iy+\zeta(\tilde{y}-y))$ has no zeros in ${\mathbb{C}}_{+}$. Factors of $q$ are of the form $(1+\alpha \zeta)$ with ${\text{Im}}\alpha \leq 0$. Since $|1+i\alpha| \geq 1$ we have $|q(0)| \leq |q(i)|$.

We can get a slightly better bound on ${\mathbb{R}}^n$ by modifying the argument of Theorem \[msz\].

Suppose $p\in {\mathbb{C}}[z_1,\dots, z_n]$ is stable. If $p(0)=1$ then for $x\in {\mathbb{R}}^n$ $$\log |p(x)| \leq {\text{Re}}(\nabla p(0)\cdot x) + \frac{1}{2}(|\nabla p(0)|^2+ \|{\text{Re}}(Hp)(0)\|) |x|^2$$ where $Hp$ is the Hessian matrix of $p$.

We can write $x \in {\mathbb{R}}^n$ as $x= x_{+} - x_{-}$ where $(x_+)_j = \begin{cases} x_j & \text{ if } x_j \geq 0 \\ 0 & \text{ if
  } x_j<0 \end{cases}$. Define $$P(z_1,z_2) = p(z_1 x_+ + z_2 x_{-})$$ which has no zeros in ${\mathbb{C}}_{+}^2$ and $P(0)=1$. Set $S_+ = \{j: x_j \geq 0\}$, $S_{-}=\{j:x_j<0\}$.

Note that $$P_1(z) = \sum_{j\in S_{+}} p_j(z_1x_+ + z_2 x_{-}) |x_j| \qquad P_2(z) = \sum_{j\in S_{-}} p_j(z_1x_+ + z_2 x_{-}) |x_j|$$ $$P_{11}(0)= \sum_{j,k\in S_{+}} p_{jk}(0)|x_j||x_k|,  \quad P_{12}(0) =
\sum_{j\in S_{+},k\in S_{-}} p_{jk}(0)|x_j||x_k|, \quad 
P_{22}(0) = \sum_{j,k\in S_{-}} p_{jk}(0) |x_j||x_k|.$$ Now, since $P(1,-1) = p(x)$ we have $$\log|p(x)| \leq {\text{Re}}( \nabla p(0) \cdot x) 
+ (1/2)(|\nabla p(0)|^2|x|^2 - {\text{Re}}( \sum_{jk} p_{jk}(0) |x_j||x_k|))$$ by Theorem \[bisz\].

Multivariable inequalities when $p(0)=0$ {#sec:msz2}
========================================

In this section we prove Theorem \[msz2\]. Write the homogeneous expansion of $p$ $$p(z) = \sum_{j=r}^{d} P_j(z).$$ Notice that $P_r(z)$ is stable itself by Hurwitz’s theorem because $$P_r(z) = \lim_{t\searrow 0} t^{-r} p(t z)$$ exhibits $P_r$ as a limit of polynomials with no zeros in ${\mathbb{C}}_+^{n}$.

We can make some of the reductions as in the previous section. We may assume $z= x+iy \in {\mathbb{C}}_{+}^n$ by Lemma \[Polya2\]. Define $m =\max\{|x_j|,y_j: 1\leq j
\leq n\}$ and $\tilde{y} = m\vec{1}$. Then, $\tilde{y}\geq \pm x$, $\tilde{y}\geq y$ and $|p(z)| \leq |p(x+i\tilde{y})|$. Define $$q(w_1,w_2) = p(w_1 (\tilde{y} +x) + w_2(\tilde{y}-x))$$ which is stable and has homogeneous expansion $$\sum_{j=r}^{d} Q_j(w) = \sum_{j=r}^{d} P_j(w_1 (\tilde{y} +x) + w_2(\tilde{y}-x)).$$ All of the terms above are homogeneous of the correct degree but it is conceivable that the first term vanishes. Setting $w_1=w_2=1$ we see the first term evaluates to $P_r(2\tilde{y}) = (2m)^r
P_r(\vec{1})$ which is non-zero.

The data we need for Theorem \[bisz2\] is: $$\begin{aligned}
Q_{j}(\vec{1}) &= (2m)^{j}P_{j}(\vec{1}) \\
\frac{\partial Q_j}{\partial w_1}(\vec{1}) &= (2m)^{j-1}\sum_{k=1}^{n}
\frac{\partial P_j}{\partial z_k}(\vec{1}) (m+x_k) \\
\frac{\partial Q_j}{\partial w_2}(\vec{1}) &= (2m)^{j-1}\sum_{k=1}^{n}
\frac{\partial P_j}{\partial z_k}(\vec{1}) (m-x_k) 
\end{aligned}$$ and so (omitting some details) $$\left|q\left(\frac{i+1}{2},\frac{i-1}{2}\right)\right| \leq (2m)^{r}|P_r(\vec{1})| e^{-r/2}
\exp({\text{Re}}( A ) + \frac{1}{2} B )$$ where $$A = \frac{1}{P_r(\vec{1})} \left( \nabla P_r(\vec{1})\cdot
  (x+i\tilde{y})\left(\frac{1}{2m}- \frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right)
  + \nabla P_{r+1}(\vec{1})\cdot (x+i\tilde{y})\right)$$ $$B = \frac{1}{2}\left( (2m)^2
  \left|\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right|^2 - 2(2m)^2
  {\text{Re}}\left(\frac{P_{r+2}(\vec{1})}{P_r(\vec{1})}\right) +r\right).$$ Note $m\leq \|z\|_{\infty}$ and $\|x+i\tilde{y}\|_{\infty} \leq
\sqrt{2} m$. We can crudely estimate $A$: $$\begin{aligned}
|A| &\leq \frac{1}{|P_r(\vec{1})|} \left( \|\nabla
  P_r(\vec{1})\|_{1}(\sqrt{2}m)\left(\frac{1}{2m} +
  \left|\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right|\right)
+ \|\nabla P_{r+1}(\vec{1})\|_{1} \sqrt{2}m\right) \\
&\leq \frac{\|\nabla P_r(\vec{1})\|_1}{\sqrt{2} |P_r(\vec{1})|}
+ \frac{\sqrt{2}}{|P_r(\vec{1})|}\left( \|\nabla P_r(\vec{1})\|_1
    \left|\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right| +
\|\nabla P_{r+1}(\vec{1})\|_{1}\right) \|z\|_{\infty}
\end{aligned}$$ and $$|B| \leq 2\|z\|^2_{\infty}\left(
\left|\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right|^2 - 2
  {\text{Re}}\left(\frac{P_{r+2}(\vec{1})}{P_r(\vec{1})}\right)\right) + r/2.$$ Here we use $\|\cdot \|_1$ for $\ell^1$ norm of a vector. Putting everything together $$|p(z)| \leq \|z\|_{\infty}^{r} |P_r(\vec{1})|
\exp(C_0+C_1\|z\|_{\infty} + C_2\|z\|_{\infty}^{2})$$ where $$C_0 = r(\log(2)-1/4) + \frac{\|\nabla P_r(\vec{1})\|_1}{\sqrt{2} |P_r(\vec{1})|}$$ $$C_1 = \frac{\sqrt{2}}{|P_r(\vec{1})|}\left( \|\nabla P_r(\vec{1})\|_1
    \left|\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right| +
\|\nabla P_{r+1}(\vec{1})\|_{1}\right)$$ $$C_2 = \left(
\left|\frac{P_{r+1}(\vec{1})}{P_r(\vec{1})}\right|^2 - 2
  {\text{Re}}\left(\frac{P_{r+2}(\vec{1})}{P_r(\vec{1})}\right)\right).$$
