Ferroelectric phase in Stockmayer fluids

Critical Casimir forces in soft matter

We review recent advances in the theoretical, numerical, and experimental studies of critical Casimir forces in soft matter, with particular emphasis on their relevance for the structures of colloidal suspensions and on their dynamics. Distinct from other interactions which act in soft matter, such as electrostatic and van der Waals forces, critical Casimir forces are effective interactions characterised by the possibility to control reversibly their strength via minute temperature changes, while their attractive or repulsive character is conveniently determined via surface treatments or by structuring the involved surfaces. These features make critical Casimir forces excellent candidates for controlling the equilibrium and dynamical properties of individual colloids or colloidal dispersions as well as for possible applications in micro-mechanical systems. In the past 25 years a number of theoretical and experimental studies have been devoted to investigating these forces primarily under thermal equilibrium conditions, while their dynamical and non-equilibrium behaviour is a largely unexplored subject open for future investigations.

Density functional approach to elastic properties of three-dimensional dipole-spring models for magnetic gels

Magnetic gels are composite materials consisting of a polymer matrix and embedded magnetic particles. Those are mechanically coupled to each other, giving rise to the magnetostrictive effects as well as to a controllable overall elasticity responsive to external magnetic fields. Due to their inherent composite and thereby multiscale nature, a theoretical framework bridging different levels of description is indispensable for understanding the magnetomechanical properties of magnetic gels. In this study, we extend a recently developed density functional approach from two spatial dimensions to more realistic three-dimensional systems. Along these lines, we connect a mesoscopic characterization resolving the discrete structure of the magnetic particles to macroscopic continuum parameters of magnetic gels. In particular, we incorporate the long-range nature of the magnetic dipole-dipole interaction and consider the approximate incompressibility of the embedding media and relative rotations with respect to an external magnetic field breaking rotational symmetry. We then probe the shape of the model system in its reference state, confirming the dependence of magnetostrictive effects on the configuration of the magnetic particles and on the shape of the considered sample. Moreover, calculating the elastic and rotational coefficients on the basis of our mesoscopic approach, we examine how the macroscopic types of behavior are related to the mesoscopic properties. Implications for real systems of random particle configurations are also discussed.

Molecular theory of the static dielectric constant of dipolar fluids

The link between the static dielectric constant and the microscopic intermolecular interactions is the Kirkwood g₁ factor, which depends on the orientational structure of the fluid. Over the years, there have been several attempts to provide an accurate description of the orientational structure of dipolar fluids using molecular theories. However, these approaches were either limited to mean-field approximations for the pair correlation function or, more recently, limited to adjusting the orientational dependence to simulation data. Here, we derive a theory for the dielectric constant of dipolar hard-sphere fluids using the augmented modified mean-field approximation. Qualitative agreement is achieved throughout all relevant thermodynamic states, as demonstrated by a comparison with simulation data from the literature. Excellent quantitative agreement can be obtained using a single empirical scaling factor, the physical origin of which is analyzed and accounted for. In order to predict the dielectric constant of the Stockmayer fluid (Lennard-Jones plus dipole potential), we use an adjusted version of the expression for the dipolar hard-sphere fluid. Comparing theoretical predictions with newly generated simulation data, we show that it is possible to obtain excellent agreement with simulation by performing the calculations at a corresponding state using the same scaling factor. Finally, we compare the theoretical orientational structure of the Stockmayer fluid with that obtained from simulations. The simulated structure is calculated following a post-processing methodology that we introduce by deriving an original expression that relates the proposed theory to the histogram of relative dipole angles.

Phase behavior of a lattice-gas model for biaxial nematics

We employ a lattice-gas extension of the Maier-Saupe model with discrete orientation states to study the phase behavior of a statistical model for biaxial nematogenic units in mean-field theory. The phase behavior of the system is investigated in terms of the strength of isotropic interaction between anisotropic objects, as well as the degree of biaxiality and the concentration of those units. We obtain phase diagrams with isotropic phases and stable biaxial and uniaxial nematic structures, various phase coexistences, many types of critical and multicritical behaviors, such as ordinary vapor-liquid critical points, critical end points, and tricritical points, and distinct Landau-like multicritical points. Our results widen the possibilities of relating the phenomenological coefficients of the Landau-de Gennes expansion to microscopic parameters, allowing an improved interpretation of theoretical fittings to experimental data.

Ferrofluids and bio-ferrofluids: looking back and stepping forward

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

Liquid-liquid phase separation in an inhomogeneous ternary colloid-polymer mixture

Suspended colloids are often considered as models for molecules, which are sufficiently big so that they can be observed directly in (light) microscopes and for which the effective interaction among each other can be tailored. The Asakura-Oosawa model of ideal colloid-polymer mixtures captures the idea of tuning the interaction between the colloids via a potential, which possesses a range set by the size of the polymers and an attractive strength characterized by the (reservoir) number density of the polymers, which plays the role of an inverse temperature. The celebrated Asakura-Oosawa depletion potential allows one to recreate the bulk phase diagram of a simple fluid by employing a colloid-polymer mixture. This has been verified in theory, by computer simulations, and via experiments. Here, we study the phase behavior of a confined colloid-polymer mixture with two polymer species. The sizes and densities are chosen such that the resulting bulk phase diagram exhibits a second stable critical point within the framework of the classical density functional theory. Our results suggest that a suitably tuned colloid-polymer mixture can be an interesting model system to study fluids with two critical points.

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.

Self-assembly of rigid magnetic rods consisting of single dipolar beads in two dimensions

Molecular dynamics simulations are used to investigate the structural properties of a two-dimensional ensemble of magnetic rods, which are modeled as aligned single dipolar beads. The obtained self-assembled configurations can be characterized as (1) clusters, (2) percolated, and (3) ordered structures, and their structural properties are investigated in detail. By increasing the aspect ratio of the magnetic rods, we show that the percolation transition is suppressed due to the reduced mobility of the rods in two dimensions. Such a behavior is opposite to the one observed in three dimensions. A magnetic bulk phase is found with local ferromagnetic order and an unusual nonmonotonic behavior of the nematic order is observed.

Smectic phases in ionic liquid crystals

Ionic liquid crystals (ILCs) are anisotropic mesogenic molecules which carry charges and therefore combine properties of liquid crystals, e.g. the formation of mesophases, and of ionic liquids, such as low melting temperatures and tiny triple-point pressures. Previous density functional calculations have revealed that the phase behavior of ILCs is strongly affected by their molecular properties, i.e. their aspect ratio, the loci of the charges, and their interaction strengths. Here, we report new findings concerning the phase behavior of ILCs as obtained by density functional theory and Monte Carlo simulations. The most important result is the occurrence of a novel, wide smectic-A phase [Formula: see text], at low temperature, the layer spacing of which is larger than that of the ordinary high-temperature smectic-A phase [Formula: see text]. Unlike the ordinary smectic S _A phase, the structure of the [Formula: see text] phase consists of alternating layers of particles oriented parallel to the layer normal and oriented perpendicular to it.

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.

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.

Ferroelectric domain formation in discotic liquid crystals: Monte Carlo study on the influence of boundary conditions

The realization of a spontaneous macroscopic ferroelectric order in fluids of anisotropic mesogens is a topic of both fundamental and technological interest. Recently we demonstrated that a system of dipolar achiral disklike ellipsoids can exhibit long-searched ferroelectric liquid crystalline phases of dipolar origin. In the present work, extensive off-lattice Monte Carlo simulations are used to investigate the phase behavior of the system under the influences of the electrostatic boundary conditions that restrict any global polarization. We find that the system develops strongly ferroelectric slablike domains periodically arranged in an antiferroelectric fashion. Exploring the phase behavior at different dipole strengths, we find existence of the ferroelectric nematic and ferroelectric columnar order inside the domains. For higher dipole strengths, a biaxial phase is also obtained with a similar periodic array of ferroelectric slabs of antiparallel polarizations. We have studied the depolarizing effects by using both the Ewald summation and the spherical cutoff techniques. We present and compare the results of the two different approaches of considering the depolarizing effects in this anisotropic system. It is explicitly shown that the domain size increases with the system size as a result of considering a longer range of dipolar interactions. The system exhibits pronounced system size effects for stronger dipolar interactions. The results provide strong evidence to the novel understanding that the dipolar interactions are indeed sufficient to produce long-range ferroelectric order in anisotropic fluids.

Linear and nonlinear magnetic properties of ferrofluids

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

Active dipole clusters: From helical motion to fission

The structure of a finite particle cluster is typically determined by total energy minimization. Here we consider the case where a cluster of soft-sphere dipoles becomes active, i.e., when the individual particles exhibit an additional self-propulsion along their dipole moments. We numerically solve the overdamped equations of motion for soft-sphere dipoles in a solvent. Starting from an initial metastable dipolar cluster, the self-propulsion generates a complex cluster dynamics. The final cluster state has in general a structure widely different to the initial one, the details depend on the model parameters and on the protocol of how the self-propulsion is turned on. The center of mass of the cluster moves on a helical path, the details of which are governed by the initial cluster magnetization. An instantaneous switch to a high self-propulsion leads to fission of the cluster. However, fission does not occur if the self-propulsion is increased slowly to high strengths. Our predictions can be verified through experiments with self-phoretic colloidal Janus particles and for macroscopic self-propelled dipoles in a highly viscous solvent.

Condensation, demixing, and orientational ordering of magnetic colloidal suspensions

In this work we study the phase behavior of magnetic particles suspended in a simple nonmagnetic solvent. Magnetic particles are modelled as spherical particles carrying a three-dimensional, classical Heisenberg spin, whereas solvent molecules are treated as spherically symmetric Lennard-Jones particles. The binary mixture of magnetic particles and solvent is studied within the framework of classical density functional theory (DFT). Within DFT pair correlations are treated at the modified mean-field level at which they are approximated by orientation dependent Mayer f functions. In the absence of an external magnetic field four generic types of phase diagrams are observed depending on the concentration of magnetic particles. In this case we observe liquid-liquid phase coexistence between an orientationally ordered (polarized) and a disordered phase characterized by slightly different concentrations of magnetic particles. Liquid-liquid phase coexistence is suppressed by an external field and vanishes completely if the strength of the field is sufficiently large.

Phase diagram of fluid phases in (3)He-(4)He mixtures

Fluid parts of the phase diagram of (3)He-(4)He mixtures are obtained from a mean-field analysis of a suitable lattice gas model for binary liquid mixtures. The proposed model takes into account the continuous rotational symmetry O(2) of the superfluid degrees of freedom associated with (4)He and includes the occurrence of vacancies. This latter degree of freedom allows the model to exhibit a vapor phase and hence can provide the theoretical framework to describe the experimental conditions for measurements of tricritical Casimir forces in (3)He-(4)He wetting films.

Molecular theory of proper ferroelectricity in bent-core liquid crystals

Antiferro- and ferro-electric ordering has been discovered in orthogonal smectic phases composed of nonchiral bent-core molecules. These systems are the only proper fluid ferroelectrics confirmed experimentally so far. We consider a molecular theory of proper ferroelectric ordering in isotropic, nematic and smectic A phases and conclude that the delicate balance between the tendencies for local parallel and antiparallel ordering of molecular electric and steric dipoles is strongly shifted in restricted geometries. This is a reason why dipolar ordering is more likely to occur within a smectic layer. We derive model interaction potentials for polar bent-core molecules and present the results of the mean-field theory of ferroelectric ordering in the orthogonal smectic phase taking into account also the molecular biaxiality. Order parameter profiles have been calculated numerically and phase diagrams are presented which enable one to analyze the relative importance of dipole-dipole interaction and intermolecular attraction modulated by polar bent-core molecular shape.

Density-functional theory and Monte Carlo simulations of the phase behavior of a simple model liquid crystal

We consider the phase behavior of a simple model of a liquid crystal by means of modified mean-field density-functional theory (MMF DFT) and Monte Carlo simulations in the grand canonical ensemble (GCEMC). The pairwise additive interactions between liquid-crystal molecules are modeled via a Lennard-Jones potential in which the attractive contribution depends on the orientation of the molecules. We derive the form of this orientation dependence through an expansion in terms of rotational invariants. Our MMF DFT predicts two topologically different phase diagrams. At weak to intermediate coupling of the orientation dependent attraction, there is a discontinuous isotropic-nematic liquid-liquid phase transition in addition to the gas-isotropic liquid one. In the limit of strong coupling, the gas-isotropic liquid critical point is suppressed in favor of a fluid- (gas- or isotropic-) nematic phase transition which is always discontinuous. By considering three representative isotherms in parallel GCEMC simulations, we confirm the general topology of the phase diagram predicted by MMF DFT at intermediate coupling strength. From the combined MMF DFT-GCEMC approach, we conclude that the isotropic-nematic phase transition is very weakly first order, thus confirming earlier computer simulation results for the same model [see M. Greschek and M. Schoen, Phys. Rev. E 83, 011704 (2011)].

Ferroelectric order in liquid crystal phases of polar disk-shaped ellipsoids

The demonstration of a spontaneous macroscopic ferroelectric order in liquid phases in the absence of any long range positional order is considered an outstanding problem of both fundamental and technological interest. Recently, we reported that a system of polar achiral disklike ellipsoids can spontaneously exhibit a long searched ferroelectric nematic phase and a ferroelectric columnar phase with strong axial polarization. The major role is played by the dipolar interactions. The model system of interest consists of attractive-repulsive Gay-Berne oblate ellipsoids embedded with two parallel point dipoles positioned symmetrically on the equatorial plane of the ellipsoids. In the present work, we investigate in detail the profound effects of changing the separation between the two symmetrically placed dipoles and the strength of the dipoles upon the existence of different ferroelectric discotic liquid crystal phases via extensive off-lattice N-P-T Monte Carlo simulations. Ferroelectric biaxial phases are exhibited in addition to the uniaxial ferroelectric fluids where the phase biaxiality results from the dipolar interactions. The structures of all the ferroelectric configurations of interest are presented in detail. Simple phase diagrams are determined which include different polar and apolar discotic fluids generated by the system.

Size effect and stability of polarized fluid phases

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

Phase behavior of an amphiphilic fluid

We invoke mean-field density functional theory (DFT) to investigate the phase behavior of an amphiphilic fluid composed of a hard-sphere core plus a superimposed anisometric Lennard-Jones perturbation. The orientation dependence of the interactions consists of a contribution analogous to the interaction potential between a pair of "spins" in the classical, three-dimensional Heisenberg fluid and another one reminiscent of the interaction between (electric or magnetic) point dipoles. At fixed orientation both contributions are short-range in nature decaying as r-6 (r being the separation between the centers of mass of a pair of amphiphiles). Based upon two mean-field-like approximations for the pair correlation function that differ in the degree of sophistication we derive expressions for the phase boundaries between various isotropic and polar phases that we solve numerically by the Newton-Raphson method. For sufficiently strong coupling between the Heisenberg "spins" both mean-field approximations generate three topologically different and generic types of phase diagrams that are observed in agreement with earlier work [see, for example, Tavares et al., Phys. Rev. E 52, 1915 (1995)]. Whereas the dipolar contribution alone is incapable of stabilizing polar phases on account of its short-range nature it is nevertheless important for details of the phase diagram such as location of the gas-isotropic liquid critical point, triple, and tricritical points. By tuning the dipolar coupling constant suitably one may, in fact, switch between topologically different phase diagrams. Employing also Monte Carlo simulations in the isothermal-isobaric ensemble the general topology of the DFT phase diagrams is confirmed.

Monte Carlo simulations of spontaneous ferroelectric order in discotic liquid crystals

The demonstration of a spontaneous macroscopic ferroelectric order in liquid phases in the absence of any long-range positional order is considered as an outstanding problem of great fundamental and technological interest. We report here off-lattice Monte Carlo simulations of a system of polar achiral disklike ellipsoids which spontaneously exhibit a novel ferroelectric nematic phase which is a liquid in three dimensions, considering attractive-repulsive pair interaction suitable for the anisotropic particles. At lower temperature, the ferroelectric nematic phase condenses to a ferroelectric hexagonal columnar fluid with an axial macroscopic polarization. A spontaneous ferroelectric order of dipolar origin is established here for the first time in columnar liquid crystals. Our study demonstrates that simple dipolar interactions are indeed sufficient to produce a class of novel ferroelectric fluids of essential interest. The present work reveals the structure-property relationship of achieving long searched ferroelectric liquid crystal phases and transitions between them, and we hope these findings will help in future development of technologically important fluid ferroelectric materials.

Isotropic-polar phase transitions in an amphiphilic fluid: density functional theory versus computer simulations

We investigate the critical line separating isotropic from polar phases in an amphiphilic bulk fluid by means of density functional theory (DFT) and Monte Carlo (MC) simulations in the isothermal-isobaric ensemble. The intermolecular interactions are described by a Lennard-Jones potential in which the attractive contribution is modified by an orientation-dependent function. The latter consists of two terms: The first one has the orientation dependence of a classical three-dimensional Heisenberg interaction, whereas, the second one has the orientation dependence of a classical dipole-dipole interaction. However, both contributions are short range. Employing DFT together with a modified mean-field (MMF) approximation for the orientation-dependent pair correlation function, we derive an analytical expression for the critical line separating isotropic from polar liquidlike phases. In parallel MC simulations, we locate the line of critical points through an analysis of Binder's second-order cumulant of the polar-order parameter. Comparison with DFT shows that the dipolelike contribution is irrelevant for the isotropic-polar phase transition. As far as the Heisenberg contribution is concerned, the MC data are in semiquantitative agreement with the DFT predictions for sufficiently strong coupling between molecular orientations. For weaker coupling, the variation in the ratio of critical density and temperature ρ(c)/T(c) with the Heisenberg coupling constant ε(H) is underestimated by the MMF treatment. The MC results suggest that this is because ρ(c) increases with decreasing ε(H) such that the assumption on which the MMF approach rests becomes less applicable in the weaker-coupling limit.

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

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

Nonequilibrium condensation and coarsening of field-driven dipolar colloids

In colloidal suspensions, self-organization processes can be easily fueled by external fields. Here we consider monolayers of particles with permanent dipole moments that are driven by rotating external fields. In recent experiments, it has been shown that the particles in such systems self-organize into two-dimensional clusters. Here we report results from a computer simulation study of these pattern forming systems. Specifically, we employ Langevin dynamics simulations, Brownian dynamics simulations that include hydrodynamic interactions, and Wang-Landau Monte Carlo simulations of soft spheres interacting via dipolar potentials. We investigate at which field strengths and frequencies clusters form and explore the influence of hydrodynamic interactions. We also examine the phase behavior of the equilibrium system resulting from a time average of the colloidal interactions in the rotating field. In this way we demonstrate that the clustering described in the driven system arises from a first-order phase transition between a vapor and a condensed phase.

Microscopic and macroscopic theories for the dynamics of polar liquid crystals

We derive and analyze the dynamic equations for polar liquid crystals in two spatial dimensions in the framework of classical dynamical density functional theory (DDFT). Translational density variations, polarization, and quadrupolar order are used as order-parameter fields. The results are critically compared with those obtained using the macroscopic approach of time-dependent Ginzburg-Landau (GL) equations for the analogous order-parameter fields. We demonstrate that, for both the microscopic DDFT and the macroscopic GL approach, the resulting dissipative dynamics can be derived from a dissipation function. We obtain microscopic expressions for all diagonal contributions and for many of the cross-coupling terms emerging from a GL approach. Thus, we establish a bridge between molecular correlations and macroscopic modeling for the dissipative dynamics of polar liquid crystals.

Magnetization of multicomponent ferrofluids

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

Nanodrop of an Ising magnetic fluid on a solid surface

The density functional theory of inhomogeneous simple fluids is extended to an Ising magnetic fluid in contact with a solid surface, which is subjected to an external uniform or nonuniform magnetic field. The system is described by two coupled integral equations regarding the magnetic moment and fluid density distributions. The dependence of the contact angle that a nanodrop makes with the solid surface on the parameters involved in the magnetic interactions between the molecules of fluid and between the molecules of fluid and an external magnetic field is calculated. For the uniform magnetic field, the contact angle increases with increasing magnetic field, approaching an asymptotic value that depends on the strength of the fluid-fluid magnetic interactions. In the nonuniform field generated by a permanent magnet, the contact angle first increases with increasing magnetic field B(M) and then decreases, with the decrease being almost linear for large values of B(M). The obtained results are in qualitative agreement with the experimental data on the contact angle of magnetic drops on a solid surface available in the literature.

Steric effects in a mean-field model for polar nematic liquid crystals

The existence of uniaxial liquid crystals comprising polar molecules, with all the dipoles aligned in a parallel pattern, is classically ruled out. Generally, there are two different avenues to a mean-field theory for liquid crystals: one is based on short-range, repulsive, steric forces, and the other is based on long-range, globally attractive, dispersion forces. Purely polar steric interactions have been shown to have the potential of inducing unexpected orientationally ordered states. In real molecules, anisotropies both in shape and in polarizability coexist; it has been shown that dispersion forces interaction can be combined with hard-core repulsion in a formal theory, based on a steric tensor. Starting from this, we build an interaction Hamiltonian featuring the average electric dipolar energy exchanged between molecules with the same excluded region. Under the assumption that the molecular shape is spheroidal, we propose a mean-field model for polar nematic liquid crystals which can exhibit both uniaxial and biaxial polar phases. By means of a numerical bifurcation analysis, we discuss the stability of the equilibrium against the choice of two model parameters, one describing the degree of molecular shape biaxiality and the other describing the relative orientation of the electric dipole within each molecule. We find only uniaxial stable phases, which are effectively characterized by a single scalar order parameter.

Phase transitions and ordering of confined dipolar fluids

We apply a modified mean-field density functional theory to determine the phase behavior of Stockmayer fluids in slit-like pores formed by two walls with identical substrate potentials. Based on the Carnahan-Starling equation of state, a fundamental-measure theory is employed to incorporate the effects of short-ranged hard-sphere-like correlations while the long-ranged contributions to the fluid interaction potential are treated perturbatively. The liquid-vapor, ferromagnetic-liquid-vapor, and ferromagnetic-liquid-isotropic-liquid first-order phase separations are investigated. The local orientational structure of the anisotropic and inhomogeneous ferromagnetic liquid phase is also studied. We discuss how the phase diagrams are shifted and distorted upon varying the pore width.

Magnetization and susceptibility of ferrofluids

A second-order Taylor series expansion of the free energy functional provides analytical expressions for the magnetic field dependence of the free energy and of the magnetization of ferrofluids, here modeled by dipolar Yukawa interaction potentials. The corresponding hard core dipolar Yukawa reference fluid is studied within the framework of the mean spherical approximation. Our findings for the magnetic and phase equilibrium properties are in quantitative agreement with previously published and new Monte Carlo simulation data.

Phase behavior of the Stockmayer fluid via molecular dynamics simulation

The gas-isotropic liquid-nematic liquid phase behavior of the Stockmayer fluid is studied using molecular dynamics simulation together with a mean field lattice model. We obtain coexistence curves of the Stockmayer fluid over a wide range of dipole strengths, temperatures, and densities, including the transition from the isotropic liquid to the ferroelectric liquid. In our simulations we do not observe the disappearance of the isotropic gas-isotropic liquid coexistence at high dipole strength contrary to earlier findings based on Monte Carlo techniques. Even though the formation of reversible dipole chains strongly affects the location of the critical point, it does not lead to its disappearance. These results are supported by a mean field lattice model which yields good qualitative, and in parts quantitative, agreement with our simulations. In addition, we also investigate the gas-isotropic liquid phase behavior for different polarizabilities.

Dynamic electric polarization of nematic liquid crystals subjected to a shear flow

The effect of coupling between the dipole moment and the orientation is explored based on a relaxation equation for the second rank alignment tensor characterizing the molecular order in liquid crystals and a corresponding equation for the electric polarization. The orientational dynamics leads to a time dependence of the electric polarization. We propose a possibility to measure these effects via the resulting magnetic fields of the magnitude |B| approximately equal to 10(-9)T . Furthermore, the presence of the electric polarization modifies the orientational dynamics as demonstrated in solution phase diagrams.

Equilibrium polymerization and gas-liquid critical behavior in the Stockmayer fluid

We develop a simple theory explaining the dependence of the gas-liquid critical point in the Stockmayer fluid on dipole strength. The theory is based on the Flory-Huggins lattice description for polymer systems in conjunction with a transfer matrix model for isolated chains of reversibly assembled dipolar particles. We find that the shift of the critical point as a function of dipole strength, which originally was found in computer simulation, strongly resembles the critical point shift as a function of chain length in ordinary linear polymer systems. In particular, the decrease of the critical density with increasing dipole strength is a consequence of the existence of reversible chains near criticality. In addition we report simulation results for gas-liquid critical points well above the limiting dipole strength found previously.

Capillary condensation and orientational ordering of confined polar fluids

The phase behavior and the orientational structure of polar model fluids confined to slit pores is investigated by means of density functional theory in a modified mean-field approximation. We focus on fluid states and further assume a uniform number density throughout the pore. Our results for spherical dipolar particles with additional van der Waals-like interactions (Stockmayer fluids) reveal complex fluid-fluid phase behavior involving condensation and first- and second-order isotropic-to-ferroelectric phase transitions, where the ferroelectric ordering occurs parallel to the confining walls. The relative importance of these phase transitions depends on two "tuning" parameters, that is the strength of the dipolar interactions (relative to the isotropic attractive ones) between fluid particles, and on the pore width. In particular, in narrow pores the condensation transition seen in bulk Stockmayer fluids is entirely suppressed. For dipolar hard spheres, on the other hand, the impact of confinement consists in a decrease of the isotropic-to-ferroelectric transition temperatures. We also demonstrate that the local orientational structure is inhomogeneous and anisotropic even in globally isotropic systems, in agreement with computer simulation results.

Phase behavior of a symmetrical binary fluid mixture

We have investigated the phase behavior of a symmetrical binary fluid mixture for the situation where the chemical potentials mu(1) and mu(2) of the two species differ. Attention is focused on the set of interparticle interaction strengths for which, when mu(1)=mu(2), the phase diagram exhibits both a liquid-vapor critical point and a tricritical point. The corresponding phase behavior for the case mu(1) not equalmu(2) is investigated via integral-equation theory calculations within the mean spherical approximation and grand canonical Monte Carlo (GCMC) simulations. We find that two possible subtypes of phase behavior can occur, these being distinguished by the relationship between the triple lines in the full phase diagram in the space of temperature, density, and concentration. We present the detailed form of the phase diagram for both subtypes and compare with the results from GCMC simulations, finding good overall agreement. The scenario via which one subtype evolves into the other is also studied, revealing interesting features.

Topography of phase diagrams in binary fluid mixtures: a mean-field lattice density functional study

We employ mean-field lattice density functional theory to calculate complete phase diagrams of binary fluid mixtures composed of molecules of equal size. We consider asymmetric binary mixtures in which the attraction strength between like molecules of either species differs as well as the attractivity between a pair of unlike molecules. Focusing on the topology of phase diagrams in the space spanned by the thermodynamic fields temperature , (mean) chemical potential , and incremental chemical potential (, are chemical potentials of pure mixture components A and B, respectively), we present an argument which precludes the existence of tricritical points (TCPs) in binary mixtures in general. This is a consequence of a purely geometrical argument based upon an analysis of the number of ways in which coexistence surfaces can be joined in the (Euclidian) space of , , and . However, we show that by the same token, TCPs may exist in cases where the mixture possesses some special symmetry. These latter results are in qualitative agreement with earlier works where, however, only special cuts through the complete phase diagrams were considered so that the important relation between existence of TCPs and symmetry properties of the mixture cannot be fully appreciated.

Type-IV phase behavior in fluids with an internal degree of freedom

We have identified a fourth archetype of phase diagram in binary symmetrical mixtures, which is encountered when the ratio of the interaction between the unlike and the like particles is sufficiently small. This type of phase diagram is characterized by the fact that the lambda line (i.e., the line of the second-order demixing transition) intersects the first-order liquid-vapor curve at densities smaller than the liquid-vapor critical density.