SEARCH · Search NASA

Results for “Computational Icing”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

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 OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗

DOE OSTI2023

Two-layer cathode architecture for high-energy density and high-power density solid state batteries

Solid state batteries with high-energy density and high-power density require the development of thick and energy dense cathodes. Structured cathode electrodes with a double-layer configuration were enabled using a freeze tape casting technique. A bottom dense layer was utilized to enhance the energy density whereas a top porous layer with vertically aligned walls was utilized to enhance the power density. The porous structure of the power layer was achieved by ice templating this layer on top of the densified energy layer of the cathode. This configuration was found to better utilize the active material of the cathodes and was optimized using numerical simulation and computer modeling. Cells with Li metal anode and LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) at approximately 5 and 20 mg/cm 2 were cycled at 70 °C at different C-rates. Poly(ethylene oxide) (PEO) with lithium bis-trifluoromethanesulfonimide (LiTFSI) was used for the catholyte and the solid-state electrolyte. The structured cathodes exhibited more than double capacity values as well as better Coulombic efficiency compared to non-structured (single-layer) thick cathodes. In conclusion, synchrotron X-ray tomography and scanning electron microscopy were used to characterize the microstructure of the cathodes.

25 ENERGY STORAGE↗

DOE OSTI2023

Energy Frontier Research Centers: Center for the Computational Design of Functional Layered Materials (CCDM) August 1, 2014 - July 31, 2018; Center for Complex Materials from First Principles (CCM) August 1, 2018 - July 31, 2021 (Final Report)

The mission of the DOE Energy Frontier Research Centers CCDM (2014-2018) and CCM (2018-2021) was to theoretically develop, computationally apply, and experimentally validate electronic structure methods for all materials, with a focus on the complex materials, especially layered and two-dimensional materials, strongly-correlated materials, and liquid water. This was achieved by over 200 published journal articles authored by about 17 senior investigators from physics and chemistry and from theory, computation, and experiment, plus their collaborators. In particular, the Centers confirmed the predictive power of the SCAN (strongly constrained and appropriately normed) density functional, which was constructed to satisfy 17 known exact constraints and several appropriate norms. Without being fitted to real bonded systems, and at a modest computational cost, SCAN correctly predicted covalent, ionic, metallic, hydrogen, and van der Waals bonds in many challenging materials. SCAN gave an improved description of defects in semiconductors, surface properties of metals, seven phases of ice, liquid water, liquid and supercooled silicon, subtle structural distortions in ferroelectrics, formation energies and structural predictions for solids, and critical pressures for structural phase transitions. Perhaps most remarkably, SCAN correctly described some strongly-correlated materials that were previously believed to be beyond the reach of density-functional approximations. SCAN is the only density functional that correctly predicts the band gap closing under chemical doping of the cuprate high-temperature superconducting materials. SCAN also predicts a landscape of competing stripe and magnetic phases in the cuprates. For some materials with some codes, SCAN has convergence problems that are greatly reduced by the CCM-developed r 2 SCAN, without loss of accuracy or rigor. SCAN and r 2 SCAN still make some self-interaction error, which is greatly reduced by the CCDM/ CCM-developed local orbital scaling correction (LOSC). These Centers further proved that the fundamental energy gaps of a solid from an orbital energy difference and from total energy differences are the same for a large class of generalized Kohn-Sham (GKS) functionals, including SCAN and standard hybrid functionals, and that symmetry breaking arises when a dynamic density fluctuation drops to zero frequency. The Centers identified new mechanisms for catalysis in layered materials with ions intercalated between the layers, investigated charge density waves both in model systems and in real layered materials, studied changes of band gap with the number of layers, and explored topological ultrathin films, bent nanoribbons, and defects.

08 HYDROGEN↗

DOE OSTI2021

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. Here, the fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII" and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗

DOE OSTI2024

Enhancing the Functionality of a Hollow Scaffold Solid State Bioreactor via Computer Aided Design Optimization

The concern over greenhouse gases, methane (CH 4 ) and carbon dioxide (CO 2 ), is increasing rapidly. There have been strides to find solutions to this global issue but there is not a clear path to a successful end goal. The concentration of CH 4 and CO 2 in the atmosphere has increased significantly over the last 60 years, methane is a great source of concern due to its ability to trap a high amount of heat in the atmosphere. These greenhouse gases contribute to global warming which has caused changes in the environment, including the melting of ice caps, and altered weather patterns. Solutions for these pressing challenges have led to different avenues of methane mitigation one of which is the development of solid-state bioreactors. These reactors harness the power of biological species that have evolved to use methane as an energy source. The development and optimization of Hollow Scaffold Solid State Bioreactors (HS-SSBR) has become readily available due to the advancements in additive manufacturing technology and accessibility of computer aided design (CAD) software. With laboratory scale experiments, time and effort are of great importance. Enhancing the design of the HS-SSBR to create a more user-friendly interface, but also increase the functionality of the reactor. The reactor's design improvements focus on better dispersion of methane and circulation of media.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Particle Track Classification Using Quantum Associative Memory (Final Technical Report)

This project explored the use of quantum-assisted algorithms for pattern matching in sub-atomic physics experiments. Pattern matching algorithms are commonly employed to prune data of random noise and to help discriminate between signals generated by particle tracks of interest and signals generated by background events. The quantum-assisted algorithms explored in this project were based on an Ising formulation of quantum associative model (QAMM) recall and quantum content-addressable memory (QCAM) recall. The recall is performed by comparing a probe pattern with those stored in a library of patterns encoded in the QAMM/QCAM model. The classification accuracy of QAMM and QCAM recall was determined as a function of detector resolution, noise, and efficiency and pattern density, where pattern density is defined as the ratio of the number of reference signal patterns encoded in the library to each pattern’s length. We found that QAMM achieved high classification accuracy when applied to datasets with low pattern density. QCAM achieved high classification accuracy for datasets with high pattern density and was found to be more robust to detector noise. The project methodology and results are described in detail in our arXiv preprint (arXiv:2011.11848) . This project was conducted by scientists at the Johns Hopkins University Applied Physics Laboratory and Oak Ridge National Laboratory from August 2018 to August 2020 and was supported by DOE grant DE-SC0019497.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Logical gates embedding in artificial spin ice

The realization and study of arrays of interacting magnetic nanoislands, such as artificial spin ices, have reached mature levels of control that allow design and demonstration of exotic, collective behaviors not seen in natural materials. Advances in the direct manipulation of their local, binary moments also suggest a use as nanopatterned, interacting memory media, for computation within a magnetic memory. Recent experimental work has demonstrated the possibility of building logic gates from clusters of interacting magnetic domains, and yet the possibility of large scale integration of such gates can prove problematic even at the theoretical level. Here we introduce theoretically complete sets of logical gates, in principle realizable in an experiment, and we study the feasibility of their integration into tree-like circuits. By evaluating the fidelity control parameter between their collective behavior and their expected logic functionality we determine conditions for integration. Also, we test our numerical results against the presence of disorder in the couplings, showing that the design gate structure is robust to small coupling perturbations, and thus possibly to small imperfections in the fabrication of the islands.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Quality Data Essential for Modeling Water Cycles Effectively

The experimental and observational data is the backbone of effective modeling of all aspects of water cycles. As the computer technology used for these data gathering efforts has improved dramatically in the recent decades, the environmental conditions these instruments are subjected to remain pretty much the same. These measuring instruments are exposed to the sun light, wind, rain, snow and ice; and they are immersed in dirt, water and mud. All of these are damaging to the electronics components used for data gathering, storage, and transmission, which creates many different kinds of data quality issues that affect the quality and reliability of modeling and prediction efforts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Hybrid Cooling and Water Treatment for Resilient, Water-Self-Sufficient Data Centers

The continued growth of data center infrastructure is intensifying demand for freshwater resources, particularly in water-stressed regions, and is increasingly limiting sustainable capacity expansion. This work investigates a conceptual hybrid system that integrates freeze desalination with ultrasonic-assisted ice separation to enable on-site production of pure water from diverse sources, including seawater, brackish groundwater, and reclaimed industrial or agricultural wastewater. Simultaneously, the system produces low-temperature cooling streams that enhance heat removal in high power density computing environments. A process-level thermodynamic analysis is performed across a range of boundary and operating conditions, including variations in feed concentration, freezing temperature, and mass flow rate. The results are presented as performance curves relating feedwater concentration and target purified water output to the corresponding intake flow requirements, enabling estimation of source water demand per unit of Information Technology Equipment (ITE) energy consumption (kWh) across varying data center scales and operational scenarios. The corresponding electrical energy consumption for integrated cooling and water treatment is also evaluated as a function of system operating parameters and target water production levels. These results provide a basis for evaluating system feasibility across different conditions and for identifying parameter ranges in which integrated water treatment and cooling improve resource efficiency, thermal performance, and operational flexibility in data centers.

Elhefny, Aly [ORNL] (ORCID:0000000284907923)↗

DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗

DOE OSTI2022

MPAS-Seaice (v1.0.0): sea-ice dynamics on unstructured Voronoi meshes

Abstract. We present MPAS-Seaice, a sea-ice model which uses the Model for Prediction Across Scales (MPAS) framework and spherical centroidal Voronoi tessellation (SCVT) unstructured meshes. As well as SCVT meshes, MPAS-Seaice can run on the traditional quadrilateral grids used by sea-ice models such as CICE. The MPAS-Seaice velocity solver uses the elastic–viscous–plastic (EVP) rheology and the variational discretization of the internal stress divergence operator used by CICE, but adapted for the polygonal cells of MPAS meshes, or alternatively an integral (“finite-volume”) formulation of the stress divergence operator. An incremental remapping advection scheme is used for mass and tracer transport. We validate these formulations with idealized test cases, both planar and on the sphere. The variational scheme displays lower errors than the finite-volume formulation for the strain rate operator but higher errors for the stress divergence operator. The variational stress divergence operator displays increased errors around the pentagonal cells of a quasi-uniform mesh, which is ameliorated with an alternate formulation for the operator. MPAS-Seaice shares the sophisticated column physics and biogeochemistry of CICE and when used with quadrilateral meshes can reproduce the results of CICE. We have used global simulations with realistic forcing to validate MPAS-Seaice against similar simulations with CICE and against observations. We find very similar results compared to CICE, with differences explained by minor differences in implementation such as with interpolation between the primary and dual meshes at coastlines. We have assessed the computational performance of the model, which, because it is unstructured, runs with 70 % of the throughput of CICE for a comparison quadrilateral simulation. The SCVT meshes used by MPAS-Seaice allow removal of equatorial model cells and flexibility in domain decomposition, improving model performance. MPAS-Seaice is the current sea-ice component of the Energy Exascale Earth System Model (E3SM).

58 GEOSCIENCES↗

DOE OSTI2021

High-Resolution Spatio-Temporal Estimation of Net Ecosystem Exchange in Ice-Wedge Polygon Tundra Using In Situ Sensors and Remote Sensing Data

Land-atmosphere carbon exchange is known to be extremely heterogeneous in arctic ice-wedge polygonal tundra regions. In this study, a Kalman filter-based method was developed to estimate the spatio-temporal dynamics of daytime average net ecosystem exchange (NEEday) at 0.5-m resolution over a 550 m by 700 m study site. We integrated multi-scale, multi-type datasets, including normalized difference vegetation indices (NDVIs) obtained from a novel automated mobile sensor system (or tram system) and a greenness index map obtained from airborne imagery. We took advantage of the significant correlations between NDVI and NEEday identified based on flux chamber measurements. The weighted average of the estimated NEEday within the flux-tower footprint agreed with the flux tower data in term of its seasonal dynamics. We then evaluated the spatial variability of the growing season average NEEday, as a function of polygon geomorphic classes; i.e., the combination of polygon types—which are known to present different degradation stages associated with permafrost thaw—and microtopographic features (i.e., troughs, centers and rims). Our study suggests the importance of considering microtopographic features and their spatial coverage in computing spatially aggregated carbon exchange.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Measurement of Photovoltaic Module Deformation Dynamics During Hail Impact Using Digital Image Correlation

Stereo high-speed video of photovoltaic modules undergoing laboratory hail tests was processed using digital image correlation to determine module surface deformation during and immediately following impact. The purpose of this work was to demonstrate a methodology for characterizing module impact response differences as a function of construction and incident hail parameters. Video capture and digital image analysis were able to capture out-of-plane module deformation to a resolution of ±0.1 mm at 11 kHz on an in-plane grid of 10 × 10 mm over the area of a 1 × 2 m commercial photovoltaic module. With lighting and optical adjustments, the technique was adaptable to arbitrary module designs, including size, backsheet color, and cell interconnection. Furthermore, impacts were observed to produce an initially localized dimple in the glass surface, with peak deflection proportional to the square root of incident energy. Subsequent deformation propagation and dissipation were also captured, along with behavior for instances when the module glass fractured. Natural frequencies of the module were identifiable by analyzing module oscillations postimpact. Limitations of the measurement technique were that the impacting ice ball obscured the data field immediately surrounding the point of contact, and both ice and glass fracture events occurred within 100 μs, which was not resolvable at the chosen frame rate. Increasing the frame rate and visualizing the back surface of the impact could be applied to avoid these issues. Applications for these data include validating computational models for hail impacts, identifying the natural frequencies of a module, and identifying damage initiation mechanisms.

14 SOLAR ENERGY↗

DOE OSTI2023

LANL Institutional Computing Close-out Report for Project t22_ocean_time_step

This year, my team made a very productive use of LANL IC time. Four papers were published that used IC resources, and two more are under review. These publications fall into three categories: 1.) Improving tide modeling in the global ocean by adding self attraction and loading (Barton et. al. 2022), ice shelf cavities (Pal et al. 2023), and local time stepping (Lilly et al. 2023). 2.) A new sea ice numerical formulation (Capodaglio 2023). 3.) Performance comparisons, methods, and test cases for ocean model development, This includes a verification suite for ocean models (Bishnu, submitted) and a comparison between Julia and Fortran (Strauss, submitted). A new time-stepping method was introduced in Calandrini et al. (2022). These publications, and the use of LANL IC resources, go hand-in-hand with our mentoring efforts to train young scientists. Two of the lead authors are graduate students who are conducting their PhD research: Kristin Barton at the University of Michigan and Jeremy Lilly at Oregon State. These are both their first publications, and they both have DOE funding and DOE mentors. In addition, Bishnu is a post-doctoral researcher at LANL; Strauss conducted his research as a senior in high school; and Capodaglio, Calandrini and Pal are all early-career scientists who were converted to staff in 2021 or 2022. Here we highlight two publications on improvements in tidal modeling: Barton et al. 2022 and Lilly et al. 2023.

58 GEOSCIENCES↗

DOE OSTI2025

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗