Heterogeneous nucleation in the low-barrier regime

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.

Influence on crystal nucleation of an order-disorder transition among the subcritical clusters

Studies of nucleation generally focus on the properties of the critical cluster, but the presence of defects within the crystal lattice means that the population of nuclei necessarily evolve through a distribution of precritical clusters with varying degrees of structural disorder on their way to forming a growing stable crystal. To investigate the role precritical clusters play in nucleation, we develop a simple thermodynamic model for crystal nucleation in terms of cluster size and the degree of cluster order that allows us to alter the work of forming the precritical clusters without affecting the properties of the critical cluster. The steady state and transient nucleation behavior of the system are then studied numerically, for different microscopic ordering kinetics. We find that the model exhibits a generic order-disorder transition in the precritical clusters. Independent of the types of ordering kinetics, increasing the accessibility of disordered precritical clusters decreases both the steady state nucleation rate and the nucleation lag time. Furthermore, the interplay between the free-energy surface and the microscopic ordering kinetics leads to three distinct nucleation pathways.

Stability of bubble-like fluxons in disk-shaped Josephson junctions in the presence of a coaxial dipole current

We investigate analytically and numerically the stability of bubble-like fluxons in disk-shaped heterogeneous Josephson junctions. Using ring solitons as a model of bubble fluxons in the two-dimensional sine-Gordon equation, we show that the insertion of coaxial dipole currents prevents their collapse. We characterize the onset of instability by introducing a single parameter that couples the radius of the bubble fluxon with the properties of the injected current. For different combinations of parameters, we report the formation of stable oscillating bubbles, the emergence of internal modes, and bubble breakup due to internal mode instability. We show that the critical germ depends on the ratio between its radius and the steepness of the wall separating the different phases in the system. If the steepness of the wall is increased (decreased), the critical radius decreases (increases). Our theoretical findings are in good agreement with numerical simulations.

Diffusion in quasi-one-dimensional channels: A small system n, p, T, transition state theory for hopping times

Particles confined to a single file, in a narrow quasi-one-dimensional channel, exhibit a dynamic crossover from single file diffusion to Fickian diffusion as the channel radius increases and the particles begin to pass each other. The long time diffusion coefficient for a system in the crossover regime can be described in terms of a hopping time, which measures the time it takes for a particle to escape the cage formed by its neighbours. In this paper, we develop a transition state theory approach to the calculation of the hopping time, using the small system isobaric-isothermal ensemble to rigorously account for the volume fluctuations associated with the size of the cage. We also describe a Monte Carlo simulation scheme that can be used to calculate the free energy barrier for particle hopping. The theory and simulation method correctly predict the hopping times for a two-dimensional confined ideal gas system and a system of confined hard discs over a range of channel radii, but the method breaks down for wide channels in the hard discs' case, underestimating the height of the hopping barrier due to the neglect of interactions between the small system and its surroundings.

CRYSTAL GROWTH. Crystallization by particle attachment in synthetic, biogenic, and geologic environments

Field and laboratory observations show that crystals commonly form by the addition and attachment of particles that range from multi-ion complexes to fully formed nanoparticles. The particles involved in these nonclassical pathways to crystallization are diverse, in contrast to classical models that consider only the addition of monomeric chemical species. We review progress toward understanding crystal growth by particle-attachment processes and show that multiple pathways result from the interplay of free-energy landscapes and reaction dynamics. Much remains unknown about the fundamental aspects, particularly the relationships between solution structure, interfacial forces, and particle motion. Developing a predictive description that connects molecular details to ensemble behavior will require revisiting long-standing interpretations of crystal formation in synthetic systems, biominerals, and patterns of mineralization in natural environments.

Crystallization of Lennard-Jones nanodroplets: From near melting to deeply supercooled

We carry out molecular dynamics (MD) and Monte Carlo (MC) simulations to characterize nucleation in liquid clusters of 600 Lennard-Jones particles over a broad range of temperatures. We use the formalism of mean first-passage times to determine the rate and find that Classical Nucleation Theory (CNT) predicts the rate quite well, even when employing simple modelling of crystallite shape, chemical potential, surface tension, and particle attachment rate, down to the temperature where the droplet loses metastability and crystallization proceeds through growth-limited nucleation in an unequilibrated liquid. Below this crossover temperature, the nucleation rate is still predicted when MC simulations are used to directly calculate quantities required by CNT. Discrepancy in critical embryo sizes obtained from MD and MC arises when twinned structures with five-fold symmetry provide a competing free energy pathway out of the critical region. We find that crystallization begins with hcp-fcc stacked precritical nuclei and differentiation to various end structures occurs when these embryos become critical. We confirm that using the largest embryo in the system as a reaction coordinate is useful in determining the onset of growth-limited nucleation and show that it gives the same free energy barriers as the full cluster size distribution once the proper reference state is identified. We find that the bulk melting temperature controls the rate, even though the solid-liquid coexistence temperature for the droplet is significantly lower. The value of surface tension that renders close agreement between CNT and direct rate determination is significantly lower than what is expected for the bulk system.

Vapor condensation onto a non-volatile liquid drop

Molecular dynamics simulations of miscible and partially miscible binary Lennard-Jones mixtures are used to study the dynamics and thermodynamics of vapor condensation onto a non-volatile liquid drop in the canonical ensemble. When the system volume is large, the driving force for condensation is low and only a submonolayer of the solvent is adsorbed onto the liquid drop. A small degree of mixing of the solvent phase into the core of the particles occurs for the miscible system. At smaller volumes, complete film formation is observed and the dynamics of film growth are dominated by cluster-cluster coalescence. Mixing into the core of the droplet is also observed for partially miscible systems below an onset volume suggesting the presence of a solubility transition. We also develop a non-volatile liquid drop model, based on the capillarity approximations, that exhibits a solubility transition between small and large drops for partially miscible mixtures and has a hysteresis loop similar to the one observed in the deliquescence of small soluble salt particles. The properties of the model are compared to our simulation results and the model is used to study the formulation of classical nucleation theory for systems with low free energy barriers.