0000000000170762
AUTHOR
Fabian Schmitz
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…
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…
Monte Carlo Tests of Nucleation Concepts in the Lattice Gas Model
The conventional theory of homogeneous and heterogeneous nucleation in a supersaturated vapor is tested by Monte Carlo simulations of the lattice gas (Ising) model with nearest-neighbor attractive interactions on the simple cubic lattice. The theory considers the nucleation process as a slow (quasi-static) cluster (droplet) growth over a free energy barrier $\Delta F^*$, constructed in terms of a balance of surface and bulk term of a "critical droplet" of radius $R^*$, implying that the rates of droplet growth and shrinking essentially balance each other for droplet radius $R=R^*$. For heterogeneous nucleation at surfaces, the barrier is reduced by a factor depending on the contact angle. U…
Determination of the origin and magnitude of logarithmic finite-size effects on interfacial tension: Role of interfacial fluctuations and domain breathing
The ensemble-switch method for computing wall excess free energies of condensed matter is extended to estimate the interface free energies between coexisting phases very accurately. By this method, system geometries with linear dimensions $L$ parallel and $L_z$ perpendicular to the interface with various boundary conditions in the canonical or grandcanonical ensemble can be studied. Using two- and three-dimensional Ising models, the nature of the occurring logarithmic finite size corrections is studied. It is found crucial to include interfacial fluctuations due to "domain breathing".
The ensemble switch method and related approaches to obtain interfacial free energies between coexisting phases from simulations: a brief review
The accurate estimation of the excess free energy due to an interface between coexisting phases of a model system by computer simulation often is a challenging task. We review here two methods, whi...
The ensemble switch method for computing interfacial tensions
We present a systematic thermodynamic integration approach to compute interfacial tensions for solid-liquid interfaces, which is based on the ensemble switch method. Applying Monte Carlo simulations and finite-size scaling techniques, we obtain results for hard spheres, which are in agreement with previous computations. The case of solid-liquid interfaces in a variant of the effective Asakura-Oosawa model and of liquid-vapor interfaces in the Lennard-Jones model are discussed as well. We demonstrate that a thorough finite-size analysis of the simulation data is required to obtain precise results for the interfacial tension.
Equilibrium between a Droplet and Surrounding Vapor: A Discussion of Finite Size Effects
In a theoretical description of homogeneous nucleation one frequently assumes an "equilibrium" coexistence of a liquid droplet with surrounding vapor of a density exceeding that of a saturated vapor at bulk vapor-liquid two-phase coexistence. Thereby one ignores the caveat that in the thermodynamic limit, for which the vapor would be called supersaturated, such states will at best be metastable with finite lifetime, and thus not be well-defined within equilibrium statistical mechanics. In contrast, in a system of finite volume stable equilibrium coexistence of droplet and supersaturated vapor at constant total density is perfectly possible, and numerical analysis of equilibrium free energie…
Logarithmic finite-size effects on interfacial free energies: Phenomenological theory and Monte Carlo studies
The computation of interfacial free energies between coexisting phases (e.g.~saturated vapor and liquid) by computer simulation methods is still a challenging problem due to the difficulty of an atomistic identification of an interface, and due to interfacial fluctuations on all length scales. The approach to estimate the interfacial tension from the free energy excess of a system with interfaces relative to corresponding single-phase systems does not suffer from the first problem but still suffers from the latter. Considering $d$-dimensional systems with interfacial area $L^{d-1}$ and linear dimension $L_z$ in the direction perpendicular to the interface, it is argued that the interfacial …