Groh and Dietrich reply

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

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.

Phase diagram of a binary symmetric hard-core Yukawa mixture

We assess the accuracy of the self-consistent Ornstein-Zernike approximation for a binary symmetric hard-core Yukawa mixture by comparison with Monte Carlo simulations of the phase diagrams obtained for different choices of the ratio alpha of the unlike-to-like interactions. In particular, from the results obtained at alpha=0.75 we find evidence for a critical endpoint in contrast to recent studies based on integral equation and hierarchical reference theories. The variation of the phase diagrams with range of the Yukawa potential is investigated.

Phase diagram of symmetric binary fluid mixtures: first-order or second-order demixing

Binary fluid mixtures of 1:1 concentration can demix in a phase transition of first order or of second order. We analyze the two scenarios in density-concentration space and relate them to the structure of the line at which the demixing coexistence surface cuts the liquid-vapor coexistence surface. These scenarios help us to decide between first and second order for a model of a symmetric Lennard-Jones mixture. An optimized reference hypernetted chain integral equation method is employed for calculating the correlation functions and from there the pressure and chemical potentials. We conclude that demixing of a 1:1 mixture is of first order in the whole range of parameters that we have investigated. We did not find a critical point in the 1:1 concentration plane.

Influence of solvent granularity on the effective interaction between charged colloidal suspensions

We study the effect of solvent granularity on the effective force between two charged colloidal particles by computer simulations of the primitive model of strongly asymmetric electrolytes with an explicitly added hard-sphere solvent. Apart from molecular oscillating forces for nearly touching colloids that arise from solvent and counterion layering, the counterions are attracted towards the colloidal surfaces by solvent depletion providing a simple statistical description of hydration. This, in turn, has an important influence on the effective forces for larger distances which are considerably reduced as compared to the prediction based on the primitive model. When these forces are repulsive, the long-distance behavior can be described by an effective Yukawa pair potential with a solvent-renormalized charge. As a function of colloidal volume fraction and added salt concentration, this solvent-renormalized charge behaves qualitatively similar to that obtained via the Poisson-Boltzmann cell model, but there are quantitative differences. For divalent counterions and nanosized colloids, on the other hand, the hydration may lead to overscreened colloids with mutual attraction while the primitive model yields repulsive forces. All these new effects can be accounted for through a solvent-averaged primitive model (SPM) which is obtained from the full model by integrating out the solvent degrees of freedom. The SPM was used to access larger colloidal particles without simulating the solvent explicitly.

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

For ferrofluids that are described by a system of hard spheres interacting via dipolar forces we evaluate the magnetization as a function of the internal magnetic field with a Born-Mayer technique and an expansion in the dipolar coupling strength. Two different approximations are presented for the magnetization, considering different contributions to a series expansion in terms of the volume fraction of the particles and the dipolar coupling strength.

Ferromagnetic liquid thin films under applied field

Theoretical calculations, computer simulations, and experiments indicate the possible existence of a ferromagnetic liquid state, although definitive experimental evidence is lacking. Should such a state exist, demagnetization effects would force a nontrivial magnetization texture. Since liquid droplets are deformable, the droplet shape is coupled with the magnetization texture. In a thin-film geometry in zero applied field, the droplet has a circular shape and a rotating magnetization texture with a point vortex at the center. We calculate the elongation and magnetization texture of such a ferromagnetic thin-film liquid droplet confined between two parallel plates under a weak applied magnetic field. The vortex stretches into a domain wall and exchange forces break the reflection symmetry. This behavior contrasts qualitatively and quantitatively with the elongation of paramagnetic thin films.

Orientational order in model dipolar fluids

Fluids of hard spheres each carrying two parallel point dipoles have been investigated using constant-volume Monte Carlo computer simulations. The results show that both ferroelectric and antiferroelectric fluid phases can be stabilized at high density and low temperature by dipolar interactions alone, if the separation between the dipoles on each sphere is sufficiently large. A simple lattice calculation provides some insight into the balance between dipole energy and orientational entropy which governs the polarization state.