  The dynamics of non-equilibrium Ising model with parallel updates is
investigated using a generalized mean field approximation that incorporates
multiple two-site correlations at any two time steps, which can be obtained
recursively. The proposed method shows significant improvement in predicting
local system properties compared to other mean field approximation techniques,
particularly in systems with symmetric interactions. Results are also evaluated
against those obtained from Monte Carlo simulations. The method is also
employed to obtain parameter values for the kinetic inverse Ising modeling
problem, where couplings and local fields values of a fully connected spin
system are inferred from data.
