Condensation, demixing, and orientational ordering of magnetic colloidal suspensions

Pair correlations of the easy magnetisation axes of superparamagnetic nanoparticles in a ferrofluid/ferrocomposite

The widespread use of magnetic nanoparticles in modern technologies and medical applications highlights the need for reliable theoretical models that can predict their physical properties. The pair correlation function of two randomly selected superparamagnetic nanoparticles in a ferrofluid/ferrocomposite is studied to depict the joint probability density of the easy magnetisation axes across the planes of parameters of major importance; these are the interaction of ferroparticles with an external magnetic field, the energy of magnetic anisotropy inside the superparamagnetic nanoparticle, and the interparticle magnetic dipole-dipole interaction. Assuming the rotational symmetry of the system, we come to the conclusion that the pair correlations of interest are dependent only on the polar angles, determining the inclinations of the ferroparticle easy axes from the direction of an external magnetic field. The dimer configuration, where two ferroparticles are in close contact along a magnetic field with their easy magnetisation axes aligned, is the most probable. This configuration becomes more pronounced with increasing anisotropy energy, dipolar coupling constant, and external magnetic field.

Magnetic properties of textured ferrocomposite consisting of immobilized superparamagnetic nanoparticles

Wide use of magnetic nanoparticles in modern technologies and biomedical applications requires reliable theoretical models capable of predicting physical properties. Solidification of a ferroparticle suspension under the action of permanent magnetic field allows us to obtain a ferrocomposite, characterized by some orientational texture of the nanoparticle easy magnetization axes. The static magnetic response of this ferrocomposite differs from that of the parent magnetic suspension due to "freezing" of nanoparticle translational and rotational degrees of freedom. Here the superparamagnetic fluctuations of the nanoparticle magnetic moments play a key role in the formation of the ferrocomposite magnetic response depending on the degree of orientational ordering, obtained during synthesis of a ferrocomposite. With the help of statistical mechanics we calculate the magnetization and the initial magnetic susceptibility of the textured ferrocomposite for various temperatures and magnetic field strengths. The easy axis texturing leads to a significant increase of the magnetic properties, and the effect intensifies with the growth of nanoparticle magnetocrystalline anisotropy. Theoretical predictions are supported by Monte Carlo simulations. The obtained results evidence that the texturing of a ferroparticle suspension and transforming it into a textured ferrocomposite are a real way to enhance the magnetic response of these magnetic soft materials.

Arrested dynamics of the dipolar hard sphere model

We report the combined results of molecular dynamics simulations and theoretical calculations concerning various dynamical arrest transitions in a model system representing a dipolar fluid, namely, N (soft core) rigid spheres interacting through a truncated dipole-dipole potential. By exploring different regimes of concentration and temperature, we find three distinct scenarios for the slowing down of the dynamics of the translational and orientational degrees of freedom: at low (η = 0.2) and intermediate (η = 0.4) volume fractions, both dynamics are strongly coupled and become simultaneously arrested upon cooling. At high concentrations (η≥ 0.6), the translational dynamics shows the features of an ordinary glass transition, either by compressing or cooling down the system, but with the orientations remaining ergodic, thus indicating the existence of partially arrested states. In this density regime, but at lower temperatures, the relaxation of the orientational dynamics also freezes. The physical scenario provided by the simulations is discussed and compared against results obtained with the self-consistent generalized Langevin equation theory, and both provide a consistent description of the dynamical arrest transitions in the system. Our results are summarized in an arrested states diagram which qualitatively organizes the simulation data and provides a generic picture of the glass transitions of a dipolar fluid.

Molecular theory of a ferromagnetic nematic liquid crystal

We employ a version of classical density functional theory to study the phase behavior of a simple model liquid crystal in an external field. The uniaxially symmetric molecules have a spherically symmetric core with superimposed orientation-dependent attractions. The interaction between the cores consists of a hard-sphere repulsion plus an isotropic square-well attraction. The anisotropic part of the interaction potential allows for the formation of a uniaxially symmetric nematic phase. The orientation of the molecules couples to an external polar field. The external field is capable of rotating the nematic director n[over ̂] in the x-z plane. The field is also capable of changing the topology of the phase diagram in that it suppresses the phase coexistence between an isotropic liquid and a nematic phase observed in the absence of the field. We study the transition from an unpolar to a polar nematic phase in terms of the orientation-distribution function (odf), nematic and polar order parameters, and components of n[over ̂]. If represented suitably the odf allows us to study orientational changes during the switching process between nonpolar and polar nematic phases. We also give a simple argument that explains why nematic order is lost whereas polar order persists up to the gas-liquid critical point along the coexistence curve. We also discuss the relevance of our theory for future experimental studies.

Mean-field density functional theory of a nanoconfined classical, three-dimensional Heisenberg fluid. II. The interplay between molecular packing and orientational order

As in Paper I of this series of papers [S. M. Cattes et al., J. Chem. Phys. 144, 194704 (2016)], we study a Heisenberg fluid confined to a nanoscopic slit pore with smooth walls. The pore walls can either energetically discriminate specific orientations of the molecules next to them or are indifferent to molecular orientations. Unlike in Paper I, we employ a version of classical density functional theory that allows us to explicitly account for the stratification of the fluid (i.e., the formation of molecular layers) as a consequence of the symmetry-breaking presence of the pore walls. We treat this stratification within the White Bear version (Mark I) of fundamental measure theory. Thus, in this work, we focus on the interplay between local packing of the molecules and orientational features. In particular, we demonstrate why a critical end point can only exist if the pore walls are not energetically discriminating specific molecular orientations. We analyze in detail the positional and orientational order of the confined fluid and show that reorienting molecules across the pore space can be a two-dimensional process. Last but not least, we propose an algorithm based upon a series expansion of Bessel functions of the first kind with which we can solve certain types of integrals in a very efficient manner.

Phase Behavior of Magnetic Nanocolloids of Different Sizes Suspended in an Apolar Solvent

We employ classical density functional theory (DFT) to investigate the phase behavior and composition of binary mixtures; each compound consists of hard spheres of different sizes with superimposed dispersion attraction. In addition to the dispersion attraction, molecules of one component carry an additional three-dimensional magnetic "spin" where the orientation-dependent spin-spin interaction is accounted for by the Heisenberg model. We are treating the excess free energy using a modified mean-field approximation (second virial coefficient) for the orientation-dependent pair correlation function. Depending on the concentration of the magnetic particles, the strength of the spin-spin coupling, and the size ratio of the particles, the model predicts the formation of ordered (polar) phases in addition to the more conventional gas and isotropic liquid phases. Key features of our model are a particle-size dependent shift of the gas-liquid critical point (critical temperature and density) and a change in the width of the phase diagram. In the near-critical region, the latter can be analyzed quantitativly in terms of an effective critical exponent β_eff that may differ from the classical critical exponent [Formula: see text]; the classical value is attained in the immediate vicinity of the critical point as it must. The deviation between β_eff and β can be linked to nontrivial composition effects along the phase boundaries.

Dynamic Design of Spatial Patterns of Colloidal Suspensions

We study the collective dynamics of colloidal suspensions in the presence of a time-dependent potential by means of dynamic density functional theory. We consider a nonlinear diffusion equation for the density and show that spatial patterns emerge from a sinusoidal external potential with a time-dependent wavelength. These patterns are characterized by a sinusoidal density with the average wavelength and a Bessel-function envelope with an induced wavelength that depends only on the amplitude of the temporal oscillations. As a generalization of this result, we propose a design strategy to obtain a family of spatial patterns using time-dependent potentials of practically arbitrary shape.

Perturbation Theory versus Thermodynamic Integration. Beyond a Mean-Field Treatment of Pair Correlations in a Nematic Model Liquid Crystal

The phase behavior and structure of a simple square-well bulk fluid with anisotropic interactions is described in detail. The orientation dependence of the intermolecular interactions allows for the formation of a nematic liquid-crystalline phase in addition to the more conventional isotropic gas and liquid phases. A version of classical density functional theory (DFT) is employed to determine the properties of the model, and comparisons are made with the corresponding data from Monte Carlo (MC) computer simulations in both the grand canonical and canonical ensembles, providing a benchmark to assess the adequacy of the DFT results. A novel element of the DFT approach is the assumption that the structure of the fluid is dominated by intermolecular interactions in the isotropic fluid. A so-called augmented modified mean-field (AMMF) approximation is employed accounting for the influence of anisotropic interactions. The AMMF approximation becomes exact in the limit of vanishing density. We discuss advantages and disadvantages of the AMMF approximation with respect to an accurate description of isotropic and nematic branches of the phase diagram, the degree of orientational order, and orientation-dependent pair correlations. The performance of the AMMF approximations is found to be good in comparison with the MC data; the AMMF approximation has clear advantages with respect to an accurate and more detailed description of the fluid structure. Possible strategies to improve the DFT are discussed.

Wrinkled labyrinths in critical demixing ferrofluid

A thin film of a critical ferrofluid mixture undergoes a sequence of transitions in a magnetic field. First the application of a field induces a critical demixing of the fluid into cylindrical droplets of the minority phase immersed in an extended majority phase. At a second critical field the cylindrical shape is destabilized and transforms into a labyrinth pattern. A third wrinkling transition occurs at even higher field if the liquid has a liquid/air interface. The wrinkling is absent if the droplet has a cover-slide on top. We explain the wrinkling by the wetting behavior of the liquid/air interface that shifts the surface region away from a critical demixing point.

Magnetic field induced modulated phases in a ferrofluid lutidine silicone oil mixture

A mixture of an ester based ferrofluid with silicone oil and 2,6-lutidine is exposed to an external magnetic field. We find a region of composition of the ternary mixture, where weak magnetic fields of the order of a few kA m^-1 induce a modulated phase with a pattern characterized by equilibrium size droplets of the minority phase immersed into the extended majority phase. While the pattern resembles in many ways the pattern of immiscible magnetic fluids, the dependence of the characteristic parameters of the pattern on the magnetic field are completely different than in immiscible fluids. We theoretically explain the pattern formation as a magnetic field induced polymerization of magnetic particles into magnetic chains that goes along with a reduction of the entropy of mixing. This entropy reduction causes the Ostwald ripening of chains into mesoscopic droplets the size of which is limited by repulsive dipolar interactions between the chains.

Mean-field density functional theory of a nanoconfined classical, three-dimensional Heisenberg fluid. I. The role of molecular anchoring

In this work, we employ classical density functional theory (DFT) to investigate for the first time equilibrium properties of a Heisenberg fluid confined to nanoscopic slit pores of variable width. Within DFT pair correlations are treated at modified mean-field level. We consider three types of walls: hard ones, where the fluid-wall potential becomes infinite upon molecular contact but vanishes otherwise, and hard walls with superimposed short-range attraction with and without explicit orientation dependence. To model the distance dependence of the attractions, we employ a Yukawa potential. The orientation dependence is realized through anchoring of molecules at the substrates, i.e., an energetic discrimination of specific molecular orientations. If the walls are hard or attractive without specific anchoring, the results are "quasi-bulk"-like in that they can be linked to a confinement-induced reduction of the bulk mean field. In these cases, the precise nature of the walls is completely irrelevant at coexistence. Only for specific anchoring nontrivial features arise, because then the fluid-wall interaction potential affects the orientation distribution function in a nontrivial way and thus appears explicitly in the Euler-Lagrange equations to be solved for minima of the grand potential of coexisting phases.