Why the Eigenvalues Matter
A quadratic form can look complicated in its original coordinates. For a Gaussian vector, however, an orthogonal change of coordinates turns it into a weighted sum of squares. The signs of those weights—the eigenvalues of the matrix—then reveal whether the form can be negative or positive, give exact bounds when the vector’s length is fixed, and identify the support of its distribution.
Throughout, let \(U\sim N_p(\mu,I_p)\), where \(I_p\) is the \(p\times p\) identity matrix, and let \(A\) be a real symmetric \(p\times p\) matrix. We study \(Q=U^\mathsf{T}AU\). The identity covariance ensures that an orthogonal change of coordinates preserves the standard-normal structure; the mean may be nonzero.
Diagonal Coordinates and Eigenvalue Bounds
By the spectral theorem, there is an orthogonal matrix \(H\) and a diagonal matrix \(\Lambda=\operatorname{diag}(\lambda_1,\ldots,\lambda_p)\) such that \(A=H\Lambda H^\mathsf{T}\). Set \(V=H^\mathsf{T}U\) and \(\nu=H^\mathsf{T}\mu\). Then \(V\sim N_p(\nu,I_p)\), and
Thus the eigenvalues are the weights on the squared coordinates. Since the squared coordinates are nonnegative, nonnegative eigenvalues give a nonnegative form; nonpositive eigenvalues give a nonpositive form. Eigenvalues of both signs permit both signs of the form.
Proof. Write \(x=Hy\), which preserves length because \(H\) is orthogonal. Then \(\lVert x\rVert^2=\sum_i y_i^2\) and \(x^\mathsf{T}Ax=\sum_i\lambda_i y_i^2\). For every \(i\), \(\lambda_{\min}\leq\lambda_i\leq\lambda_{\max}\); multiplying by \(y_i^2\geq 0\) and summing gives the two inequalities. If \(y\) is a unit coordinate vector for an eigenvalue \(\lambda_{\min}\) or \(\lambda_{\max}\), the corresponding bound is an equality. Scaling that eigenvector to length \(r\) proves the sphere statement. \(\square\)
Worked Example: A Positive Semidefinite Form With a Zero Eigenvalue
Let \(A=\begin{pmatrix}2&0\\0&0\end{pmatrix}\), and let \(U\sim N_2((1,-2),I_2)\). The eigenvalues are \(2\) and \(0\), so
For every vector \(u=(u_1,u_2)^\mathsf{T}\), \(u^\mathsf{T}Au=2u_1^2\geq 0\). The eigenvalue bounds give \(0\leq u^\mathsf{T}Au\leq 2\lVert u\rVert^2\); on the circle \(\lVert u\rVert=3\), the exact range is \([0,18]\). The lower endpoint occurs at \(u=(0,3)^\mathsf{T}\), and the upper endpoint at \(u=(3,0)^\mathsf{T}\).
For the random vector \(U\), \(Q\geq 0\) always, and \(Q=0\) only if \(U_1=0\), an event of probability zero. Nevertheless, values of \(U_1\) arbitrarily close to zero have positive probability, so values of \(Q\) arbitrarily close to zero are possible with positive probability.
Worked Example: An Indefinite Form in Non-diagonal Coordinates
Take \(A=\begin{pmatrix}1&2\\2&1\end{pmatrix}\). Its eigenvectors are proportional to \((1,1)^\mathsf{T}\) and \((1,-1)^\mathsf{T}\), with eigenvalues \(3\) and \(-1\). For \(u=(u_1,u_2)^\mathsf{T}\), the form is
At \(u=(1,1)^\mathsf{T}\), the value is \(6>0\); at \(u=(1,-1)^\mathsf{T}\), it is \(-2<0\). Hence \(A\) is indefinite. On the unit circle, the theorem gives the sharp bounds \(-1\leq u^\mathsf{T}Au\leq 3\). On a circle of radius \(2\), they become \(-4\leq u^\mathsf{T}Au\leq 12\).
Sign and Support for a Gaussian Vector
The matrix determines the sign for every possible input vector. For a Gaussian vector with identity covariance, it also determines which signs occur with positive probability. A Gaussian density with positive-definite covariance is positive everywhere, so every nonempty open ball has positive probability. In particular, if the form is negative at even one vector, continuity makes it negative in a neighborhood of that vector, and the Gaussian vector enters that neighborhood with positive probability.
- \(Q=U^\mathsf{T}AU\geq 0\) almost surely if and only if \(A\) is positive semidefinite.
- \(Q\leq 0\) almost surely if and only if \(A\) is negative semidefinite.
- If \(A\) is indefinite, then \(P(Q>0)>0\) and \(P(Q<0)>0\), and the support of \(Q\) is all of \(\mathbb{R}\).
- If \(A\) is nonzero positive semidefinite, the support of \(Q\) is \([0,\infty)\), although \(Q=0\) has probability zero. For nonzero negative semidefinite \(A\), the support is \((-\infty,0]\), and \(Q=0\) has probability zero.
- If \(A=0\), then \(Q=0\) almost surely and its support is \(\{0\}\).
Proof. First, a positive semidefinite \(A\) makes \(u^\mathsf{T}Au\geq 0\) for every \(u\), so \(Q\geq 0\) almost surely. Conversely, if \(A\) is not positive semidefinite, there is a vector \(x\) with \(x^\mathsf{T}Ax<0\). Continuity gives an open ball around \(x\) on which the form remains negative. Since \(U\) has a density positive everywhere, that ball has positive probability. Thus \(Q\geq 0\) almost surely is impossible. Applying the same argument to \(-A\) proves the negative-semidefinite assertion and shows that an indefinite form has both signs with positive probability.
To identify support, let \(g(x)=x^\mathsf{T}Ax\). Continuity and full support of \(U\) imply that the support of \(g(U)\) is the closure of \(g(\mathbb{R}^p)\): if \(t=g(x)\), every neighborhood of \(t\) contains values of \(g\) on an open neighborhood of \(x\), which has positive Gaussian probability. Conversely, if \(t\) lies outside the closure of the image, some neighborhood of \(t\) cannot be reached by \(g(U)\).
If \(A\) is nonzero positive semidefinite, its image under \(g\) is exactly \([0,\infty)\): zero is attained at \(x=0\), and for a unit eigenvector \(v\) with eigenvalue \(\lambda>0\), \(g(cv)=c^2\lambda\) ranges over all nonnegative values as \(c\) varies. The negative-semidefinite case gives \((-\infty,0]\). If \(A\) is indefinite, choose vectors with positive and negative form values. Scaling those vectors produces arbitrarily large values of both signs, and the continuous path between them takes every intermediate value. Its image is therefore \(\mathbb{R}\). Finally, for nonzero positive semidefinite \(A\), \(g(U)=0\) requires \(U\) to lie in the proper subspace \(\ker(A)\), which has probability zero under a density; the negative-semidefinite case follows by considering \(-A\). If \(A=0\), the form is identically zero. \(\square\)
Worked Example: Positive Definite Does Not Mean Support Excludes Zero
Let \(A=\begin{pmatrix}1&0\\0&4\end{pmatrix}\) and \(U\sim N_2((2,-1),I_2)\). Then \(Q=U_1^2+4U_2^2\). Both eigenvalues are positive, so \(Q>0\) unless \(U=(0,0)^\mathsf{T}\), an event of probability zero. Yet the support is \([0,\infty)\): for every \(t\geq 0\), the vector \(u=(\sqrt{t},0)^\mathsf{T}\) gives \(u^\mathsf{T}Au=t\), and the Gaussian density is positive near that vector. In particular, every interval \([0,\varepsilon)\) with \(\varepsilon>0\) has positive probability.
Worked Example: A Negative Semidefinite Form
Let \(A=\begin{pmatrix}-4&0\\0&-1\end{pmatrix}\) and \(U\sim N_2((0,1),I_2)\). Then \(Q=-4U_1^2-U_2^2\), with eigenvalues \(-4\) and \(-1\). For every \(u\),
On the circle \(\lVert u\rVert=3\), the sharp range is \([-36,-9]\). Without a restriction on length, the support of \(Q\) is \((-\infty,0]\). The value zero is in the support because vectors near the origin have positive probability, but it is not attained with positive probability: \(Q=0\) requires \(U=(0,0)^\mathsf{T}\).
Common Pitfall and Practical Use
A common mistake is to say that a positive-definite form has support \((0,\infty)\) because its value is strictly positive almost surely. That reasoning confuses values attained with positive probability and points in the support. A neighborhood of zero still has positive probability, however small it may be; therefore zero belongs to the support. The same distinction applies at the upper or lower endpoints of supports more generally.
The eigenvalue check is useful before attempting a more detailed distributional calculation. It establishes whether negative values are possible, supplies sharp bounds when the vector length is fixed, and reveals whether a proposed support is plausible. It also shows why the mean vector does not change these sign and support conclusions in this setting: a nondegenerate Gaussian vector still has positive density throughout \(\mathbb{R}^p\). The mean changes probabilities, but not the range available to the quadratic map.
Check Your Understanding
Use eigenvalues and the full-support property of the Gaussian vector to answer these questions.
- A symmetric matrix has eigenvalues \(0,2,5\). What can be said about the sign of its quadratic form, and what is the support of the corresponding noncentral Gaussian form?
- A symmetric matrix has eigenvalues \(-3,1\). What are the sharp bounds for its quadratic form on the circle of radius \(2\)?
- For a nonzero positive-definite matrix, can the associated Gaussian quadratic form equal zero with positive probability? Is zero in its support?
- Why does an indefinite matrix imply positive probability of both a positive and a negative quadratic-form value?
- If all eigenvalues are zero, what are the value and support of the quadratic form?