The covariance of a stationary process X is diagonalized by Fourier transform. It does not take into account the complex phase and defines Gaussian maximum entropy models. We introduce general family harmonic moments, which rely on phases to capture non-Gaussian properties. They are defined as Hˆ(LX), where L linear operator Hˆ non-linear multiplies each coefficient integers. can also be calcul...