Crystal nucleation in binary hard sphere mixtures: a Monte Carlo simulation study

We present calculations of the nucleation barrier during crystallization in binary hard sphere mixtures under moderate degrees of supercooling using Monte Carlo simulations in the isothermal-isobaric semigrand ensemble in conjunction with an umbrella sampling technique. We study both additive and negatively nonadditive binary hard sphere systems. The solid-fluid phase diagrams of such systems show a rich variety of behavior, ranging from simple spindle shapes to the appearance of azeotropes and eutectics to the appearance of substitutionally ordered solid phase compounds. We investigate the effect of these types of phase behavior upon the nucleation barrier and the structure of the critical nucleus. We find that the underlying phase diagram has a significant effect on the mechanism of crystal nucleation. Our calculations indicate that fractionation of the species upon crystallization increases the difficulty of crystallization of fluid mixtures and in the absence of fractionation (azeotropic conditions) the nucleation barrier is comparable to pure fluids. We also calculate the barrier to nucleation of a substitutionally ordered compound solid. In such systems, which also show solid-solid phase separation, we find that the phase that nucleates is the one whose equilibrium composition is closer to the composition of the fluid phase.     

Thermodynamics of dipolar hard sphere mixtures

The thermodynamic properties of mixtures of hard spheres with imbedded point dipoles are investigated by an extension of the perturbation theory used by Rushbrooke et al. for pure fluids. Equations are presented for the general multicomponent mixture in which all the hard spheres have the same diameter. Numerical calculations are presented of the phase behaviour and excess thermodynamic mixing functions for the special case of the binary mixture in which only one species is polar. A brief discussion is given of the relationship of this model to experimental results for real fluid mixtures and of possible extensions of this work.     

Alternative fundamental measure theory for additive hard sphere mixtures

The purpose of this short paper is to present an alternative fundamental measure theory (FMT) for hard sphere mixtures. Keeping the main features of the original Rosenfeld's FMT [Phys. Rev. Lett. 63, 980 (1989)] and using the dimensional and the low-density limit conditions a new functional is derived incorporating Boublik's multicomponent extension [Mol. Phys. 59, 371 (1986)] of highly accurate Kolafa's equation of state for pure hard spheres. We test the theory for pure hard spheres and hard sphere mixtures near a planar hard wall and compare the results with the original Rosenfeld's FMT and one of its modifications and with new very accurate simulation data. The test reveals an excellent agreement between the results based on the alternative FMT and simulation data for density profile near a contact and some improvement over the original Rosenfeld's FMT and its modification at the contact region.     

Integral equation theory for symmetric nonadditive hard sphere mixtures

An integral equation theory is presented for the pair correlation functions and phase behavior of symmetric nonadditive hard sphere mixtures with hard sphere diameters given by sigma(A)(A)() = sigma(BB) = lambdad and sigma(AB) = d. This mixture exhibits a fluid-fluid phase separation into an A-rich phase and a B-rich phase at high densities. The theory incorporates, into the closure approximation, all terms that can be calculated exactly in the density expansion of the direct correlation functions. We find that the closure approximation developed in this work is accurate for the structure and phase behavior over the entire range of lambda, when compared to computer simulations, and is significantly more accurate than the previous theories.     

Thermodynamics of heterogeneous crystal nucleation in contact and immersion modes

One of the most intriguing problems of heterogeneous crystal nucleation in droplets is its strong enhancement in the contact mode (when the foreign particle is presumably in some kind of contact with the droplet surface) compared to the immersion mode (particle immersed in the droplet). Heterogeneous centers can have different nucleation thresholds when they act in contact or immersion modes. The underlying physical reasons for this enhancement have remained largely unclear. In this paper we present a model for the thermodynamic enhancement of heterogeneous crystal nucleation in the contact mode compared to the immersion one. To determine if and how the surface of a liquid droplet can thermodynamically stimulate its heterogeneous crystallization, we examine crystal nucleation in the immersion and contact modes by deriving and comparing with each other the reversible works of formation of crystal nuclei in these cases. The line tension of a three-phase contact gives rise to additional terms in the formation free energy of a crystal cluster and affects its Wulff (equilibrium) shape. As an illustration, the proposed model is applied to the heterogeneous nucleation of hexagonal ice crystals on generic macroscopic foreign particles in water droplets at T = 253 K. Our results show that the droplet surface does thermodynamically favor the contact mode over the immersion one. Surprisingly, the numerical evaluations suggest that the line tension contribution (from the contact of three water phases (vapor-liquid-crystal)) to this enhancement may be of the same order of magnitude as or even larger than the surface tension contribution.     

Poisson's ratio of the fcc hard sphere crystal at high densities

Elastic constants and the Poisson ratio of the fcc hard-sphere crystalline phases, free of defects and with vacancies, are determined by two Monte Carlo methods: (i) the analysis of the box fluctuations in the constant pressure ensemble with variable box shape (N-P-T) and (ii) by the free-energy differentiation with respect to deformation in the fixed box ensemble (N-V-T). Very good agreement is observed for the extrapolated to the infinitely large system limit results of both the methods. The coefficients of the leading singularities of the elastic constants near close packing are estimated; they are well described by the free volume approximation. Two mechanisms influencing the Poisson ratio are studied. (i) It is shown that at high densities particle motions decrease the Poisson ratio with respect to the static case which corresponds to zero temperature. Simulations performed for systems of soft spheres, interacting through n-inverse-power potentials, r-n, show that the elastic constants of the hard spheres can be obtained in the limit n-->infinity. When T-->0 the elastic constants of the soft spheres tend to those of the static model. (ii) It is also shown that vacancies decrease C11 and C44 and increase C12 and, hence, increase the Poisson ratio with respect to the defect-free state of the system.     

Application of hard convex body and hard sphere equations of state to the critical properties of binary mixtures

The critical properties of hydrocarbon mixtures, perfluorocarbon + hydrocarbon, perfluoromethylcyclohexane + siloxane, acetone + hydrocarbon and polydimethyl siloxane mixtures have been calculated from an equation of state for hard convex bodies and from Guggenheim's equation of state for hard spheres. In general, the results of both equations agree well with experimental data.It appears, however, that taking shape factors into account (by using the hard convex body equation) does not lead to a significant improvement in the agreement between theory and experiment for the critical properties.     

Viscosity of hard sphere suspensions

A simple functional representation of the concentration dependence of the low-shear viscosity eta of hard sphere suspensions is proposed. The representation, which agrees with published literature at all volume fractions phi, has a hitherto-unremarked transition in its functional form at phi approximately 0.42 identical with phi(t). phi(t) is definitely less than the volume fraction 0.49 of the hard sphere melting transition.     

Molecular dynamics simulations of a hard sphere crystal and reaction-like mechanism for homogeneous melting

Molecular dynamics simulations of a hard sphere crystal are performed for volume fractions ranging from solidification point to melting point. A local bond order parameter is chosen to assign a nature, liquid or solid, to a particle. The probability for a liquid or solid particle to change state presents a typical sigmoid shape as the nature of its neighbors changes. Using this property, I propose a reaction-like mechanism and introduce a small number of rate constants. A mean-field approach to melting and a kinetic Monte Carlo algorithm on a lattice are derived from these chemical processes. The results of these models successfully compare with molecular dynamics simulations, proving that the main properties of melting can be captured by a small number of dynamical parameters.     

Crystal nucleation of colloidal hard dumbbells

Using computer simulations, we investigate the homogeneous crystal nucleation in suspensions of colloidal hard dumbbells. The free energy barriers are determined by Monte Carlo simulations using the umbrella sampling technique. We calculate the nucleation rates for the plastic crystal and the aperiodic crystal phase using the kinetic prefactor as determined from event driven molecular dynamics simulations. We find good agreement with the nucleation rates determined from spontaneous nucleation events observed in event driven molecular dynamics simulations within error bars of one order of magnitude. We study the effect of aspect ratio of the dumbbells on the nucleation of plastic and aperiodic crystal phases, and we also determine the structure of the critical nuclei. Moreover, we find that the nucleation of the aligned close-packed crystal structure is strongly suppressed by a high free energy barrier at low supersaturations and slow dynamics at high supersaturations.     

Theory of correlated two-particle activated glassy dynamics: general formulation and heterogeneous structural relaxation in hard sphere fluids

We generalize the nonlinear Langevin equation theory of activated single particle dynamics to describe the correlated motion of two tagged spherical particles in a glass- or gel-forming fluid as a function of their initial separation. The theory is built on the concept of a two-dimensional dynamic free energy surface which quantifies the forces on two particles moving in a cooperative manner. For the hard sphere fluid, above a threshold volume fraction we generically find two relaxation channels corresponding largely, but not exclusively, to a center-of-mass-like displacement and a radial separation of the two tagged particles. The entropic barriers and mean first passage times are computed and found to systematically vary with volume fraction and initial particle separation; both oscillate as a function of the latter in a manner related to the equilibrium pair correlation function. A dynamic correlation length is estimated as the length scale beyond which the two-particle activated dynamics becomes uncorrelated in space and time, and is found to modestly grow with increasing mean relaxation time. The theory is also applied to a simplified model of cage escape, the elementary step of structural relaxation. Predictions for characteristic relaxation times, translation-relaxation decoupling, and stretched-exponential decay of time correlation functions are obtained. A novel mechanism for understanding why strong decoupling emerges in the activated regime, but stretched nonexponential time correlation functions do not change shape as the mean relaxation time grows, is presented and favorably compared with experiment. The theory may serve as a starting point for constructing a predictive model of multiple correlated caging and hopping (forward and backward) events of a pair of tagged particles.     

Microscopic investigations of homogeneous nucleation in charged sphere suspensions

We studied the homogeneous nucleation kinetics of an aqueous suspension of charged colloidal spheres under de-ionized conditions. Samples of equilibrium crystalline structure were shear molten and the metastable melt left to solidify after cessation of shear. At low particle number densities n, corresponding to low metastability of the melt, nucleation was monitored directly via video microscopy. We determined the nucleation rates gamma(t) by counting the number of newly appearing crystals in the observation volume per unit time. Using a suitable discrete adaptation of Avrami's [J. Chem. Phys. 7, 1003 (1939); ibid.8, 212 (1940); ibid.9, 177 (1941)] model for solidification via homogeneous nucleation and subsequent growth, we calculate the remaining free volume VF(t) to obtain the rate densities J(t) = gamma(t)/VF(t). We observe J(t) to rise steeply, display a plateau at a maximum rate density Jmax, and to decrease again. With increased n the plateau duration shrinks while Jmax increases. At low to moderate number densities fully solidified samples were analyzed by microscopy to obtain the grain-size distribution and the average crystallite size angle brackets(L). Under the assumption of stationarity, we obtained the nucleation rate density J(Avr), which increased strongly with increasing n. Interestingly, J(Avr) agrees quantitatively to Jmax and to J(Avr) as obtained previously from scattering data taken on the same sample at large n. Thus, by combination of different methods, reliable nucleation rate densities are now available over roughly one order of magnitude in n and eight orders of magnitude in J.     

Fractal aggregates in protein crystal nucleation

Monte Carlo simulations of homogeneous nucleation for a protein model with an exceedingly short-ranged attractive potential yielded a nonconventional crystal nucleation mechanism, which proceeds by the formation of fractal, low-dimensional aggregates followed by a concurrent collapse and increase of the crystallinity of these aggregates to become compact ordered nuclei. This result corroborates a recently proposed two-step mechanism for protein crystal nucleation from solution.     

Suppression of binary nucleation in amorphous La-Fe-Sb mixtures

The nucleation energy of a series of La(x)Fe(y)Sb(z) modulated elemental reactants was measured as a function of the Fe/Sb ratio over a large composition range while holding the La content constant. The nucleation energy of the ternary compound La(0.5)Fe(4)Sb(12) with the skutterudite crystal structure was found to depend very strongly on the Fe/Sb ratio in the modulated elemental reactant, with a higher nucleation energy as the Fe/Sb ratio is moved away from the 1:3 stoichiometric value. When the results of this study are compared with those from Fe(y)Sb(z) modulated reactants, the addition of lanthanum was found to suppress the nucleation of FeSb(2), thereby broadening the Fe/Sb composition range in which the ternary skutterudite compound La(x)Fe(4)Sb(12) nucleates. This suppression of nucleation of a binary phase on addition of a ternary component to an amorphous intermediate is in agreement with theoretical arguments. The observed suppression of nucleation also provides rational for the observed nucleation of metastable ternary and higher-order compounds from homogeneous amorphous reactants.     

In-plane stacking disorder in polydisperse hard sphere crystals

We demonstrate that in random-stacking hard-sphere colloidal crystals the stacking disorder not only exists in the direction perpendicular to the close-packed layers, but also extends in the lateral direction. The existence of such in-plane stacking disorder is suggested by a recent observation of lateral broadening of the Bragg scattering rods in microradian X-ray diffraction and is further confirmed here by real-space confocal microscopy in two hard-sphere colloidal systems with different relative gravity effects. Due to the in-plane stacking disorder, the hexagonal planes consist of islands with different lateral A, B, and C positions with characteristic line defects in between them. The real-space layer-by-layer stacks of images also reveal the 3-D structure of the defects. The chance zeta to find another line-defect above a line-defect in the layer below turns out to be close to 1/2--independent of relative gravity--which can be explained by the two different stacking options above a defect. The stacking of a few sets of several line defects situated on top of each other turns out to be predominantly FCC-like.     

Stacking in sediments of colloidal hard spheres

We use computer simulations to investigate the crystallization dynamics of sedimenting hard spheres in large systems (hundreds of thousands of particles). We show that slow sedimentation results primarily in face-centered cubic (fcc) stacked crystals, instead of random hexagonal close packed or hexagonal close packed (hcp) crystals. We also find slanted stacking faults, in the fcc regions. However, we attribute the formation of fcc to the free energy difference between fcc and hcp and not to the presence of these slanted stacking faults. Although the free energy difference between hcp and fcc per particle is small (only 10(-3) times the thermal energy), it can become considerable, when multiplied by the number of particles in each domain. The ratio of fcc to hcp obtained from dynamic simulations is in excellent agreement with well-equilibrated Monte Carlo simulations, in which no slanted stacking faults were found. Our results explain a range of experiments on colloids, in which the amount of fcc increases upon lowering the sedimentation rate or decreasing the initial volume fraction.     

Equation of state of the hard sphere liquid

Using the results of Monte Carlo simulation, equations of state of hard sphere liquids are calculated for 106 values of the fill factor η= 0.005–0.530 (step of 0.005). In the region of liquid phase stability the absolute accuracy of about 0.00001–0.00008 is reached. Correctness of the accuracy estimate is discussed. The results obtained are compared with reported equations of state of the hard sphere liquid.     

Homopolymer micelles in heterogeneous solvent mixtures

Amphiphilic homopolymers containing a hydrophilic and a hydrophobic functionality in each monomer unit have been shown to form polar or apolar containers depending on the solvent environment. When presented with a mixture of solvents, these polymeric containers are capable of releasing certain guest molecules. The fundamental mechanism behind these properties is investigated, and the utility of these assemblies in separations has been demonstrated with an example.     

Phase diagram of the adhesive hard sphere fluid

The phase behavior of the Baxter adhesive hard sphere fluid has been determined using specialized Monte Carlo simulations. We give a detailed account of the techniques used and present data for the fluid-fluid coexistence curve as well as parametrized fits for the supercritical equation of state and the percolation threshold. These properties are compared with the existing results of Percus-Yevick theory for this system.