Magnetization and susceptibility of ferrofluids

Effects of interactions on magnetization relaxation dynamics in ferrofluids

The dynamics of magnetization relaxation in ferrofluids are studied with statistical-mechanical theory and Brownian dynamics simulations. The particle dipole moments are initially perfectly aligned, and the magnetization is equal to its saturation value. The magnetization is then allowed to decay under zero-field conditions toward its equilibrium value of zero. The time dependence is predicted by solving the Fokker-Planck equation for the one-particle orientational distribution function. Interactions between particles are included by introducing an effective magnetic field acting on a given particle and arising from all of the other particles. Two different approximations are proposed and tested against simulations: a first-order modified mean-field theory and a modified Weiss model. The theory predicts that the short-time decay is characterized by the Brownian rotation time τ_{B}, independent of the interaction strength. At times much longer than τ_{B}, the asymptotic decay time is predicted to grow with increasing interaction strength. These predictions are borne out by the simulations. The modified Weiss model gives the best agreement with simulation, and its range of validity is limited to moderate, but realistic, values of the dipolar coupling constant.

Electric Double Layer and Orientational Ordering of Water Dipoles in Narrow Channels within a Modified Langevin Poisson-Boltzmann Model

The electric double layer (EDL) is an important phenomenon that arises in systems where a charged surface comes into contact with an electrolyte solution. In this work we describe the generalization of classic Poisson-Boltzmann (PB) theory for point-like ions by taking into account orientational ordering of water molecules. The modified Langevin Poisson-Boltzmann (LPB) model of EDL is derived by minimizing the corresponding Helmholtz free energy functional, which includes also orientational entropy contribution of water dipoles. The formation of EDL is important in many artificial and biological systems bound by a cylindrical geometry. We therefore numerically solve the modified LPB equation in cylindrical coordinates, determining the spatial dependencies of electric potential, relative permittivity and average orientations of water dipoles within charged tubes of different radii. Results show that for tubes of a large radius, macroscopic (net) volume charge density of coions and counterions is zero at the geometrical axis. This is attributed to effective electrolyte charge screening in the vicinity of the inner charged surface of the tube. For tubes of small radii, the screening region extends into the whole inner space of the tube, leading to non-zero net volume charge density and non-zero orientational ordering of water dipoles near the axis.

Effect of an external field on capillary waves in a dipolar fluid

The role of an external field on capillary waves at the liquid-vapor interface of a dipolar fluid is investigated using molecular dynamics simulations. For fields parallel to the interface, the interfacial width squared increases linearly with respect to the logarithm of the size of the interface across all field strengths tested. The value of the slope decreases with increasing field strength, indicating that the field dampens the capillary waves. With the inclusion of the parallel field, the surface stiffness increases with increasing field strength faster than the surface tension. For fields perpendicular to the interface, the interfacial width squared is linear with respect to the logarithm of the size of the interface for small field strengths, and the surface stiffness is less than the surface tension. Above a critical field strength that decreases as the size of the interface increases, the interface becomes unstable due to the increased amplitude of the capillary waves.

Modified mean-field theory of the magnetic properties of concentrated, high-susceptibility, polydisperse ferrofluids

The effects of particle-size polydispersity on the magnetostatic properties of concentrated ferrofluids are studied using theory and computer simulation. The second-order modified mean-field (MMF2) theory of Ivanov and Kuznetsova [Phys. Rev. E 64, 041405 (2001)1063-651X10.1103/PhysRevE.64.041405] has been extended by calculating additional terms of higher order in the dipolar coupling constant in the expansions of the initial magnetic susceptibility and the magnetization curve. The theoretical predictions have been tested rigorously against results from Monte Carlo simulations of model monodisperse, bidisperse, and highly polydisperse ferrofluids. Comparisons have been made between systems with the same Langevin susceptibility and the same saturation magnetization. In all cases, the new theoretical magnetization curve shows better agreement with simulation data than does the MMF2 theory. As for the initial susceptibility, MMF2 theory is most accurate for the monodisperse model, while the new theory works best for polydisperse systems with a significant proportion of large particles. These results are important for the analysis and characterization of recently synthesized polydisperse ferrofluids with record-breaking values of the initial magnetic susceptibility.

Linear and nonlinear magnetic properties of ferrofluids

Within a high-magnetic-field approximation, employing Ruelle's algebraic perturbation theory, a field-dependent free-energy expression is proposed which allows one to determine the magnetic properties of ferrofluids modeled as dipolar hard-sphere systems. We compare the ensuing magnetization curves, following from this free energy, with those obtained by Ivanov and Kuznetsova [Phys. Rev. E 64, 041405 (2001)] as well as with new corresponding Monte Carlo simulation data. Based on the power-series expansion of the magnetization, a closed expression for the magnetization is also proposed, which is a high-density extension of the corresponding equation of Ivanov and Kuznetsova. From both magnetization equations the zero-field susceptibility expression due to Tani et al. [Mol. Phys. 48, 863 (1983)] can be obtained, which is in good agreement with our MC simulation results. From the closed expression for the magnetization the second-order nonlinear magnetic susceptibility is also derived, which shows fair agreement with the corresponding MC simulation data.

Liquid-vapor interface of the Stockmayer fluid in a uniform external field

The effect of a uniform (nonspatially varying) external field on the liquid-vapor interface of the Stockmayer fluid (Lennard-Jones particles embedded with a point dipole) has been investigated by molecular-dynamics simulations. The long-ranged parts of both the dipole and Lennard-Jones interactions are treated using an Ewald summation, which removes the effects of the cutoff. The direction of the field shifts the critical point and interfacial properties in different directions. For an external field parallel to the interface, the critical temperature increases, while for a field applied perpendicular to the interface, it decreases. The effects of the field on surface tension and interfacial width are also investigated. For zero field, dipoles near the liquid-vapor interface show a weak orientation parallel to the interface. For fields parallel to the interface, ordering in the liquid phase is greater than the vapor, while for fields perpendicular to the interface, the opposite is true.

The effects of polydispersity on the initial susceptibilities of ferrofluids

The effects of particle-size polydispersity on the initial susceptibilities of concentrated ferrofluids are analyzed using a combination of theory and computer simulation. The study is focused on a model ferrofluid with a prescribed magnetic-core diameter distribution, a fixed non-magnetic surface layer (corresponding to a demagnetized layer and adsorbed surfactant) and a combination of dipolar and hard-core interactions. The non-trivial effects of polydispersity are identified by comparing the initial susceptibilities of monodisperse and polydisperse ferrofluids with the same Langevin susceptibility. The theory is based on a correction to the second-order modified mean-field theory arising from a formal Mayer-type cluster expansion; this correction is dependent on a parameter similar to the normal dipolar coupling constant, except that it contains a complicated double average over the particle-size distribution, which means that the initial susceptibility should depend significantly on polydispersity. Specifically, the theory predicts that the initial susceptibility is enhanced significantly by polydispersity. This prediction is tested rigorously against results from Monte Carlo simulations and is found to be robust. The qualitative agreement between theory and simulation is already satisfactory, but the quantitative agreement could be improved by a systematic extension of the cluster expansion. The overall conclusion is that polydispersity should be accounted for carefully in magnetogranulometric analyses of real ferrofluids.

Comparison between theory and simulations for the magnetization and the susceptibility of polydisperse ferrofluids

The influence of polydispersity on the magnetization of ferrofluids is studied based on a previously published magnetization equation of state (Szalai and Dietrich, 2011 J. Phys.: Condens. Matter 23 326004) and computer simulations. The polydispersity of the particle diameter is described by the gamma distribution function. Canonical ensemble Monte Carlo simulations have been performed in order to test these theoretical results for the initial susceptibility and the magnetization. The results for the magnetic properties of the polydisperse systems turn out to be in quantitative agreement with our present simulation data. In addition, we find good agreement between our theory and experimental data for magnetite-based ferrofluids.

An analytic solution for the magnetization of two-dimensional ferrofluids based on the mean spherical approximation

An analytic formula is derived for the magnetization of a two-dimensional dipolar hard disk fluid using a variational functional series expansion of the free energy as a function of the orientational distribution function. The excess term expressing the effect of the intermolecular forces is calculated on the basis of the mean spherical approximation. Comparison with our own Monte Carlo simulation data shows excellent agreement for large external fields and for the zero-field susceptibility. At intermediate field strengths, the agreement is satisfactory for moderate dipole moments and densities.

Decrease of permittivity of an electrolyte solution near a charged surface due to saturation and excluded volume effects

The dipole moment of a water molecule in liquid water differs from that of an isolated one because each molecule is further polarized by the electric field of its neighbours. In this work a formula for the spatial dependence of the relative permittivity of an electrolyte near a highly charged surface is obtained in which the mutual influence of the water molecules is taken into account by means of the cavity field. The orientational ordering of water dipoles is considered in the saturation regime. It is predicted that the relative permittivity of an electrolyte solution near the highly charged surface (i.e. in saturation regime) may be substantially decreased due to orientational ordering of water (saturation effect) and depletion of water molecules (excluded volume effect) due to accumulation of counterions.

Magnetization of multicomponent ferrofluids

The solution of the mean spherical approximation (MSA) integral equation for isotropic multicomponent dipolar hard sphere fluids without external fields is used to construct a density functional theory (DFT), which includes external fields, in order to obtain an analytical expression for the external field dependence of the magnetization of ferrofluidic mixtures. This DFT is based on a second-order Taylor series expansion of the free energy density functional of the anisotropic system around the corresponding isotropic MSA reference system. The ensuing results for the magnetic properties are in quantitative agreement with our canonical ensemble Monte Carlo simulation data presented here.

Behavior of bulky ferrofluids in the diluted low-coupling regime: theory and simulation

A theoretical formalism to predict the structure factors observed in dipolar soft-sphere fluids based on a virial expansion of the radial distribution function is presented. The theory is able to account for cases with and without externally applied magnetic fields. A thorough comparison of the theoretical results to molecular-dynamics simulations shows a good agreement between theory and numerical simulations when the fraction of particles involved in clustering is low; i.e., the dipolar coupling parameter is lambda less, similar 2, and the volume fraction is phi less, similar 0.25. When magnetic fields are applied to the system, special attention is paid to the study of the anisotropy of the structure factor. The theory reasonably accounts for the structure factors when the Langevin parameter is smaller than 5.

Dipolar particles in an external field: Molecular dynamics simulation and mean field theory

Using molecular dynamics computer simulation we compute gas-liquid phase coexistence curves for the Stockmayer fluid in an external electric field. We observe a field-induced shift of the critical temperature DeltaTc. The sign of DeltaTc depends on whether the potential or the surface charge density is held constant assuming that the dielectric material fills the space between capacitor plates. Our own as well as previous literature data for DeltaTc are compared to and interpreted in terms of a simple mean field theory. Despite considerable errors in the simulation results, we find consistency between the simulation results obtained by different groups including our own and the mean field description. The latter ties the sign of DeltaTc to the outside constraints via the electric field dependence of the orientation part of the mean field free energy.