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

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.

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.

Coarse-grained treatment of the self-assembly of colloids suspended in a nematic host phase

The complex interplay of molecular scale effects, nonlinearities in the orientational field and long-range elastic forces makes liquid-crystal physics very challenging. A consistent way to extract information from the microscopic, molecular scale up to the meso- and macroscopic scale is still missing. Here, we develop a hybrid procedure that bridges this gap by combining extensive Monte Carlo (MC) simulations, a local Landau-de Gennes theory, classical density functional theory, and finite-size scaling theory. As a test case to demonstrate the power and validity of our novel approach we study the effective interaction among colloids with Boojum defect topology immersed in a nematic liquid crystal. In particular, at sufficiently small separations colloids attract each other if the angle between their center-of-mass distance vector and the far-field nematic director is about 30°. Using the effective potential in coarse-grained two-dimensional MC simulations we show that self-assembled structures formed by the colloids are in excellent agreement with experimental data.

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.

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)].

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.