Expansion of the Glauber equation of motion in terms of cumulants powers

We present an expansion of the Glauber equation of motion in terms of cumulants of the effective fields. We show that we recover a part of the fluctuations with a simple third order expansion. We illustrate the method with an application to the 2D square lattice Ising model with nearest neighbors. For the special case of perturbed non gaussian effective fields, we show that only the first cumulants are significative. Furthermore it can be shown that the resulting distribution is closer to the normal distribution than the unperturbed fields one. It leads to the conclusion that we can consistently truncate the expansion of the exact Glauber equation of motion.