The fluid-fluid interface in a model colloid-polymer mixture: Application of grand canonical Monte Carlo to asymmetric binary mixtures

We present a Monte Carlo method to simulate asymmetric binary mixtures in the grand canonical ensemble. The method is used to study the colloid-polymer model of Asakura and Oosawa. We determine the phase diagram of the fluid-fluid unmixing transition and the interfacial tension, both at high polymer density and close to the critical point. We also present density profiles in the two-phase region. The results are compared to predictions of a recent density functional theory.Comment: 4 pages, 4 figure

Distribution and Implications of Sponge Spicules in Surficial Deposits in Ohio

Author Institution: Agronomy Department, The Ohio State University, Columbus, OhioMicroscopic examination of biogenic opal isolated from the 0.05-0.02-mm total mineral fraction of 12 upland soil profiles indicates that fragments of sponge spicules are minor but ubiquitous constituents of Ohio soils, with major concentrations in the upper 10 to 15 inches of the profile. Quantities range from about 30 to 2000 parts per million biogenic opal or 1 to 65 parts per 10 million parts soil. Spicules are absent or extremely rare in calcareous Wisconsin-age till deposits. Their correlation with horizons high in silt content (50-75%), and their size and depth distribution in landscape positions which preclude an authigenic origin, indicate their aeolian transport from aquatic source areas with other loessial materials. Identification of spicules thus provides direct evidence that these horizons have been derived from loess or loess-till admixtures. This microscopic technique may serve useful for the identification of loess when field or laboratory particle-size analysis yields inconclusive evidence

The Lennard-Jones-Devonshire cell model revisited

We reanalyse the cell theory of Lennard-Jones and Devonshire and find that in addition to the critical point originally reported for the 12-6 potential (and widely quoted in standard textbooks), the model exhibits a further critical point. We show that the latter is actually a more appropriate candidate for liquid-gas criticality than the original critical point.Comment: 5 pages, 3 figures, submitted to Mol. Phy

Free energies of crystalline solids: a lattice-switch Monte Carlo method

We present a method for the direct evaluation of the difference between the free energies of two crystalline structures, of different symmetry. The method rests on a Monte Carlo procedure which allows one to sample along a path, through atomic-displacement-space, leading from one structure to the other by way of an intervening transformation that switches one set of lattice vectors for another. The configurations of both structures can thus be sampled within a single Monte Carlo process, and the difference between their free energies evaluated directly from the ratio of the measured probabilities of each. The method is used to determine the difference between the free energies of the fcc and hcp crystalline phases of a system of hard spheres.Comment: 5 pages Revtex, 3 figure

Wetting of a symmetrical binary fluid mixture on a wall

We study the wetting behaviour of a symmetrical binary fluid below the demixing temperature at a non-selective attractive wall. Although it demixes in the bulk, a sufficiently thin liquid film remains mixed. On approaching liquid/vapour coexistence, however, the thickness of the liquid film increases and it may demix and then wet the substrate. We show that the wetting properties are determined by an interplay of the two length scales related to the density and the composition fluctuations. The problem is analysed within the framework of a generic two component Ginzburg-Landau functional (appropriate for systems with short-ranged interactions). This functional is minimized both numerically and analytically within a piecewise parabolic potential approximation. A number of novel surface transitions are found, including first order demixing and prewetting, continuous demixing, a tricritical point connecting the two regimes, or a critical end point beyond which the prewetting line separates a strongly and a weakly demixed film. Our results are supported by detailed Monte Carlo simulations of a symmetrical binary Lennard-Jones fluid at an attractive wall.Comment: submitted to Phys. Rev.

Self-trapping at the liquid vapor critical point

Experiments suggest that localization via self-trapping plays a central role in the behavior of equilibrated low mass particles in both liquids and in supercritical fluids. In the latter case, the behavior is dominated by the liquid-vapor critical point which is difficult to probe, both experimentally and theoretically. Here, for the first time, we present the results of path-integral computations of the characteristics of a self-trapped particle at the critical point of a Lennard-Jones fluid for a positive particle-atom scattering length. We investigate the influence of the range of the particle-atom interaction on trapping properties, and the pick-off decay rate for the case where the particle is ortho-positronium.Comment: 12 pages, 3 figures, revtex4 preprin

Interfacial tension of the isotropic--nematic interface in suspensions of soft spherocylinders

The isotropic to nematic transition in a system of soft spherocylinders is studied by means of grand canonical Monte Carlo simulations. The probability distribution of the particle density is used to determine the coexistence density of the isotropic and the nematic phases. The distributions are also used to compute the interfacial tension of the isotropic--nematic interface, including an analysis of finite size effects. Our results confirm that the Onsager limit is not recovered until for very large elongation, exceeding at least L/D=40, with L the spherocylinder length and D the diameter. For smaller elongation, we find that the interfacial tension increases with increasing L/D, in agreement with theoretical predictions.Comment: 8 pages, 7 figures, and also 1 tabl

Are critical finite-size scaling functions calculable from knowledge of an appropriate critical exponent?

Critical finite-size scaling functions for the order parameter distribution of the two and three dimensional Ising model are investigated. Within a recently introduced classification theory of phase transitions, the universal part of the critical finite-size scaling functions has been derived by employing a scaling limit that differs from the traditional finite-size scaling limit. In this paper the analytical predictions are compared with Monte Carlo simulations. We find good agreement between the analytical expression and the simulation results. The agreement is consistent with the possibility that the functional form of the critical finite-size scaling function for the order parameter distribution is determined uniquely by only a few universal parameters, most notably the equation of state exponent.Comment: 11 pages postscript, plus 2 separate postscript figures, all as uuencoded gzipped tar file. To appear in J. Phys. A