Sign up ×
Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. It's 100% free.

In the classic experimental physics text "Statistical Theory of Signal Detection" by Carl. W. Helstrom, Chapter II, section 4 concerns Gaussian Stochastic Processes. Such a process is observed at times $t_1, t_2, t_3, ... t_n$ to obtain $n$ random variables $x_1, x_2, x_3, ... x_n$, then the probability density for the $n$ variables is of the form:
$$p_n(x_1,t_1; x_2,t_2; ...; x_n,t_n) = M_n\;\exp(-0.5\Sigma_j\Sigma_k\mu_{jk}x_jx_k).$$
where $\mu_{jk}$ is a positive definite matrix and $M_n$ is a normalization constant to give unit probability when you integrate over all possible values:
$$M_n = (2\pi)^{-n/2}\;|\;\textrm{det}\;\mu|^{0.5},$$ where $\mu$ is the determinant of the matrix $\mu_{jk}$. We assume that the expected values are all zero: $E(x_k)=0$.

He notes that the expected value of the product of any odd number of $x_j$ is zero (which seems to follow from symmetry), and gives a formula for the expected value of a product of even numbers of variables. We define $\phi_{jk}$ as:
$$\phi_{jk} = E(x_jx_k).$$ He then notes that:
$$E(x_1x_2x_3x_4) = \phi_{12}\phi_{34} + \phi_{13}\phi_{24} + \phi_{14}\phi_{23}.$$

Is there a simple proof? And is there a simple proof that relates to the methods of quantum mechanics / quantum field theory?

share|cite|improve this question
I'm asking this because Gaussian stochastic processes are entirely classical physics, but the methods used here somehow seem reminiscent of Feynman diagrams to me. – Carl Brannen Aug 11 '11 at 18:36

1 Answer 1

up vote 4 down vote accepted

Yes, it is simple to prove using moment generating functions. And yes, the mathematics is very closely related to that of quantum field theory.

You compute $G(j) = <exp(\sum j_i x_i)>$ where each $j_i$ is a "source" for the corresponding $x_i$. This is easily shown to be something like $G(j) = exp(\sum j_i \mu_{ij}^{-1} j_j)$ To get expectation values you then take $ <x_i x_j ...> = \frac{\partial}{\partial j_i} \frac{\partial}{\partial j_j} ... G(j)|_{j=0}$ . The rest follows simply. In particular, you can see how the variables must be grouped in pairs to get a nonzero results when you set $j=0$ after taking derivatives. You essentially have a Feynman expansion of a non-interacting 0-dimensional field theory.

This is covered nicely in Zee's "Quantum Field Theory in a Nutshell", where it is a simple application of what he calls the "central identity of quantum field theory"

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.