Monte Carlo tests of nucleation concepts in the lattice gas model

Low-temperature nucleation rate calculations using the N-Fold way

We present a numerical study to determine nucleation rates for magnetization reversal within the Ising model (lattice gas model) in the low-temperature regime, a domain less explored in previous research. To achieve this, we implemented the N-Fold way algorithm, a well-established method for low-temperature simulations, alongside a novel, highly efficient cluster identification algorithm. Our method can access nucleation rates up to 50 orders of magnitude lower than previously reported results. We examine three cases: homogeneous pure system, system with static impurities, and system with mobile impurities, where impurities are defined as sites with zero interactions with neighboring spins (the spin value of impurities is set to 0). Classical nucleation theory holds across the entire temperature range studied in the paper, for both the homogeneous system and the static impurity case. However, in the case of mobile impurities, the umbrella sampling technique appears ineffective at low mobility values. These findings provide valuable insights into nucleation phenomena at low temperatures, contributing to theoretical and experimental understanding.

Mapping the influence of impurity interaction energy on nucleation in a lattice-gas model of solute precipitation

We study nucleation in the two dimensional Ising lattice-gas model of solute precipitation in the presence of randomly placed static and dynamic impurities. Impurity-solute and impurity-solvent interaction energies are varied whilst keeping other interaction energies fixed. In the case of static impurities, we observe a monotonic decrease in the nucleation rate when the difference between impurity-solute and impurity-solvent interaction energies is increased. The nucleation rate saturates to a minimum value with increasing interaction energy difference when the impurity density is low. However the nucleation rate does not saturate for high impurity densities. Similar behaviour is observed with dynamic impurities both at low and high densities. We explore a broad range of both symmetric and anti-symmetric interactions with impurities and map the regime for which the impurities act as a surfactant, decreasing the surface energy of the nucleating phase. We also characterise different nucleation regimes observed at different values of interaction energy. These include additional regimes where impurities play the role of inert-spectators, bulk-stabilizers or cluster together to create heterogeneous nucleation sites for solute clusters to form.

Two-step nucleation in confined geometry: Phase diagram of finite particles on a lattice gas model

We use a degenerated Ising model to describe nucleation and crystallization from solution in a confined two-component system. The free energy is calculated using metadynamics simulation with coordination numbers as the reaction coordinates. We deploy nudged elastic band simulation to determine the minimum energy path and give properties of the crystallization path. In this confined system, depletion effects, which could also be caused by slow material transport in the solution, prevent the post-critical cluster from further growth, and the crystalline state would only be stable at larger cluster sizes. Fluctuation of the higher coupling strength of the crystalline state enables further growth until the crystalline cluster is in equilibrium with the solvent, and this way, a second barrier is crossed. From the parameters and setup, we find necessary conditions for the occurrence of two-step nucleation in our system. These findings can be adapted to real systems as biomineralization, colloidal crystallization, and the solidification of metals.

Condensation and Crystal Nucleation in a Lattice Gas with a Realistic Phase Diagram

We reconsider model II of Orban et al. (J. Chem. Phys. 1968, 49, 1778−1783), a two-dimensional lattice-gas system featuring a crystalline phase and two distinct fluid phases (liquid and vapor). In this system, a particle prevents other particles from occupying sites up to third neighbors on the square lattice, while attracting (with decreasing strength) particles sitting at fourth- or fifth-neighbor sites. To make the model more realistic, we assume a finite repulsion at third-neighbor distance, with the result that a second crystalline phase appears at higher pressures. However, the similarity with real-world substances is only partial: Upon closer inspection, the alleged liquid−vapor transition turns out to be a continuous (albeit sharp) crossover, even near the putative triple point. Closer to the standard picture is instead the freezing transition, as we show by computing the free-energy barrier relative to crystal nucleation from the “liquid”.

Nucleation rate in the two dimensional Ising model in the presence of random impurities

Nucleation phenomena are ubiquitous in nature and the presence of impurities in every real and experimental system is unavoidable. Yet numerical studies of nucleation are nearly always conducted for entirely pure systems. We have studied the behaviour of the droplet free energy in two dimensional Ising model in the presence of randomly positioned static and dynamic impurities. We have shown that both the free energy barrier height and critical nucleus size monotonically decreases with increasing the impurity density for the static case. We have compared the nucleation rates obtained from the Classical Nucleation Theory and the Forward Flux Sampling method for different densities of the static impurities. The results show good agreement. In the case of dynamic impurities, we observe preferential occupancy of the impurities at the boundary positions of the nucleus when the temperature is low. This further boosts enhancement of the nucleation rate due to lowering of the effective interfacial free energy.

Droplet condensation in the lattice gas with density functional theory

A density functional for the lattice gas with next-neighbor attractions (Ising model) from fundamental measure theory is applied to the problem of droplet states in three-dimensional, finite systems. The density functional is constructed via an auxiliary model with hard lattice gas particles and lattice polymers to incorporate the attractions. Similar to previous simulation studies, the sequence of droplets changing to cylinders and to planar slabs is found upon increasing the average density ρ[over ¯] in the system. Owing to the discreteness of the lattice, additional effects in the state curve for the chemical potential μ(ρ[over ¯]) are seen upon lowering the temperature away from the critical temperature [oscillations in μ(ρ[over ¯]) in the slab portion and spiky undulations in μ(ρ[over ¯]) in the cylinder portion as well as an undulatory behavior of the radius of the surface of tension R_{s} in the droplet region]. This behavior in the cylinder and droplet region is related to washed-out layering transitions at the surface of liquid cylinders and droplets. The analysis of the large-radius behavior of the surface tension γ(R_{s}) gave a dominant contribution ∝1/R_{s}^{2}, although the consistency of γ(R_{s}) with the asymptotic behavior of the radius-dependent Tolman length seems to suggest a weak logarithmic contribution ∝lnR_{s}/R_{s}^{2} in γ(R_{s}). The coefficient of this logarithmic term is smaller than a universal value derived with field-theoretic methods.

Weak scaling of the contact distance between two fluctuating interfaces with system size

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ^{*} depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ^{*} scales weakly with system size, i.e., the dependence of ℓ^{*} on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim_{L→∞}(ℓ^{*}/L)=0. In particular, for (2+1)-dimensional interfaces ℓ^{*} is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. This weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.

Nucleation work, surface tension, and Gibbs-Tolman length for nucleus of any size

In the framework of the Gibbs approach to nucleation thermodynamics, expressions are derived for the nucleation work, nucleus size, surface tension, and Gibbs-Tolman length in homogeneous single-component nucleation at a fixed temperature. These expressions are in terms of the experimentally controlled overpressure of the nucleating phase and are valid for the entire overpressure range, i.e., for nucleus of any size. Analysis of available data for bubble and droplet nucleation in Lennard-Jones fluid shows that the theory describes well the data by means of a single free parameter, the Gibbs-Tolman length of the planar liquid/vapor interface. It is found that this length is about one-tenth of the Lennard-Jones molecular-diameter parameter and that it is positive for the bubble nucleus and negative for the droplet nucleus. In a sufficiently narrow temperature range, the nucleation work, nucleus radius, scaled surface tension, and Gibbs-Tolman length are apparently universal functions of scaled overpressure.

Do the contact angle and line tension of surface-attached droplets depend on the radius of curvature?

Results from Monte Carlo simulations of wall-attached droplets in the three-dimensional Ising lattice gas model and in a symmetric binary Lennard-Jones fluid, confined by antisymmetric walls, are analyzed, with the aim to estimate the dependence of the contact angle [Formula: see text] on the droplet radius [Formula: see text] of curvature. Sphere-cap shape of the wall-attached droplets is assumed throughout. An approach, based purely on 'thermodynamic' observables, e.g. chemical potential, excess density due to the droplet, etc, is used, to avoid ambiguities in the decision which particles belong (or do not belong, respectively) to the droplet. It is found that the results are compatible with a variation [Formula: see text], [Formula: see text] being the contact angle in the thermodynamic limit ([Formula: see text]). The possibility to use such results to estimate the excess free energy related to the contact line of the droplet, namely the line tension, at the wall, is discussed. Various problems that hamper this approach and were not fully recognized in previous attempts to extract the line tension are identified. It is also found that the dependence of wall tensions on the difference of chemical potential of the droplet from that at the bulk coexistence provides effectively a change of the contact angle of similar magnitude. The simulation approach yields precise estimates for the excess density due to wall-attached droplets and the corresponding free energy excess, relative to a system without a droplet at the same chemical potential. It is shown that this information suffices to estimate nucleation barriers, not affected by ambiguities on droplet shape, contact angle and line tension.

Heterogeneous nucleation of a droplet pinned at a chemically inhomogeneous substrate: A simulation study of the two-dimensional Ising case

Heterogeneous nucleation is studied by Monte Carlo simulations and phenomenological theory, using the two-dimensional lattice gas model with suitable boundary fields. A chemical inhomogeneity of length b at one boundary favors the liquid phase, while elsewhere the vapor is favored. Switching on the bulk field H_b favoring the liquid, nucleation and growth of the liquid phase starting from the region of the chemical inhomogeneity are analyzed. Three regimes occur: for small fields, H_b<H_b^crit, the critical droplet radius is so large that a critical droplet having the contact angle θ_c required by Young's equation in the region of the chemical inhomogeneity does not yet "fit" there since the baseline length of the circle-cut sphere droplet would exceed b. For H_b^crit<H_b<H_b^*, such droplets fit inside the inhomogeneity and are indeed found in simulations with large enough observation times, but these droplets remain pinned to the chemical inhomogeneity when their baseline has grown to the length b. Assuming that these pinned droplets have a circle cut shape and effective contact angles θ_eff in the regime θ_c < θ_eff < π/2, the density excess due to these droplets can be predicted and is found to be in reasonable agreement with the simulation results. On general grounds, one can predict that the effective contact angle θ_eff and the excess density of the droplets, scaled by b, are functions of the product bH_b but do not depend on both variables separately. Since the free energy barrier for the "depinning" of the droplet (i.e., growth of θ_eff to π - θ_c) vanishes when θ_eff approaches π/2, in practice only angles θ_eff up to about θ_eff^max≃70° were observed. For larger fields (H_b>H_b^*), the droplets nucleated at the chemical inhomogeneity grow to the full system size. While the relaxation time for the growth scales as τ_G∝H_b^-1, the nucleation time τ_N scales as lnτ_N∝H_b^-1. However, the prefactor in the latter relation, as evaluated for our simulations results, is not in accord with an extension of the Volmer-Turnbull theory to two-dimensions, when the theoretical contact angle θ_c is used.

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 us to extract the surface free energy of the nucleus. Neither the knowledge of the (in general nonspherical) 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 colloid-polymer mixtures, which form face-centered cubic colloidal crystals. For a polymer to colloid size ratio of 0.15, the colloid packing fraction in the fluid phase can be varied from melt values to zero by the variation of an effective attractive potential between the colloids. It is found that the approximation of spherical crystal nuclei often underestimates actual nucleation barriers significantly. Nucleation barriers are found to scale as ΔF^{*}=(4π/3)^{1/3}γ[over ¯](V^{*})^{2/3}+const with the nucleus volume V^{*}, and the effective surface tension γ[over ¯] that accounts implicitly for the nonspherical shape can be precisely estimated.

Thermodynamic formulation of the barrier for heterogeneous pinned nucleation: Implication to the crossover scenarios associated with barrierless and homogeneous nucleation

The effect of contact line pinning on nucleation is reported using continuum thermodynamics. Based on the principle of the free-energy maximization, closed-form expressions in the dimensionless form for the free-energy of the three-phase metastable system and the thermodynamic barrier are formulated with respect to the system geometry and the substrate wettability. The condition of maximality limits the dynamic contact angle within the cluster-phase-phobic regime. The dimensionless nucleation barrier or the potency factor can be divided into two components related to the system geometry and the pinning effect. Depending on the relative value of the equilibrium and the critical dynamic contact angle, the contact line pinning can either have favorable or adverse effects. Associated pinning-depinning transition can also lead to the crossovers related to barrierless and homogeneous nucleation. Contact line tension is found to have a considerable effect during these transitional scenarios. Complete wetting transition associated with barrierless nucleation can take place due to the presence of tensile (negative) line tension. On the other hand, complete drying transition related to homogeneous nucleation can occur when line tension is compressive (positive) in nature. The pinning has a favorable effect only when the substrate wettability is within the cluster-phase-philic regime. There can be favorable, adverse, or no pinning effects when the substrate wettability is within the cluster-phase-phobic regime. Although the contact line is pinned, the minimum value of the potency factor is obtained when equilibrium and dynamic contact angles are equal.

Exploring different regimes in finite-size scaling of the droplet condensation-evaporation transition

We present a finite-size scaling analysis of the droplet condensation-evaporation transition of a lattice gas (in two and three dimensions) and a Lennard-Jones gas (in three dimensions) at fixed density. Parallel multicanonical simulations allow sampling of the required system sizes with precise equilibrium estimates. In the limit of large systems, we verify the theoretical leading-order scaling prediction for both the transition temperature and the finite-size rounding. In addition, we present an emerging intermediate scaling regime, consistent in all considered cases and with similar recent observations for polymer aggregation. While the intermediate regime locally may show a different effective scaling, we show that it is a gradual crossover to the large-system scaling behavior by including empirical higher-order corrections. This implies that care has to be taken when considering scaling ranges, possibly leading to completely wrong predictions for the thermodynamic limit. In this study, we consider a crossing of the phase boundary orthogonal to the usual fixed temperature studies. We show that this is an equivalent approach and, under certain conditions, may show smaller finite-size corrections.

A design equation for low dosage additives that accelerate nucleation

Additives are used to control nucleation in many natural and industrial environments. However, the mechanisms by which additives inhibit or accelerate solute precipitate nucleation are not well understood. We propose an equation that predicts changes in nucleation barriers based on the adsorption properties and concentrations of trace additives. The equation shows that nucleant efficacy depends on the product of an adsorption equilibrium constant and the reduction in interfacial tension. Moreover, the two factors that determine the potency of additives are related to each other, suggesting that assays of just one property might facilitate additive design. We test the design equation for a Potts lattice gas model with surfactant-like additives in addition to solutes and solvents.

Anisotropic interfacial tension, contact angles, and line tensions: a graphics-processing-unit-based Monte Carlo study of the Ising model

As a generic example for crystals where the crystal-fluid interface tension depends on the orientation of the interface relative to the crystal lattice axes, the nearest-neighbor Ising model on the simple cubic lattice is studied over a wide temperature range, both above and below the roughening transition temperature. Using a thin-film geometry L(x)×L(y)×L(z) with periodic boundary conditions along the z axis and two free L(x)×L(y) surfaces at which opposing surface fields ±H(1) act, under conditions of partial wetting, a single planar interface inclined under a contact angle θ<π/2 relative to the yz plane is stabilized. In the y direction, a generalization of the antiperiodic boundary condition is used that maintains the translational invariance in the y direction despite the inhomogeneity of the magnetization distribution in this system. This geometry allows a simultaneous study of the angle-dependent interface tension, the contact angle, and the line tension (which depends on the contact angle, and on temperature). All these quantities are extracted from suitable thermodynamic integration procedures. In order to keep finite-size effects as well as statistical errors small enough, rather large lattice sizes (of the order of 46 million sites) were found to be necessary, and the availability of very efficient code implementation of graphics processing units was crucial for the feasibility of this study.

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 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 fluctuations cause logarithmic finite-size effects of order ln(L)/L(d-1) and order ln(L(z))/L(d-1), in addition to regular corrections (with leading-order const/L(d-1)). A phenomenological theory predicts that the prefactors of the logarithmic terms are universal (but depend on the applied boundary conditions and the considered statistical ensemble). The physical origin of these corrections are the translational entropy of the interface as a whole, "domain breathing" (coupling of interfacial fluctuations to the bulk order parameter fluctuations of the coexisting domains), and capillary waves. Using a new variant of the ensemble switch method, interfacial tensions are found from Monte Carlo simulations of d = 2 and d = 3 Ising models and a Lennard-Jones fluid. The simulation results are fully consistent with the theoretical predictions.

Dynamic metastability in the two-dimensional Potts ferromagnet

We investigate the nonequilibrium dynamics of the two-dimensional (2D) Potts model on the square lattice after a quench below the discontinuous transition point. By means of numerical simulations of systems with q=12, 24, and 48, we observe the onset of a stationary regime below the temperature-driven transition, in a temperature interval decreasing with the system size and increasing with q. These results obtained dynamically agree with those obtained from the analytical continuation of the free energy [J. L. Meunier and A. Morel, Eur. Phys. J. B 13, 341 (2000)], from which metastability in the 2D Potts model results to be a finite-size effect.