Antonia Statt

Finite-size effects on liquid-solid phase coexistence and the estimation of crystal nucleation barriers.

A fluid in equilibrium in a finite volume V with particle number N at a density ρ=N/V exceeding the onset density ρ_{f} of freezing may exhibit phase coexistence between a crystalline nucleus and surrounding fluid. Using a method suitable for the estimation of the chemical potential of dense fluids, we obtain the excess free energy due to the surface of the crystalline nucleus. There is neither a need to precisely locate the interface nor to compute the (anisotropic) interfacial tension. As a test case, a soft version of the Asakura-Oosawa model for colloid-polymer mixtures is treated. While our analysis is appropriate for crystal nuclei of arbitrary shape, we find the nucleation barrier to be compatible with a spherical shape and consistent with classical nucleation theory.

Crystal nuclei in melts: a Monte Carlo simulation of a model for attractive colloids

As a model for a suspension of hard-sphere-like colloidal particles where small non-adsorbing dissolved polymers create a depletion attraction, we introduce an effective colloid–colloid potential closely related to the Asakura–Oosawa model, but that does not have any discontinuities. In simulations, this model straightforwardly allows the calculation of the pressure from the virial formula, and the phase transition in the bulk from the liquid to crystalline solid can be accurately located from a study where a stable coexistence of a crystalline slab with a surrounding liquid phase occurs. For this model, crystalline nuclei surrounded by fluid are studied both by identifying the crystal–fluid interface on the particle level (using suitable bond orientational order parameters to distinguish the phases) and by ‘thermodynamic’ means, i.e. the latter method amounts to compute the enhancement of chemical potential and pressure relative to their coexistence values. We show that the chemical potential can be obtained from simulating thick films, where one wall with a rather long-range repulsion is present, since near this wall, the Widom particle insertion method works, exploiting the fact that the chemical potential in the system is homogeneous. Finally, the surface excess free energy of the nucleus is obtained, for a wide range of nuclei volumes. From this method, it is established that classical nucleation theory works, showing that for the present model, the anisotropy of the interface excess free energy of crystals and their resulting non-spherical shape has only a very small effect on the barrier.

Direct observation in 3d of structural crossover in binary hard sphere mixtures

For binary fluid mixtures of spherical particles in which the two species are sufficiently different in size, the dominant wavelength of oscillations of the pair correlation functions is predicted to change from roughly the diameter of the large species to that of the small species along a sharp crossover line in the phase diagram [C. Grodon et al., J. Chem. Phys. 121, 7869 (2004)]. Using particle-resolved colloid experiments in 3d we demonstrate that crossover exists and that its location in the phase diagram is in quantitative agreement with the results of both theory and our Monte-Carlo simulations. In contrast with previous work [J. Baumgartl et al., Phys. Rev. Lett. 98, 198303 (2007)], where a correspondence was drawn between crossover and percolation of both species, in our 3d study we find that structural crossover is unrelated to percolation.

Monte Carlo Simulation of Crystal-Liquid Phase Coexistence

When a crystal nucleus is surrounded by coexisting fluid in a finite volume in thermal equilibrium, the thermodynamic properties of the fluid (density, pressure, chemical potential) are uniquely related to the surface excess free energy of the nucleus. Using a model for weakly attractive soft colloidal particles, it is shown that this surface excess free energy can be determined accurately from Monte Carlo simulations over a wide range of nucleus volumes, and the resulting nucleation barriers are completely independent from the size of the total volume of the system. A necessary ingredient of the analysis, the pressure at phase coexistence in the thermodynamic limit, is obtained from the interface velocity method. Computing the solid-liquid interface excess free energy via the ensemble switch method, a detailed test of classical nucleation theory is possible (assuming the hypothetical spherical nucleus shape, which is not realized when the crystal nucleus is facetted). Consequences for the interpretation of experiments will be briefly discussed.

Estimation of Nucleation Barriers from Simulations of Crystal Nuclei Surrounded by Fluid in Equilibrium

Nucleation rates for homogeneous nucleation are commonly estimated in terms of an Arrhenius law involving the nucleation barrier, written in terms of a competition of the contribution in surface free energy of the nucleus and the free energy gain proportional to the nucleus volume. For crystal nuclei this “classical nucleation theory” is hampered by the problem that the nucleus in general is non spherical, since the interfacial excess free energy depends on the orientation of the interface relative to the crystal axes. This problem can be avoided by analyzing the equilibrium of a crystal nucleus surrounded by fluid in a small simulation box in thermal equilibrium. Estimating the fluid pressure and the chemical potential, as well as the volume of the nucleus, suffices to obtain the nucleation barrier, if the equation of state of the pure phases as well as the coexistence pressure are known. This method is demonstrated to work using a coarse-grained model for colloids with an effective attraction due to small polymers, comparing two choices of the attraction strength.

Investigation of Finite-Size Effects in the Determination of Interfacial Tensions

The interfacial tension between coexisting phases of a material is an important parameter in the description of many phenomena such as crystallization, and even today its accurate measurement remains difficult. We have studied logarithmic finite-size corrections in the determination of the interfacial tension with large scale Monte Carlo simulations, and have identified several novel contributions which not only depend on the ensemble, but also on the type of the applied boundary conditions. We present results for the Lennard-Jones system and the Ising model, as well as for hard spheres, which are particularly challenging. In the future, these findings will contribute to the understanding and determination of highly accurate interfacial properties with computer simulations, and will be used in the study of nucleation of colloidal crystals. As a first application, we compare the Laplace pressure of a crystalline nucleus surrounded by liquid as obtained from simulations with classical nucleation theory.

Phase Separation of Colloid Polymer Mixtures Under Confinement

Colloid polymer mixtures exhibit vapor-liquid like and liquid-solid like phase transitions in bulk suspensions, and are well-suited model systems to explore confinement effects on these phase transitions. Static aspects of these phenomena are studied by large-scale Monte Carlo simulations, including novel “ensemble switch” methods to estimate excess free energies due to confining walls. The kinetics of phase separation is investigated by a Molecular Dynamics method, where hydrodynamic effects due to the solvent are included via the multiparticle collision dynamics method.