Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

Ricardo Peredo Ortíz, Martin Hernández Contreras and Raquel Hernández Gómez

## Abstract

In this chapter, we review the experimental and theoretical modeling of structural and dynamical properties of colloidal magnetic fluids at equilibrium. Presently, several prototype experimental systems are very well characterized. We survey the different models, which help to reach a comprehensive knowledge of these complex magnetic fluids. One prime example is the ongoing investigation of the realistic interparticle potentials that drive the formation of the different phase states observed experimentally. Further, a stochastic equation approach for the description of tracer diffusion, viscoelasticity, and dielectric relaxation at equilibrium in colloidal ferrofluids is discussed.

Keywords: ferrofluids, viscoelasticity, colloidal magnetic fluids, magnetohydrodynamics, dielectric relaxation, diffusion stochastic dynamics

## 1. Introduction

Colloidal magnetic fluids are made up of nanometer and micron-size ferromagnetic particles dispersed either in electrolytes or in organic solvents [1]. Due to their easy manipulation with magnetic or electric fields, they are having a significant impact on diverse technological applications ranging from biomedicine [2] and photonic devices [3, 4] up to fundamental studies that motivate the development of modern theories of nonequilibrium condensed matter [5–10]. In this chapter, we review the studies of static structural and dynamical properties in colloidal magnetic fluids, which provide a comprehensive description of their bulk phase behavior at thermal equilibrium. Presently, model magnetic colloidal systems can be prepared with tunable interaction among particles of the hard sphere and longrange dipolar types [11–14], especially the colloidal system constituted of micronsize polystyrene spheres that are electrically charged and dispersed in organic solvents [11]. In such a system, particles acquire an induced electric dipole moment under the external static electric field. Confocal microscopy has allowed the determination of its three-dimensional bulk phase diagram as a function of applied electric field and by taking into account the presence of a 1:1 salt. The electric and magnetic dipole pair interaction potential is symmetric. Consequently, a corresponding change of physical units leads to an equivalent description of the magnetic fluid phase diagrams. Novel Monte Carlo simulation methods have

confirmed the different experimental phases [15, 16]. Integral equations and density functional theories have assisted in the understanding of the phases formed by molecular liquids [17], and they have been used to gain qualitative insight into the expected phases as a function of volume fraction and dipolar strength of colloidal magnetic fluids [18–23]. Likewise, the structure factor and diffusion coefficient of magnetic suspensions made of maghemite nanoparticles dispersed in water and equilibrium with electrolyte solutions have also been reported [24–31]. The measured structure factor provides the microscopic arrangements of particles accurately. Here, we provide a Langevin stochastic approach that allows the determination of the translational and rotational diffusion of the particles in ferrofluids [32, 33]. These dynamic properties, in turn, allow quantifying the viscoelastic and dielectric moduli of the structural evolution of the fluid toward their equilibrium states [34]. This approach may help to interpret recent experiments of passive microrheology [35, 36] and Altern Current spectroscopic techniques that measure the viscoelastic dynamics and dielectric relaxation of colloidal magnetic fluids. It is expected these theoretical methods be extended to help to understand the corresponding experimental observations of similar dynamics in two-dimensional paramagnetic colloids [37] and under external fields.

volume fraction of particles at the position r

DOI: http://dx.doi.org/10.5772/intechopen.81091

!

dr<sup>2</sup> ∇ϕ r ! h i � � <sup>2</sup>

<sup>4</sup> � � 3rij

dipole m^ <sup>i</sup> with field H

volume is 1ð Þ <sup>=</sup><sup>2</sup> <sup>Ð</sup>

Ð d<sup>2</sup> k !

43

a δm k � � � �! �

� � � 2

<sup>þ</sup> <sup>b</sup>δm k � �!

δϕ � k � �!

material parameters of ferrofluid external field m and ϕ. The uniform phase becomes unstable to the formation of modulated nonuniform phases if the determinant of the integrand gets negative at the minimum of wave number [38]. It has been recognized that ferromagnetic particles of 10 -nm size remain dispersed [39]; however, when particles have larger sizes up to the micron scale, experiments find they do assemble into chains, rings, and nets were several chains bound together into an amorphous structure [40]. A mean field theory that takes into account such topological structures was developed to explain the liquid-gas phase transition by the formation of Y-like bounds of chains in the network [41]. There remains a quantitative verification of this theory with its experimental

þ c δϕ k � � � �! �

<sup>2</sup> � � with a, b, c related to the

� � �

given by Vdd ¼ � <sup>m</sup><sup>2</sup>=rij

of the pair dipole-dipole interactions 1ð Þ <sup>=</sup><sup>2</sup> <sup>∑</sup>i,j ið Þ 6¼<sup>j</sup> Vdd <sup>r</sup>

Their model is a thin film of ferrofluid with a constant magnetic field acting perpendicular to the fluid layer. Since external fields produce an inhomogeneous distribution of particles, the volume fraction becomes a local function of position. FH is constructed at mean field level to have the following terms: the internal energy

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

energy Frot of an ideal gas of paramagnetic dipoles with polar angles as the only degrees of freedom under a constant magnetic field. Frot is obtained from the partition function in the canonical ensemble calculated with the total energy of

to all possible spatial configurations of the finite size particles in the volume available �kBT½ � ϕ ln ϕ þ ð Þ 1 � ϕ ln 1ð Þ � ϕ as given by the lattice gas. kB is the Boltzmann constant and the room temperature T. The energy of particles that arise from gradients in their concentration and therefore setting their average available local

lattice gas model. The ferroparticles are all equal, and their energy of pairwise interaction for two dipoles separated with the center-to-center vector distance r

dipole moment, and u^<sup>i</sup> is the unitary vector orientation of particle i. Thus, the resulting mean Gibbs energy G is mapped from the lattice gas model into a continuous position dependence representation. The free energy is a function of the local dipole moment and concentration m rð Þ, ϕð Þr . A ferrofluid in a continuous single phase such as a dispersed homogenous ferrofluid is an equilibrium state with an average concentration ϕ and bulk dipole moment m. This state is characterized by the minimum of the free energy, that is, <sup>∂</sup>G=∂<sup>m</sup> <sup>¼</sup> <sup>0</sup>, <sup>∂</sup>G=∂<sup>ϕ</sup> <sup>¼</sup> 0, which provides the average values of bulk dipole moment and concentration above. However, it is observed that for specific magnitudes of the external field and ferrofluid concentration, there appears chaining of particles that coalescence into sheets, which at higher field strength evolve to labyrinthine structures. The nucleation of these structures starts with the formation of magnetization domains associated with small clustering of particles with sizes on the order of a few microns. Short-wavelength thermal fluctuations of this size are capable of describing the formation of nonuniform phases. One way to analyze the formation of these phases uses the expansion of the free energy difference [38] ΔGð Þ¼ φ; m Gð Þ� φ; m Gð Þ φ; m up to second order in the fluctuations of dipole magnetization δm ¼ m rð Þ� m, and concentration δϕ ¼ ϕð Þ� r ϕ about their mean values. In the two-dimensional Fourier space representation, the free energy difference is a quadratic form <sup>Δ</sup><sup>G</sup> <sup>¼</sup> <sup>1</sup>=2π<sup>2</sup> ð Þ

! � u^<sup>j</sup> � u^<sup>i</sup> � u^<sup>i</sup>

!

interactions � ∑<sup>i</sup> m^ <sup>i</sup>� H

! � u^irij

! to have a total fixed value Ð

!

, which is provided by a continuous version of the

h i, where <sup>m</sup> is the particle's constant

dr<sup>2</sup> ϕ r !� �.

! is

ij � � and the rotational free

, the entropic energy contribution due

## 2. Experimental observations of structural properties of magnetic fluids

The importance of knowing the phase diagram of a ferrofluid resides on the information it provides about the undetermined underlying effective interaction of the ferromagnetic particles [13]. For maghemite nanoparticles in aqueous solutions, there have been attempts to determine such electrostatic potentials by using the measured parameters such as dipole strength and fluid density as inputs in Brownian dynamics simulations to reproduce the observed bulk structure factor [30, 31].

## 2.1 Importance of an experimental phase diagram of maghemite colloid

There is a set of experiments in a well-characterized system of electrically stabilized maghemite nanometer-size particles, which are dispersed in water [24–29]. Its stabilization is reached with the citrate electrolyte. By fine-tuning the salt concentration, the particle's interaction was shifted from repulsive, where they form a solid glass phase, into long-range attractive interactions that yield a fluid or gas state behavior. These experiments have prompted the determination of the corresponding pairwise interaction among particles [13]. The proposed interaction consists of a Yukawa repulsive part, a Van der Waals short-range attraction and includes an angular averaged attractive long-range pair dipole potential. There have not been attempts to predict the phase diagram with these potentials. However, the modeled interactions were used to determine with Brownian dynamics the observed structure factor [30, 31]. By using a generic model of pair dipole interactions together with an effective attraction associated with density gradients, Lacoste, et al. [38] considered a quasi-two-dimensional layer of a dilute ferrofluid subjected to a perpendicular magnetic field, and with the help of a mean field energy approach, they found modulated phases given by stripes and hexagonal formations. The Gibbs free energy they considered for inhomogeneous systems has the general shape G ¼ FH � μϕ r ! , where FH is the Helmholtz energy; the chemical potential μ constrains the average local composition ϕ r ! or space-dependent

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids DOI: http://dx.doi.org/10.5772/intechopen.81091

volume fraction of particles at the position r ! to have a total fixed value Ð dr<sup>2</sup> ϕ r !� �. Their model is a thin film of ferrofluid with a constant magnetic field acting perpendicular to the fluid layer. Since external fields produce an inhomogeneous distribution of particles, the volume fraction becomes a local function of position. FH is constructed at mean field level to have the following terms: the internal energy of the pair dipole-dipole interactions 1ð Þ <sup>=</sup><sup>2</sup> <sup>∑</sup>i,j ið Þ 6¼<sup>j</sup> Vdd <sup>r</sup> ! ij � � and the rotational free energy Frot of an ideal gas of paramagnetic dipoles with polar angles as the only degrees of freedom under a constant magnetic field. Frot is obtained from the partition function in the canonical ensemble calculated with the total energy of dipole m^ <sup>i</sup> with field H ! interactions � ∑<sup>i</sup> m^ <sup>i</sup>� H ! , the entropic energy contribution due to all possible spatial configurations of the finite size particles in the volume available �kBT½ � ϕ ln ϕ þ ð Þ 1 � ϕ ln 1ð Þ � ϕ as given by the lattice gas. kB is the Boltzmann constant and the room temperature T. The energy of particles that arise from gradients in their concentration and therefore setting their average available local volume is 1ð Þ <sup>=</sup><sup>2</sup> <sup>Ð</sup> dr<sup>2</sup> ∇ϕ r ! h i � � <sup>2</sup> , which is provided by a continuous version of the lattice gas model. The ferroparticles are all equal, and their energy of pairwise interaction for two dipoles separated with the center-to-center vector distance r ! is given by Vdd ¼ � <sup>m</sup><sup>2</sup>=rij <sup>4</sup> � � 3rij ! � u^irij ! � u^<sup>j</sup> � u^<sup>i</sup> � u^<sup>i</sup> h i, where <sup>m</sup> is the particle's constant dipole moment, and u^<sup>i</sup> is the unitary vector orientation of particle i. Thus, the resulting mean Gibbs energy G is mapped from the lattice gas model into a continuous position dependence representation. The free energy is a function of the local dipole moment and concentration m rð Þ, ϕð Þr . A ferrofluid in a continuous single phase such as a dispersed homogenous ferrofluid is an equilibrium state with an average concentration ϕ and bulk dipole moment m. This state is characterized by the minimum of the free energy, that is, <sup>∂</sup>G=∂<sup>m</sup> <sup>¼</sup> <sup>0</sup>, <sup>∂</sup>G=∂<sup>ϕ</sup> <sup>¼</sup> 0, which provides the average values of bulk dipole moment and concentration above. However, it is observed that for specific magnitudes of the external field and ferrofluid concentration, there appears chaining of particles that coalescence into sheets, which at higher field strength evolve to labyrinthine structures. The nucleation of these structures starts with the formation of magnetization domains associated with small clustering of particles with sizes on the order of a few microns. Short-wavelength thermal fluctuations of this size are capable of describing the formation of nonuniform phases. One way to analyze the formation of these phases uses the expansion of the free energy difference [38] ΔGð Þ¼ φ; m Gð Þ� φ; m Gð Þ φ; m up to second order in the fluctuations of dipole magnetization δm ¼ m rð Þ� m, and concentration δϕ ¼ ϕð Þ� r ϕ about their mean values. In the two-dimensional Fourier space representation, the free energy difference is a quadratic form <sup>Δ</sup><sup>G</sup> <sup>¼</sup> <sup>1</sup>=2π<sup>2</sup> ð Þ Ð d<sup>2</sup> k ! a δm k � � � �! � � � � 2 <sup>þ</sup> <sup>b</sup>δm k � �! δϕ � k � �! þ c δϕ k � � � �! � � � � <sup>2</sup> � � with a, b, c related to the

material parameters of ferrofluid external field m and ϕ. The uniform phase becomes unstable to the formation of modulated nonuniform phases if the determinant of the integrand gets negative at the minimum of wave number [38].

It has been recognized that ferromagnetic particles of 10 -nm size remain dispersed [39]; however, when particles have larger sizes up to the micron scale, experiments find they do assemble into chains, rings, and nets were several chains bound together into an amorphous structure [40]. A mean field theory that takes into account such topological structures was developed to explain the liquid-gas phase transition by the formation of Y-like bounds of chains in the network [41]. There remains a quantitative verification of this theory with its experimental

confirmed the different experimental phases [15, 16]. Integral equations and density functional theories have assisted in the understanding of the phases formed by molecular liquids [17], and they have been used to gain qualitative insight into the expected phases as a function of volume fraction and dipolar strength of colloidal magnetic fluids [18–23]. Likewise, the structure factor and diffusion coefficient of magnetic suspensions made of maghemite nanoparticles dispersed in water and equilibrium with electrolyte solutions have also been reported [24–31]. The measured structure factor provides the microscopic arrangements of particles accurately. Here, we provide a Langevin stochastic approach that allows the determination of the translational and rotational diffusion of the particles in ferrofluids [32, 33]. These dynamic properties, in turn, allow quantifying the viscoelastic and dielectric moduli of the structural evolution of the fluid toward their equilibrium states [34]. This approach may help to interpret recent experiments of passive microrheology [35, 36] and Altern Current spectroscopic techniques that measure the viscoelastic dynamics and dielectric relaxation of colloidal magnetic fluids. It is

expected these theoretical methods be extended to help to understand the

paramagnetic colloids [37] and under external fields.

Pattern Formation and Stability in Magnetic Colloids

[30, 31].

42

corresponding experimental observations of similar dynamics in two-dimensional

2. Experimental observations of structural properties of magnetic fluids

The importance of knowing the phase diagram of a ferrofluid resides on the information it provides about the undetermined underlying effective interaction of the ferromagnetic particles [13]. For maghemite nanoparticles in aqueous solutions, there have been attempts to determine such electrostatic potentials by using the measured parameters such as dipole strength and fluid density as inputs in Brownian dynamics simulations to reproduce the observed bulk structure factor

2.1 Importance of an experimental phase diagram of maghemite colloid

modeled interactions were used to determine with Brownian dynamics the

!

cal potential μ constrains the average local composition ϕ r

the general shape G ¼ FH � μϕ r

observed structure factor [30, 31]. By using a generic model of pair dipole interactions together with an effective attraction associated with density gradients, Lacoste, et al. [38] considered a quasi-two-dimensional layer of a dilute ferrofluid subjected to a perpendicular magnetic field, and with the help of a mean field energy approach, they found modulated phases given by stripes and hexagonal formations. The Gibbs free energy they considered for inhomogeneous systems has

, where FH is the Helmholtz energy; the chemi-

!

or space-dependent

There is a set of experiments in a well-characterized system of electrically stabilized maghemite nanometer-size particles, which are dispersed in water [24–29]. Its stabilization is reached with the citrate electrolyte. By fine-tuning the salt concentration, the particle's interaction was shifted from repulsive, where they form a solid glass phase, into long-range attractive interactions that yield a fluid or gas state behavior. These experiments have prompted the determination of the corresponding pairwise interaction among particles [13]. The proposed interaction consists of a Yukawa repulsive part, a Van der Waals short-range attraction and includes an angular averaged attractive long-range pair dipole potential. There have not been attempts to predict the phase diagram with these potentials. However, the counterpart and with computer simulations. Another method to validate the proposed model of particle's interaction is based on the parameters that appear in the Lennard-Jones potential that represents the Van der Waals attractions together with the pair dipole-dipole potential, and by fitting the predicted magnetization curve as a function of the applied magnetic field to the experimental one [42]. Even though there is not yet a realistic model interaction potential that represents the interaction of ferromagnetic particles, the use of fitting parameters as mentioned above has led to the prediction of the structure factor and the birefringence as a function of the applied field, which shows qualitative agreement with the experimental values [22, 29–31].

Lennard-Jones interaction energy VLJð Þ¼ r 4ε<sup>0</sup> ð Þ d=r

DOI: http://dx.doi.org/10.5772/intechopen.81091

simulation and optimization techniques

tial depth that can be estimated experimentally. For ferrofluid in an aqueous solvent of dielectric constant ε and monodisperse spherical particles, due to the presence of

2.3 Theoretical models of phase diagrams of dipolar colloidal fluids: computer

Finite-temperature calculations of ground state free energies with Yukawa repulsion and dipolar-dipolar interactions of varying strength were performed by Hynninen and Dijkstra [15]. In their calculations, they used canonical Monte Carlos simulations to obtain the free energy minimum as a function of volume fraction and dipolar strength in three-dimensional model systems. They found a phase diagram that contains the same phases as was observed by Yethiraj et al. [11]. However, additionally, they predicted a new hexagonal close-packed phase at the hard sphere repulsion limit and body-centered orthorhombic when the repulsion becomes long range by lowering the content of electrolyte in solution. Their method allows quicker and reliable determinations of the three-dimensional phase diagram than the evolutionary algorithm method that searches for the ground state at zero temperature for two-dimensional systems [45]. The bulk phase diagram has not been searched yet with genetic algorithm (GA) at finite temperature. A particular detail of GA is that to look for the minimum of the free energy, the derivative of the energy function at each evolution step is required for finding the best-adapted crystal structures (individuals) in the population, which is made up of several crystal structures produced by the algorithm. The search of the best-adapted individuals (crystal phase) is done in the energy landscape. For fixed pressure and finite temperature, the Gibbs free energy as a function of pairwise interaction strength and particles concentration in the fluid yields a three-dimensional plot of the free energy known as energy landscape, which has local minima and maxima. The global minimum of an equilibrium crystal structure is reached with the assistance of a local gradient algorithm that requires the derivative of the energy. Such a procedure makes an indirect search of the crystal structure by sampling the energy landscape. Thus, it becomes technically more involved to implement than the Hynninen et al.

method of MC ground state energy searches. In Figure 1, we provide two

algorithm approach at zero temperature and fixed pressure.

2.4 Experimental structure factors and diffusion coefficients

at a finite temperature also.

45

nonequilibrium states of local minima that result from the application of a genetic

Recently, Spiteri and Messina [16] have proposed an efficient nonlinear optimization technique that allows predicting the crystal phases of dipolar colloids at zero temperature without the need to use derivatives of the free energies. It was found that for a monodisperse repulsive hard sphere plus dipole-dipole interaction, the predicted phase diagram has the same phases as was previously found by Hynninen et al. [15] and that a knew so-called clinohexagonal prism span all known ground state structures at any density. There remains to be verified if this prediction fulfills

The microstructural order inside of the ferrofluid is determined by the potential interaction among particles. The measured structure factor yields the details of the particles'spatial geometric arrangement. The water-based maghemite ferrofluid is

salt at density <sup>ρ</sup>i, the Yukawa potential Vyukð Þ¼ <sup>r</sup> <sup>Q</sup><sup>∗</sup> exp ½ � �κð Þ <sup>r</sup> � <sup>d</sup> <sup>=</sup>r, <sup>κ</sup><sup>2</sup> <sup>≔</sup> 4πe2ρi=εkBT is used instead, where e is the electron charge. Q<sup>∗</sup> is the effective electric charge of a particle that can be measured with electrophoresis experiments.

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

<sup>12</sup> � ð Þ <sup>d</sup>=<sup>r</sup>

<sup>6</sup> h i with <sup>ε</sup><sup>0</sup> the poten-

## 2.2 Experimental phase diagram of a charged and sterically stabilized electrorheological colloidal suspension under an electric field

Yethiraj et al. [11] made a colloidal system of charged and sterically stabilized polymethyl methacrylate spheres of micron size in an organic solvent to have dipole moments on particles induced by an electric field. This monodisperse suspension shows both long-range repulsion, and attractive anisotropic interaction potential that is fixed with the addition of salt, whereas the dipolar interaction is controlled by the external electric field. For low Reynolds numbers and in an infinite fluid, a solid spherical particle that sinks in the fluid reaches a terminal velocity due to a balance of the friction force exerted by the fluid that opposes the gravitational force on the particle. If the fluid has a dynamic viscosity η and mass density ρ<sup>f</sup> , and the spherical particle possesses a mass density ρm, then the steady velocity is V ¼ ρ<sup>f</sup> � ρ<sup>m</sup> � �gd<sup>2</sup> =18η, where g is the gravitational acceleration and d the particle diameter. Thus, the effect of gravity is diminished by matching the density of the particles to that of the solvent. The Van der Waals attraction that otherwise would produce aggregation was reduced by matching the index of refraction to the visible. The Van der Waals attractive energy of two spherical particles is Vvdwð Þ¼� r ð Þ A=12 d2 <sup>=</sup> <sup>r</sup><sup>2</sup> � <sup>d</sup><sup>2</sup> � � <sup>þ</sup> <sup>d</sup><sup>2</sup> <sup>=</sup>r<sup>2</sup> <sup>þ</sup> 2 ln 1 � <sup>d</sup><sup>2</sup> =r<sup>2</sup> � � � � , where r is the separation distance of the spheres' centers [43]. Therefore, the Hamaker constant [44] A ¼ 3hν np � ns � �<sup>2</sup> np þ ns � �<sup>2</sup> =16 ffiffi 2 <sup>p</sup> np <sup>2</sup> <sup>þ</sup> ns <sup>2</sup> � �<sup>3</sup>=<sup>2</sup> can be made zero if the optical refractive index of the colloidal particles np is equal to that of the solvent ns for an appropriated selection of the frequency ν of the visible light. Here, h is the Planck constant. With confocal microscopy studies, they observed in real-space representation crystalline structures controlled by the repulsive part of the potential, which induces a structure of random hexagonal close packing. For longer-range repulsion, it becomes facecentered cubic (fcc), and one can also find a body-centered cubic phase (bcc). As a function of the electric field, a dipolar phase diagram was established. For low field, there is a sequence of fluid-bcc-fcc (phase-centered cubic) phases. For larger values of the field, there is string formation that coalesces into sheets that in turn coarsen to form body-centered tetragonal crystallites. Moreover, for increasing volume fractions, the string fluid crystallizes to a space-filling tetragonal phase. Such a rich variety of phases can be generated controllably by the concentration of the electrolyte. High salt concentration screens the electrostatic repulsion becoming shortdistance hard sphere interaction whereas the small content of salt increases the spatial range of the potential. The dipolar interaction is controlled by the external electric field. Consequently, due to the symmetric form of the electric and magnetic dipolar potential, a simple change of physical units leads to the similar conclusions to be expected for a colloidal magnetic ferrofluid. For oily ferrofluids, the shortrange Van der Waals attraction and hardcore repulsion are modeled by the

Lennard-Jones interaction energy VLJð Þ¼ r 4ε<sup>0</sup> ð Þ d=r <sup>12</sup> � ð Þ <sup>d</sup>=<sup>r</sup> <sup>6</sup> h i with <sup>ε</sup><sup>0</sup> the potential depth that can be estimated experimentally. For ferrofluid in an aqueous solvent of dielectric constant ε and monodisperse spherical particles, due to the presence of salt at density <sup>ρ</sup>i, the Yukawa potential Vyukð Þ¼ <sup>r</sup> <sup>Q</sup><sup>∗</sup> exp ½ � �κð Þ <sup>r</sup> � <sup>d</sup> <sup>=</sup>r, <sup>κ</sup><sup>2</sup> <sup>≔</sup> 4πe2ρi=εkBT is used instead, where e is the electron charge. Q<sup>∗</sup> is the effective electric charge of a particle that can be measured with electrophoresis experiments.

## 2.3 Theoretical models of phase diagrams of dipolar colloidal fluids: computer simulation and optimization techniques

Finite-temperature calculations of ground state free energies with Yukawa repulsion and dipolar-dipolar interactions of varying strength were performed by Hynninen and Dijkstra [15]. In their calculations, they used canonical Monte Carlos simulations to obtain the free energy minimum as a function of volume fraction and dipolar strength in three-dimensional model systems. They found a phase diagram that contains the same phases as was observed by Yethiraj et al. [11]. However, additionally, they predicted a new hexagonal close-packed phase at the hard sphere repulsion limit and body-centered orthorhombic when the repulsion becomes long range by lowering the content of electrolyte in solution. Their method allows quicker and reliable determinations of the three-dimensional phase diagram than the evolutionary algorithm method that searches for the ground state at zero temperature for two-dimensional systems [45]. The bulk phase diagram has not been searched yet with genetic algorithm (GA) at finite temperature. A particular detail of GA is that to look for the minimum of the free energy, the derivative of the energy function at each evolution step is required for finding the best-adapted crystal structures (individuals) in the population, which is made up of several crystal structures produced by the algorithm. The search of the best-adapted individuals (crystal phase) is done in the energy landscape. For fixed pressure and finite temperature, the Gibbs free energy as a function of pairwise interaction strength and particles concentration in the fluid yields a three-dimensional plot of the free energy known as energy landscape, which has local minima and maxima. The global minimum of an equilibrium crystal structure is reached with the assistance of a local gradient algorithm that requires the derivative of the energy. Such a procedure makes an indirect search of the crystal structure by sampling the energy landscape. Thus, it becomes technically more involved to implement than the Hynninen et al. method of MC ground state energy searches. In Figure 1, we provide two nonequilibrium states of local minima that result from the application of a genetic algorithm approach at zero temperature and fixed pressure.

Recently, Spiteri and Messina [16] have proposed an efficient nonlinear optimization technique that allows predicting the crystal phases of dipolar colloids at zero temperature without the need to use derivatives of the free energies. It was found that for a monodisperse repulsive hard sphere plus dipole-dipole interaction, the predicted phase diagram has the same phases as was previously found by Hynninen et al. [15] and that a knew so-called clinohexagonal prism span all known ground state structures at any density. There remains to be verified if this prediction fulfills at a finite temperature also.

## 2.4 Experimental structure factors and diffusion coefficients

The microstructural order inside of the ferrofluid is determined by the potential interaction among particles. The measured structure factor yields the details of the particles'spatial geometric arrangement. The water-based maghemite ferrofluid is

counterpart and with computer simulations. Another method to validate the proposed model of particle's interaction is based on the parameters that appear in the Lennard-Jones potential that represents the Van der Waals attractions together with the pair dipole-dipole potential, and by fitting the predicted magnetization curve as a function of the applied magnetic field to the experimental one [42]. Even though there is not yet a realistic model interaction potential that represents the interaction of ferromagnetic particles, the use of fitting parameters as mentioned above has led to the prediction of the structure factor and the birefringence as a function of the applied field, which shows qualitative agreement with the experimental values

Pattern Formation and Stability in Magnetic Colloids

2.2 Experimental phase diagram of a charged and sterically stabilized electrorheological colloidal suspension under an electric field

Yethiraj et al. [11] made a colloidal system of charged and sterically stabilized polymethyl methacrylate spheres of micron size in an organic solvent to have dipole moments on particles induced by an electric field. This monodisperse suspension shows both long-range repulsion, and attractive anisotropic interaction potential that is fixed with the addition of salt, whereas the dipolar interaction is controlled by the external electric field. For low Reynolds numbers and in an infinite fluid, a solid spherical particle that sinks in the fluid reaches a terminal velocity due to a balance of the friction force exerted by the fluid that opposes the gravitational force on the particle. If the fluid has a dynamic viscosity η and mass density ρ<sup>f</sup> , and the spherical particle possesses a mass density ρm, then the steady velocity is

diameter. Thus, the effect of gravity is diminished by matching the density of the particles to that of the solvent. The Van der Waals attraction that otherwise would produce aggregation was reduced by matching the index of refraction to the visible. The Van der Waals attractive energy of two spherical particles is Vvdwð Þ¼� r ð Þ A=12

=r<sup>2</sup> � � � � , where r is the separation distance of the

colloidal particles np is equal to that of the solvent ns for an appropriated selection of the frequency ν of the visible light. Here, h is the Planck constant. With confocal microscopy studies, they observed in real-space representation crystalline structures controlled by the repulsive part of the potential, which induces a structure of random hexagonal close packing. For longer-range repulsion, it becomes facecentered cubic (fcc), and one can also find a body-centered cubic phase (bcc). As a function of the electric field, a dipolar phase diagram was established. For low field, there is a sequence of fluid-bcc-fcc (phase-centered cubic) phases. For larger values of the field, there is string formation that coalesces into sheets that in turn coarsen to form body-centered tetragonal crystallites. Moreover, for increasing volume fractions, the string fluid crystallizes to a space-filling tetragonal phase. Such a rich variety of phases can be generated controllably by the concentration of the electrolyte. High salt concentration screens the electrostatic repulsion becoming shortdistance hard sphere interaction whereas the small content of salt increases the spatial range of the potential. The dipolar interaction is controlled by the external electric field. Consequently, due to the symmetric form of the electric and magnetic dipolar potential, a simple change of physical units leads to the similar conclusions to be expected for a colloidal magnetic ferrofluid. For oily ferrofluids, the shortrange Van der Waals attraction and hardcore repulsion are modeled by the

spheres' centers [43]. Therefore, the Hamaker constant [44] A ¼ 3hν np � ns

=18η, where g is the gravitational acceleration and d the particle

<sup>2</sup> � �<sup>3</sup>=<sup>2</sup> can be made zero if the optical refractive index of the

� �<sup>2</sup>

[22, 29–31].

V ¼ ρ<sup>f</sup> � ρ<sup>m</sup> � �

<sup>=</sup> <sup>r</sup><sup>2</sup> � <sup>d</sup><sup>2</sup> � � <sup>þ</sup> <sup>d</sup><sup>2</sup>

=16 ffiffi 2 <sup>p</sup> np

d2

44

np þ ns � �<sup>2</sup> gd<sup>2</sup>

<sup>=</sup>r<sup>2</sup> <sup>þ</sup> 2 ln 1 � <sup>d</sup><sup>2</sup>

<sup>2</sup> <sup>þ</sup> ns

Figure 3 is a plot of the collective diffusion coefficient DC normalized to the free diffusion of a single particle D<sup>0</sup> in the hydrodynamic limit of zero wave number. It is given as Dc <sup>¼</sup> <sup>D</sup>0=S kð Þ <sup>¼</sup> <sup>0</sup> . In the left-hand side plot of Figure 3(a) with symbol <sup>∘</sup> the collective diffusion coefficient increases as the magnitude of dipole moment increases at high-density limit of <sup>ρ</sup><sup>∗</sup> <sup>¼</sup> <sup>0</sup>:9. However, a reverse behavior is shown by this collective diffusion property when the density is very low <sup>ρ</sup><sup>∗</sup> <sup>¼</sup> <sup>ρ</sup>d<sup>3</sup> <sup>¼</sup> <sup>0</sup>:1 (plot with a symbol ●, ρ ¼ N=V, N the number of particles, and V the volume of the fluid). We should note that in the former case, there occurs the formation of chains,

The researchers in Ref. [46] also measured the wave vector-dependent collective translational diffusion coefficient in the absence of an external field using X-ray

Simulation results of structure factor versus wave number. The contact values at k=0 yield the compressibility modulus. The material parameters of maghemite colloid were used [40]. At low densities and magnetic

Calculated collective diffusion coefficient Dc normalized to the free particle diffusion constant at the hydrodynamic limit of low wavenumber k = 0. Notice that at high density, <sup>ρ</sup><sup>∗</sup> <sup>¼</sup> <sup>0</sup>:9, this property increases

moment, there is liquid order. For higher dipole moments formation of chain appears.

and at the highest concentration, the particles are dispersed.

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

DOI: http://dx.doi.org/10.5772/intechopen.81091

Figure 2.

Figure 3.

47

(left, black square symbol).

#### Figure 1.

Genetic algorithm prediction of minimum local structures at two densities of equally sized ferroparticles in the aqueous colloidal magnetic fluid at T ¼ 0K and constant pressure. X, Y, and Z are Cartesian coordinates. Both structures are simple cubic cells. The bottom structure has the highest density and the lowest enthalpy energy. Black arrows depict the direction of the particle's dipole moment.

electrostatically stabilized. The interparticle interaction was modeled in Ref. [13] with a Yukawa repulsion part together with the anisotropic dipolar potential. The highest peak of the structure factor regarding wave number k, as a function of applied field, shows the same qualitative tendency as the experimental peak. For the same type of pair potential, a calculation of the structure factor using the rescaled mean spherical approximation of liquid theory was developed by Wagner et al. [46]. Their calculated structure factor compares well with the measured one in a cobalt-based aqueous colloid. Their fabricated ferrofluid has a polydispersity in particle sizes, which were successfully taken into account via a Shultz distribution. The microstructure in the suspension has also been modeled by Pyanzina et al. [47] with an analytical pair correlation function obtained through a thermodynamic expansion of the density and dipole strength. They find excellent agreement between the theoretical calculation and molecular dynamics data for the scattering vector dependence of the pair correlation at the dilute limit. The agreement improves at low volume fraction whereas at higher concentrations, the coincidence occurs only about the main first peak. Recently, we determined with Langevin dynamics simulations the structure factor in model monodisperse ferrofluids of spherical particles that interact through Lennard-Jones plus dipolar pair potentials. As previously found [17, 42], we also find that at the low colloid density and dipole strength, the particles are well dispersed; however, for higher dipole moment, there is chaining of particles. If the concentration rises, then the system shows a small cluster with few particles as can be seen in Figure 2 of the structure factor S kð Þ.

The contact values at k = 0 yield the compressibility modulus. The material parameters of maghemite colloid were used [40]. At low densities and magnetic moment, there is liquid order. For higher dipole moments, formation of chain appears. Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids DOI: http://dx.doi.org/10.5772/intechopen.81091

Figure 3 is a plot of the collective diffusion coefficient DC normalized to the free diffusion of a single particle D<sup>0</sup> in the hydrodynamic limit of zero wave number. It is given as Dc <sup>¼</sup> <sup>D</sup>0=S kð Þ <sup>¼</sup> <sup>0</sup> . In the left-hand side plot of Figure 3(a) with symbol <sup>∘</sup> the collective diffusion coefficient increases as the magnitude of dipole moment increases at high-density limit of <sup>ρ</sup><sup>∗</sup> <sup>¼</sup> <sup>0</sup>:9. However, a reverse behavior is shown by this collective diffusion property when the density is very low <sup>ρ</sup><sup>∗</sup> <sup>¼</sup> <sup>ρ</sup>d<sup>3</sup> <sup>¼</sup> <sup>0</sup>:1 (plot with a symbol ●, ρ ¼ N=V, N the number of particles, and V the volume of the fluid). We should note that in the former case, there occurs the formation of chains, and at the highest concentration, the particles are dispersed.

The researchers in Ref. [46] also measured the wave vector-dependent collective translational diffusion coefficient in the absence of an external field using X-ray

#### Figure 2.

electrostatically stabilized. The interparticle interaction was modeled in Ref. [13] with a Yukawa repulsion part together with the anisotropic dipolar potential. The highest peak of the structure factor regarding wave number k, as a function of applied field, shows the same qualitative tendency as the experimental peak. For the same type of pair potential, a calculation of the structure factor using the rescaled mean spherical approximation of liquid theory was developed by Wagner et al. [46]. Their calculated structure factor compares well with the measured one in a cobalt-based aqueous colloid. Their fabricated ferrofluid has a polydispersity in particle sizes, which were successfully taken into account via a Shultz distribution. The microstructure in the suspension has also been modeled by Pyanzina et al. [47] with an analytical pair correlation function obtained through a thermodynamic expansion of the density and dipole strength. They find excellent agreement between the theoretical calculation and molecular dynamics data for the scattering vector dependence of the pair correlation at the dilute limit. The agreement improves at low volume fraction whereas at higher concentrations, the coincidence occurs only about the main first peak.

Genetic algorithm prediction of minimum local structures at two densities of equally sized ferroparticles in the aqueous colloidal magnetic fluid at T ¼ 0K and constant pressure. X, Y, and Z are Cartesian coordinates. Both structures are simple cubic cells. The bottom structure has the highest density and the lowest enthalpy energy.

Black arrows depict the direction of the particle's dipole moment.

Pattern Formation and Stability in Magnetic Colloids

Recently, we determined with Langevin dynamics simulations the structure factor in model monodisperse ferrofluids of spherical particles that interact through Lennard-Jones plus dipolar pair potentials. As previously found [17, 42], we also find that at the low colloid density and dipole strength, the particles are well dispersed; however, for higher dipole moment, there is chaining of particles. If the concentration rises, then the system shows a small cluster with few particles as can be seen in Figure 2 of

The contact values at k = 0 yield the compressibility modulus. The material parameters of maghemite colloid were used [40]. At low densities and magnetic moment, there is liquid order. For higher dipole moments, formation of chain appears.

the structure factor S kð Þ.

46

Figure 1.

Simulation results of structure factor versus wave number. The contact values at k=0 yield the compressibility modulus. The material parameters of maghemite colloid were used [40]. At low densities and magnetic moment, there is liquid order. For higher dipole moments formation of chain appears.

#### Figure 3.

Calculated collective diffusion coefficient Dc normalized to the free particle diffusion constant at the hydrodynamic limit of low wavenumber k = 0. Notice that at high density, <sup>ρ</sup><sup>∗</sup> <sup>¼</sup> <sup>0</sup>:9, this property increases (left, black square symbol).

correlation spectroscopy. Such a dynamical property has a known linear dependence on the particles' hydrodynamic interactions (HIs). The proposed theoretical hydrodynamic function that captures the HI among particles overestimates the experimental data. On the other hand, Meriguet et al. [29] used small-angle neutron scattering to measure the incoherent scattering function at the fixed magnetic field. The incoherent part gives in the long-time overdamped regime, the single translational diffusion coefficient for long wave vectors k, and the coherent part provides the collective diffusion as a function of k. Using liquid theory definitions of these two dynamical functions, their comparison with the experimental values yields good agreement at all wave scattering vectors k. The same collective diffusion coefficient can be measured by forced Rayleigh scattering also [26]. In this technique, the fluctuations in concentrations of the particles are considered a dynamic variable. We note that similar studies have been undertaken in the case of rodlike macromolecules in suspension such as Tobacco mosaic virus [48]. Explicit statistical mechanics derivations for the dynamical collective correlation functions (intermediate scattering function) have been given that fit well Brownian dynamics simulations just at long wave numbers but fail at intermediate and short k values. In Refs. [32–34], we proposed a Langevin equation theory for a tracer ferroparticle whose motion is coupled to the cloud of the other colloidal magnetic particles with which it interacts.

$$\begin{split} M\frac{d\mathcal{V}(t)}{dt} &= -\zeta^{0} \cdot \mathcal{V}(t) - \zeta^{0}\_{TR} \cdot \mathcal{W}(t) + f^{0}(t) + F(t), \\ I \cdot \frac{d\mathcal{W}(t)}{dt} &= -\zeta^{0}\_{R} \cdot \mathcal{W}(t) - \zeta^{0}\_{RT} \cdot \mathcal{V}(t) + t^{0}(t) + T(t). \end{split} \tag{1}$$

equation of the concentration of particles around the tracer through the total force on

ψ r !; Ω h i � � n r!; <sup>Ω</sup>; <sup>t</sup>

� � around the local equilibrium concentration of particles

D E � � . Thus, in the center of mass of the tracer, and using general

<sup>0</sup> � � <sup>þ</sup> <sup>∇</sup>

� j r!; Ω; t � �� <sup>∇</sup>

� �: and the diffusion equation of the con-

n tð Þ h i can be linearized altogether with Eq. (3)

dr00dΩ<sup>00</sup>L r!; r

\$

\$ � j † r !; Ω; t D E � � is related to

> t � t <sup>0</sup> ð Þ� V \$ ð Þþt F \$

\$ ð Þt F \$†ð Þ <sup>0</sup> D E <sup>¼</sup> kBT<sup>Δ</sup> <sup>ζ</sup>

� � that embodies the microstructure of the particles inside of

!0 ; Ω; Ω<sup>0</sup>

� j r!; Ω; t � �:

� �. † means transpose and

� �, and the result is

<sup>0</sup> � �

; t � t

ð Þt (5)

\$ ð Þt .

:

(6)

(4)

� � (3)

the tracer, which is linear in the instantaneous concentration of particles.

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

! �∇ þ ∇<sup>Ω</sup>

\$

results of linear irreversible thermodynamics, the diffusion equation reads

; <sup>Ω</sup><sup>00</sup> � �δn r!00; <sup>Ω</sup>00; <sup>t</sup>

\$

complex conjugate (adjoint operator). The inverse of the static correlation function σ fulfills a relationship with the potential and concentration function [32–34]. Using

further substituted into the linearized version in the concentration of Eq. (3) yield-

are fluctuating forces arising from the spontaneous departure from zero of the

dΩ<sup>00</sup> ∇ \$ neq r !; Ω h i � �

; Ω00; t � � <sup>∇</sup>

> dΩ<sup>00</sup> ∇ \$ ψ r !; Ω h i � �

\$00neq r !00; <sup>Ω</sup><sup>00</sup> h i � � †

\$00neq r !00; <sup>Ω</sup><sup>00</sup> h i � � †

net force and torque on the tracer. It satisfies the fluctuation-dissipation relation-

T ¼ 300K is the room temperature. As a result, explicit expressions for the tracer translational and rotational friction coefficients are given regarding the equilibrium

!00dΩdΩ<sup>0</sup>

!0 ; r !00; Ω<sup>0</sup>

!0 ; r !00; Ω<sup>0</sup> ; <sup>Ω</sup><sup>00</sup> � � <sup>∇</sup>

!00dΩdΩ<sup>0</sup>

\$ � ðt 0 dt<sup>0</sup> ð d<sup>3</sup> r !0 dΩ<sup>0</sup> ð

\$0� <sup>V</sup> \$ ð Þþt f \$ ðÞþt ð d<sup>3</sup> r !dΩ∇ \$

\$

!; Ω

\$ neq r !; Ω h i � � � <sup>V</sup>

�σ�<sup>1</sup> <sup>r</sup> !0 ; r !00; Ω<sup>0</sup>

where the Onsager coefficient L ¼ ∇

the spatial variation of the stochastic current j r!; Ω; t

methods of Green functions Eq. (4), one solves for δn r!; Ω; t

\$<sup>0</sup>� <sup>V</sup> \$

ship with the tracer-effective friction coefficient F

!; Ω

the liquid, and the direct pair potential ψ,

ð d<sup>3</sup> r !d<sup>3</sup> r !0 d<sup>3</sup> r

!; r !0 ; Ω; Ω<sup>0</sup> � �<sup>χ</sup> <sup>r</sup>

> ð d<sup>3</sup> r !d<sup>3</sup> r !0 d<sup>3</sup> r

kBT ð Þ <sup>4</sup><sup>π</sup> <sup>2</sup>

� <sup>σ</sup>�<sup>1</sup> <sup>r</sup>

� χ r !; r !0 ; Ω; Ω<sup>0</sup> ; t � �<sup>σ</sup> <sup>r</sup>

1 kBTð Þ <sup>4</sup><sup>π</sup> <sup>2</sup> ð Þþt f \$ ðÞ�t ðt 0 dt0 Δ ζ \$

� �.∂<sup>t</sup> ¼ �<sup>∇</sup> � <sup>V</sup>

¼ ∇; r

¼ ð Þ V;W , ∇

DOI: http://dx.doi.org/10.5772/intechopen.81091

M \$ � dV \$ ð Þt dt ¼ �<sup>ζ</sup>

Here, V \$

neq r !; Ω

centration ∂n r!; Ω;t

in <sup>δ</sup>n tðÞ¼ n tð Þ‐neq <sup>r</sup>

∂δn r!; Ω; t � �

� � <sup>≔</sup> n r!; <sup>Ω</sup>; <sup>t</sup>

<sup>∂</sup><sup>t</sup> <sup>¼</sup> <sup>∇</sup>

ing the tracer Langevin equation

M \$ � dV \$ ð Þt dt ¼ �<sup>ζ</sup>

concentration neq r

Δζ \$ ðÞ¼ t

Δζ \$ ðÞ¼ t

49

F \$

V,W are the translational and angular velocities of the tracer particle and referred to a space-fixed frame at the center of mass of the tracer and following the orientation of the main tracer axis of symmetry. M and I are the mass and particle's matrix of the moment of inertia, respectively. There are neither hydrodynamic interactions among particles nor external magnetic fields. The first two terms in Eq. (1) are the solvent friction force and torque. The short-time free particle diagonal friction tensors ζ<sup>0</sup>, ζ<sup>0</sup> TR, ζ<sup>0</sup> RT, ζ<sup>0</sup> <sup>R</sup> represent hydrodynamic drag forces and torques. These friction and random forces are the only quantities that convey information on the nature of the solvent. They need to be measured experimentally or provided by a free parameter. Moreover, they ignore the molecular degrees of freedom of position and orientation coordinates and momenta. They are coupled to the thermally driven solvent random forces f <sup>0</sup> and torques t <sup>0</sup> by fluctuationdissipation theorems [32–34]. The total force and torque F, T on the tracer by the rest of the colloidal particles, which are at the concentration n r!; Ω; t � � <sup>¼</sup>

$$\sum\_{i=1}^{N} \delta \left( \overrightarrow{r} - \overrightarrow{r}\_i(t) \right) \delta(\Omega - \Omega\_i(t)), \text{ are given by } $$

$$\begin{aligned} F(t) &= \int d^3 \vec{r} \, d\Omega \nabla \varphi \left( \vec{r}, \Omega \right) n \left( \vec{r}, \Omega, t \right), \\ T(t) &= \int d^3 \vec{r} \, d\Omega \left[ \vec{r} \times \nabla + \nabla\_{\Omega} \right] \varphi \left( \vec{r}, \Omega \right) n \left( \vec{r}, \Omega, t \right) \end{aligned} \tag{2}$$

with Ω ¼ ð Þ θ; φ being the polar angles, ψ r !; Ω � � <sup>¼</sup> VLJ <sup>þ</sup> Vddð Þ <sup>μ</sup>0=4<sup>π</sup> is the pair potential with μ<sup>0</sup> the magnetic permeability of vacuum, and ∇<sup>Ω</sup> ¼ u^ � d=du^ the angular gradient operator. The total force and torque are calculated using the probe's dipole located along the Z axis direction of its local frame, that is Ω ¼ ð Þ 0; 0 . The unitary Cartesian vector u^ of the orientation of any other particle is in the direction of that particle's axis of symmetry. The Langevin equation is coupled to the diffusion

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids DOI: http://dx.doi.org/10.5772/intechopen.81091

equation of the concentration of particles around the tracer through the total force on the tracer, which is linear in the instantaneous concentration of particles.

$$
\overrightarrow{M} \cdot \frac{d\overrightarrow{V}(t)}{dt} = -\overrightarrow{\zeta}^{0} \cdot \overrightarrow{V}(t) + \overrightarrow{f}(t) + \int d^{3} \overrightarrow{r} d\Omega \overrightarrow{\nabla} \left[ \psi \left( \overrightarrow{r}, \Omega \right) \right] n \left( \overrightarrow{r}, \Omega, t \right) \tag{3}
$$

Here, V \$ ¼ ð Þ V;W , ∇ \$ ¼ ∇; r ! �∇ þ ∇<sup>Ω</sup> � �: and the diffusion equation of the concentration ∂n r!; Ω;t � �.∂<sup>t</sup> ¼ �<sup>∇</sup> � <sup>V</sup> \$ n tð Þ h i can be linearized altogether with Eq. (3) in <sup>δ</sup>n tðÞ¼ n tð Þ‐neq <sup>r</sup> !; Ω � � around the local equilibrium concentration of particles neq r !; Ω � � <sup>≔</sup> n r!; <sup>Ω</sup>; <sup>t</sup> D E � � . Thus, in the center of mass of the tracer, and using general results of linear irreversible thermodynamics, the diffusion equation reads

$$
\begin{split}
\frac{\partial \delta n\left(\overrightarrow{r},\Omega,t\right)}{\partial t} &= \left[\overrightarrow{\nabla}n^{eq}\left(\overrightarrow{r},\Omega\right)\right] \cdot \overrightarrow{V} - \int\_{0}^{t} dt' \int d^{3}\overrightarrow{r}' d\Omega' \int dr' d\Omega'' L\left(\overrightarrow{r},\overrightarrow{r}',\Omega,\Omega';t-t'\right) \\ &\times \sigma^{-1}\left(\overrightarrow{r}',\overrightarrow{r}'',\Omega',\Omega''\right) \delta n\left(\overrightarrow{r}'',\Omega'',t'\right) + \overrightarrow{\nabla} \cdot j\left(\overrightarrow{r},\Omega,t\right).
\end{split} \tag{4}$$

where the Onsager coefficient L ¼ ∇ \$ � j r!; Ω; t � �� <sup>∇</sup> \$ � j † r !; Ω; t D E � � is related to the spatial variation of the stochastic current j r!; Ω; t � �. † means transpose and complex conjugate (adjoint operator). The inverse of the static correlation function σ fulfills a relationship with the potential and concentration function [32–34]. Using methods of Green functions Eq. (4), one solves for δn r!; Ω; t � �, and the result is further substituted into the linearized version in the concentration of Eq. (3) yielding the tracer Langevin equation

$$
\overrightarrow{\dot{M}} \cdot \frac{d\overrightarrow{V}(t)}{dt} = -\overrightarrow{\zeta}^{0} \cdot \overrightarrow{V}\ (t) + \overrightarrow{f}\ (t) - \int\_{0}^{t} dt' \Delta \ \overrightarrow{\zeta}\ (t - t') \cdot \overrightarrow{V}(t) + \overrightarrow{F}(t) \tag{5}
$$

F \$ are fluctuating forces arising from the spontaneous departure from zero of the net force and torque on the tracer. It satisfies the fluctuation-dissipation relationship with the tracer-effective friction coefficient F \$ ð Þt F \$†ð Þ <sup>0</sup> D E <sup>¼</sup> kBT<sup>Δ</sup> <sup>ζ</sup> \$ ð Þt . T ¼ 300K is the room temperature. As a result, explicit expressions for the tracer translational and rotational friction coefficients are given regarding the equilibrium concentration neq r !; Ω � � that embodies the microstructure of the particles inside of the liquid, and the direct pair potential ψ,

$$\begin{split} \Delta \overrightarrow{\zeta}(t) &= \frac{k\_B T}{\left(4\pi\right)^2} \int d^3 \vec{r} d^3 \vec{r}' d^3 \vec{r}'' d\Omega d\Omega' d\Omega'' \left[\overleftrightarrow{\nabla} n'^q \left(\overrightarrow{r}, \Omega\right)\right] \\ &\times \sigma^{-1} \left(\overrightarrow{r}, \overrightarrow{r}', \Omega, \Delta'\right) \chi \left(\overrightarrow{r}', \overrightarrow{r}'', \Omega', \Omega'', t\right) \left[\overleftrightarrow{\nabla}'' n'^q \left(\overrightarrow{r}'', \Omega''\right)\right]^\dagger. \\ \Delta \overleftrightarrow{\zeta}'(t) &= \frac{1}{k\_B T \left(4\pi\right)^2} \int d^3 \vec{r} d^3 \vec{r}' d^3 \vec{r}'' \eta d\Omega d\Omega' d\Omega'' \left[\overleftrightarrow{\nabla}\nu \left(\overrightarrow{r}, \Omega\right)\right] \\ &\times \chi \left(\overrightarrow{r}', \overrightarrow{r}', \Omega, \Delta', t\right) \sigma \left(\overrightarrow{r}', \overrightarrow{r}'', \Delta', \Omega'\right) \left[\overleftrightarrow{\nabla}'' n'^q \left(\overrightarrow{r}'', \Omega''\right)\right]^\dagger \end{split} \tag{6}$$

correlation spectroscopy. Such a dynamical property has a known linear dependence on the particles' hydrodynamic interactions (HIs). The proposed theoretical hydrodynamic function that captures the HI among particles overestimates the experimental data. On the other hand, Meriguet et al. [29] used small-angle neutron scattering to measure the incoherent scattering function at the fixed magnetic field. The incoherent part gives in the long-time overdamped regime, the single translational diffusion coefficient for long wave vectors k, and the coherent part provides the collective diffusion as a function of k. Using liquid theory definitions of these two dynamical functions, their comparison with the experimental values yields good agreement at all wave scattering vectors k. The same collective diffusion coefficient can be measured by forced Rayleigh scattering also [26]. In this technique, the fluctuations in concentrations of the particles are considered a dynamic variable. We note that similar studies have been undertaken in the case of rodlike macromolecules in suspension such as Tobacco mosaic virus [48]. Explicit statistical mechanics derivations for the dynamical collective correlation functions (intermediate scattering function) have been given that fit well Brownian dynamics simulations just at long wave numbers but fail at intermediate and short k values. In Refs. [32–34], we proposed a Langevin equation theory for a tracer ferroparticle whose motion is coupled to the cloud of the other colloidal magnetic particles with which it interacts.

<sup>M</sup> dV tð Þ

dW tð Þ

Pattern Formation and Stability in Magnetic Colloids

the thermally driven solvent random forces f

F tðÞ¼ <sup>Ð</sup>

T tðÞ¼ <sup>Ð</sup>

dt ¼ �ζ<sup>0</sup>

TR, ζ<sup>0</sup> RT, ζ<sup>0</sup>

I �

onal friction tensors ζ<sup>0</sup>, ζ<sup>0</sup>

∑<sup>N</sup> <sup>i</sup>¼<sup>1</sup><sup>δ</sup> <sup>r</sup> ! �r ! <sup>i</sup>ð Þt � �

48

dt ¼ �ζ<sup>0</sup> � V tðÞ� <sup>ζ</sup><sup>0</sup>

<sup>R</sup> � W tðÞ� <sup>ζ</sup><sup>0</sup>

V,W are the translational and angular velocities of the tracer particle and referred to a space-fixed frame at the center of mass of the tracer and following the orientation of the main tracer axis of symmetry. M and I are the mass and particle's matrix of the moment of inertia, respectively. There are neither hydrodynamic interactions among particles nor external magnetic fields. The first two terms in Eq. (1) are the solvent friction force and torque. The short-time free particle diag-

torques. These friction and random forces are the only quantities that convey information on the nature of the solvent. They need to be measured experimentally or provided by a free parameter. Moreover, they ignore the molecular degrees of freedom of position and orientation coordinates and momenta. They are coupled to

dissipation theorems [32–34]. The total force and torque F, T on the tracer by the

!; Ω � �

! �∇ þ ∇<sup>Ω</sup> h i

potential with μ<sup>0</sup> the magnetic permeability of vacuum, and ∇<sup>Ω</sup> ¼ u^ � d=du^ the angular gradient operator. The total force and torque are calculated using the probe's dipole located along the Z axis direction of its local frame, that is Ω ¼ ð Þ 0; 0 . The unitary Cartesian vector u^ of the orientation of any other particle is in the direction of that particle's axis of symmetry. The Langevin equation is coupled to the diffusion

rest of the colloidal particles, which are at the concentration n r!; Ω; t

! dΩ∇ψ r

δð Þ Ω � Ωið Þt , are given by

d<sup>3</sup> r

d<sup>3</sup> r ! dΩ r

with Ω ¼ ð Þ θ; φ being the polar angles, ψ r

TR � W tðÞþ f

RT � V tðÞþ t

<sup>R</sup> represent hydrodynamic drag forces and

<sup>0</sup> and torques t

,

n r!; Ω; t

n r!; Ω; t � �

!; Ω � �

ψ r !; Ω � �

<sup>0</sup>ðÞþ<sup>t</sup> F tð Þ,

<sup>0</sup>ð Þþ<sup>t</sup> T tð Þ:

<sup>0</sup> by fluctuation-

� � (2)

¼ VLJ þ Vddð Þ μ0=4π is the pair

¼

� �

(1)

The propagator χ is obtained from the diffusion Eq. (4), and it is related to the van Hove function C r!; Ω; r !0 ; Ω<sup>0</sup> ; t � � <sup>≔</sup> <sup>δ</sup>n r!; <sup>Ω</sup>; <sup>t</sup> � �δn r<sup>0</sup> Ω0 0 D E � � <sup>¼</sup> <sup>Ð</sup> d<sup>3</sup> r !00dΩ<sup>00</sup> χ r !; r !00; Ω; Ω00; t � �<sup>σ</sup> <sup>r</sup> !00; r !0 ; Ω00; Ω<sup>0</sup> � � whose dynamical evolution is attained from Eq. (4) after multiplying both sides by <sup>δ</sup>n r!; <sup>Ω</sup>; <sup>t</sup> <sup>¼</sup> <sup>0</sup> � � and taking further an ensemble average. The resulting dynamical equation of the van Hove function has at t = 0 the initial condition Ð d<sup>3</sup> r !00dΩ00σ r !; r !00; <sup>Ω</sup>; <sup>Ω</sup><sup>00</sup> � �<sup>χ</sup> <sup>r</sup> !00; r !0 ; Ω00; Ω<sup>0</sup> ; t ¼ 0 � � <sup>¼</sup> σ r !; r !0 ; Ω; Ω<sup>0</sup> � �. The properties in Eq. (6) depend on the knowledge of the structure factor of the ferrofluid, which takes into account the micro-structural order of particles inside the ferrofluid as dictated by the underlying interaction potential. It is valid at high concentrations but ignores hydrodynamic interactions among particles. For monodisperse colloids, the tracer is equal to the other particles and the fluid is homogeneous. Thus, the distribution of particles is <sup>n</sup>eqð Þ¼ <sup>r</sup>; <sup>Ω</sup> ρg r ¼ r ! � � � � � �; <sup>Ω</sup>; <sup>Ω</sup><sup>0</sup> ð Þ 0; 0 � �. Notice that the pair correlation function is defined between the tracer orientation Ω<sup>0</sup> ð Þ 0; 0 , and another particle with Ωð Þ θ; φ , separated by the mean distance r ¼ r ! �<sup>r</sup> !0 � � � � � �. Due to the spatial symmetry of the potential <sup>ψ</sup>, the pair correlation function has the rotationally invariant expansion g r; Ω; Ω<sup>0</sup> ð Þ¼ g rð Þþ <sup>h</sup>Δð Þ<sup>r</sup> <sup>Δ</sup> <sup>þ</sup> hDð Þ<sup>r</sup> <sup>D</sup>1, with <sup>Δ</sup> <sup>≔</sup> <sup>u</sup>^<sup>1</sup> � <sup>u</sup>^2, D<sup>1</sup> <sup>≔</sup> <sup>3</sup>^ r ! <sup>12</sup> � u^<sup>1</sup> ^ r ! <sup>12</sup> � u^<sup>2</sup> � u^<sup>2</sup> � u^<sup>1</sup> h i, where ^ r ! ¼r !=r, u^ is the unitary dipole vector with orientation u^ ¼ u^ð Þ θ; φ . Using the Fick approximation for the Onsager coefficient, L r;r<sup>0</sup> ; <sup>Ω</sup>; <sup>Ω</sup><sup>0</sup> ð Þ¼ <sup>ρ</sup> <sup>D</sup><sup>0</sup>∇<sup>2</sup> <sup>þ</sup> <sup>D</sup><sup>0</sup> <sup>R</sup>∇<sup>Ω</sup> � � <sup>δ</sup> <sup>r</sup> � <sup>r</sup><sup>0</sup> ð Þ<sup>δ</sup> <sup>Ω</sup> � <sup>Ω</sup><sup>0</sup> ð Þ, the unknown propagator in Eq. (4) <sup>χ</sup> <sup>¼</sup> exp �tL∘σ�<sup>1</sup> ð Þ,∘ ≔ <sup>Ð</sup> d<sup>3</sup> r ! dΩ can be fully determined. For a spherical particle, the free translational and rotational diffusion coefficients are <sup>D</sup><sup>0</sup> <sup>¼</sup> kBT=ζ<sup>0</sup>, D<sup>0</sup> <sup>R</sup> <sup>¼</sup> kBT=ζ<sup>0</sup> <sup>R</sup>, with <sup>ζ</sup><sup>0</sup> <sup>¼</sup> <sup>3</sup>πηsold, ζ0 <sup>R</sup> <sup>¼</sup> πηsold<sup>3</sup> , and ηsol the solvent viscosity. Taking the Laplace transform Δζ ωð Þ¼ Ð <sup>∞</sup> <sup>0</sup> dte<sup>i</sup>ω<sup>t</sup> Δζð Þt of the friction tensor in Eq. (6), it yields the real part of the translational and rotational diffusion coefficients of the tracer at the longtime overdamped regime <sup>D</sup> <sup>¼</sup> kBT<sup>=</sup> <sup>ζ</sup><sup>0</sup> <sup>þ</sup> ReΔζ ωð Þ <sup>¼</sup> <sup>0</sup> � �, DR <sup>¼</sup> kBT<sup>=</sup> <sup>ζ</sup><sup>0</sup> <sup>R</sup> <sup>þ</sup> ReΔζRð Þ <sup>ω</sup> <sup>¼</sup> <sup>0</sup> � � [34], where

These correlation function components are obtained from the Fourier transform

<sup>∑</sup>i,j ið Þ 6¼<sup>j</sup> <sup>δ</sup> <sup>r</sup>� r

<sup>∑</sup>i,j ið Þ 6¼<sup>j</sup> <sup>δ</sup> <sup>r</sup> � <sup>r</sup>

In Figure 4, Ref. [34] is provided a plot of the translational and rotational selfdiffusion coefficients of such a theory. In this picture, we compare the diffusion coefficients with Langevin dynamics simulation results for the same diffusion coef-

made of Fe2O3. μ<sup>0</sup> is the magnetic permeability of vacuum. The material data are:

2.5 Diffusion microrheology of colloidal magnetic fluids under no external

strain rates γ\_ð Þt is represented by a constitutive equation for the shear stress

Ferrofluids are complex fluids whose viscoelastic response to weekly applied

is G tð Þ. In the representation of frequency space ω, the shear is σ ωð Þ¼ Gð Þ ω iωγ ωð Þ. Mason and Weitz [49] assumed that the macroscopic temporal relaxation of G tð Þ is the same time scale of the viscosity that affects a tracer particle that diffuses in the viscoelastic fluid. Because the shear stress divided by the strain rate has dimensions of viscosity, Mason et al. [49] defined the viscosity as ηðÞ¼ t σð Þt =γ\_ð Þt . Thus, they proposed that for spherical particles undergoing translational diffusion and in the frequency space, it is Gð Þ¼ ω iωη ωð Þ¼ iωζ ωð Þ=3πd for a particle experiencing

Kg=ms, respectively.

<sup>0</sup> ð Þγ\_ð Þt where the dot means time derivative and the shear modulus

R

strength (b) ((d)). Simulation calculations are depicted with a black symbol and theory of Eq. (6) with ∘.

� � versus density (a) ((c)), and dipole

[40] <sup>ε</sup><sup>0</sup> <sup>¼</sup> <sup>5</sup>:<sup>45311</sup> � <sup>10</sup>�23J, d <sup>¼</sup> <sup>10</sup>�8m, particle mass and the viscosity are

!ij � � � � � � � � <sup>3</sup>

! ij � � � � � � � �u^<sup>i</sup> � <sup>u</sup>^<sup>j</sup> D E

^ r !ij � u^<sup>i</sup> ^ r

D E h i

!ij � u^j�u^<sup>i</sup> � u^<sup>j</sup>

kBT � � for the real ferrofluid

<sup>N</sup>4πρ∗r<sup>2</sup> ,

<sup>N</sup>4πρ<sup>∗</sup>r<sup>2</sup> (9)

2

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

of their definitions

<sup>∑</sup>i,j ið Þ 6¼<sup>j</sup> <sup>δ</sup> <sup>r</sup>� r

D E � �

!ij � � � � � �

DOI: http://dx.doi.org/10.5772/intechopen.81091

<sup>m</sup><sup>0</sup> <sup>¼</sup> <sup>2</sup>:70710Kg, <sup>η</sup>sol <sup>¼</sup> <sup>0</sup>:<sup>852</sup> � <sup>10</sup>�<sup>3</sup>

electric or magnetic fields

G t � t

<sup>N</sup>4πρ∗r<sup>2</sup> , hDð Þ¼ <sup>r</sup> <sup>3</sup>

hΔð Þ¼ r 3

ficients by using typical parameters of <sup>μ</sup><sup>∗</sup><sup>2</sup> <sup>¼</sup> <sup>μ</sup>0m2<sup>=</sup> <sup>4</sup>πd<sup>3</sup>

g rð Þ¼

<sup>σ</sup>ðÞ¼ <sup>t</sup> <sup>Ð</sup><sup>t</sup>

Figure 4.

51

Translational (rotational) self-diffusion coefficients D=D<sup>0</sup> DR=D<sup>0</sup>

<sup>0</sup> dt<sup>0</sup>

$$\mathrm{Re}\Delta\zeta = \frac{\zeta^0}{6\pi^2\rho^\*} \int\_0^\infty d\mathbf{x} \mathbf{x}^4 \frac{\left(\mathcal{S}\_{,0}^{00}(\mathbf{x}) - \mathbf{1}\right)^2 \tau\_{s,i\alpha}^\*}{\mathcal{S}\_{,0}^{00}(\mathbf{x})}$$

$$+ \frac{128\pi\zeta^0\rho^\*\mu^{\*4}}{3} \int\_0^\infty d\mathbf{x} j\_2(\mathbf{x})^2 \left(\mathbf{6} + \mathbf{x}^2\right) \left[2\tau\_{s,1}^{\*2} + 3\tau\_{s,0}^{\*2}\right],\tag{7}$$

$$\mathrm{Re}\Delta\zeta\_R = \frac{128\pi\rho^\*\mu^{\*4}\zeta\_R^0}{105} \int\_0^\infty d\mathbf{x} j\_1(\mathbf{x})^2 \left(\mathbf{6} + \mathbf{x}^2\right) \left[67\tau\_{s,1}^{\*2} + 38\tau\_{s,0}^{\*2}\right],$$

j <sup>l</sup> is the spherical Bessel function of order l, τ<sup>∗</sup> s,iso <sup>¼</sup> <sup>1</sup><sup>=</sup> <sup>x</sup><sup>2</sup> <sup>1</sup> <sup>þ</sup> <sup>1</sup>=S<sup>00</sup> ,<sup>0</sup> ð Þ <sup>x</sup> � � � � , τ∗ s,<sup>α</sup> <sup>¼</sup> <sup>S</sup><sup>11</sup> ,αð Þ <sup>x</sup> <sup>=</sup> <sup>8</sup><sup>π</sup> <sup>x</sup> ½ � ð Þ <sup>2</sup> <sup>þ</sup> <sup>6</sup> , x <sup>≔</sup> kd. The structure factor components of the dipolar liquid are

$$S\_{,0}^{00}(\mathbf{x}) = \mathbf{1} + \frac{\rho^\*}{d^3}[\mathbf{g}(\mathbf{x}) - \mathbf{1}],\\ S\_{,0}^{11}(\mathbf{x}) = \mathbf{1} + \frac{\rho^\*}{\mathbf{3}d^3}[h\_\Delta(\mathbf{x}) + 2h\_D(\mathbf{x})],$$

$$S\_{,1}^{11}(\mathbf{x}) = \mathbf{1} + \frac{\rho^\*}{\mathbf{3}d^3}[h\_\Delta(\mathbf{x}) - h\_D(\mathbf{x})].\tag{8}$$

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids DOI: http://dx.doi.org/10.5772/intechopen.81091

These correlation function components are obtained from the Fourier transform of their definitions

$$\mathbf{g}(r) = \frac{\left\langle \sum\_{\hat{\nu}, \hat{\nu} \neq \hat{\nu}} \delta\left(r - \left| \stackrel{\scriptstyle \to}{\left| r \right|} \right|\right) \right\rangle}{N 4 \pi \rho^\* r^t},\\h\_D(r) = \frac{\left\langle \sum\_{\hat{\nu}, \hat{\nu} \neq \hat{\nu}} \delta\left(r - \left| \stackrel{\scriptstyle \to}{\left| r \right|} \right|\right) \left[ \stackrel{\scriptstyle \to}{\left| r \right|} \hat{\nu}, \hat{\stackrel{\scriptstyle \to}{\left| r \right|} \hat{\nu}} \stackrel{\scriptstyle \to}{\left| \hat{u}\_i \right|} \hat{u}\_j \cdot \hat{u}\_i \right]}\right\rangle},$$

$$h\_\Delta(r) = 3 \frac{\left\langle \sum\_{\hat{u}\_i \, j \left(i \neq j\right)} \delta\left(r - \left| \stackrel{\scriptstyle \to}{\left| r \right|} \right|\right) \hat{u}\_i \cdot \hat{u}\_j \right\rangle}{N 4 \pi \rho^\* r^2} \tag{9}$$

In Figure 4, Ref. [34] is provided a plot of the translational and rotational selfdiffusion coefficients of such a theory. In this picture, we compare the diffusion coefficients with Langevin dynamics simulation results for the same diffusion coefficients by using typical parameters of <sup>μ</sup><sup>∗</sup><sup>2</sup> <sup>¼</sup> <sup>μ</sup>0m2<sup>=</sup> <sup>4</sup>πd<sup>3</sup> kBT � � for the real ferrofluid made of Fe2O3. μ<sup>0</sup> is the magnetic permeability of vacuum. The material data are: [40] <sup>ε</sup><sup>0</sup> <sup>¼</sup> <sup>5</sup>:<sup>45311</sup> � <sup>10</sup>�23J, d <sup>¼</sup> <sup>10</sup>�8m, particle mass and the viscosity are <sup>m</sup><sup>0</sup> <sup>¼</sup> <sup>2</sup>:70710Kg, <sup>η</sup>sol <sup>¼</sup> <sup>0</sup>:<sup>852</sup> � <sup>10</sup>�<sup>3</sup> Kg=ms, respectively.

## 2.5 Diffusion microrheology of colloidal magnetic fluids under no external electric or magnetic fields

Ferrofluids are complex fluids whose viscoelastic response to weekly applied strain rates γ\_ð Þt is represented by a constitutive equation for the shear stress <sup>σ</sup>ðÞ¼ <sup>t</sup> <sup>Ð</sup><sup>t</sup> <sup>0</sup> dt<sup>0</sup> G t � t <sup>0</sup> ð Þγ\_ð Þt where the dot means time derivative and the shear modulus is G tð Þ. In the representation of frequency space ω, the shear is σ ωð Þ¼ Gð Þ ω iωγ ωð Þ. Mason and Weitz [49] assumed that the macroscopic temporal relaxation of G tð Þ is the same time scale of the viscosity that affects a tracer particle that diffuses in the viscoelastic fluid. Because the shear stress divided by the strain rate has dimensions of viscosity, Mason et al. [49] defined the viscosity as ηðÞ¼ t σð Þt =γ\_ð Þt . Thus, they proposed that for spherical particles undergoing translational diffusion and in the frequency space, it is Gð Þ¼ ω iωη ωð Þ¼ iωζ ωð Þ=3πd for a particle experiencing

Figure 4.

Translational (rotational) self-diffusion coefficients D=D<sup>0</sup> DR=D<sup>0</sup> R � � versus density (a) ((c)), and dipole strength (b) ((d)). Simulation calculations are depicted with a black symbol and theory of Eq. (6) with ∘.

The propagator χ is obtained from the diffusion Eq. (4), and it is related to the

≔ δn r!; Ω; t � �

ensemble average. The resulting dynamical equation of the van Hove function has

!; r

!00dΩ00σ r

fluid is homogeneous. Thus, the distribution of particles is <sup>n</sup>eqð Þ¼ <sup>r</sup>; <sup>Ω</sup>

� �

! �<sup>r</sup> !0

� � �

g rð Þþ <sup>h</sup>Δð Þ<sup>r</sup> <sup>Δ</sup> <sup>þ</sup> hDð Þ<sup>r</sup> <sup>D</sup>1, with <sup>Δ</sup> <sup>≔</sup> <sup>u</sup>^<sup>1</sup> � <sup>u</sup>^2, D<sup>1</sup> <sup>≔</sup> <sup>3</sup>^

approximation for the Onsager coefficient, L r;r<sup>0</sup>

rotational diffusion coefficients are <sup>D</sup><sup>0</sup> <sup>¼</sup> kBT=ζ<sup>0</sup>, D<sup>0</sup>

regime <sup>D</sup> <sup>¼</sup> kBT<sup>=</sup> <sup>ζ</sup><sup>0</sup> <sup>þ</sup> ReΔζ ωð Þ <sup>¼</sup> <sup>0</sup> � �, DR <sup>¼</sup> kBT<sup>=</sup> <sup>ζ</sup><sup>0</sup>

ð<sup>∞</sup> 0

128πζ<sup>0</sup>ρ<sup>∗</sup>μ<sup>∗</sup><sup>4</sup> 3

105

<sup>l</sup> is the spherical Bessel function of order l, τ<sup>∗</sup>

S11

dxx<sup>4</sup> <sup>S</sup><sup>00</sup>

R

<sup>d</sup><sup>3</sup> ½ � g xð Þ� <sup>1</sup> , S<sup>11</sup>

, <sup>1</sup>ð Þ¼ <sup>x</sup> <sup>1</sup> <sup>þ</sup> <sup>ρ</sup><sup>∗</sup>

ð<sup>∞</sup> 0

ð<sup>∞</sup> 0

6π<sup>2</sup>ρ<sup>∗</sup>

þ

ReΔζ<sup>R</sup> <sup>¼</sup> <sup>128</sup>πρ<sup>∗</sup>μ<sup>∗</sup><sup>4</sup>ζ<sup>0</sup>

ReΔ<sup>ζ</sup> <sup>¼</sup> <sup>ζ</sup><sup>0</sup>

factor of the ferrofluid, which takes into account the micro-structural order of particles inside the ferrofluid as dictated by the underlying interaction potential. It is valid at high concentrations but ignores hydrodynamic interactions among particles. For monodisperse colloids, the tracer is equal to the other particles and the

the pair correlation function has the rotationally invariant expansion g r; Ω; Ω<sup>0</sup> ð Þ¼

<sup>δ</sup> <sup>r</sup> � <sup>r</sup><sup>0</sup> ð Þ<sup>δ</sup> <sup>Ω</sup> � <sup>Ω</sup><sup>0</sup> ð Þ, the unknown propagator in Eq. (4) <sup>χ</sup> <sup>¼</sup> exp �tL∘σ�<sup>1</sup> ð Þ,∘ ≔ <sup>Ð</sup>

!=r, u^ is the unitary dipole vector with orientation u^ ¼ u^ð Þ θ; φ . Using the Fick

! dΩ can be fully determined. For a spherical particle, the free translational and

, and ηsol the solvent viscosity. Taking the Laplace transform Δζ ωð Þ¼

τ∗ s,iso

<sup>2</sup> <sup>6</sup> <sup>þ</sup> <sup>x</sup><sup>2</sup> � � <sup>2</sup>τ<sup>∗</sup><sup>2</sup>

<sup>2</sup> <sup>6</sup> <sup>þ</sup> <sup>x</sup><sup>2</sup> � � <sup>67</sup>τ<sup>∗</sup><sup>2</sup>

Δζð Þt of the friction tensor in Eq. (6), it yields the real part of the translational and rotational diffusion coefficients of the tracer at the longtime overdamped

,<sup>0</sup> ð Þ� <sup>x</sup> <sup>1</sup> � �<sup>2</sup>

S<sup>00</sup> ,<sup>0</sup> ð Þ x

dxj2ð Þ x

dxj1ð Þ x

,αð Þ <sup>x</sup> <sup>=</sup> <sup>8</sup><sup>π</sup> <sup>x</sup> ½ � ð Þ <sup>2</sup> <sup>þ</sup> <sup>6</sup> , x <sup>≔</sup> kd. The structure factor components of the dipolar

,0ð Þ¼ <sup>x</sup> <sup>1</sup> <sup>þ</sup> <sup>ρ</sup><sup>∗</sup>

δn r 0 Ω0 0

whose dynamical evolution is attained from

χ r !00; r !0 ; Ω00; Ω<sup>0</sup>

ð Þ 0; 0 , and another particle with Ωð Þ θ; φ , separated

; <sup>Ω</sup>; <sup>Ω</sup><sup>0</sup> ð Þ¼ <sup>ρ</sup> <sup>D</sup><sup>0</sup>∇<sup>2</sup> <sup>þ</sup> <sup>D</sup><sup>0</sup>

<sup>R</sup> <sup>¼</sup> kBT=ζ<sup>0</sup>

!<sup>12</sup> � u^<sup>2</sup> � u^<sup>2</sup> � u^<sup>1</sup> h i

� �

<sup>R</sup> <sup>þ</sup> ReΔζRð Þ <sup>ω</sup> <sup>¼</sup> <sup>0</sup> � � [34], where

s, <sup>1</sup> <sup>þ</sup> <sup>3</sup>τ<sup>∗</sup><sup>2</sup> s,0 � �,

> s, <sup>1</sup> <sup>þ</sup> <sup>38</sup>τ<sup>∗</sup><sup>2</sup> s,0

,<sup>0</sup> ð Þ <sup>x</sup> � � � � ,

� �,

s,iso <sup>¼</sup> <sup>1</sup><sup>=</sup> <sup>x</sup><sup>2</sup> <sup>1</sup> <sup>þ</sup> <sup>1</sup>=S<sup>00</sup>

<sup>3</sup>d<sup>3</sup> ½ � <sup>h</sup>Δð Þþ <sup>x</sup> <sup>2</sup>hDð Þ <sup>x</sup> ,

<sup>3</sup>d<sup>3</sup> ½ � <sup>h</sup>Δð Þ� <sup>x</sup> hDð Þ <sup>x</sup> : (8)

�. Due to the spatial symmetry of the potential <sup>ψ</sup>,

¼ Ð d<sup>3</sup> r !00dΩ<sup>00</sup>

and taking further an

� �

; t ¼ 0

¼

, where

(7)

<sup>R</sup>∇<sup>Ω</sup>

<sup>R</sup>, with <sup>ζ</sup><sup>0</sup> <sup>¼</sup> <sup>3</sup>πηsold,

D E � �

� �

!00; Ω; Ω<sup>00</sup> � �

. The properties in Eq. (6) depend on the knowledge of the structure

. Notice that the pair correlation function is defined

r ! <sup>12</sup> � u^<sup>1</sup> ^ r

van Hove function C r!; Ω; r

at t = 0 the initial condition Ð

� �

by the mean distance r ¼ r

σ r !00; r !0 ; Ω00; Ω<sup>0</sup> � �

Pattern Formation and Stability in Magnetic Colloids

ð Þ 0; 0

between the tracer orientation Ω<sup>0</sup>

!00; Ω; Ω00; t � �

χ r !; r

σ r !; r !0 ; Ω; Ω<sup>0</sup> � �

ρg r ¼ r ! � � � � � �; <sup>Ω</sup>; <sup>Ω</sup><sup>0</sup>

^ r ! ¼r

ζ0

Ð <sup>∞</sup> <sup>0</sup> dte<sup>i</sup>ω<sup>t</sup>

d<sup>3</sup> r

j

liquid are

S<sup>00</sup>

,<sup>0</sup> ð Þ¼ <sup>x</sup> <sup>1</sup> <sup>þ</sup> <sup>ρ</sup><sup>∗</sup>

τ∗ s,<sup>α</sup> <sup>¼</sup> <sup>S</sup><sup>11</sup>

50

<sup>R</sup> <sup>¼</sup> πηsold<sup>3</sup>

!0 ; Ω<sup>0</sup> ; t

d<sup>3</sup> r

Eq. (4) after multiplying both sides by <sup>δ</sup>n r!; <sup>Ω</sup>; <sup>t</sup> <sup>¼</sup> <sup>0</sup>

� �

effective friction ζ ωð Þ. Here, we propose under the same conditions that the complex shear modulus of a ferrofluid is given by <sup>G</sup>ð Þ¼ <sup>ω</sup> <sup>G</sup><sup>0</sup> <sup>þ</sup> iG<sup>00</sup> <sup>¼</sup> <sup>i</sup>ω ζ<sup>0</sup> <sup>þ</sup> <sup>Δ</sup>ζ ωð Þ � �=3π<sup>d</sup> <sup>¼</sup> iωη ωð Þ. Since the complex viscosity of a viscoelastic media is defined by η ¼ η<sup>0</sup> � iη<sup>0</sup> , we find for the elastic and dissipative moduli, respectively [34]

2.6 Dynamic susceptibility response of a ferrofluid

DOI: http://dx.doi.org/10.5772/intechopen.81091

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

magnetic tensor function

where C<sup>11</sup>

function

Figure 6.

53

factor, and k ¼ 0.

h i Mkð Þ<sup>t</sup> <sup>M</sup>�kð Þ <sup>t</sup> <sup>¼</sup> <sup>0</sup> ¼ � <sup>m</sup><sup>2</sup>

� ∑ N i¼1

<sup>∂</sup>C kð Þ ; <sup>Ω</sup>; <sup>t</sup> ∂t ¼

by <sup>χ</sup>ð Þ¼ <sup>k</sup>;<sup>ω</sup> <sup>4</sup>πρ∗μ∗<sup>2</sup> ð Þ <sup>=</sup><sup>3</sup> <sup>S</sup><sup>11</sup>

ð d<sup>3</sup> r !e ik�r !i

> ðt 0 dt0 dΩ<sup>0</sup>

,0ð Þ� <sup>x</sup> <sup>i</sup>ωC<sup>11</sup>

The time-dependent magnetic susceptibility χð Þt of a ferrofluid can be measured either under a time-oscillating external magnetic field or in its absence with Altern Current dielectric spectroscopic techniques [51]. Our approach to calculating this property is based on the diffusion equation Eq. (4) that leads us to an expression for the collective dynamical property χð Þt regarding the friction Eq. (6). For this purpose, we multiplied Eq. (4) by δn tð Þ ¼ 0 ∑i,j u^ið Þt u^jð Þt , then we made a statistical average and Fourier transformed to write the dynamical equation of the collective

μ0N=12π � �C<sup>11</sup>

the hydrodynamic approach of Hess and Klein [52] to include rotational Brownian motion of the particles and derived the dynamical equation of the Van Hove

M k; Ω; Ω<sup>0</sup>

This equation is valid for colloidal suspensions of particles with axially symmetric potentials, where the memory function M k; Ω; Ω<sup>0</sup> ð Þ ; t is related to the friction Eq. (6), which does not require such an approximation on L [52]. The solution of this equation for the transversal component of the van Hove function C kð Þ ; Ω; t leads to the complex susceptibility function in the frequency domain, and it is given

Real χ<sup>0</sup> and imaginary χ<sup>00</sup> magnetic susceptibilities of a ferrofluid normalized to the initial value of susceptibility χ0. Picture (a) is a magnetic Cole-Cole plot and (b) depicts the susceptibility moduli versus frequency ω. The susceptibilities were evaluated at kmax the value of the wave number at the highest first peak of the structure

, <sup>1</sup>ð Þ k; t , Mkð Þt ≔ μ ffiffiffiffiffi

, <sup>1</sup> is the transversal component of the van Hove function. We extended

; t � t <sup>0</sup> ð ÞC k; Ω<sup>0</sup>

, <sup>1</sup>ð Þ <sup>x</sup>;<sup>ω</sup> � � <sup>þ</sup> 1. This equation reproduces

ð Þt dΩu^xð Þt δn rð Þ ; Ω; t , (11)

; t

<sup>0</sup> ð Þ (12)

μ0 <sup>p</sup> <sup>=</sup>ð Þ <sup>4</sup><sup>π</sup> <sup>3</sup>=<sup>2</sup> � �

G0 ð Þ¼ w ωη00ð Þ¼ ω ð<sup>∞</sup> 0 dx P<sup>00</sup> <sup>0</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup> s,iso � �<sup>2</sup> 1 þ ω<sup>∗</sup>τ<sup>∗</sup> s,iso � �<sup>2</sup> <sup>þ</sup> P<sup>11</sup> <sup>1</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup> s, 1 � �<sup>2</sup> 1 þ ω<sup>∗</sup>τ<sup>∗</sup> s, 1 � �<sup>2</sup> þ P<sup>11</sup> <sup>0</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup> s,0 � �<sup>2</sup> 1 þ ω<sup>∗</sup>τ<sup>∗</sup> s,0 � �<sup>2</sup> 2 6 4 3 7 5, G00ð Þ¼ w ωη<sup>0</sup> ð Þ¼ <sup>ω</sup> kBTω<sup>∗</sup> <sup>3</sup>πd<sup>3</sup> <sup>þ</sup> ð<sup>∞</sup> 0 dx P<sup>00</sup> <sup>0</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup> s,iso 1 þ ω<sup>∗</sup>τ<sup>∗</sup> s,iso � �<sup>2</sup> <sup>þ</sup> P<sup>11</sup> <sup>1</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup> s, 1 1 þ ω<sup>∗</sup>τ<sup>∗</sup> s, 1 � �<sup>2</sup> þ P<sup>11</sup> <sup>0</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup> s,0 1 þ ω<sup>∗</sup>τ<sup>∗</sup> s,0 � �<sup>2</sup> 2 6 4 3 7 5, P<sup>00</sup> <sup>0</sup> ð Þ¼ x kBTx<sup>4</sup> S<sup>00</sup> ,<sup>0</sup> ð Þ� <sup>x</sup> <sup>1</sup> � �<sup>2</sup> 18π<sup>3</sup>d<sup>3</sup> ρ∗S<sup>00</sup> ,<sup>0</sup> ð Þ x , P<sup>11</sup> <sup>1</sup> ð Þ¼� x 32kBTρ∗μ∗<sup>4</sup>j <sup>2</sup>ð Þ x 2 S11 , <sup>1</sup>ð Þ x <sup>9</sup>πd<sup>3</sup> , P<sup>11</sup> <sup>0</sup> ð Þ¼ x 48kBTρ∗μ∗<sup>4</sup>j <sup>2</sup>ð Þ x 2 S11 ,0ð Þ x <sup>9</sup>πd<sup>3</sup> , (10)

where <sup>ω</sup><sup>∗</sup> <sup>¼</sup> <sup>ω</sup>tB, tB <sup>¼</sup> <sup>d</sup><sup>2</sup> =D0. Figure 5 provides a comparison of the viscoelastic moduli of Eq. (10) with Langevin dynamics simulations for the same parameter as those of Figure 4.

We notice that the effective viscosity at the overdamped regime ω ¼ 0, which results from Eq. (10), is η<sup>0</sup> ¼ tBGo, where the storage modulus of the ferrofluid is Go <sup>¼</sup> kBT=3πd<sup>3</sup> � �D<sup>0</sup>=D. At low volume fraction <sup>ϕ</sup> <sup>¼</sup> πρ∗=6 of colloidal particles in the magnetic fluid, we found the following expression of the viscosity η<sup>0</sup> ≈ ηsolð Þ 1 þ 5ϕ=2 , which is the known exact Einstein result for low volume fractions if we assume that the free parameter ζ<sup>0</sup> can be replaced by the approximated form proposed by Mazur and Geigenmüller <sup>ζ</sup><sup>0</sup>ð Þ <sup>1</sup> <sup>þ</sup> <sup>1</sup>:5<sup>ϕ</sup> <sup>=</sup>ð Þ <sup>1</sup> � <sup>ϕ</sup> [50], which describes very well the free particle diffusion coefficient for all volume fractions.

#### Figure 5.

Logarithmic plot of elastic G<sup>0</sup> and loss G<sup>00</sup> moduli as a function of dimensionless frequency ω<sup>∗</sup> at three different concentrations and fixed dipole moment. Theory for G<sup>0</sup> (G00) is denoted by ∘(∗) and that with simulation by •, ∗.

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids DOI: http://dx.doi.org/10.5772/intechopen.81091

## 2.6 Dynamic susceptibility response of a ferrofluid

effective friction ζ ωð Þ. Here, we propose under the same conditions that the complex shear modulus of a ferrofluid is given by <sup>G</sup>ð Þ¼ <sup>ω</sup> <sup>G</sup><sup>0</sup> <sup>þ</sup> iG<sup>00</sup> <sup>¼</sup> <sup>i</sup>ω ζ<sup>0</sup> <sup>þ</sup> <sup>Δ</sup>ζ ωð Þ � �=3π<sup>d</sup> <sup>¼</sup> iωη ωð Þ. Since the complex viscosity of a viscoelastic media is defined by η ¼ η<sup>0</sup> � iη<sup>0</sup>

we find for the elastic and dissipative moduli, respectively [34]

<sup>0</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup>

1 þ ω<sup>∗</sup>τ<sup>∗</sup>

ð<sup>∞</sup> 0 dx

,<sup>0</sup> ð Þ� <sup>x</sup> <sup>1</sup> � �<sup>2</sup>

ρ∗S<sup>00</sup> ,<sup>0</sup> ð Þ x

s,iso � �<sup>2</sup> P<sup>11</sup>

s,iso

s,iso � �<sup>2</sup> þ

<sup>1</sup> ð Þ¼� x

<sup>2</sup>ð Þ x 2 S11 ,0ð Þ x

<sup>1</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup>

1 þ ω<sup>∗</sup>τ<sup>∗</sup>

s, 1 � �<sup>2</sup> P<sup>11</sup>

s, 1

<sup>2</sup>ð Þ x 2 S11 , <sup>1</sup>ð Þ x <sup>9</sup>πd<sup>3</sup> ,

s, 1 � �<sup>2</sup> þ

<sup>9</sup>πd<sup>3</sup> , (10)

(G00) is denoted by ∘(∗) and that with simulation by •, ∗.

<sup>0</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup>

1 þ ω<sup>∗</sup>τ<sup>∗</sup>

s,0 � �<sup>2</sup>

3 7 5,

s,0

3 7 5,

s,0 � �<sup>2</sup>

s,0 � �<sup>2</sup>

1 þ ω<sup>∗</sup>τ<sup>∗</sup>

P<sup>11</sup> <sup>0</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup>

s, 1 � �<sup>2</sup> þ

> P<sup>11</sup> <sup>1</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup>

1 þ ω<sup>∗</sup>τ<sup>∗</sup>

32kBTρ∗μ∗<sup>4</sup>j

=D0. Figure 5 provides a comparison of the viscoelastic

s,iso � �<sup>2</sup> þ

> P<sup>00</sup> <sup>0</sup> ð Þ <sup>x</sup> <sup>ω</sup><sup>∗</sup>τ<sup>∗</sup>

48kBTρ∗μ∗<sup>4</sup>j

the magnetic fluid, we found the following expression of the viscosity

2 6 4

1 þ ω<sup>∗</sup>τ<sup>∗</sup>

, P<sup>11</sup>

moduli of Eq. (10) with Langevin dynamics simulations for the same parameter as

We notice that the effective viscosity at the overdamped regime ω ¼ 0, which results from Eq. (10), is η<sup>0</sup> ¼ tBGo, where the storage modulus of the ferrofluid is Go <sup>¼</sup> kBT=3πd<sup>3</sup> � �D<sup>0</sup>=D. At low volume fraction <sup>ϕ</sup> <sup>¼</sup> πρ∗=6 of colloidal particles in

η<sup>0</sup> ≈ ηsolð Þ 1 þ 5ϕ=2 , which is the known exact Einstein result for low volume fractions if we assume that the free parameter ζ<sup>0</sup> can be replaced by the approximated form proposed by Mazur and Geigenmüller <sup>ζ</sup><sup>0</sup>ð Þ <sup>1</sup> <sup>þ</sup> <sup>1</sup>:5<sup>ϕ</sup> <sup>=</sup>ð Þ <sup>1</sup> � <sup>ϕ</sup> [50], which describes very well the free particle diffusion coefficient for all volume fractions.

Logarithmic plot of elastic G<sup>0</sup> and loss G<sup>00</sup> moduli as a function of dimensionless frequency ω<sup>∗</sup> at three different

P<sup>00</sup>

2 6 4

G0

ð Þ¼ w ωη00ð Þ¼ ω

P<sup>00</sup> <sup>0</sup> ð Þ¼ x

those of Figure 4.

Figure 5.

52

concentrations and fixed dipole moment. Theory for G<sup>0</sup>

where <sup>ω</sup><sup>∗</sup> <sup>¼</sup> <sup>ω</sup>tB, tB <sup>¼</sup> <sup>d</sup><sup>2</sup>

G00ð Þ¼ w ωη<sup>0</sup>

ð<sup>∞</sup> 0 dx

Pattern Formation and Stability in Magnetic Colloids

ð Þ¼ <sup>ω</sup> kBTω<sup>∗</sup>

<sup>3</sup>πd<sup>3</sup> <sup>þ</sup>

kBTx<sup>4</sup> S<sup>00</sup>

18π<sup>3</sup>d<sup>3</sup>

P<sup>11</sup> <sup>0</sup> ð Þ¼ x ,

The time-dependent magnetic susceptibility χð Þt of a ferrofluid can be measured either under a time-oscillating external magnetic field or in its absence with Altern Current dielectric spectroscopic techniques [51]. Our approach to calculating this property is based on the diffusion equation Eq. (4) that leads us to an expression for the collective dynamical property χð Þt regarding the friction Eq. (6). For this purpose, we multiplied Eq. (4) by δn tð Þ ¼ 0 ∑i,j u^ið Þt u^jð Þt , then we made a statistical average and Fourier transformed to write the dynamical equation of the collective magnetic tensor function

$$
\langle \mathsf{M}\_k(t)\mathsf{M}\_{-k}(t=0)\rangle = -\left(m^2\mu\_0 N/12\pi\right)\mathsf{C}\_{,1}^{11}(k,t),\\\mathsf{M}\_k(t) \coloneqq \left(\mu\sqrt{\mu\_0}/(4\pi)^{3/2}\right)
$$

$$
\times \sum\_{i=1}^N \int d^3 \vec{r} \,\vec{\epsilon}^{k\cdot\overrightarrow{r}\_i}(t) d\Omega \hat{u}\_\mathbf{x}(t) \delta n(r,\Omega,t),\tag{11}
$$

where C<sup>11</sup> , <sup>1</sup> is the transversal component of the van Hove function. We extended the hydrodynamic approach of Hess and Klein [52] to include rotational Brownian motion of the particles and derived the dynamical equation of the Van Hove function

$$\frac{\partial \mathcal{C}(k,\Omega,t)}{\partial t} = \int\_0^t dt' d\Omega' M(k,\Omega,\Omega',t-t') \mathcal{C}(k,\Omega',t') \tag{12}$$

This equation is valid for colloidal suspensions of particles with axially symmetric potentials, where the memory function M k; Ω; Ω<sup>0</sup> ð Þ ; t is related to the friction Eq. (6), which does not require such an approximation on L [52]. The solution of this equation for the transversal component of the van Hove function C kð Þ ; Ω; t leads to the complex susceptibility function in the frequency domain, and it is given by <sup>χ</sup>ð Þ¼ <sup>k</sup>;<sup>ω</sup> <sup>4</sup>πρ∗μ∗<sup>2</sup> ð Þ <sup>=</sup><sup>3</sup> <sup>S</sup><sup>11</sup> ,0ð Þ� <sup>x</sup> <sup>i</sup>ωC<sup>11</sup> , <sup>1</sup>ð Þ <sup>x</sup>;<sup>ω</sup> � � <sup>þ</sup> 1. This equation reproduces

#### Figure 6.

Real χ<sup>0</sup> and imaginary χ<sup>00</sup> magnetic susceptibilities of a ferrofluid normalized to the initial value of susceptibility χ0. Picture (a) is a magnetic Cole-Cole plot and (b) depicts the susceptibility moduli versus frequency ω. The susceptibilities were evaluated at kmax the value of the wave number at the highest first peak of the structure factor, and k ¼ 0.

similar results of a model molecular liquid of Wei and Patey [53] who, however, used a Kerr dynamical equation for C. Since χ ¼ χ<sup>0</sup> þ iχ<sup>00</sup> in Figure 6, we provide typical plots of the real and imaginary components of the magnetic susceptibility of the ferrofluid Fe2O3.

Conflict of interest

Author details

55

Ricardo Peredo Ortíz<sup>1</sup>

Raquel Hernández Gómez<sup>2</sup>

Polytechnic Institute, CD Mexico, Mexico

provided the original work is properly cited.

National Polytechnic Institute, CD Mexico, Mexico

\*Address all correspondence to: marther@fis.cinvestav.mx

, Martin Hernández Contreras<sup>1</sup>

1 Physics Department, Center for Research and Advanced Studies of the National

2 Computation Department, Center for Research and Advanced Studies of the

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

\* and

There are no conflicts to declare.

DOI: http://dx.doi.org/10.5772/intechopen.81091

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

Alternative kinetics methods for low-density ferrofluids are given in [40].

## 3. Conclusion

This review highlights several active research lines on the structure and dynamics in ferrofluids. These investigations are motivated by many experimental observations that colloidal ferromagnetic particles both in bulk or in thin film aggregate, building several macroscopic structures of practical and scientific interest. As for all applications where colloidal materials are being processed, it is necessary to understand the interparticle interactions to have better control of the resulting structures and their properties. The interactions are known experimentally under specific conditions through the measured phase diagrams [11]; it has been found [29] that a useful potential that represents well the structure factor of ferrofluids is a Lennard-Jones plus dipolar interaction between pairs of particles. Therefore, our use of this potential in an evolutionary genetic algorithm leads to Figure 1, which shows transient metastable crystalline structures of a monodisperse ferrofluid at T ¼ 0K. On the other hand, Figure 2 shows the predicted structure factors using material parameters of typical Fe2O<sup>3</sup> ferrofluid. There we observe that for low particle density and dipole moment, the particles are dispersed in the fluid, whereas for the highest dipole moment system, they form chains. As a consequence, there are noticeable effects on the collective diffusion coefficient as it is shown in Figure 3, where we notice that as the strength of the dipole moment is increased for lowdensity suspensions, this property drops because of particles' association in larger structures (chains). However, for higher-density suspensions, particles re-disperse again, therefore, increasing their collective diffusion. In Section 2.4, we presented the stochastic Langevin approach to describe the diffusion of a tracer particle in colloidal magnetic fluids. Figure 4 provides theoretical predictions of this theory for the longtime translational and rotational self-diffusion coefficients Eq. (7) of a ferroparticle for the same material parameters of Fe2O3. We notice that the theory prediction has an agreement with Langevin dynamics simulations. In Section 2.5, we described our extension of a hydrodynamic theory of colloid dynamics [52] to include the rotational Brownian motions of ferrofluid's particles. As a result, we derived an equation of the complex susceptibility χð Þ k; w that depends on the diffusion constants Dð Þ ω , DRð Þ ω of a tracer particle in the colloid, which were given in Section 2.4. Thus, in Figure 6, we show the predicted magnetic Cole-Cole plot and real and imaginary components of the susceptibility moduli for Fe2O<sup>3</sup> ferrofluid. There remain further studies to find the successful comparison of this method with well-characterized magnetic colloidal fluids and computer simulations in model ferrofluid systems.

## Acknowledgements

The authors acknowledge the General Coordination of Information and Communications Technologies (CGSTIC) at CINVESTAV for providing HPC resources on the Hybrid Supercomputer <sup>00</sup>Xiuhcoatl<sup>00</sup> that has contributed to the research results reported within this chapter.

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids DOI: http://dx.doi.org/10.5772/intechopen.81091

## Conflict of interest

similar results of a model molecular liquid of Wei and Patey [53] who, however, used a Kerr dynamical equation for C. Since χ ¼ χ<sup>0</sup> þ iχ<sup>00</sup> in Figure 6, we provide typical plots of the real and imaginary components of the magnetic susceptibility of

Alternative kinetics methods for low-density ferrofluids are given in [40].

This review highlights several active research lines on the structure and dynamics in ferrofluids. These investigations are motivated by many experimental observations that colloidal ferromagnetic particles both in bulk or in thin film aggregate, building several macroscopic structures of practical and scientific interest. As for all applications where colloidal materials are being processed, it is necessary to understand the interparticle interactions to have better control of the resulting structures and their properties. The interactions are known experimentally under specific conditions through the measured phase diagrams [11]; it has been found [29] that a useful potential that represents well the structure factor of ferrofluids is a Lennard-Jones plus dipolar interaction between pairs of particles. Therefore, our use of this potential in an evolutionary genetic algorithm leads to Figure 1, which shows transient metastable crystalline structures of a monodisperse ferrofluid at T ¼ 0K. On the other hand, Figure 2 shows the predicted structure factors using material parameters of typical Fe2O<sup>3</sup> ferrofluid. There we observe that for low particle density and dipole moment, the particles are dispersed in the fluid, whereas for the highest dipole moment system, they form chains. As a consequence, there are noticeable effects on the collective diffusion coefficient as it is shown in Figure 3, where we notice that as the strength of the dipole moment is increased for lowdensity suspensions, this property drops because of particles' association in larger structures (chains). However, for higher-density suspensions, particles re-disperse again, therefore, increasing their collective diffusion. In Section 2.4, we presented the stochastic Langevin approach to describe the diffusion of a tracer particle in colloidal magnetic fluids. Figure 4 provides theoretical predictions of this theory for the longtime translational and rotational self-diffusion coefficients Eq. (7) of a ferroparticle for the same material parameters of Fe2O3. We notice that the theory prediction has an agreement with Langevin dynamics simulations. In Section 2.5, we described our extension of a hydrodynamic theory of colloid dynamics [52] to include the rotational Brownian motions of ferrofluid's particles. As a result, we derived an equation of the complex susceptibility χð Þ k; w that depends on the diffusion constants Dð Þ ω , DRð Þ ω of a tracer particle in the colloid, which were given in Section 2.4. Thus, in Figure 6, we show the predicted magnetic Cole-Cole plot

and real and imaginary components of the susceptibility moduli for Fe2O<sup>3</sup> ferrofluid. There remain further studies to find the successful comparison of this method with well-characterized magnetic colloidal fluids and computer simulations

The authors acknowledge the General Coordination of Information and Communications Technologies (CGSTIC) at CINVESTAV for providing HPC resources on the Hybrid Supercomputer <sup>00</sup>Xiuhcoatl<sup>00</sup> that has contributed to the research

in model ferrofluid systems.

results reported within this chapter.

Acknowledgements

54

the ferrofluid Fe2O3.

Pattern Formation and Stability in Magnetic Colloids

3. Conclusion

There are no conflicts to declare.

## Author details

Ricardo Peredo Ortíz<sup>1</sup> , Martin Hernández Contreras<sup>1</sup> \* and Raquel Hernández Gómez<sup>2</sup>

1 Physics Department, Center for Research and Advanced Studies of the National Polytechnic Institute, CD Mexico, Mexico

2 Computation Department, Center for Research and Advanced Studies of the National Polytechnic Institute, CD Mexico, Mexico

\*Address all correspondence to: marther@fis.cinvestav.mx

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## References

[1] Rosensweig RE. Ferrohydrodynamics. Cambridge: Cambridge University Press; 1985

[2] Jordan A, Scholz R, Maier-Hauff K, Johannsen M, Wust P, Nadobny J, et al. Presentation of a new magnetic field therapy system for the treatment of human solid tumors with magnetic fluid hyperthermia. Journal of Magnetism and Magnetic Materials. 2001;225:118-126

[3] Kim SH, Lee SY, Yang SM, Yi GR. Self-assembled colloidal structures for photonics. NPG Asia Materials. 2011;3: 25-33

[4] He L, Wang H, Ge J, Yin Y. Magnetic assembly route to colloidal responsive photonic nanostructures. Accounts of Chemical Research. 2012;45:1431-1440

[5] Andelman D, Rosensweig RE. Modulated Phases: Review and Recent Results. The Journal of Physical Chemistry. B. 2009;113:3785-3798

[6] Klapp SHL. Collective dynamics of dipolar and multipolar colloids: from passive to active systems. Current Opinion in Colloid & Interface Science. 2016;21:76-85

[7] Tierno P. Recent advances in anisotropic magnetic colloids: Realization, assembly and applications. Physical Chemistry Chemical Physics. 2014;16:23515-23528

[8] Schilling R, Scheidsteger T. Mode coupling approach to the ideal glass transition of molecular liquids: Linear molecules. Physical Review E. 1997;56: 2932-2949

[9] Snezhko A. Non-equilibrium magnetic colloidal dispersions at liquid– air interfaces: dynamic patterns, magnetic order and self-assembled swimmers. Journal of Physics: Condensed Matter. 2011;15:153101 (1-21) [10] Ramírez-González PE, Medina Noyola M. Glass transition in soft-sphere dispersions. Journal of Physics: Condensed Matter. 2009;21:075101 (1-13) The Journal of Chemical Physics. 2006;

DOI: http://dx.doi.org/10.5772/intechopen.81091

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids

[27] Meriguet G, Dubois E, Bourdon A, Demouchy G, Dupuis V, Perzynski R. Forced Rayleigh scattering experiments in concentrated magnetic fluids: effect of interparticle interactions on the diffusion coefficient. Journal of Magnetism and Magnetic Materials.

[28] Lal J, Abernathy D, Auvray L, Diat O, Grubel G. Dynamics and correlations in magnetic colloidal systems studied by X-ray photon correlation spectroscopy. European Physical Journal E: Soft Matter and Biological. Physics. 2001;4:

[29] Meriguet G, Dubois E, Jardat M, Bourdon A, Demouchy G, Dupuis V, et al. Understanding the structure and the dynamics of magnetic fluids: coupling of experiment and simulation. Journal of Physics: Condensed Matter.

[30] Meriguet G, Jardat M, Turq P. Structural properties of chargestabilized ferrofluids under a magnetic field: A Brownian dynamics study. The Journal of Chemical Physics. 2004;121:

[31] Meriguet G, Jardat M, Turq P. Brownian dynamics investigation of magnetization and birefringence relaxations in ferrofluids. The Journal of Chemical Physics. 2005;123:144915 (1-8)

[32] Hernández-Contreras M, Ruíz Estrada H. Transport properties of ferrofluids. Physical Review E. 2003;68:

[33] Peredo-Ortíz R, Hernández Contreras M, Hernández Gómez R. Diffusion in monodisperse ferrofluids. Journal of Physics: Conference Series.

[34] Peredo-Ortíz R, Hernández Contreras M. Diffusion microrheology of ferrofluids. Revista Mexicana de

2017;792:012099 (1-11)

Fisica. 2018;64:82-93

2006;18:S2685-S2696

6078-6085

031202 (1-7)

2005;289:39-42

263-271

[19] Wagner J, Fischer B, Autenrieth T,

dynamics of charged magnetic colloids. Journal of Physics: Condensed Matter.

[20] Martin GAR, Bradbury A, Chantrell RW. A comparison between a threedimensional Monte Carlo model and a hypernnetted chain approximation of a ferrofluid. Journal of Magnetism and Magnetic Materials. 1987;65:177-180

[21] Bradbury A, Martin GAR, Chantrell RW. Zero field particle correlations in a

simulations and the hypernetted chain approximation. Journal of Magnetism and Magnetic Materials. 1987;69:5-11

[22] Hayter JB, Pynn R. Structure factor of a magnetically saturated ferrofluid. Physical Review Letters. 1982;49:1103-

[23] Pynn R, Hayter JB, Charles SW. Determination of ferrofluid structure by neutron polarization analysis. Physical Review Letters. 1983;51:710-713

[24] Bacri JC, Cebers A, Bourdon A, Demouchy G, Heegaard BM, Perzynski R. Forced rayleigh experiment in a magnetic fluid. Physical Review Letters.

[25] Gazeau F, Dubois E, Bacri JC, Boue F, Cebers A, Perzynski R. Anisotropy of the structure factor of magnetic fluids under a field probed by small-angle neutron scattering. Physical Review E.

[26] Bacri JC, Cebers A, Bourdon A, Demouchy G, Heegaard BM, Kashevsky B, et al. Transient grating in a ferrofluid under magnetic field: Effect of magnetic interactions on the diffusion coefficient of translation. Physical Review E. 1995;

1995;74:5032-5935

2002;65:031403 (1-15)

52:3936-3942

57

ferrofluid using Monte Carlo

1106

Hempelmann R. Structure and

124:114901 (1-10)

2006;18:S2697-S2711

[11] Yethiraj A, van Blaaderen A. A colloidal model system with an interaction tunable from hard sphere to soft and dipolar. Nature. 2003;421: 513-517

[12] Klokkenburg M, Erne BH, Meeldijk JD, Wiedenmann A, Pethukov AV, Dullens RPA, et al. In situ imaging of field-induced hexagonal columns in magnetite ferrofluids. Physical Review 36 Letters. 2006;97:185702 (1-4)

[13] Cousin F, Dubois E, Cabuil V. Tuning the interactions of a magnetic colloid suspension. Physical Review E. 2003;68:021405 (1-9)

[14] Pham AT, Zhuang Y, Detwiler P, Socolar JES, Charbonneau P, Yellen BB. Phase diagram and aggregation dynamics of a monolayer of paramagnetic colloids. Physical Review E. 2017;95:052607 (1-12)

[15] Hynninen AP, Dijkstra M. Phase Diagram of dipolar hard and soft spheres: Manipulation of colloidal crystal structures by an external field. Physical Review Letters. 2005;94: 138303 (1-4)

[16] Spiteri L, Messina R. Dipolar Crystals: The crucial role of the clinohexagonal prism phase. Physical Review Letters. 2017;119:155501 (1-4)

[17] Weis JJ, Levesque D. Advanced Computer Simulation Approaches for Soft Matter Sciences II. Berlin Heidelberg: Springer; 2005. pp. 163-225

[18] Wagner J, Fischer B, Autenrieth T. Field induced anisotropy of charged magnetic colloids: A rescaled mean spherical approximation study.

Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids DOI: http://dx.doi.org/10.5772/intechopen.81091

The Journal of Chemical Physics. 2006; 124:114901 (1-10)

References

25-33

2016;21:76-85

2932-2949

56

2014;16:23515-23528

[1] Rosensweig RE.

Ferrohydrodynamics. Cambridge: Cambridge University Press; 1985

[2] Jordan A, Scholz R, Maier-Hauff K, Johannsen M, Wust P, Nadobny J, et al. Presentation of a new magnetic field therapy system for the treatment of human solid tumors with magnetic fluid hyperthermia. Journal of Magnetism and Magnetic Materials. 2001;225:118-126

Pattern Formation and Stability in Magnetic Colloids

[10] Ramírez-González PE, Medina Noyola M. Glass transition in soft-sphere

[11] Yethiraj A, van Blaaderen A. A colloidal model system with an

Condensed Matter. 2009;21:075101 (1-13)

interaction tunable from hard sphere to soft and dipolar. Nature. 2003;421:

[12] Klokkenburg M, Erne BH, Meeldijk JD, Wiedenmann A, Pethukov AV, Dullens RPA, et al. In situ imaging of field-induced hexagonal columns in magnetite ferrofluids. Physical Review 36 Letters. 2006;97:185702 (1-4)

[13] Cousin F, Dubois E, Cabuil V. Tuning the interactions of a magnetic colloid suspension. Physical Review E.

[14] Pham AT, Zhuang Y, Detwiler P, Socolar JES, Charbonneau P, Yellen BB.

paramagnetic colloids. Physical Review

[15] Hynninen AP, Dijkstra M. Phase Diagram of dipolar hard and soft spheres: Manipulation of colloidal crystal structures by an external field. Physical Review Letters. 2005;94:

[16] Spiteri L, Messina R. Dipolar Crystals: The crucial role of the clinohexagonal prism phase. Physical Review Letters. 2017;119:155501 (1-4)

[17] Weis JJ, Levesque D. Advanced Computer Simulation Approaches for

Heidelberg: Springer; 2005. pp. 163-225

[18] Wagner J, Fischer B, Autenrieth T. Field induced anisotropy of charged magnetic colloids: A rescaled mean spherical approximation study.

Soft Matter Sciences II. Berlin

Phase diagram and aggregation dynamics of a monolayer of

E. 2017;95:052607 (1-12)

138303 (1-4)

2003;68:021405 (1-9)

dispersions. Journal of Physics:

513-517

[3] Kim SH, Lee SY, Yang SM, Yi GR. Self-assembled colloidal structures for photonics. NPG Asia Materials. 2011;3:

[4] He L, Wang H, Ge J, Yin Y. Magnetic assembly route to colloidal responsive photonic nanostructures. Accounts of Chemical Research. 2012;45:1431-1440

[5] Andelman D, Rosensweig RE. Modulated Phases: Review and Recent Results. The Journal of Physical Chemistry. B. 2009;113:3785-3798

[6] Klapp SHL. Collective dynamics of dipolar and multipolar colloids: from passive to active systems. Current Opinion in Colloid & Interface Science.

[7] Tierno P. Recent advances in anisotropic magnetic colloids:

Realization, assembly and applications. Physical Chemistry Chemical Physics.

[8] Schilling R, Scheidsteger T. Mode coupling approach to the ideal glass transition of molecular liquids: Linear molecules. Physical Review E. 1997;56:

[9] Snezhko A. Non-equilibrium

air interfaces: dynamic patterns, magnetic order and self-assembled swimmers. Journal of Physics:

magnetic colloidal dispersions at liquid–

Condensed Matter. 2011;15:153101 (1-21)

[19] Wagner J, Fischer B, Autenrieth T, Hempelmann R. Structure and dynamics of charged magnetic colloids. Journal of Physics: Condensed Matter. 2006;18:S2697-S2711

[20] Martin GAR, Bradbury A, Chantrell RW. A comparison between a threedimensional Monte Carlo model and a hypernnetted chain approximation of a ferrofluid. Journal of Magnetism and Magnetic Materials. 1987;65:177-180

[21] Bradbury A, Martin GAR, Chantrell RW. Zero field particle correlations in a ferrofluid using Monte Carlo simulations and the hypernetted chain approximation. Journal of Magnetism and Magnetic Materials. 1987;69:5-11

[22] Hayter JB, Pynn R. Structure factor of a magnetically saturated ferrofluid. Physical Review Letters. 1982;49:1103- 1106

[23] Pynn R, Hayter JB, Charles SW. Determination of ferrofluid structure by neutron polarization analysis. Physical Review Letters. 1983;51:710-713

[24] Bacri JC, Cebers A, Bourdon A, Demouchy G, Heegaard BM, Perzynski R. Forced rayleigh experiment in a magnetic fluid. Physical Review Letters. 1995;74:5032-5935

[25] Gazeau F, Dubois E, Bacri JC, Boue F, Cebers A, Perzynski R. Anisotropy of the structure factor of magnetic fluids under a field probed by small-angle neutron scattering. Physical Review E. 2002;65:031403 (1-15)

[26] Bacri JC, Cebers A, Bourdon A, Demouchy G, Heegaard BM, Kashevsky B, et al. Transient grating in a ferrofluid under magnetic field: Effect of magnetic interactions on the diffusion coefficient of translation. Physical Review E. 1995; 52:3936-3942

[27] Meriguet G, Dubois E, Bourdon A, Demouchy G, Dupuis V, Perzynski R. Forced Rayleigh scattering experiments in concentrated magnetic fluids: effect of interparticle interactions on the diffusion coefficient. Journal of Magnetism and Magnetic Materials. 2005;289:39-42

[28] Lal J, Abernathy D, Auvray L, Diat O, Grubel G. Dynamics and correlations in magnetic colloidal systems studied by X-ray photon correlation spectroscopy. European Physical Journal E: Soft Matter and Biological. Physics. 2001;4: 263-271

[29] Meriguet G, Dubois E, Jardat M, Bourdon A, Demouchy G, Dupuis V, et al. Understanding the structure and the dynamics of magnetic fluids: coupling of experiment and simulation. Journal of Physics: Condensed Matter. 2006;18:S2685-S2696

[30] Meriguet G, Jardat M, Turq P. Structural properties of chargestabilized ferrofluids under a magnetic field: A Brownian dynamics study. The Journal of Chemical Physics. 2004;121: 6078-6085

[31] Meriguet G, Jardat M, Turq P. Brownian dynamics investigation of magnetization and birefringence relaxations in ferrofluids. The Journal of Chemical Physics. 2005;123:144915 (1-8)

[32] Hernández-Contreras M, Ruíz Estrada H. Transport properties of ferrofluids. Physical Review E. 2003;68: 031202 (1-7)

[33] Peredo-Ortíz R, Hernández Contreras M, Hernández Gómez R. Diffusion in monodisperse ferrofluids. Journal of Physics: Conference Series. 2017;792:012099 (1-11)

[34] Peredo-Ortíz R, Hernández Contreras M. Diffusion microrheology of ferrofluids. Revista Mexicana de Fisica. 2018;64:82-93

[35] Mertelj A, Rešetič A, Gyergyek S, Markovec D, Čopič M. Anisotropic microrheological properties of chainforming magnetic fluids. Soft Matter. 2011;7:125-131

[36] Yendeti B, Thirupathi G, Vudaygiri A, Singh R. Field-dependent anisotropic microrheological and microstructural properties of dilute ferrofluids. European Physical Journal E: Soft Matter and Biological Physics. 2014;37: 10 (1-8)

[37] Mazoyer S, Ebert F, Maret G, Keim P. Dynamics of particles and cages in an experimental 2D glass former. Europhysics Letters. 2009;88:66004 (1-6)

[38] Lacoste D, Lubensky TC. Phase transitions in a ferrofluid at magneticfield-induced microphase separation. Physical Review E. 2001;64:041506 (1-8)

[39] de Gennes PG, Pincus PA. Pair correlations in a ferromagnetic colloid. Physics Kondens Matter. 1970;11:189-194

[40] Odenbach S. Colloidal Magnetic Fluids. 1st ed. Berlin: Springer; 2009

[41] Tlusty T, Safran SA. Defect-induced phase separation in dipolar fluids. Science. 2000;290:1328-1331

[42] Sreekumari A, Ilg P. Slow relaxation in structure-forming ferrofluids. Physical Review E. 2013;88:042315 (1-9)

[43] Russel WB, Saville DA, Showalter WR. Colloidal Dispersions. New York NY: Cambridge University Press; 1989

[44] Israelachvili JN. van der Waals dispersion force contribution to works of adhesion and contact angles on the basis of macroscopic theory. Journal of the Chemical Society, Faraday Transactions 2. 1973;69:1729-1738

[45] Fornleitner J, Loverso F, Kahl G, Likos CN. Genetic algorithms predict formation of exotic ordered

configurations for two-component dipolar monolayers. Soft Matter. 2008; 4:480-484

Chapter 4

Abstract

Raoudha Chaabane

the streamlines widely.

free convection, LBM

1. Introduction

59

The Study of Magneto-Convection

A free convection heat transfer in a sinusoidally heated enclosure filled with conducting fluid is presented in this chapter by Lattice Boltzmann Method (LBM). The horizontal walls in the enclosures are insulated and there is an opening part on the right wall. The right non-open parts of the vertical wall of the square cavity are maintained at constant cold temperature and the left wall of the cavity is sinusoidally heated. The cavity is get under a uniform in-plane magnetic field. The main aim of this study is to highlight the effectiveness of the LBM mesoscopic approach to predict the effects of pertinent parameters such as the Hartmann number varying from 0 to 150 where Rayleigh number is fixed at moderate value of 10<sup>5</sup> on flow patterns. This in-house numerical code used in this chapter is ascertained and a good agreement with literature is highlighted. The appropriate validation with previous numerical investigations demonstrated that this attitude is a suitable method and a powerful approach for engineering MHD problems. Findings and results show the alterations of Hartmann number that influence the isotherms and

Keywords: sinusoidal thermal boundary condition, MHD, partially open cavity,

Convective flow and heat transfer in an open cavity has been studied due to the extensive range of applications in engineering science and technology that consider various combinations of imposed temperature gradients and enclosure sketches [1–10]. Open cavity with a modified linear or sinusoidal thermal boundary condition is encountered in many practical engineering and industrial applications, such as solar energy collection, cooling of electronic devices, material processing, grain storage, flow and heat transfer in solar ponds, high-performance insulation for buildings, dynamics of lakes, reservoirs and cooling ponds, crystal growing, float glass production, metal casting, food processing, galvanizing, metal coating, and so on. Besides, pertinent useful numerical research works had been conducted to simulate the MHD free convection under nonuniform thermic boundaries where recent attention has been intensively focused on the cases of mixed boundary conditions [11–14]. However, few results have been reported for free convection caused

Heat Transfer in a Partially Open

Cavity Based on LBM

[46] Autenrieth T, Robert A, Wagner J, Grubel G. The dynamic behavior of magnetic colloids in suspension. Journal of Applied Crystallography. 2007;40: 5250-5253

[47] Pyanzina E, Kantorovich S, Cerda JJ, Ivanov A, Holm C. How to analyze the structure factor in ferrofluids with strong magnetic interactions: a combined analytic and simulation approach. Molecular Physics. 2009;107: 571-590

[48] Kirchoff T, Lowen H, Klein R. Dynamical correlations in suspensions of charged rodlike macromolecules. Physical Review E. 1996;53:5011-5022

[49] Mason TG, Weitz DA. Optical measurements of frequency-dependent viscoelastic moduli of complex fluids. Physical Review Letters. 1995;74:1250- 1253

[50] Mazur P, Geigenmüller U. A simple formula for the short-time self-diffusion coefficient in concentrated suspensions. Physica A. 1987;146:657-661

[51] Roeben E, Roeder L, Teusch S, Effertz M, Deiters UK, Schmidt AM. Magnetic particle nanorheology. Colloid & Polymer Science. 2014;292:2013-2023

[52] Hess W, Klein R. Generalized hydrodynamics of systems of Brownian particles. Advances in Physics. 1983;32: 172-283

[53] Wei D, Patey GN. Rotational motion in molecular liquids. The Journal of Chemical Physics. 1989;91:7113-7129

## Chapter 4

[35] Mertelj A, Rešetič A, Gyergyek S, Markovec D, Čopič M. Anisotropic microrheological properties of chainforming magnetic fluids. Soft Matter.

Pattern Formation and Stability in Magnetic Colloids

configurations for two-component dipolar monolayers. Soft Matter. 2008;

[46] Autenrieth T, Robert A, Wagner J, Grubel G. The dynamic behavior of magnetic colloids in suspension. Journal of Applied Crystallography. 2007;40:

[47] Pyanzina E, Kantorovich S, Cerda JJ, Ivanov A, Holm C. How to analyze the structure factor in ferrofluids with strong magnetic interactions: a combined analytic and simulation approach. Molecular Physics. 2009;107:

[48] Kirchoff T, Lowen H, Klein R. Dynamical correlations in suspensions of charged rodlike macromolecules. Physical Review E. 1996;53:5011-5022

[49] Mason TG, Weitz DA. Optical measurements of frequency-dependent viscoelastic moduli of complex fluids. Physical Review Letters. 1995;74:1250-

[50] Mazur P, Geigenmüller U. A simple formula for the short-time self-diffusion coefficient in concentrated suspensions.

Physica A. 1987;146:657-661

[51] Roeben E, Roeder L, Teusch S, Effertz M, Deiters UK, Schmidt AM. Magnetic particle nanorheology. Colloid & Polymer Science. 2014;292:2013-2023

[52] Hess W, Klein R. Generalized hydrodynamics of systems of Brownian particles. Advances in Physics. 1983;32:

[53] Wei D, Patey GN. Rotational motion in molecular liquids. The Journal of Chemical Physics. 1989;91:7113-7129

4:480-484

5250-5253

571-590

1253

172-283

[36] Yendeti B, Thirupathi G, Vudaygiri A, Singh R. Field-dependent anisotropic microrheological and microstructural properties of dilute ferrofluids. European Physical Journal E: Soft Matter and Biological Physics. 2014;37:

[37] Mazoyer S, Ebert F, Maret G, Keim P. Dynamics of particles and cages in an experimental 2D glass former.

Europhysics Letters. 2009;88:66004 (1-6)

[38] Lacoste D, Lubensky TC. Phase transitions in a ferrofluid at magneticfield-induced microphase separation. Physical Review E. 2001;64:041506 (1-8)

[39] de Gennes PG, Pincus PA. Pair correlations in a ferromagnetic colloid. Physics Kondens Matter. 1970;11:189-194

[40] Odenbach S. Colloidal Magnetic Fluids. 1st ed. Berlin: Springer; 2009

[41] Tlusty T, Safran SA. Defect-induced phase separation in dipolar fluids. Science. 2000;290:1328-1331

[42] Sreekumari A, Ilg P. Slow relaxation

Physical Review E. 2013;88:042315 (1-9)

[43] Russel WB, Saville DA, Showalter WR. Colloidal Dispersions. New York NY: Cambridge University Press; 1989

[44] Israelachvili JN. van der Waals dispersion force contribution to works of adhesion and contact angles on the basis of macroscopic theory. Journal of

[45] Fornleitner J, Loverso F, Kahl G, Likos CN. Genetic algorithms predict

the Chemical Society, Faraday Transactions 2. 1973;69:1729-1738

formation of exotic ordered

58

in structure-forming ferrofluids.

2011;7:125-131

10 (1-8)
