SEARCH · Search NASA

Results for “Supercooled water”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Tethered Balloon System (TBS) Instrument Handbook

The Tethered Balloon System (TBS) is an unmanned aerial system composed of a helium-filled balloon, tether, winch, and sensors. Individual components of the system may change with each flight based on the desired measurements, atmospheric conditions, and flight strategy. The TBS operates within the U.S. Department of Energy (DOE)’s R-2204 Restricted Airspace, which encompasses a two-nautical-mile radius centered on Oliktok Point, Alaska that is segmented by altitude into R-2204 Low (0-1,499’ MSL or 0-457 m MSL) and R-2204 High (1,500-6,999’ MSL or 457-2133 m MSL). The TBS may operate outside of R-2204 if allowed to do so under a Certificate of Authorization from the Federal Aviation Administration.

42 ENGINEERING↗

DOE OSTI2024

Isotope effects in supercooled H 2 O and D 2 O and a corresponding-states-like rescaling of the temperature and pressure

Water shows anomalous properties that are enhanced upon supercooling. The unusual behavior is observed in both H 2 O and D 2 O, however, with different temperature dependences for the two isotopes. It is often noted that comparing the properties of the isotopes at two different temperatures (i.e., a temperature shift) approximately accounts for many of the observations—with a temperature shift of 7.2 K in the temperature of maximum density being the most well-known example. However, the physical justification for such a shift is unclear. Motivated by recent work demonstrating a “corresponding-states-like” rescaling for water properties in three classical water models that all exhibit a liquid–liquid transition and critical point, the applicability of this approach for reconciling the differences in the temperature- and pressure-dependent thermodynamic properties of H 2 O and D 2 O is investigated here. Utilizing previously published data and equations-of-state for H 2 O and D 2 O, we show that the available data and models for these isotopes are consistent with such a low temperature correspondence. These observations provide support for the hypothesis that a liquid–liquid critical point, which is predicted to occur at low temperatures and high pressures, is the origin of many of water’s anomalies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Standard Characterization Techniques for Inorganic Phase Change Materials

Salt hydrates are promising candidates for phase change materials (PCMs) because of their low cost and high energy storage capacity. However, practical applications of salt hydrate based PCMs are limited due to their poor thermal cycling performance caused by phase separation and supercooling. Various techniques, including differential scanning calorimetry (DSC) and transient plane source (TPS) have been used to characterize salt hydrate based PCMs but may not provide full understanding of the degradation mechanism. In this manuscript, salt hydrate and graphite composites are prepared. Temperature history (T-History) method and X-ray diffraction (XRD) are utilized for measuring melting temperature, energy storage, and latent heat degradation of the PCM composites. Highlighting the importance of each technique and how each technique plays an important role in analyzing salt hydrates performance accurately. Temperature history, and X-ray diffraction are also used in conjunction to precisely predict the cause of degradation in salt hydrate PCMs with thermal cycling, as in most literatures the cause of degradation in salt hydrate are attributed to phase segregation, supercooling, and change in stichometry water concentration without any scientific validation.

Kumar, Navin↗

DOE OSTI2026

Observationally-based relations among water cycle parameters during snowfall events in the Upper Colorado River Basin

Cold-season precipitation predictability within the complex terrain of the Upper Colorado River Basin is vital for water resource management across many western states, as the Colorado River serves as a primary source of water for over 40 million people in the Western United States and Mexico. This study uses the remote sensing measurements of such water cycle parameters as the liquid equivalent snowfall rate, S, and accumulation, A, vertically integrated amounts of supercooled cloud liquid and ice expressed as liquid water path(LWP) and ice water path (IWP), respectively, and vertically integrated water vapor IWV.

Matrosov, Sergey [CIRES University of Colorado Bou↗

DOE OSTI2021

The prevalence of precipitation from polar supercooled clouds

Abstract. Supercooled clouds substantially impact polar surface energy budgets, but large-scale models often underestimate their occurrence, which motivates accurately establishing metrics of basic processes. An analysis of long-term measurements at Utqiaġvik, Alaska, and McMurdo Station, Antarctica, combines lidar-validated use of soundings to identify supercooled cloud layers and colocated ground-based profiling radar measurements to quantify cloud base precipitation. We find that more than 85 % (75 %) of sampled supercooled layers are precipitating over the Arctic (Antarctic) site, with more than 75 % (50 %) precipitating continuously to the surface. Such high frequencies can be reconciled with substantially lesser spaceborne estimates by considering differences in radar hydrometeor detection sensitivity. While ice precipitation into supercooled clouds from aloft is common, we also find that the great majority of supercooled cloud layers without ice falling into them are themselves continuously generating precipitation. Such sustained primary ice formation is consistent with continuous activation of immersion-mode ice-nucleating particles (INPs), suggesting that supercooled cloud formation is a principal gateway to ice formation at temperatures greater than ∼-38 ∘C over polar regions. The prevalence of weak precipitation fluxes is also consistent with supercooled cloud longevity and with well-observed and widely simulated case studies. An analysis of colocated microwave radiometer retrievals suggests that weak precipitation fluxes can be nonetheless consequential to moisture budgets for supercooled clouds owing to small liquid water paths. The results here also demonstrate that the observed abundance of mixed-phase clouds can vary substantially with instrument sensitivity and methodology. Finally, we suggest that these ground-based precipitation rate statistics offer valuable guidance for improving the representation of polar cloud processes in large-scale models.

54 ENVIRONMENTAL SCIENCES↗

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 OSTI2020

Thermodynamic and Kinetic Transitions of Liquids in Nanoconfinement

Core principles of chemistry are ubiquitously invoked to shed light on the nature of molecular level interactions in nanoconfined fluids, which play a pivotal role in a wide range of processes in geochemistry, biology, and engineering. A detailed understanding of the physicochemical processes involved in the flow, structural transitions, and freezing or melting behavior of fluids confined within nanometer-sized pores of solid materials is thus of enormous importance for both basic research and technological applications. This Account provides a perspective on new insights into the thermodynamic and kinetic transitions of nanoconfined fluids in their stable and metastable forms. After briefly introducing the unique properties of mesoporous silicas from the SBA, MCM, and FDU families that serve as the confinement matrices, combining highly ordered single and bimodal mesopore architectures with tunable pore sizes in the ~2-15 nm range and narrow size distributions, recent studies on melting/freezing behavior of water confined in these host matrices are reviewed. While differential scanning calorimetry (DSC) reveals a linear relationship between melting point depression and pore size (independent of the pore shape), as predicted by the Gibbs-Thomson relation, variable temperature 2 H wide-line nuclear magnetic resonance (NMR) spectroscopy studies confirm the core-shell model of water and give evidence for a layer-by-layer freezing mechanism, which gives rise to an apparent fragile-to-strong transition in the solidification dynamics. In contrast to the freezing/melting behavior of water, the effect of nanoconfinement on the glass transition of supercooled liquids is nonuniversal and the glass transition temperature T g can either increase or decrease with the dimensionality and extent of confinement. This nonuniversal behavior is exemplified by the two glass-forming molecular liquids, glycerol and ortho-terphenyl (OTP). While glycerol shows an increase in T g and a pronounced slowdown of the rotational dynamics of the constituent molecules due to a change in the molecular packing between the bulk and the confined liquid, OTP displays a linear and confining-media-dependent depression of T g with increased confinement that is strongly influenced by the pore-liquid interface characteristics. This Account concludes with a focus on recent experimental evidence of extreme spatial and dynamical heterogeneity in both freezing and glass transition processes. This discovery was enabled by the unique mesoporous structures of SBA-16 and FDU-5, possessing bimodal architectures with two interconnected pore types of different size and shape (spherical and cylindrical). For the very first time, two melting points for water and two glass transitions for supercooled OTP, corresponding to a specific pore type, were observed. Collectively, these observations strongly suggest a close mechanistic connection between the local fluctuations in the structure and dynamics of nanoconfined liquids. Finally, while the findings reviewed in this Account provide new insights into thermodynamic and kinetic transitions of fluids, there remain many unanswered questions regarding the effects of nanoconfinement on the fundamental properties of fluids, which offer exciting future opportunities in chemical research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Understanding supercooling mechanism in sodium sulfate decahydrate phase-change material

Salt hydrate-based phase-change materials are considered promising for future heat storage applications in residential heating/cooling systems. Smooth phase transition from the liquid to solid phase and vice versa is essential for effective heat exchanger; however, supercooling in salt hydrates delays the onset of liquid–solid phase transition. We investigate the molecular level mechanism of supercooling in sodium sulfate decahydrate (SSD). SSD is a complex salt hydrate whose properties are governed by electrostatic forces that include pure Coulombic interactions as well as hydrogen bonds. Experimentally, we examine the importance of a nucleator in reducing supercooling temperatures. We investigated the effect of various mass concentrations of a borax nucleator on a decrease of supercooling temperatures. Molecular dynamics simulation techniques are used to obtain a basic understanding of supercooling in SSD. We observe that by introducing borax as a nucleator, there is a decrease in the supercooling temperature before nucleation. Our molecular dynamics simulations show that long-range electrostatics between sodium and sulfate ion pairs and that with polar water molecules is responsible for delayed nucleation in SSD that results in supercooling, and also, dynamics of charged molecules slows down. The lack of crystallization leads to amorphous structures in supercooled SSD.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Hydrogel-Salt Hydrate Composite for Highly Stable Heat Energy Storage with Reduced Supercooling

Phase change materials (PCM) have potential for use in thermal energy storage in buildings, medical devices, and water heat pumps. Sodium sulfate hydrate (SSD) is appealing due to its high energy storage capability and affordability. However, SSD has issues including high supercooling (> 15°C) and low thermal cyclic stability. In this study, we introduced an ionic molecular nucleating agent that decreases the supercooling temperature to under 2°C. When this SSD was combined with a hydrogel, it maintained its thermal energy storage capacity for over 100 cycles without any decline. The success is attributed to the polymer confining the SSD crystals, preventing large-scale phase separation and the nucleating agent which resulted in nucleation of many small SSD crystals at small undercoolings rather than a small number of larger crystals. As a proof-of-application, we synthesized this composite at a kg scale and demonstrated its properties in a close to real-world demonstration.

chemical composition, thermodynamics↗

DOE OSTI2025

Supercooling of Alaskan Beetle Larvae as a Winter Survival Strategy

Insects are able to survive subfreezing temperatures by either limiting ice crystal formation in their bodies or through freeze avoidance. Beetle larvae are able to avoid freezing in winter by dehydrating in the fall months and replacing their body water content with high concentrations of glycerol. This enables the body fluid of the insect to supercool, and even vitrify, recovering unharmed when the temperature warms in the spring. Using nondestructive, high‐energy X‐ray synchrotron diffraction experiments, direct insight into how cryopreservation occurs at the atomic level within the beetle larvae has been obtained. The results shed light on the molecular‐level interactions associated with the mechanism responsible for surviving freezing temperatures. The molecular models of severely dehydrated Alaskan beetle larvae, based on glycerol‐water mixtures, yield a total of 4.2 ± 1.2 intermolecular hydrogen bonds per glycerol molecule at 275 K, in good agreement with existing molecular dynamics simulations. Most importantly, they show that if just over half the body fluid content is water, the water clusters are too small to form ice crystals that cause cellular damage.

Benmore, Chris J. [X‐Ray Science Division Advanced↗

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 OSTI2020

A new method for operating a continuous flow diffusion chamber to investigate immersion freezing: assessment and performance study

Glaciation in mixed-phase clouds predominately occurs through the immersion freezing mode where ice nucleating particles (INPs) immersed within supercooled droplets induce nucleation of ice. Currently, model representations of this process are a large source of uncertainty in simulating cloud radiative properties, and to constrain these estimates, continuous flow diffusion chamber (CFDC)-style INP devices are commonly used to assess the immersion freezing efficiencies of INPs. In this study, a new approach was explored to operating such an ice chamber that provides maximum activation of particles without droplet breakthrough and correction factor ambiguity to obtain high-quality INP measurements in a manner that has not been demonstrated as possible previously. The conditioning section of the chamber was maintained at ~ –20 °C and water relative humidity (RH w ) ~ 113 % conditions to maximize the droplet activation, and the droplets were supercooled with an independently temperature-controlled nucleation section at a steady cooling rate (0.5 °C min –1 ) to induce the freezing of droplets and evaporation of unfrozen droplets. The performance of the modified ice chamber was evaluated using four INP species: K-feldspar, illite-NX, Argentinian soil dust, and airborne arable dust that had shown ice nucleation over a wide span of supercooled temperatures. Dry dispersed and size-selected K-feldspar particles were generated in the laboratory. Illite-NX and soil dust particles were sampled during the second phase of the Fifth International Ice Nucleation Workshop (FIN-02) campaign, and airborne arable dust particles were sampled from the aerosol inlet located on the rooftop of the laboratory. The measured ice nucleation efficiencies of model aerosols with a surface active site density (n s ) metric were higher, but mostly agreed within one order of magnitude compared to literature results.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A new method for operating a continuous-flow diffusion chamber to investigate immersion freezing: assessment and performance study

Abstract. Glaciation in mixed-phase clouds predominantly occurs through the immersion-freezing mode where ice-nucleating particles (INPs) immersed within supercooled droplets induce the nucleation of ice. Model representations of this process currently are a large source of uncertainty in simulating cloud radiative properties, so to constrain these estimates, continuous-flow diffusion chamber (CFDC)-style INP devices are commonly used to assess the immersion-freezing efficiencies of INPs. This study explored a new approach to operating such an ice chamber that provides maximum activation of particles without droplet breakthrough and correction factor ambiguity to obtain high-quality INP measurements in a manner that previously had not been demonstrated to be possible. The conditioning section of the chamber was maintained at −20 ∘C and water relative humidity (RHw) conditions of 113 % to maximize the droplet activation, and the droplets were supercooled with an independently temperature-controlled nucleation section at a steady cooling rate (0.5 ∘C min−1) to induce the freezing of droplets and evaporation of unfrozen droplets. The performance of the modified compact ice chamber (MCIC) was evaluated using four INP species: K-feldspar, illite-NX, Argentinian soil dust, and airborne soil dusts from an arable region that had shown ice nucleation over a wide span of supercooled temperatures. Dry-dispersed and size-selected K-feldspar particles were generated in the laboratory. Illite-NX and soil dust particles were sampled during the second phase of the Fifth International Ice Nucleation Workshop (FIN-02) campaign, and airborne soil dust particles were sampled from an ambient aerosol inlet. The measured ice nucleation efficiencies of model aerosols that had a surface active site density (ns) metric were higher but mostly agreed within 1 order of magnitude compared to results reported in the literature.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Inherent-state melting and the onset of glassy dynamics in two-dimensional supercooled liquids

Below the onset temperature T o , the equilibrium relaxation time of most glass-forming liquids exhibits glassy dynamics characterized by a super-Arrhenius temperature dependence. In this supercooled regime, the relaxation dynamics also proceeds through localized elastic excitations corresponding to hopping events between inherent states, i.e., potential-energy-minimizing configurations of the liquid. Despite its importance in distinguishing the supercooled regime from the high-temperature regime, the microscopic origin of T o is not yet known. Here, we construct a theory for the onset temperature in two dimensions and find that an inherent-state melting transition, described by the binding–unbinding transition of dipolar elastic excitations, delineates the supercooled regime from the high-temperature regime. The corresponding melting transition temperature is in good agreement with the onset temperature found in various two-dimensional (2D) atomistic models of glass formers and an experimental binary colloidal system confined to a water–air interface. Additionally, we find the predictions for the renormalized elastic moduli to agree with the experimentally observed values for the latter 2D colloidal system. We further discuss the predictions of our theory on the displacement and density correlations at supercooled conditions, which are consistent with observations of the Mermin–Wagner fluctuations in experiments and molecular simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗

DOE OSTI2025

High-Density, Low-Hysteresis Storage Using Hydrated Salts in Surface-Functionalized Hydrogels (Final Technical Report)

Nearly 70 years ago, Glauber’s salt was identified as a leading phase change material (PCM) in terms of its heat storage density (~2x paraffin), thermal conductivity (~1W/m·K), safety, availability and cost (~$\$$100/ton). However, the complex issues of supercooling and incongruent melting due to phase separation have prevented realization of the promise. The addition of thickeners and nucleating agents such as borax solve these issues but only over few cycles. This work aims to (a) resolve long-standing challenges with Glauber’s salt as a thermal storage material through a unique materials approach, (b) to characterize the new material’s properties that are relevant to performance and (c) to explore its incorporation into commercial water heaters. The materials concept involves encapsulating the salt in custom-designed, large-mesh hydrogels that enable breakthrough advances. Specifically, (1) the choice of mesh size and polymer chemistry control diffusion of salt/water and help to eliminate phase segregation. With the hydrogel itself occupying <10% volume, there is little loss in storage density compared to another encapsulation. (2) Specific nucleation centers that covalently tether to the hydrogel trigger heterogeneous nucleation, eliminating supercooling-associated hysteresis losses. The fact that they are spatially tethered, prevents the loss in performance over multiple freeze/thaw cycles (>100). We report extensive characterization of the hydrogel complex in terms of its storage density, freezing/melting temperature, cycling losses, rheological properties, aging and thermal conductivity. The novel material developed in this work is a significant advancement over the state-of-art. Finally, we investigate its potential as a thermal storage material for commercial/residential water heating and identify scenarios in which its deployment is advantageous.

25 ENERGY STORAGE↗

DOE OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗