Siegert relation ($g_{12}^{(2)}(\tau )=1+|g_{12}^{(1)}(\tau ){|}^2$)

For one selected polarization of waves, write the detected complex fields as $E_1(t),E_2(t+\tau )$ and their light intensities as $I_j=|E_j{|}^2$, absorbing a common conversion constant. With stationary ensemble averages and nonzero ${\bar{I}}_j=⟨I_j⟩$, defineHere second-order optical coherence concerns distinct detector channels, with the corresponding normal ordering in the quantum formulation. Independent detection has $g^{(2)}=1$. Subtracting the means gives the normalized covarianceThus multiplying raw photocurrents contains a large mean-product contribution; multiplying their fluctuations isolates the small excess. Independent detector shot noise does not supply a mean cross-channel excess, although it increases measurement uncertainty.

The displayed Siegert relation requires a zero-mean, jointly circular complex Gaussian random field, the usual chaotic thermal light model. Circular means invariance under a common arbitrary complex phase rotation; in particular $⟨E_i{E}_j⟩=0$. Apply the Gaussian moment pairing theorem to the real and imaginary parts:The remaining pairing vanishes by circularity. Division proves the Siegert relation. Gaussianity is an assumption about fourth moments, not a consequence of merely measuring first-order optical coherence; a coherent state or photon antibunching cannot be interpreted by this thermal formula.

For equal-intensity beams combined in a common optical mode with matched wave polarizations, interferometric visibility is $V=|g_{12}^{(1)}|$. Consequently ideal thermal stellar intensity interferometry measures $g^{(2)}-1=V^2$, while the full $g^{(2)}$ includes the independent-count reference one. Unpolarized thermal light with two independent equally populated wave polarizations has half this excess: adding their covariances doubles the numerator while multiplying total means quadruples the denominator. Finite detector response further suppresses the peak and must be calibrated.^[ref]