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 in…
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.
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 a…
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 $\rho = N/V$ exceeding the onset density $\rho_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 nucleat…
Perspective: The Asakura Oosawa model: A colloid prototype for bulk and interfacial phase behavior
In many colloidal suspensions, the micrometer-sized particles behave like hard spheres, but when non-adsorbing polymers are added to the solution a depletion attraction (of entropic origin) is created. Since 60 years the Asakura-Oosawa model, which simply describes the polymers as ideal soft spheres, is an archetypical description for the statistical thermodynamics of such systems, accounting for many features of real colloid-polymer mixtures very well. While the fugacity of the polymers (which controls their concentration in the solution) plays a role like inverse temperature, the size ratio of polymer versus colloid radii acts as a control parameter to modify the phase diagram: when this …
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 press…
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, M. Dijkstra, R. Evans & R. Roth, 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, R. Dullens, M. …
The phase coexistence method to obtain surface free energies and nucleation barriers: a brief review
A recently developed method where one analyses the finite size effects associated with liquid–solid phase equilibria including vapour–crystal coexistence is briefly reviewed. It is shown that the e...
Free-energy barriers for crystal nucleation from fluid phases.
Monte Carlo simulations of crystal nuclei coexisting with the fluid phase in thermal equilibrium in finite volumes are presented and analyzed, for fluid densities from dense melts to the vapor. Generalizing the lever-rule for two-phase coexistence in the canonical ensemble to finite volume, "measurements" of the nucleus volume together with the pressure and chemical potential of the surrounding fluid allows to extract the surface free energy of the nucleus. Neither the knowledge of the (in general non-spherical) nucleus shape nor of the angle-dependent interface tension is required for this task. The feasibility of the approach is demonstrated for a variant of the Asakura-Oosawa model for c…
Molecular dynamics simulations in hybrid particle-continuum schemes: Pitfalls and caveats
Heterogeneous multiscale methods (HMM) combine molecular accuracy of particle-based simulations with the computational efficiency of continuum descriptions to model flow in soft matter liquids. In these schemes, molecular simulations typically pose a computational bottleneck, which we investigate in detail in this study. We find that it is preferable to simulate many small systems as opposed to a few large systems, and that a choice of a simple isokinetic thermostat is typically sufficient while thermostats such as Lowe-Andersen allow for simulations at elevated viscosity. We discuss suitable choices for time steps and finite-size effects which arise in the limit of very small simulation bo…
Controlling the wetting properties of the Asakura-Oosawa model and applications to spherical confinement.
We demonstrate for the Asakura-Oosawa model and an extension of this model that uses continuous rather than hard potentials, how wetting properties at walls can be easily controlled. By increasing the interaction range of the repulsive wall potential acting on the colloids (while keeping the polymer-wall interactions constant) polymers begin to substitute colloids at walls and the system can be driven from complete wetting of colloids via partial wetting to complete wetting of polymers. As an application, we discuss the morphology and wetting behavior of colloid-polymer mixtures in spherical confinement. We apply the recently developed 'ensemble switch method' where the Hamiltonian is exten…
Phase transitions and phase equilibria in spherical confinement
Phase transitions in finite systems are rounded and shifted and affected by boundary effects due to the surface of the system. This interplay of finite size and surface effects for fluids confined inside of a sphere of radius $R$ is studied by a phenomenological theory and Monte Carlo simulations of a model for colloid-polymer mixtures. For this system the phase separation in a colloid-rich phase and a polymer-rich phase has been previously studied extensively in the bulk. It is shown that spherical confinement can strongly enhance the miscibility of the mixture. Depending on the wall potentials at the confining surface, the wetting properties of the wall can be controlled, and this interpl…
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 nonadsorbing 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 …