Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Supercooling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗
DOE OSTI2026

Explainable machine learning reveals that local structural motifs encode the thermodynamic state across the CuZr metallic glass-forming range

Metallic glasses derive their properties from the statistics of local atomic motifs rather than from long-range order, yet a quantitative, chemistry-specific link between motif populations and the underlying glassy state has remained elusive. In this work we combine large-scale molecular dynamics, Voronoi tessellation, deep neural networks, and SHapley Additive exPlanations (SHAP) to identify which local structural motifs define the glassy state of Cu—Zr metallic glasses. A dataset of 17,180 atomistic configurations spanning ten compositions (Cu 20 Zr 80 –Cu 80 Zr 20 ) and four quench rates (10 9 –10 12 K/s) is used to train a feed-forward neural network that regresses temperature across the 50–2000 K liquid–supercooled–glass range, achieving a mean absolute error of 19.89 K and R 2 = 0.9974, confirming that the local structural state is faithfully encoded in motif-level structure. SHAP analysis then reveals that a tightly coupled near-icosahedral family of motifs (coordination numbers (CN) 11–13, including the full icosahedron 001200 and its single-atom-perturbation sibling 10930) collectively encodes the thermodynamic state of the system across the full glass-forming range. The CN = 11–13 ordered members carry negative SHAP values at high populations, tracking the most deeply-quenched configurations, while 10930 shows the reversed signature consistent with its role as a soft-spot host whose population shrinks as the icosahedral network deepens. The analysis demonstrates that explainable machine learning can isolate the minimal motif vocabulary defining the glassy state and recovers the near-icosahedral building blocks previously identified by data-driven analyses of Cu—Zr. The approach provides a general, chemistry-specific route for characterizing the structural state of disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermal property characterization of phase change materials in building applications: A systematic review of fundamentals, recent progress, and future directions

Phase change materials (PCMs) can reduce building peak loads and enable demand-responsive thermal energy storage (TES), but their deployment depends on reliable measurement and interpretation of thermal properties across laboratory, intermediate, and application scales. Here, this review systematically examines characterization methods, testing protocols, and recent advances for neat PCMs and PCM composites, emphasizing thermal conductivity, enthalpy-related properties (phase change temperature, latent heat, specific heat), and cycling stability. For thermal conductivity, we compare steady-state and transient techniques and note limitations when phase transition and contact resistance affect measurements. For enthalpy–temperature characterization, we discuss differential scanning calorimetry together with intermediate- and bulk-scale methods, including T-history, heat flow meter testing, and three-layer calorimetry (3LC), to generate application-relevant enthalpy–temperature profiles. Cycling stability is organized into four experimental families: thermoelectric–air, fully thermoelectric, water-bath, and in situ chamber approaches, with attention to separating reversible supercooling from true degradation such as phase segregation. We highlight emerging noncontact diagnostics, including infrared thermography and embedded sensing, for spatially resolved validation and multiscale interpretation. Finally, we review the growing use of AI and machine learning for property prediction, inverse characterization from experimental signals, and real-time state estimation in building-integrated TES. Key needs include harmonized protocols, interlaboratory benchmarking, uncertainty reporting, and metadata-rich datasets to accelerate reproducible PCM qualification for grid-flexible buildings.

AI↗
DOE OSTI2025

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Engineering Antifreeze Proteins to Optimally Resist Engulfment by Ice

Antifreeze proteins (AFPs) facilitate the survival of organisms in cold climates by inhibiting the growth and/or recrystallization of ice. To function, AFPs must first bind to ice crystals; bound AFPs must then resist engulfment by using their nonbinding side (NBS) to pin the ice–water interface. Here, we seek to understand how the molecular characteristics of an NBS, such as its ice-phobicity or shape, influence its ability to resist engulfment. By characterizing the free energy barriers that impede the engulfment of model AFPs, we find that the critical supercooling ΔT*, above which an AFP is engulfed, is dictated by an optimal pinning site on the NBS. We further find that the optimal pinning site is determined by an interplay between the contact line perimeter P and a pinning efficiency η, with ΔT* ∝ ηP at the optimal pinning site. For a hemispherical AFP, which displays progressively inward tapering, we find that η increases during engulfment, whereas P decreases; conversely, an NBS with outward tapering can achieve high P, but it suffers from low η. Because the product of η and P determines ΔT*, the inverse correlation between them limits ΔT. To circumvent such limiting behavior, we propose an NBS shape with an outward bulge; by initially tapering outward, a bulged NBS permits higher P, and by subsequently tapering inward, it promotes high η as well. Importantly, we find that ΔT* is enhanced by more than a factor of 2 with an outward bulge of only 1 nm. We also find that the more ice-phobic an NBS is, the more efficiently it pins the ice–water interface, resulting in a higher ΔT. Furthermore, our findings shed light on how the NBS molecular characteristics influence ΔT* and suggest strategies for engineering the NBS to optimally resist engulfment by ice.

Antifreeze↗
DOE OSTI2025

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗
DOE OSTI2025

Atomic Structure, Dynamics, Changes in Chemical Bonding and Semiconductor-Metal Transition in Sb 2 Se 3 : A Remarkable Material for Quantum Networks and Energy Applications

Antimony sesquiselenide has become an outstanding functional material for photovoltaics, energy storage and transformation, memory and photonic applications. Sb 2 Se 3 is one of the most successful emerging solar light absorbers and has also been identified as a highly promising ultralow-loss phase-change material (PCM) for next-generation coherent nanophotonic processors, photonic tensor cores, quantum and neuromorphic networks. Unlike benchmark telluride PCMs, Sb 2 Se 3 features a quasi-one-dimensional (1D) crystalline structure consisting of (Sb 4 Se 6 ) ∞ ribbons, lacks the typical PCM chemical bonding, and undergoes an extended semiconductor-metal transition above the melting point. Consequently, the origin of high optical contrast between crystalline (SET) and amorphous (RESET) logic states remains elusive and presents a significant challenge. Using high-energy X-ray diffraction and Raman spectroscopy over a wide temperature range, supported by first-principles simulations and complemented by thermal, optical and electrical measurements, as well as by 121 Sb-Mossbauer spectroscopy, the quasi-1D network of orthorhombic antimony sesquiselenide was found to undergo significant evolution in amorphous and supercooled Sb 2 Se 3 , leading to lower coordination, shorter interatomic distances and a higher p-electron density on antimony, indicating changes in chemical bonding. The observed novel Sb 2 Se 3 nanocrystalline polymorph, characterized by trigonal antimony coordination and more isolated Sb-Se ribbons, could help reduce multiple trapping defect states in the bandgap, which are typical of orthorhombic Sb 2 Se 3 , thereby enhancing the power-conversion efficiency of photovoltaic devices. Semimetallic and metallic liquid Sb 2 Se 3 exhibit a gradual transformation into a denser 2D and/or 3D network with higher antimony coordination. Localized electron states in the pseudogap are becoming extended, leading to an increase in electronic conductivity σ following the relationship σ ∝ N(E F ) 2 . Liquid Sb 2 Se 3 also appears to be strongly fragile, with a nonmonotonic change in viscosity and higher atomic mobility in the metallic liquid. Furthermore, these results explain extraordinary functionalities of Sb 2 Se 3 for photonic and energy applications.

antimony↗
DOE OSTI2025

Investigating the Relative Roles of INPs and CCN in a Simulated Thunderstorm Using a New Immersion Freezing Algorithm

Microphysical processes in deep convective clouds are sensitive to the number concentrations of cloud condensation nuclei (CCN) and ice nucleating particles (INPs), but the effects of INPs are less studied. Modeling studies investigating the effects of INPs and/or CCN on deep convection typically retain a volume‐dependent raindrop freezing relation. The resulting neglect of aerosol accumulation in raindrops via drop collisions has likely produced unrealistic storm responses to INPs in past studies. To address this deficiency, a new immersion freezing algorithm was developed and embedded in a bulk microphysics scheme that freezes both cloud drops and raindrops using the same immersion freezing INP (IF‐INP) activity spectrum based on measurements. Multiple idealized simulations of a single case of deep convection observed during the Clouds, Aerosols, and Complex Terrain Interactions (CACTI) field campaign were conducted, with microphysical differences produced by independently altering IF‐INP temperature dependencies and CCN number concentrations from their observed values. Surface precipitation in all simulations resulted almost exclusively from riming graupel that melted upon descending to the surface. Rainfall and cold pools were substantially and systematically weakened with increased CCN due to decreased graupel riming rates but were relatively insensitive to variations in the magnitude and slope of IF‐INP spectra due to compensating depletion of supercooled liquid water. These compensating processes were a consequence of the accumulation of IF‐INPs in raindrops, encouraging caution in studying IF‐INP effects upon thunderstorms using traditional volume‐dependent drop freezing relationships.

CCN↗
DOE OSTI2025

Hemispheric Asymmetry of Phase Partition in Mixed‐Phase Clouds Based on Near Global‐Scale Airborne Observations

Mixed-phase clouds contribute to substantial uncertainties in global climate models due to their complex microphysical properties. Former model evaluations almost exclusively rely on satellite observations to assess cloud phase distributions globally. This study investigated mixed-phase cloud properties using near global-scale in situ observation data sets from 14 flight campaigns in combination with collocated output from a global climate model. The Southern Hemisphere (SH) shows significantly higher occurrence frequencies and higher mass fractions of supercooled liquid water than Northern Hemisphere (NH) based on observations at 0.2 and 100 km horizontal scales. Such hemispheric asymmetry is not captured by the model. The model also consistently overestimates liquid water content (LWC) in all cloud phases but shows ice water content (IWC) biases that vary with phase. Key processes contributing to model biases in phase partition can be identified through the combination of evaluation of phase frequency, liquid mass fraction, LWC and IWC.

Yang, Ching An [San Jose State University, CA (Uni↗
DOE OSTI2026

First EarthCARE CPR Observations of Dynamics–Microphysics Coupling in Marine Cold‐Air Outbreak Clouds

Marine cold-air outbreak (MCAO) clouds begin as shallow cloud streets near the ice edge and evolve into broken open-cell convection downstream. Their rapid transitions in cloud structure and microphysics are closely linked to vertical motions, but this interplay has not previously been observed from space. The EarthCARE satellite, carrying the first spaceborne Doppler radar (Cloud Profiling Radar, CPR), now measures hydrometeor sedimentation velocity and vertical air motion. Using CPR observations over the Norwegian and Barents Seas (December 2024–May 2025), we analyze MCAO and non-MCAO clouds. We find that strong updrafts in MCAO clouds coincide with larger supercooled liquid water (SLW) path (LWP). These conditions are accompanied by increases in sedimentation velocity and its vertical gradient, reflecting rapid riming-driven ice growth. In contrast, non-MCAO clouds show lower riming efficiency at similar LWP due to relatively weaker dynamical support. This emphasizes the importance of dynamics in shaping polar cloud structure and microphysics.

Kim, Jiseob [McGill Univ., Montreal, QC (Canada)] ↗
DOE OSTI2026

Quasi‐Steady State Supersaturation: Do High Values Derived From ESCAPE Represent Real High Supersaturations and the Potential for Condensational Invigoration?

Deep convective clouds were intensively sampled during the Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE) with coordinated flights of the NRC Convair-580 and SPEC Learjet. A total of 219 updraft core segments were sampled over coastal Texas and Louisiana under diverse meteorological conditions. Median updraft properties included widths of ∼1 km, velocities of 4.8 m s −1 , droplet number concentrations of ∼400 cm −3 , and liquid water contents of 0.9 g m −3 . The limitations of using the quasi-steady state approximation to derive supersaturations were explored. Supersaturation (S QSS ) estimated from in situ observations under a quasi-steady state assumption averaged 0.4% but occasionally exceeded 2%, with values >1% (high supersaturations) identified as statistical outliers. Two case studies illustrated the conditions linked to high supersaturations. In a storm over the Gulf, median core S QSS reached 2.46% in the developing stage compared to 2.17% in the mature stage under similar thermodynamic conditions. In a storm over coastal Louisiana, S QSS peaked near 11% within a 13.7-m s −1 updraft, accompanied by predominantly supercooled liquid droplets at −13°C and exceptionally low diameter concentrations of 0.29 mm cm −3 . Bootstrap analysis of all sampled cores showed that high supersaturations are most probable in cold and mixed-phase regimes with moderate to strong updrafts and are strongly influenced by vertical velocity and droplet number concentrations. In conclusion, while extreme supersaturations (∼10%) were rare, their occurrence underscores the need for targeted multiplatform observations to resolve their spatiotemporal variability and assess their potential role in deep convective invigoration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Metastability and Ostwald step rule in the crystallisation of diamond and graphite from molten carbon

Abstract Experimental challenges in determining the phase diagram of carbon at temperatures and pressures near the graphite-diamond-liquid triple point are often related to the persistence of metastable crystalline or glassy phases, superheated crystals, or supercooled liquids. A deeper understanding of the crystallisation kinetics of diamond and graphite is crucial for effectively interpreting the outcomes of these experiments. Here, we reveal the microscopic mechanisms of diamond and graphite nucleation from liquid carbon through molecular simulations with first-principles machine learning potentials. Our simulations accurately reproduce the experimental phase diagram of carbon near the triple point and show that liquid carbon crystallises spontaneously upon cooling. Metastable graphite crystallises in the domain of diamond thermodynamic stability at pressures above the triple point. Furthermore, whereas diamond crystallises through a classical nucleation pathway, graphite follows a two-step process in which low-density fluctuations forego ordering. Calculations of the nucleation rates of the two competing phases confirm this result and reveal a manifestation of Ostwald’s step rule, where the strong metastability of graphite hinders the transformation to the stable diamond phase. Our results provide a key to interpreting melting and recrystallisation experiments and shed light on nucleation kinetics in polymorphic materials with deep metastable states.

Science & Technology - Other Topics↗
DOE OSTI2025

Biological enhancement of cloud droplet concentrations observed off East Antarctica

Abstract The impact that biogenic emissions have on aerosol-cloud interactions across the Southern Ocean is poorly quantified. Here we use satellite and ship observations during austral summer to study these interactions. We present observational evidence that biogenic aerosols increase cloud condensation nuclei and cloud droplet number concentrations over the Southern Ocean off East Antarctica, coinciding with very low concentrations of ice-nucleating particles and higher occurrences of supercooled liquid-containing low-level clouds.

Mallet, Marc D.↗
DOE OSTI2025

Ultrathin liquid sheets: water gets in shape for VUV absorption

We present absorption spectra of thin, free-flowing liquid sheets in the vacuum ultraviolet energy range using a gas-squeezed liquid jet. Compared to liquid flow cells, operation without transmission windows eliminates restrictions on the energy range. The temperature of the water sheet is estimated at 0 ± 3 °C, at the verge of the supercooled regime. By adjusting flow conditions in situ, we recorded absorption spectra at water sheet thicknesses ranging from 20 to 50 nm. We show that the absorption spectra of thin jets contain significant contributions from interference effects that need to be deconvoluted from spectral contributions due to the electronic structure. We employ a Fresnel propagation model to model the spectral changes and understand the impact of thickness variations and thin film interference. This opens the door for the investigation of solvation, interface, and similar effects by recording valence band spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

X-ray photoelectron spectroscopy of surfactants on sub-micron aqueous aerosols

Organic molecules at the aqueous interfaces of droplets and aerosols can influence multiphase chemistry and cloud condensation nuclei activity. While surface tension measurements on flat surfaces are commonly used to quantify organic partitioning and structure, extending these measurements to more realistic aerosols that might exist in metastable, supercooled and supersaturated states remains a challenge. Here, we use aerosol velocity map imaging X-ray photoelectron spectroscopy (A-VMI-XPS), a surface-sensitive and in situ technique, to study the partitioning and structure of a model surfactant (n-octyl β-D-thioglucopyranoside, OTG) in submicron aqueous aerosols with a mean radius of 116 nm. After accounting for finite size effects and evaporative cooling, we find the bulk concentrations of OTG in the aerosol are depleted 600-fold relative to macroscopic solutions. At the low temperatures of the experiment (ca. ∼188 K), we obtain a Langmuir equilibrium constant of 700 m3 mol-1, which corresponds to 4.9 m3 mol-1 at room temperature; a value that is consistent with prior literature reports. Analysis of photoelectron angular distributions, peak areas, and secondary electron escape barriers indicates a transition from a disordered surfactant layer at low [OTG] to a more structured layer at high [OTG]. These results link macroscopic surface tension measurements to nanoparticle surfactant behavior and demonstrate the utility of A-VMI-XPS in probing surfactant coating structure on free aerosols, with implications for understanding the impact of organic coatings on cloud formation and atmospheric chemistry.

Deal, Alexandra M↗
DOE OSTI2025

Structure of molten ytterbium aluminum garnet

Rare earth aluminum garnets are important materials in optical, dielectric, and thermal barrier applications. To advance the understanding of their melt processing and glass forming ability, we report the atomic structure of molten Yb 3 Al 5 O 12 over 1770–2630 K, which spans the equilibrium and supercooled liquid regimes. The melt density at T m = 2283 K is 5.50 g cm –3 , measured via silhouette imaging of electrostatically levitated drops over 1010–2420 K. Four separate structure measurements were made with aerodynamically levitated melts using x-ray and neutron diffraction with isotope substitution of Yb ( 172 Yb, 174 Yb, or nat Yb). Empirical potential structure refinement models were developed, which are in excellent agreement with the experiments. Coordination environments for Al–O are predominantly 4- and 5-coordinate, with a mean coordination of n AlO = 4.43(8), while Yb–O environments mostly range from 5- to 8-coordinate, with n YbO = 6.26(8). The cation–oxygen polyhedra are connected primarily by corner-sharing, with edge-sharing constituting up to ~1/3 of the connectivity among polyhedra with Yb or higher-coordinated Al–O. Structurally, the –Al–O– network in molten Yb 3 Al 5 O 12 appears conducive to glass formation: n OAl = 1.85(3), there are 1.86 AlO x –AlO x connections per Al atom (e.g., a mixture of Q 3 and Q 4 units), and the modal ring size is six cations. These characterize a network that is somewhat less constrained compared to SiO 2 glass, yet Yb 3 Al 5 O 12 cannot be quenched into crystal-free glass. Here, aluminum garnet compositions with larger rare earth cations do form glass, so these characterizations help reveal the structural characteristics corresponding to the limit of glass forming ability in rare earth aluminates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The origin of the Stokes–Einstein relation in simple dense liquids

Here, we investigate the origin of the universal relation between structural relaxation and diffusion in simple dense liquids, known as the Stokes–Einstein (SE) relation. The fact that this relation, originally derived from a hydrodynamic model of a macroscopic particle in a viscous medium, can describe the microscopic-scale liquid dynamics still eludes understanding. We introduce a new universal measure of structural relaxation in a system of N identical particles based on an explicit decomposition of the configuration space into N! congruent convex polyhedra. This measure makes it possible to quantify the correlation between two distinct particle configurations in terms of their minimal Euclidean distance, optimized with respect to particle permutations. Using this measure alongside a model of independent random walkers under the single-occupancy constraint, we derive a master equation that quantifies the SE relation. It allows us to demonstrate that the universal relation between structural relaxation and diffusion in simple dense liquids is caused by two conditions: (a) the confinement of the dominant density fluctuations to the first coordination shell, manifested by de Gennes narrowing, and (b) Gaussianity of the diffusion process; the former is shown to be violated in low-density fluids, and the latter is known to be violated in supercooled liquids.

Physics - Condensed matter physics↗
DOE OSTI2026

Deriving reliable nucleation rates from metadynamics simulations: Application to Yukawa fluids

In order to solidify the usefulness of metadynamics in studying nucleation of crystals from supercooled liquids, this work provides a specific procedure to calculate nucleation free energy barriers. After a pedagogical review of the important elements of classical nucleation theory and how metadynamics is used to find nucleation free energy barriers, we explain the benefits of local collective variables over more common global collective variables. We show how a metadynamics free energy barrier must be carefully post-processed so that classical nucleation theory can be applied to calculate nucleation rates. We apply our procedure to a Yukawa plasma and show that a particular physically motivated fit to metadynamics data reproduces low-temperature reference data, justifying the usefulness of metadynamics to predict nucleation rates.

42 ENGINEERING↗