Search NASA⌕ Search

SEARCH · Search NASA

Results for “Observing System Simulation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Observed declines in upper ocean phosphate-to-nitrate availability

Climate warming is increasing ocean stratification, which in turn should decrease the nutrient flux to the upper ocean. This may slow marine primary productivity, causing cascading effects throughout food webs. However, observing changes in upper ocean nutrients is challenging because surface concentrations are often below detection limits. We show that the nutricline depth, where nutrient concentrations reach well-detected levels, is tied to productivity and upper ocean nutrient availability. Next, we quantify nutricline depths from a global database of observed vertical nitrate and phosphate profiles to assess contemporary trends in global nutrient availability (1972–2022). We find strong evidence that the P-nutricline (phosphacline) is mostly deepening, especially throughout the southern hemisphere, but the N-nutricline (nitracline) remains mostly stable. Earth System Model (ESM) simulations support the hypothesis that reduced iron stress and increased nitrogen fixation buffer the nitracline, but not phosphacline, against increasing stratification. These contemporary trends are expected to continue in the coming decades, leading to increasing phosphorus but not nitrogen stress for marine phytoplankton, with important ramifications for ocean biogeochemistry and food web dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Universality of scaled particle spectra in ultrarelativistic heavy-ion collisions

We study the transverse momentum spectra of identified particles in ultrarelativistic collisions of large and small collision systems. In order to isolate information contained in the momentum dependence, we propose to scale the spectra by the total particle number and mean transverse momentum, global quantities which are already well studied. We observe an interesting, nearly universal, centrality-independent shape in the scaled spectra, similar to scalings that have been studied previously. This scaling behavior breaks down at large transverse momentum and for very small systems, such as those produced in p − p collisions. We perform hybrid hydrodynamic simulations and show that, in these simulations, a centrality-independent shape is a consequence of an event-by-event independence. Here, our results motivate further theoretical and experimental investigations of the regime of validity of this scaling phenomenon and their physical interpretation at different collision energies and systems.

Collective flow↗
DOE OSTI2025

Simulating Marine Ecosystem Dynamics and Biogeochemical Cycling With Multiple Plankton Functional Types

Current representations of marine ecosystems in Earth System Models are greatly simplified, neglecting key interactions between dynamic food webs, biogeochemistry, and climate change. We use the Marine Biogeochemistry Library code base within the Community Earth System Model 2.2.2 to create an expanded ecosystem model with eight phytoplankton groups and four zooplankton size classes (MARBL-8P4Z). Incorporating more specific plankton types and size classes has the potential to capture a wider range of possible behaviors of the ecosystem, its complex interactions with biogeochemistry, and its feedback to climate change. It also permits stronger observational constraints, including in situ group-specific biomass and various observational estimates of plankton community composition. MARBL-8P4Z broadly captures observed global-scale patterns in biomass and community composition for both phytoplankton and zooplankton, with a good performance in simulating broad biogeochemistry fields. The model shows comparable spatial patterns and magnitudes to the observed picophytoplankton biomass (Prochlorococcus, Synechococcus, picoeukaryotes), and captures the seasonal cycle of mesozooplankton biomass. Picophytoplankton groups and microzooplankton dominate biomass and production in oligotrophic, subtropical regions, while nano-phytoplankton, diatoms and the larger zooplankton groups prevail at higher latitudes and within upwelling zones. The model simulates reasonable energy transfer efficiency through the food web, with tight linkages between the phytoplankton community composition, zooplankton grazing, and carbon export, with the potential to link to fisheries models. Thus, MARBL-8P4Z has the potential to account for key climate-driven ecological shifts in the plankton that will modify ocean biogeochemistry in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Applying Gaussian Process Machine Learning and Modern Probabilistic Programming to Satellite Data to Infer CO 2 Emissions

Satellite data provides essential insights into the spatiotemporal distribution of CO 2 concentrations. However, many atmospheric inverse models fail to adequately incorporate the spatial and temporal correlations inherent in satellite observations and often lack rigorous methods for estimating parameters like spatial length scales. We introduce an inference model that processes the spatiotemporal covariance in satellite data and estimates hyperparameters such as covariance length scales. Our approach uses the Gaussian process (GP) machine learning (ML) and modern probabilistic programming languages (PPLs) to perform atmospheric inversions of emissions from satellite data. We develop a GP ML inversion system based on modern PPLs and the GEOS-Chem chemical transport model, simulating atmospheric CO 2 concentrations corresponding to the Orbiting Carbon Observatory-2/3 (OCO-2/3) data for July 2020. In our supervised learning framework, we treat the GEOS-Chem simulated data set as the target, with predictors derived by scaling the target with sector-specific factors hidden from the GP machine. Our results show that the GP model, combined with GPU-enabled PPLs, effectively retrieves true emission scaling factors and infers noise levels concealed within the data. This suggests that our method could be applied over larger areas with more complex covariance structures, enabling comprehensive analysis of the spatiotemporal patterns observed in OCO-2/3 and similar satellite data sets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unraveling electronic correlations in warm dense quantum plasmas

The study of matter at extreme densities and temperatures has emerged as a highly active frontier at the interface of plasma physics, material science and quantum chemistry with relevance for planetary modeling and inertial confinement fusion. A particular feature of such warm dense matter is the complex interplay of Coulomb interactions, quantum effects, and thermal excitations, making its rigorous theoretical description challenging. Here, we demonstrate how ab initio path integral Monte Carlo simulations allow us to unravel this intricate interplay for the example of strongly compressed beryllium, focusing on two X-ray Thomson scattering data sets obtained at the National Ignition Facility. We find excellent agreement between simulation and experiment with a very high level of consistency between independent observations without the need for any empirical input parameters. Our results call into question previously used chemical models, with important implications for the interpretation of scattering experiments and radiation hydrodynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Specific ion effects on ion transport in charged polymer membranes

Selective separation of like-charged ions is a central challenge in applications such as critical mineral recovery. Electrochemical membrane–based separations offer promising pathways to address this need, but a limited fundamental understanding of ion transport in charged polymer membranes hampers the development of highly selective materials. This study elucidates the role of specific ion effects (SIEs) in ion transport through charged polymer membranes by integrating experimental measurements of ion mobility and in situ ion/ion and ion/water interactions in model charged polymer membranes with results from molecular dynamics simulations on analogous systems. We demonstrate that solvent-mediated ion interactions drive pronounced SIEs within the studied membranes, with the ion softness, i.e., the malleability of ion hydration shells, from hard/soft acid/base (HSAB) theory emerging as a key predictor of transport properties. HSAB theory also explains the observed mechanism of solvent-mediated ion interactions. Our findings offer a mechanistic framework for designing membranes with tailored ion selectivity, potentially enabling efficient separations of chemically similar ions.

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

CLimate Impact: Determining Etiology thRough pAthways (CLDERA)

Climate impacts have broad economic, health, political, and national security ramifications. Societally relevant impacts are typically farther downstream, are the product of multiple interacting processes, and can arise over small regions and timeframes because their sources are short-term and localized. Short-term forcings (as can be seen in volcanic eruptions, climatic tipping points (e.g., the collapse of rainforests or the disappearance of sea ice), or in increasingly plausible climate interventions) fundamentally possess low signal-to-noise and could benefit from accounting for the multiple conditional processes through which a downstream impact arises. Under the Grand Challenge LDRD CLDERA (CLimate impacts: Discovering Etiology thRough pAthways), we have developed tools to enable downstream impact attribution from geographically and temporally localized source forcings in the climate. CLDERA developed methods that can distinguish how a localized source drives the climate system to respond with particular impacts. The how is embodied in pathways – the spatio-temporally evolving chain of physical processes that connects a source to a series of increasingly distant impacts. Novel analytic methods in pursuit of downstream impact attribution were developed and demonstrated on simulations and observations of the 1991 eruption of Mt. Pinatubo in the Philippines. As described within this report we have • developed stratospheric expertise and aerosol modeling capabilities in E3SM, • created original methods to detect and model pathways from source-to-impact, and • advanced climate attribution through novel methods, cases, and approaches. Further, CLDERA developed a tiered verification process consisting of controlled datasets to prototype, verify, and refine the original method development. CLDERA increased Sandia’s footprint in the climate analytics community and developed new climate collaborations whilst also creating a cadre of climate analysts at Sandia. The products from CLDERA have been extensive with a total of 9 journal articles published, 12 articles submitted and under review, and an additional 8 articles in preparation. We have produced 1750 simulated years and developed 9 code-bases. This report details these accomplishments and serves as a summary of the work completed during the CLDERA Grand Challenge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pulse duration dependent effects of ultrafast laser induced damage on a 1030 nm multi-layer dielectric mirror for high repetition rate, high average power laser systems

High repetition rate, high peak, and average power laser systems are crucial for next-generation particle accelerators, inertial confinement fusion, and secondary particle sources. These applications demand durable laser optics, particularly interference coatings on optics lasting millions of shots at high fluence. This study focuses on designing, testing, and simulating multi-layer dielectric (MLD) mirrors for pulse durations of 260 fs, 77 fs, and 25 fs at 1030 nm wavelength and 45-degree incidence angle with p -polarization. S-on-1 laser-induced damage thresholds (LIDT) for varying pulse numbers were determined, with single-shot LIDT values of 0.98 Jcm -2 , 1.63 Jcm -2 , and 2.3 Jcm -2 for 25 fs, 77 fs, and 260 fs respectively. A strong correlation between blister shape and local fluence was observed, implying that the layer expansion in a blister depends on local fluence. We have also examined mechanisms responsible for laser-induced stress generation and energy release rates in blister formation. Damage mechanisms are further explored by finite-difference time-domain (FDTD) simulations, incorporating Keldysh strong field ionization, whose predictions were in excellent agreement with the onset of damage determined experimentally. These findings offer insights for enhancing MLD coating technology, promising more efficient and resilient laser systems for diverse scientific and industrial applications.

Noor, Mohamed Yaseen (ORCID:000000021036644X)↗
DOE OSTI2024

Development of the ARM Lagrangian Large-Scale Forcing Data (ARMLAGTRAJ) Value-Added Product Based on the lagtraj Framework

The Atmospheric Radiation Measurement (ARM) large-scale forcing data developed based on the constrained variational analysis (VARANAL) value-added product (VAP) (Zhang and Lin 1997, Zhang et al. 2001, Xie et al. 2004, Tang et al. 2019) has been widely used for single-column models (SCMs), cloud-resolving models (CRMs), and large-eddy simulation models (LESs) to understand and improve physical processes in models. Recently, the U.S. Department of Energy (DOE) ARM user facility conducted several major field campaigns using ship-based moving observational platforms. For example, the Marine ARM GPCI Investigation of Clouds (MAGIC) field campaign focused on the role of subtropical marine-boundary layer (MBL) clouds, and the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign aimed to improve understanding of the coupled climate systems in the Arctic. Observations from moving platforms are critical to provide a comprehensive characterization of coupled-system processes associated with all stages of the cloud and/or sea-ice life cycle. Traditional ARM large-scale forcing data have been developed at fixed locations. They need to be extended to include these moving platforms to address data needs for ship-based field campaigns or to support LES modeling in a Lagrangian framework. With these considerations in mind, we develop ARM-type Lagrangian large-scale forcing data sets based on the lagtraj framework (Boeing et al. 2020) with notable enhancements in generating forcings that are more suitable for ARM field campaigns. The lagtraj is a novel tool that generates forcings for LES and SCM simulation in both Lagrangian and Eulerian perspective. This technical report focuses on the major changes we performed on the lagtraj algorithm and provides an overview of the ARM Lagrangian Large-Scale Forcing Data (ARMLAGTRAJ) value-added products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Performance of SK-Gd’s Upgraded Real-time Supernova Monitoring System

Among multimessenger observations of the next Galactic core-collapse supernova, Super-Kamiokande (SK) plays a critical role in detecting the emitted supernova neutrinos, determining the direction to the supernova (SN), and notifying the astronomical community of these observations in advance of the optical signal. In 2022, SK has increased the gadolinium dissolved in its water target (SK-Gd) and has achieved a Gd concentration of 0.033%, resulting in enhanced neutron detection capability, which in turn enables more accurate determination of the supernova direction. Accordingly, SK-Gd’s real-time supernova monitoring system has been upgraded. SK_SN Notice, a warning system that works together with this monitoring system, was released on 2021 December 13, and is available through GCN Notices. When the monitoring system detects an SN-like burst of events, SK_SN Notice will automatically distribute an alarm with the reconstructed direction to the supernova candidate within a few minutes. In this paper, we present a systematic study of SK-Gd’s response to a simulated Galactic SN. Assuming a supernova situated at 10 kpc, neutrino fluxes from six supernova models are used to characterize SK-Gd’s pointing accuracy using the same tools as the online monitoring system. The pointing accuracy is found to vary from 3° to 7° depending on the models. However, if the supernova is closer than 10 kpc, SK_SN Notice can issue an alarm with three-degree accuracy, which will benefit follow-up observations by optical telescopes with large fields of view.

Core-collapse supernovae↗
DOE OSTI2025

Decoupled few-femtosecond phase transitions in vanadium dioxide

Abstract The nature of the insulator-to-metal phase transition in vanadium dioxide (VO 2 ) is one of the longest-standing problems in condensed-matter physics. Ultrafast spectroscopy has long promised to determine whether the transition is primarily driven by the electronic or structural degree of freedom, but measurements to date have been stymied by their sensitivity to only one of these components and/or their limited temporal resolution. Here we use ultra-broadband few-femtosecond pump-probe spectroscopy to resolve the electronic and structural phase transitions in VO 2 at their fundamental time scales. Our experiments show that the system transforms into a bad-metallic phase within 10 fs after photoexcitation, but requires another 100 fs to complete the transition, during which we observe electronic oscillations and a partial re-opening of the bandgap, signalling a transient semi-metallic state. Comparisons with tensor-network simulations and density-functional theory calculations show these features result from an unexpectedly fast structural transition, in which the vanadium dimers separate and untwist with two different timescales. Our results resolve the structural and electronic nature of the light-induced phase transition in VO 2 and establish ultra-broadband few-femtosecond spectroscopy as a powerful tool for studying quantum materials out of equilibrium.

Science & Technology - Other Topics↗
DOE OSTI2024

eReaxFF force field development for BaZr 0.8 Y 0.2 O 3-δ solid oxide electrolysis cells applications

The use of solid-oxide materials in electrocatalysis applications, especially in hydrogen-evolution reactions, is promising. However, further improvements are warranted to overcome the fundamental bottlenecks to enhancing the performance of solid-oxide electrolysis cells (SOECs), which is directly linked to the more-refined fundamental understanding of complex physical and chemical phenomena and mass exchanges that take place at the surfaces and in the bulk of electrocatalysis materials. Here, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium, BaZr 0.8 Y 0.2 O 3-δ (BZY20) to enable a systematic, large-length-scale, and longer-timescale atomistic simulation of solid-oxide electrocatalysis for hydrogen generation. All parameters for the eReaxFF were optimized to reproduce quantum-mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, electron localization, water adsorption, water splitting, and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed both zero-voltage (excess electrons absent) and non-zero-voltage (excess electrons present) molecular dynamics simulations to observe water adsorption, water splitting, proton migration, oxygen-vacancy migrations, and eventual hydrogen-production reactions. Based on investigations offered in the present study, we conclude that the eReaxFF force field-based approach can enable computationally efficient simulations for electron conductivity, electron leakage, and other non-zero-voltage effects on the solid oxide materials using the explicit-electron concept. Moreover, we demonstrate how the eReaxFF force field-based atomistic-simulation approach can enhance our understanding of processes in SOEC applications and potentially other renewable-energy applications.

08 HYDROGEN↗
DOE OSTI2026

Neural-Network Inverse Design of SRF Cavities and Transmons for Bosonic Quantum Computation

Three-dimensional superconducting radio-frequency (SRF) cavities provide exceptionally long-lived electromagnetic modes and, when coupled to nonlinear elements such as transmon qubits, become promising architectures for bosonic quantum information processing. The inverse design of such systems, i.e., recovering device geometries that produce specified electromagnetic and coupling targets, is generally a one-to-many problem. The qubit-cavity coupling strength depends sensitively on both the transmon geometry and its position within the cavity's electromagnetic field. As these systems scale up and their design parameter spaces grow, the cost of conventional iterative simulation becomes prohibitive. We present two deep neural network (DNN) approaches that address this inverse-design problem at complementary levels of the design stack. The first proposes SRF cavity geometries that produce target cavity observables. The second proposes transmon qubit designs that produce target qubit-cavity parameters - the coupling rate, qubit frequency, and anharmonicity $(g, ν_q, α)$. The recovered candidate designs match the targets to within ~5% (cavity) and ~2% (transmon), confirmed by end-to-end re-simulation. Both approaches map desired device behavior directly to candidate designs, a fast alternative to the iterative simulation studies usually required.

Yaker, Joseph [Fermilab; Northwestern U.]↗
DOE OSTI2025

Environmental Conditions Affecting Global Mesoscale Convective System Occurrence

Abstract The ERA5 environments of mesoscale convective systems (MCSs), tracked from satellite observations, are assessed over a 20-yr period. The use of a large set of MCS tracks allows us to robustly test the sensitivity of the results to factors such as region, latitude, and diurnal cycle. We aim to provide novel information on environments of observed MCSs for assessments of global atmospheric models and to improve their ability to simulate MCSs. Statistical analysis of all tracked MCSs is performed in two complementary ways. First, we investigate the environments when an MCS has occurred at different spatial scales before and after MCS formation. Several environmental variables are found to show marked changes before MCS initiation, particularly over land. The vertically integrated moisture flux convergence shows a robust signal across different regions and when considering MCS initiation diurnal cycle. We also found spatial scale dependence of the environments between 200 and 500 km, providing new evidence of a natural length scale for use with MCS parameterization. In the second analysis, the likelihood of MCS occurrence for given environmental conditions is evaluated, by considering all environments and determining the probability of being in an MCS core or shield region. These are compared to analogous non-MCS environments, allowing discrimination between conditions suitable for MCS and non-MCS occurrence. Three environmental variables are found to be useful predictors of MCS occurrence: total column water vapor, midlevel relative humidity, and total column moisture flux convergence. Such relations could be used as trigger conditions for the parameterization of MCSs, thereby strengthening the dependence of the MCS scheme on the environment. Significance Statement Large storm systems called mesoscale convective systems form across Earth. These are collections of thunderstorms, with associated high-level clouds that produce substantial, lighter rainfall and modulate Earth’s energy balance. They produce hazardous weather conditions, such as floods and high winds, and are responsible for a high percentage of rainfall in many regions globally. We investigate the environmental conditions under which they form, so that we can understand the spatial extent of the environment which is important for their formation, and also where and when the effects of these storms might be felt. The novel information generated here should help improve the representation of these storms in weather and climate models, improving the prediction of rainfall, thunderclouds, and high-level clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Facile Tensile Testing Platform for In Situ Transmission Electron Microscopy of Nanomaterials

In situ tensile testing using transmission electron microscopy (TEM) is a powerful technique to probe structure-property relationships of materials at the atomic scale. In this work, a facile tensile testing platform for in situ characterization of materials inside a transmission electron microscope is demonstrated. The platform consists of: 1) a commercially available, flexible, electron-transparent substrate (e.g., TEM grid) integrated with a conventional tensile testing holder, and 2) a finite element simulation providing quantification of specimen-applied strain. The flexible substrate (carbon support film of the TEM grid) mitigates strain concentrations usually found in free-standing films and enables in situ straining experiments to be performed on materials that cannot undergo localized thinning or focused ion beam lift-out. The finite element simulation enables direct correlation of holder displacement with sample strain, providing upper and lower bounds of expected strain across the substrate. The tensile testing platform is validated for three disparate material systems: sputtered gold-palladium, few-layer transferred tungsten disulfide, and electrodeposited lithium, by measuring lattice strain from experimentally recorded electron diffraction data. The results show good agreement between experiment and simulation, providing confidence in the ability to transfer strain from holder to sample and relate TEM crystal structural observations with material mechanical properties.

2D materials↗
DOE OSTI2024

Interfacial Ice Density Fluctuations Inform Surface Ice-Philicity

The propensity of a surface to nucleate ice or bind to ice is governed by its ice-philicity─its relative preference for ice over liquid water. However, the relationship between the features of a surface and its ice-philicity is not well understood, and for surfaces with chemical or topographical heterogeneity, such as proteins, their ice-philicity is not even well-defined. In the analogous problem of surface hydrophobicity, it has been shown that hydrophobic surfaces display enhanced low water-density (vapor-like) fluctuations in their vicinity. To interrogate whether enhanced ice-like fluctuations are similarly observed near ice-philic surfaces, here we use molecular simulations and enhanced sampling techniques. Using a family of model surfaces for which the wetting coefficient, k , has previously been characterized, we show that the free energy of observing rare interfacial ice-density fluctuations decreases monotonically with increasing k . By utilizing this connection, we investigate a set of fcc systems and find that the (110) surface is more ice-philic than the (111) or (100) surfaces. By additionally analyzing the structure of interfacial ice, we find that all surfaces prefer to bind to the basal plane of ice, and the topographical complementarity of the (110) surface to the basal plane explains its higher ice-philicity. Using enhanced interfacial ice-like fluctuations as a measure of surface ice-philicity, we then characterize the ice-philicity of chemically heterogeneous and topologically complex systems. In particular, we study the spruce budworm antifreeze protein (sbwAFP), which binds to ice using a known ice-binding site (IBS) and resists engulfment using nonbinding sites of the protein (NBSs). We find that the IBS displays enhanced interfacial ice-density fluctuations and is therefore more ice-philic than the two NBSs studied. We also find the two NBSs are similarly ice-phobic. By establishing a connection between interfacial ice-like fluctuations and surface ice-philicity, our findings thus provide a way to characterize the ice-philicity of heterogeneous surfaces.

crystals↗