A molecular dynamics study on the equilibrium magnetization properties and structure of ferrofluids

We investigate in detail the initial susceptibility, magnetization curves, and microstructure of ferrofluids in various concentration and particle dipole moment ranges by means of molecular dynamics simulations. We use the Ewald summation for the long-range dipolar interactions, take explicitly into account the translational and rotational degrees of freedom, coupled to a Langevin thermostat. When the dipolar interaction energy is comparable with the thermal energy, the simulation results on the magnetization properties agree with the theoretical predictions very well. For stronger dipolar couplings, however, we find systematic deviations from the theoretical curves. We analyze in detail the observed microstructure of the fluids under different conditions. The formation of clusters is found to enhance the magnetization at weak fields and thus leads to a larger initial susceptibility. The influence of the particle aggregation is isolated by studying ferro-solids, which consist of magnetic dipoles frozen in at random locations but which are free to rotate. Due to the artificial suppression of clusters in ferro-solids the observed susceptibility is considerably lowered when compared to ferrofluids.Comment: 33 pages including 12 figures, requires RevTex

Cluster structure and the first-order phase transition in dipolar systems

The Monte Carlo technique is used to simulate a 3D dipolar hard-sphere system. The spatial and magnetic structure of clusters formed by magnetic dipolar interactions in zero applied field is investigated. It is shown that the many-particle clusters are characterized by a quasi-spherical shape, extremely small magnetic moments, and a fractal dimension close to three. These clusters are regarded as nuclei of a new concentrated isotropic phase. The numerical simulation of the first-order phase transition has been realized which allows us to find the interface between two coexisting phases. It has been found that the dipole-dipole and steric interactions are sufficient to separate the system into two phases with low and high concentrations of particles. The introduction of any additional attraction potential is not required. The phase diagram of dipolar system in zero applied field has been obtained. The simulation results are in qualitative agreement with the predictions of some analytical models