Magnetization of ferrofluids with dipolar interactions: A born-mayer expansion

Correlating the Dipolar Interactions Induced Magneto-Viscoelasticity and Thermal Conductivity Enhancements in Nanomagnetic Fluids

The effective thermal management of electronic system holds the key to maximize their performance. The recent miniaturization trends require a cooling system with high heat flux capacity, localized cooling, and active control. Nanomagnetic fluids (NMFs) based cooling systems have the ability to meet the current demand of the cooling system for the miniaturized electronic system. However, the thermal characteristics of NMFs have a long way to go before the internal mechanisms are well understood. This review mainly focuses on the three aspects to establish a correlation between the thermal and rheological properties of the NMFs. First, the background, stability, and factors affecting the properties of the NMFs are discussed. Second, the ferrohydrodynamic equations are introduced for the NMFs to explain the rheological behavior and relaxation mechanism. Finally, different theoretical and experimental models are summarized that explain the thermal characteristics of the NMFs. Thermal characteristics of the NMFs are significantly affected by the morphology and composition of the magnetic nanoparticles (MNPs) in NMFs as well as the type of carrier liquids and surface functionalization that also influences the rheological properties. Thus, understanding the correlation between the thermal characteristics of the NMFs and rheological properties helps develop cooling systems with improved performance.

Ferrofluids and bio-ferrofluids: looking back and stepping forward

Ferrofluids investigated along for about five decades are ultrastable colloidal suspensions of magnetic nanoparticles, which manifest simultaneously fluid and magnetic properties. Their magnetically controllable and tunable feature proved to be from the beginning an extremely fertile ground for a wide range of engineering applications. More recently, biocompatible ferrofluids attracted huge interest and produced a considerable increase of the applicative potential in nanomedicine, biotechnology and environmental protection. This paper offers a brief overview of the most relevant early results and a comprehensive description of recent achievements in ferrofluid synthesis, advanced characterization, as well as the governing equations of ferrohydrodynamics, the most important interfacial phenomena and the flow properties. Finally, it provides an overview of recent advances in tunable and adaptive multifunctional materials derived from ferrofluids and a detailed presentation of the recent progress of applications in the field of sensors and actuators, ferrofluid-driven assembly and manipulation, droplet technology, including droplet generation and control, mechanical actuation, liquid computing and robotics.

Generic theory of the dynamic magnetic response of ferrofluids

Ferrofluids belong to an important class of highly functional soft matter, benefiting from their magnetically controllable physical properties. Therefore, it is of central importance to quantitatively predict the dynamic magnetic response of ferrofluids. Traditional dynamic theories, however, are often restricted to the near-equilibrium regime and/or only apply to nearly ideal ferrofluids that are monodisperse, dilute enough, and weakly interacting. In this paper I develop a self-consistent and nonperturbative dynamical mean field theory for typical ferrofluids which are often polydisperse, concentrated, and strongly interacting, possibly driven far from equilibrium. I obtain a general nonperturbative expression for the dynamic magnetic susceptibility, quantitatively agreeing with the spectra obtained from Brownian Dynamics simulations on both mono- and bidisperse samples. Furthermore, I derive a generic magnetization relaxation equation (MRE) for both mono- and polydisperse ferrofluids by employing the projection operator technique in nonequlibrium statistical mechanics. This MRE is in simple closed form and independent of which model is employed to approximate the equilibrium magnetization curve. Existing models can be recovered as low-order approximations of my generic and nonperturbative MRE. My theory can play a key role in studying the dynamics of ferrofluids and other polar fluids. It may also have substantial and immediate consequences to various ferrofluid applications.

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.

Equilibrium properties of assembly of interacting superparamagnetic nanoparticles

The stochastic Landau-Lifshitz equation is used to investigate the relaxation process and equilibrium magnetization of interacting assembly of superparamagnetic nanoparticles (SPMNPs) uniformly distributed in a nonmagnetic matrix. For weakly interacting assembly, the equilibrium magnetization is shown to deviate significantly from the Langevin law at moderate and large magnetic fields under the influence of their magnetic anisotropies. For dense assemblies with noticeable influence of the magneto-dipole interaction, a significant dependence of the initial susceptibility on the assembly density is revealed. The difference between the initial susceptibility and the corresponding Langevin susceptibility can serve as an indication of appreciable influence of the magneto-dipole interaction on the assembly properties. A new self-consistent approach is developed to explain the effect of mutual magneto-dipole interaction on the behavior of dense assembly of SPMNPs. The probability densities of the components of random magnetic field acting on magnetic NPs are calculated at thermodynamic equilibrium. The self-consistent probability densities of these components are found to be close to Gaussian distribution. A decreasing equilibrium assembly magnetization as a function of its density can be explained as a disorienting effect of the random magnetic field on the NPs magnetic moments.

Dynamics of interacting magnetic nanoparticles: effective behavior from competition between Brownian and Néel relaxation

The intriguing properties of magnetic nanoparticles have sparked a growing number of theoretical studies as well as practical applications. Here, we provide the first comprehensive study of the influence of interactions on the two main relaxation mechanisms: internal (Néel) and Brownian relaxation. While non-interacting magnetic nanoparticles show Debye behavior with an effective relaxation time, many authors use this model also for the interacting case. Since Néel relaxation is typically a thermally activated process on times scales that are many orders of magnitude larger than the underlying micromagnetic times, we use extensive computer simulations employing a Brownian dynamics/Monte-Carlo algorithm to show that dipolar interactions lead to significant deviations from the Debye behavior. We find that Néel and Brownian relaxation can be considered as independent processes for short enough times until dipolar interactions lead to a coupling of these mechanisms, making the interpretation more difficult. We provide mean-field arguments that describe these short and long-time, effective relaxation times well for weak up to moderate interaction strengths. Our findings about the coupling of Brownian and Néel process and the effective relaxation time provide an important theoretical insight that will have also important consequences for the interpretation of magnetic susceptibility measurements and magnetorelaxometry analysis.

Dynamic susceptibility of a concentrated ferrofluid: The role of interparticle interactions

The dynamic susceptibility of concentrated ferrofluids of magnetite-kerosene type is studied experimentally to clarify the effect of interparticle interactions on the magnetization reversal dynamics and the ferrofluid relaxation time spectrum. We synthesize six ferrofluid samples, four of which have the same wide particle size distribution with a high (more than 2kT) average energy of magnetic dipole interactions. These samples differ in particle concentration and dynamic viscosity. The two remaining samples have a lower content of large particles and a moderate energy of magnetic dipole interactions. For all samples, we measure the dynamic susceptibility in the weak probing field at frequencies up to 160 kHz and the field amplitude dependence of the susceptibility at a frequency of 27 kHz. The results show that the susceptibility dispersion at frequencies up to 10 kHz is due to the rotational diffusion of colloidal particles and aggregates. Steric and hydrodynamic interparticle interactions are the main reason for the strong concentration dependence of the viscosity and so they also strongly influence the frequency dependence of the susceptibility. The influence of van der Waals and magnetic dipole interactions on the susceptibility is manifested indirectly, through the formation of multiparticle clusters. The contribution of clusters to the low-frequency susceptibility reaches 80%. Their large sizes (about 100 nm) shift the dispersion region to frequencies of 1-100 Hz, depending on the temperature and particle concentration. Experiments at 27 kHz demonstrate the increase in the dynamic susceptibility with increasing field amplitude. This growth is unexpected since all spectral amplitudes in the Debye function expansion of the dynamic susceptibility decrease monotonically with increasing field. To clarify the situation, the auxiliary problem of the magnetodynamics of a uniaxial particle in the alternating field is solved numerically. The Fokker-Planck-Brown rotational diffusion equation is used. It is shown that an increase in the field amplitude reduces the anisotropy barrier and the Néel relaxation time of particles and increases the dynamic susceptibility by one to two orders of magnitude compared to the weak-field limit. The calculation results are in qualitative agreement with the experimental data and allow us to propose a consistent interpretation of these data. We find that the increase in dynamic susceptibility with increasing amplitude is observed when two necessary conditions are met: (i) The suspension viscosity and the field frequency are high enough to cause the blocking of the rotational degrees of freedom of particles and aggregates and (ii) particles with a large magnetic anisotropy are present in the ferrofluid.

Dimension-shifting multifunctional biocompatible nanocomposites

A facile method to prepare dimension-shifting biocompatible multifunctional nanocomposites is described. The design is based on magnetic - and electrostatic - induced transitions from the dispersed state to the assembled state of zero-dimensional nanoparticles, resulting in dimension conversion.

Polymer architecture of magnetic gels: a review

In this review article, we provide an introduction to ferrogels, i.e. polymeric gels with embedded magnetic particles. Due to the interplay between magnetic and elastic properties of these materials, they are promising candidates for engineering and biomedical applications such as actuation and controlled drug release. Particular emphasis will be put on the polymer architecture of magnetic gels since it controls the degrees of freedom of the magnetic particles in the gel, and it is important for the particle-polymer coupling determining the mechanisms available for the gel deformation in magnetic fields. We report on the different polymer architectures that have been realized so far, and provide an overview of synthesis strategies and experimental techniques for the characterization of these materials. We further focus on theoretical and simulational studies carried out on magnetic gels, and highlight their contributions towards understanding the influence of the gels' polymer architecture.

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.

Influence of dipolar interactions on the magnetic susceptibility spectra of ferrofluids

The frequency-dependent magnetic susceptibility of a ferrofluid is calculated under the assumption that the constituent particles undergo Brownian relaxation only. Brownian-dynamics simulations are carried out in order to test the predictions of a recent theory [A. O. Ivanov, V. S. Zverev, and S. S. Kantorovich, Soft Matter 12, 3507 (2016)1744-683X10.1039/C5SM02679B] that includes the effects of interparticle dipole-dipole interactions. The theory is based on the so-called modified mean-field approach and possesses the following important characteristics: in the low-concentration, noninteracting regime, it gives the correct single-particle Debye-theory results; it yields the exact leading-order results in the zero-frequency limit; it includes particle polydispersity correctly from the outset; and it is based on firm theoretical foundations allowing, in principle, systematic extensions to treat stronger interactions and/or higher concentrations. The theory and simulations are compared in the case of a model monodisperse ferrofluid, where the effects of interactions are predicted to be more pronounced than in a polydisperse ferrofluid. The susceptibility spectra are analyzed in detail in terms of the low-frequency behavior, the position of the peak in the imaginary (out-of-phase) part, and the characteristic decay time of the magnetization autocorrelation function. It is demonstrated that the theory correctly predicts the trends in all of these properties with increasing concentration and dipolar coupling constant, the product of which is proportional to the Langevin susceptibility χ_{L}. The theory is in quantitative agreement with the simulation results as long as χ_{L}≲1.

Temperature-dependent dynamic correlations in suspensions of magnetic nanoparticles in a broad range of concentrations: a combined experimental and theoretical study

We study the effects of temperature and concentration on the dynamic spectra of polydisperse magnetic nanoparticle suspensions.

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.

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.

Size effect and stability of polarized fluid phases

The existence of a ferroelectric fluid phase for systems of 1000-2000 dipolar hard or soft spheres is well established by numerical simulations. Theoretical approaches proposed to determine the stability of such a phase are either in qualitative agreement with the simulation results or disagree with them. Experimental results for systems of molecules or particles with large electric or magnetic dipole moments are also inconclusive. As a contribution to the question of existence and stability of a fluid ferroelectric phase this simulation work considers system sizes of the order of 10 000 particles, thus an order of magnitude larger than those used in previous studies. It shows that although ferroelectricity is not affected by an increase of system size, different spatial arrangements of the dipolar hard spheres in such a phase are possible whose free energies seem to differ only marginally.

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.

Characterization of the Grain Size in Ferromagnetic Colloids: Comparing Torsional-Pendulum Measurements with Standard Complementary Methods

A recently introduced shear-flow-free method for measuring the rotational viscosity of a resonantly forced torsional pendulum is used to determine the transverse magnetic relaxation time in magnetite and cobalt-based ferrofluids. From these data the average size of the ferromagnetic grains and their hydrodynamic diameter (core plus surfactant coating) are deduced under in-situ conditions, i.e. without diluting the sample. The reliability of the method is demonstrated by comparing the results with those of the complementary techniques of magneto-granulometry, X-ray diffraction, electron microscopy, and photon-correlation spectroscopy.

Equilibrium and Nonequilibrium Behaviour of Ferrofluids – Experiments and Theory

Selected results on the spatiotemporal behaviour of equilibrium and nonequilibrium properties of ferrofluids in different magnetic fields are reviewed. They have been obtained in the project B13 Transport, response and instability behaviour of ferrofluids of the SFB 277 by experiments and by various theoretical methods ranging from purely analytical calculations to full numerical approaches.

Magnetic properties of polydisperse ferrofluids: a critical comparison between experiment, theory, and computer simulation

Experimental magnetization curves for a polydisperse ferrofluid at various concentrations are examined using analytical theories and computer simulations with the aim of establishing a robust method for obtaining the magnetic-core diameter distribution function p(x). Theoretical expressions are fitted to the experimental data to yield the parameters of p(x). It is shown that the majority of available theories yield results that depend strongly on the ferrofluid concentration, even though the magnetic composition should be fixed. The sole exception is the second-order modified mean-field (MMF2) theory of Ivanov and Kuznetsova [Phys. Rev. E 64, 041405 (2001)] which yields consistent results over the full experimental range of ferrofluid concentration. To check for consistency, extensive molecular dynamics and Monte Carlo simulations are performed on systems with discretized versions of p(x) corresponding as closely as possible to that of the real ferrofluid. Essentially perfect agreement between experiment, theory, and computer simulation is demonstrated. In addition, the MMF2 theory provides excellent predictions for the initial susceptibility measured in simulations.

The dielectric virial expansion and the models of dipolar hard-sphere fluid

The virial expansion technique to determine the dielectric constant epsilon of dipolar hard-sphere fluid is developed. It is shown that the formalism allows to bring into agreement the results of Debye's, Onsager's, and Langevin's to the problem. The third virial coefficient of epsilon is considered as a series over dipolar parameter lambda=m(2)d(3)kT. The terms up to O(lambda(11)) are calculated analytically providing a correct description of the third virial coefficient for small and intermediate values of lambda (0<or=lambda<or=4). The results of the dielectric virial series are compared with the Monte Carlo data for epsilon found by Matyushov and Ladanyi [J. Chem. Phys. 110, 994 (1999)]. The theory is in agreement with simulations only at small values of lambda<or=2. At higher polarities, the virial series diverges. Realization of the renormalization procedure permits to enlarge the range of applicability of the virial series. In this way, the new expression for the dielectric constant as a function of two dipolar parameters, lambda and y=4 pi nm(2)9kT, has been found explicitly. The expression gives a perfect upper bound of the dielectric constant and is more reliable for determination of epsilon than the previously known ones.

Rigid, superparamagnetic chains of permanently linked beads coated with magnetic nanoparticles. Synthesis and rotational dynamics under applied magnetic fields

An inexpensive and versatile approach is reported for the synthesis of monodisperse magnetoresponsive rods of desired diameter, length, and magnetic susceptibility based on the confined alignment of magnetic beads in microchannels of selected channel height, followed by localized hydrolysis of sol-gel precursors within polyelectrolyte shells adsorbed on the beads. The layer-by-layer technique was used to coat the polystyrene beads with polyelectrolytes of alternating charge and with charged magnetic nanoparticles, and the polystyrene cores could be removed either by solvent dissolution or by calcination to form hollow-shelled chains. The reorientation dynamics of single and clustered chains following the application of an external magnetic field was evaluated theoretically, with favorable comparisons with the experimental data.

Magnetoviscosity of semidilute ferrofluids and the role of dipolar interactions: comparison of molecular simulations and dynamical mean-field theory

Extensive molecular simulations on a model ferrofluid are performed in order to study magnetoviscous and viscoelastic phenomena in semidilute ferrofluids. Simulation results of the nonequilibrium magnetization, shear viscosity, and normal stress differences are presented. Rotational and configurational contributions to the shear viscosity are analyzed and their influence on the magnetoviscous effect is discussed. The simplified model of noninteracting magnetic dipoles describes the nonequilibrium magnetization and the rotational viscosity, but does not account for configurational viscosity contributions and normal stress differences. Improved mean-field models that overcome these limitations show good agreement with the simulation results for weak dipolar interactions where the models should apply. Comparisons to simulation results for various interaction strengths allow us to determine the range of validity of the mean-field models.

Magnetization of polydisperse colloidal ferrofluids: effect of magnetostriction

We exploit magnetostriction in polydisperse ferrofluids in order to generate nonlinear responses and apply a thermodynamical method to derive the desired nonlinear magnetic susceptibility. For an ideal gas, this method has been demonstrated to be in excellent agreement with a statistical method. In the presence of a sinusoidal ac magnetic field, the magnetization of the polydisperse ferrofluid contains higher-order harmonics, which can be extracted analytically by using a perturbation approach. We find that the harmonics are sensitive to the particle distribution and the degree of field-induced anisotropy of the system. In addition, we find that the magnetization is higher in the polydisperse system than in the monodisperse one, as also found by a recent Monte Carlo simulation. Thus, it seems possible to detect the size distribution in a polydisperse ferrofluid by measuring the harmonics of the magnetization under the influence of magnetostriction.

Low-temperature susceptibility of concentrated magnetic fluids

The initial susceptibility of concentrated magnetic fluids (ferrocolloids) has been experimentally investigated at low temperatures. The results obtained indicate that the interparticle dipole-dipole interactions can increase the susceptibility by several times as compared to the Langevin value. It is shown that good agreement between recent theoretical models and experimental observations can be achieved by introducing a correction for coefficients in the series expansion of susceptibility in powers of density and aggregation parameter. A modified equation for equilibrium susceptibility is offered to sum over corrections made by Kalikmanov (Statistical Physics of Fluids, Springer-Verlag, Berlin, 2001) and by B. Huke and M. Lucke (Phys. Rev. E 67, 051403, 2003). The equation gives good quantitative agreement with the experimental data in the wide range of temperature and magnetic particles concentration. It has been found that in some cases the magnetic fluid solidification occurs at temperature several tens of kelvins higher than the crystallization temperature of the carrier liquid. The solidification temperature of magnetic fluids is independent of particle concentration (i.e., magneto-dipole interparticle interactions) and dependent on the surfactant type and carrier liquid. This finding allows us to suggest that molecular interactions and generation of some large-scale structure from colloidal particles in magnetic fluids are responsible for magnetic fluid solidification. If the magnetic fluid contains the particles with the Brownian relaxation mechanism of the magnetic moment, the solidification manifests itself as the peak on the "susceptibility-temperature" curve. This fact proves the dynamic nature of the observed peak: it arises from blocking the Brownian mechanism of the magnetization relaxation.

Applying the chain formation model to magnetic properties of aggregated ferrofluids

The magnetization properties of aggregated ferrofluids are calculated by combining the chain formation model developed by Zubarev with the modified mean-field theory. Using moderate assumptions for the inter- and intrachain interactions we obtain expressions for the magnetization and initial susceptibility. When comparing the results of our theory to molecular dynamics simulations of the same model we find that at large dipolar couplings (lambda>3) the chain formation model appears to give better predictions than other analytical approaches. This supports the idea that chain formation is an important structural ingredient of strongly interacting dipolar particles.

Nonlinear theory of pattern formation in ferrofluid films at high field strengths

When a magnetic field is applied to a thin layer of a suspension of magnetic nanoparticles (ferrofluid), the formation of labyrinthine and hexagonal patterns is observed. We introduce a theory to describe ferrofluid patterns at high field, where a nonlinear relationship between field and magnetization is expected. The computational difficulties due to the use of a nonlinear magnetization curve are solved by a reformulation of the magnetic energy equation. The evolution of the pattern size at intermediate and very high fields can be understood by an analysis of limiting cases of the magnetization curve. In particular, at a very high field the pattern size reaches a constant saturation value which has been recently confirmed by experiments. The field for the onset of a nonlinear behavior is shifted to higher field strength due to a demagnetization effect. This can partially explain the ability of linear approaches to reproduce experimental data even at a high field. Finally, the impact of the nonlinearity of the magnetization curve on the transition between hexagonal and labyrinthine patterns is discussed.

Structure and magnetic properties of polydisperse ferrofluids: a molecular dynamics study

We study by Langevin molecular dynamics simulations systematically the influence of polydispersity in the particle size, and subsequently in the dipole moment, on the physical properties of ferrofluids. The polydispersity is in a first approximation modeled by a bidisperse system that consists of small and large particles at different ratios of their volume fractions. In the first part of our investigations the total volume fraction of the system is fixed, and the volume fraction phiL of the large particles is varied. The initial susceptibility chi and magnetization curve of the systems show a strong dependence on the value of phiL. With the increase of phiL, the magnetization M of the system has a much faster increment at weak fields, and thus leads to a larger chi. We performed a cluster analysis that indicates that this is due to the aggregation of the large particles in the systems. The average size of these clusters increases with increasing phiL. In the second part of our investigations, we fixed the volume fraction of the large particles, and increased the volume fraction phiS of the small particles in order to study their influence on the chain formation of the large ones. We found that the average aggregate size formed by large particles decreases when phiS is increased, demonstrating a significant effect of the small particles on the structural properties of the system. A topological analysis of the structure reveals that the majority of the small particles remain nonaggregated. Only a small number of them are attracted to the ends of the chains formed by large particles.

Spontaneous ferromagnetic ordering in magnetic fluids

This paper is devoted to the theoretical justification of spontaneous orientational order in magnetic fluids. We study the self-consistent solutions of the Bogoliubov-Born-Green-Kirkwood-Yvon equation connecting the one-particle distribution function with the pair correlation function. This self-consistent approach is used in the specific density functional method and proves to be equivalent to the mean field theory. On the basis of the second-order perturbation method over the intensity of dipole-dipole interparticle interaction the following effect is discovered: the self-consistent density functional approach leads to the spontaneous "ferrimagnetic" state of the magnetic fluid induced by the dipole-dipole interaction. This strange result seems to be physically meaningless and prejudices the validity of the density functional methods and mean field theories applied to orientational microstructure in ferrofluids.

Magnetization of concentrated polydisperse ferrofluids: cluster expansion

The equilibrium magnetization of concentrated ferrofluids described by a system of polydisperse dipolar hard spheres is calculated as a function of the internal magnetic field using the Born-Mayer or cluster expansion technique. This paper extends the results of Huke and Lücke [Phys. Rev. E 62, 6875 (2000)] obtained for monodisperse ferrofluids. The magnetization is given as a power series expansion in two parameters related to the volume fraction and the coupling strength of the dipolar interaction, respectively.

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 ferrosolids the observed susceptibility is considerably lowered when compared to ferrofluids.