Assuming a distribution on $\mathbb{R}^n$ has a density function $f(x) = \frac{1}{\alpha} \exp(g(x))$
I wonder what is the condition for the distribution to be a multivariate normal distribution?
Some books says the distribution is a multivariate normal distribution, if and only if $g(x)$ is a quadratic form of $x$.
But as far as I know, a quadratic form is defined as a homogeneous polynomial of degree $2$. For a multivariate normal distribution, its density function, when exists (note that a multivariate normal distribution has a density function, if and only if its covariance matrix is nonsingular), may allow $g(x) = \sum_{ij} a_{ij} x_i x_j + \sum_{i} b_ix_i +c$, which is not a homogeneous polynomial in $x_i$'s.
when $g(x) = \sum_{ij} a_{ij} x_i x_j + \sum_{i} b_ix_i +c$, it can always (?) be written as $g(x) = (x-\mu)^T A (x-\mu) + d$. Since the distribution is a multivariate normal distribution and has a density function, doesn't $A$ need to be nonsingular? So the distribution is multivariate normal, if and only if $A$ is nonsingular?
Thanks!