Statistical calculation for a model ferromagnetic liquid

A theoretical method is developed for treating model ferromagnetic liquids. In this method the spins are treated in the mean-field approximation while spatial correlations are retained in a manner similar to conventional theories of classical liquids. Coupled self-consistency equations for the magnetization and an effective radial distribution function are derived and solved. The method is applied to a primitive model characterized by a simple spin-dependent interparticle potential. Conditions under which such a simple liquid can exhibit ferromagnetic order are determined. Numerical examples show how the transition temperature and the temperature dependence of the magnetization can be calculated.