Search NASA⌕ Search

SEARCH · Search NASA

Results for “density”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Imposing correct jellium response is key to predict the density response by orbital-free DFT

Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Skyrme pseudopotentials at next-to-next-to-leading order: Construction of local densities and first symmetry-breaking calculations

There is an ongoing quest to improve on the spectroscopic quality of nuclear energy density functionals (EDFs) of the Skyrme type through extensions of its traditional form. One direction for such activities is the inclusion of terms of higher order in gradients in the EDF. We report on exploratory symmetry-breaking calculations performed for an extension of the Skyrme EDF that includes central terms with four gradients at next-to-next-to-leading order (N2LO) and for which the high-quality parametrization SN2LO1 has been constructed recently. Up to now, the investigation of such functionals with higher-order terms was limited to infinite matter and spherically symmetric configurations of singly and doubly magic nuclei. We address here nuclei and phenomena that require us to consider axial and nonaxial deformation, both for reflection-symmetric and also reflection-asymmetric shapes, as well as the breaking of time-reversal invariance. Achieving these calculations demanded a number of formal developments. These all resulted from the formulation of the N2LO EDF requiring the introduction of new local densities with additional gradients that are not present in the EDF at NLO. Their choice is not unique, but can differ in the way the gradients are coupled. While designing a numerical implementation of N2LO EDFs in Cartesian three-dimensional coordinate-space representation, we have developed a novel definition and a new unifying notation for normal and pair densities that contain gradients at arbitrary order. Besides having mnemonic advantages, the new notation allows for the easy identification of redundancies and reducibilities in a given set of local densities, and the new definition makes it straightforward to construct densities that automatically adopt the symmetries of the many-body state they are constructed from. The resulting scheme resolves several issues with some of the choices that have been made for local densities in the past, in particular when breaking time-reversal symmetry. Guided by general practical considerations, we propose an alternative form of the N2LO contribution to the Skyrme EDF that is built from a different set of densities. It has exactly the same physics content, but is much more efficient to handle in formal discussions and, compared to the original formulation, leads to a substantial reduction of computational cost and memory requirements in deformed codes. As representative examples for the performance of SN2LO1, we have chosen the ground states of even-even Kr and Nd isotopes, the fission barrier of 240 Pu as well as the superdeformed rotational band of 194 Hg. Overall, for the nuclei and phenomena studied here, the SN2LO1 parametrization does not yet present a systematic improvement over standard NLO parametrizations. This finding calls for improved fit protocols that better discriminate between NLO and N2LO terms and better exploit the unique features of the additional degrees of freedom offered by the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

eFESTA: Ensemble Feature Exploration with Surface Density Estimates

We propose surface density estimate (SDE) to model the spatial distribution of surface features-isosurfaces, ridge surfaces, and streamsurfaces-in 3D ensemble simulation data. The inputs of SDE computation are surface features represented as polygon meshes, and no field datasets are required (e.g., scalar fields or vector fields). The SDE is defined as the kernel density estimate of the infinite set of points on the input surfaces and is approximated by accumulating the surface densities of triangular patches. We also propose an algorithm to guide the selection of a proper kernel bandwidth for SDE computation. An ensemble Feature Exploration method based on Surface densiTy EstimAtes (eFESTA) is then proposed to extract and visualize the major trends of ensemble surface features. For an ensemble of surface features, each surface is first transformed into a density field based on its contribution to the SDE, and the resulting density fields are organized into a hierarchical representation based on the pairwise distances between them. The hierarchical representation is then used to guide visual exploration of the density fields as well as the underlying surface features. We demonstrate the application of our method using isosurface in ensemble scalar fields, Lagrangian coherent structures in uncertain unsteady flows, and streamsurfaces in ensemble fluid flows.

Computational modeling↗
DOE OSTI2024

Multifactorial genetic control and magnesium levels govern the production of a Streptomyces antibiotic with unusual cell density dependence

Streptomyces bacteria are renowned both for their antibiotic production capabilities and for their cryptic metabolic potential. Their metabolic repertoire is subject to stringent genetic control, with many of the associated biosynthetic gene clusters being repressed by the conserved nucleoid-associated protein Lsr2. In an effort to stimulate new antibiotic production in wild Streptomyces isolates, we leveraged the activity of an Lsr2 knockdown construct and successfully enhanced antibiotic production in the wild Streptomyces isolate WAC07094. We determined that this new activity stemmed from increased levels of the angucycline-like family member saquayamycin. Saquayamycin has both antibiotic and anti-cancer activities, and intriguingly, beyond Lsr2-mediated repression, we found saquayamycin production was also suppressed at high density on solid or in liquid growth media; its levels were greatest in low-density cultures. This density-dependent control was exerted at the level of the cluster-situated regulatory gene sqnR and was mediated in part through the activity of the PhoRP two-component regulatory system, where deleting phoRP led to both constitutive antibiotic production and sqnR expression. This suggests that PhoP functions to repress the expression of sqnR at high cell density. We further discovered that magnesium supplementation could alleviate this density dependence, although its action was independent of PhoP. Finally, we revealed that the nitrogen-responsive regulators GlnR and AfsQ1 could relieve the repression exerted by Lsr2 and PhoP. Intriguingly, we found that this low density-dependent production of saquayamycin was not unique to WAC07094; saquayamycin production by another wild isolate also exhibited low-density activation, suggesting that this spatial control may serve an important ecological function in their native environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Fermions at finite density in the path integral approach

Abstract We study relativistic fermionic systems in 3 + 1 spacetime dimensions at finite chemical potential and zero temperature, from a path-integral point of view. We show how to properly account for theiεterm that projects on the finite density ground state, and compute the path integral analytically for free fermions in homogeneous external backgrounds, using complex analysis techniques. As an application, we show that the U(1) symmetry is always linearly realized for free fermions at finite charge density, differently from scalars. We study various aspects of finite density QED in a homogeneous magnetic background. We compute the free energy density, non-perturbatively in the electromagnetic coupling and the external magnetic field, obtaining the finite density generalization of classic results of Euler-Heisenberg and Schwinger. We also obtain analytically the magnetic susceptibility of a relativistic Fermi gas at finite density, reproducing the de Haas-van Alphen effect. Finally, we consider a (generalized) Gross-Neveu model forNinteracting fermions at finite density. We compute its non-perturbative effective potential in the large-Nlimit, and discuss the fate of the U(1) vector and$$ {\mathbb{Z}}_2^A $$ ℤ 2 A axial symmetries.

Physics↗
DOE OSTI2026

Hot Hydrogen Exposure of U x Zr 1-x C y Nuclear Fuel: The Influence of Composition and Density

Refractory carbide nuclear fuel has been one of the most promising fuel candidates for space nuclear propulsion due to its high melting point, temperature stability, and compatibility in a hot hydrogen environment. In this study, U x Zr 1-x C y fuel was produced by means of a carbothermic reduction process in different UC compositions including 5,10, 20, and 30 at.% UC in the fuel compound. The powder feedstock was consolidated via direct current sintering with densities up to 97% of the theoretical density. The samples with different U x Zr 1-x C y compositions were exposed to hot hydrogen at 2600 K for a cumulative time of 300 min. The samples were characterized by SEM, XRD, density, and measured for mass losses. The high-density samples displayed improved performance in hot hydrogen by minimizing porous sites and reducing areas of direct contact with hydrogen gas, leading to reduced mass losses. Variations in sample density proved to induce large changes in mass loss rates, increasing them up to 90%. The compositions with higher UC content reported the largest mass losses in the study. The loss of uranium occurred primarily at the surfaces exposed to the hot hydrogen where changes in the lattice constant confirmed losses exceeding 50% of the initial UC content in higher compositions, specifically to U 0.3 Zr 0.7 C y . XRD analyses revealed the presence of UH 3 in U 0.3 Zr 0.7 C y suggesting that metallic uranium formed inside the sample as a product of carbon losses. High-density U x Zr x-1 C y fuel with UC concentrations at or below 20 at.% UC exhibited stability and negligible density changes in a high temperature hydrogen environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Dislocation-density evolution and pileups in bicrystalline systems

Here, a dislocation-density crystalline plasticity (DCP) framework based on total and partial dislocation densities interactions was used to investigate the behavior of Cu/Pb bicrystals with a focus on GB effects. The modeling predictions were validated with bicrystal compression micropillar experiments. A key new aspect of the modeling approach is to account for partial dislocation-densities. A GB formulation that is directly linked to GB energies was used to monitor GB transmission and blockages, such that pileups can be monitored and predicted at the GB interfaces for misorientations. The predictions indicate that pileups can form due to fully and partially blocked slip-rates and perfect and partial dislocation-densities. As the nominal strain increases from five to fifteen percent, dislocation-densities and pileups significantly increase by almost an order of magnitude. The proposed validated approach provides a microstructural scale predictive framework that accounts for a myriad of defects related to the interactions of partial and perfect dislocation densities that interact at highly misoriented GBs; it is these interactions that are critical to the formation and evolution of dislocation-density pileups that can lead to physically limiting stress accumulations in bicrystals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

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

Forest Canopy Density Effects on Snowpack Across the Climate Gradients of the Western United States Mountain Ranges

Limitations inherent to field studies such as limited spatial and temporal coverage has made it very challenging, if not impossible, to explicitly disentangle the effect of canopy density from climate on snowpack dynamics. Hence, our understanding of how variations in canopy density can affect snow processes under diverse climate conditions has been limited. To address this knowledge gap, this study uses a physics-based modeling approach to represent most key processes that drive snow accumulation and melt in the open and forest with varied forest density, across the climate gradients of the Western U.S. mountain ranges, as represented by 228 Snow Telemetry (SNOTEL) locations. Simulations suggest that, under most winter climate conditions, canopy density exerts a strong influence on under-canopy snowpack dynamics and their controlling processes. Changing canopy density can shift the direction of the modeled differential peak snow accumulation and snow disappearance date (SDD) between the open and forest. For example, while simulated snowpack lasts longer in the open than the high-density forest at most wet/warm and dry/cold locations by up to 55 days, about 90% of all locations show a longer snow duration in the low-density forest by up to 27 days as compared to the open. Interannual variability in climate can also change the magnitude and direction of canopy impact on snow processes. By comparing model simulations to published field observations, and we also identify conditions that require additional field observations to improve process-scale understanding of forest-snow interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Inhomogeneous high temperature melting and decoupling of charge density waves in spin-triplet superconductor UTe 2

Charge, spin and Cooper-pair density waves have now been widely detected in exotic superconductors. Understanding how these density waves emerge — and become suppressed by external parameters — is a key research direction in condensed matter physics. Here we study the temperature and magnetic-field evolution of charge density waves in the rare spin-triplet superconductor candidate UTe 2 using scanning tunneling microscopy/spectroscopy. We reveal that charge modulations composed of three different wave vectors gradually weaken in a spatially inhomogeneous manner, while persisting to surprisingly high temperatures of 10–12 K. We also reveal an unexpected decoupling of the three-component charge density wave state. Our observations match closely to the temperature scale potentially related to short-range magnetic correlations, providing a possible connection between density waves observed by surface probes and intrinsic bulk features. Importantly, charge density wave modulations become suppressed with magnetic field both below and above superconducting T c in a comparable manner. Our work points towards an intimate connection between hidden magnetic correlations and the origin of the unusual charge density waves in UTe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Charge density wave in a band insulator

Charge-density-wave implies a periodic modulation of the charge density. Typically observed in metallic systems, charge-density-wave arise from Fermi surface instabilities, resulting in the total or partial gapping of the Fermi surface. Here, we present experimental evidence for a charge-density wave state emerging in a band insulator which has no Fermi surface. The bulk and surface of our material platform, 𝛂-Bi 4 Br 4 , is gapped over the entire Brillouin zone. Through topographic and spectroscopic imaging at low temperatures, we unveil an unexpected unidirectional charge modulation in 𝛂-Bi 4 Br 4 , breaking the lattice translation symmetry. The charge-density-wave develops at temperatures below 40 K and adds an energy gap atop the existing insulating gap of 𝛂-Bi 4 Br 4 . Furthermore, our transport measurements reveal nonlinear electrical conduction, a phenomenon conventionally associated with the sliding or phason mode of incommensurate charge-density-waves. These unusual observations represent a different type of charge-density-wave

Phase transitions and critical phenomena↗
DOE OSTI2022

A representation-independent electronic charge density database for crystalline materials

Abstract In addition to being the core quantity in density-functional theory, the charge density can be used in many tertiary analyses in materials sciences from bonding to assigning charge to specific atoms. The charge density is data-rich since it contains information about all the electrons in the system. With the increasing prevalence of machine-learning tools in materials sciences, a data-rich object like the charge density can be utilized in a wide range of applications. The database presented here provides a modern and user-friendly interface for a large and continuously updated collection of charge densities as part of the Materials Project. In addition to the charge density data, we provide the theory and code for changing the representation of the charge density which should enable more advanced machine-learning studies for the broader community.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Carbon ablator areal density at fusion burn: Observations and trends at the National Ignition Facility

For inertial confinement fusion experiments, the pusher is composed of a high-density deuterium tritium cyrogenic fuel layer and an ablator, often made of carbon. In an ideal, no-mix implosion, increasing the areal density of the pusher transfers more pressure to the hot spot and increases the hot spot confinement time. There has been a lack of knowledge about the final compressed state of the ablator for implosions at the National Ignition Facility. 14 MeV fusion neutrons inelastically scattering on the remaining carbon ablator excites a nuclear metastable state that emits a prompt 4.4 MeV gamma ray. The gamma reaction history diagnostic data, when reduced by a new data analysis technique, can isolate and measure the carbon gamma rays, which are proportional to the areal density of the ablator during fusion burn. The trends over many National Ignition Facility campaigns show that the ablator areal density is weakly sensitive to the maximum shell velocity, the cold fuel radius, the ablator mass remaining, or the laser picket intensity. Controlled parameter scans reveal that, for specific campaigns, ablator compression has a strong dependence on laser coast time, high Z dopants, and the laser drive foot duration. Using a model of the compressed ablator density profile reveals that the greatest variation of the ablator areal density comes from its thickness, with highly compressed, thin layers having high areal density values. The compression and thickness of the ablator are other metrics that designers should understand to differentiate the types of capsule degradation and maximize the inertial confinement fusion performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Forward and backward stimulated Raman scattering in multi-speckled beams: Density dependence and effects on cross-beam energy transfer

Particle-in-cell simulations of cross-beam energy transfer (CBET) between multi-speckled beams reveal a potentially prominent energy loss from forward stimulated Raman scattering (FSRS) in the amplified seed beam [Yin et al., Phys. Plasmas 26, 082708 (2019)]. In this work, an intensity scan is used to examine this energy loss in the CBET-amplified seed beam at intensities of a few times 1014 W/cm2. Electron tracers demonstrate both the trapping in the electron plasma waves that is sustained over multiple speckles and the subsequent acceleration to velocities well beyond the phase velocity. Motivated by expanding our understanding of this energy loss mechanism, we isolate the FSRS density dependence by using a 2×1014 W/cm2 pulse to scan densities between 0.06 ncr and 0.12 ncr; a spectral analysis characterizes the energy losses from scattering with the changing densities. At saturation, both greater beam attenuation is measured at the laser frequency and larger scattering FSRS angles are observed as density increases, along with a rapid transition to primarily backscatter SRS between 0.10 ncr and 0.12 ncr, earlier than predicted by linear theory. Density gradients are shown to effectively mitigate oblique FSRS and beam energy losses in the two lower density cases, whereas 0° forward scattering is primarily suppressed in the higher density cases. Finally returning to the two-beam CBET system of interest, the interplay and mitigation of the secondary instabilities at larger domain and beam spatial scales are discussed, here also examining backward stimulated Brillouin scattering.

Stark, D. J. (ORCID:0000000203642733)↗
DOE OSTI2022

Quantitative spatiotemporal density evolution of aluminum heated purely by monochromatic electrons

A spatially resolved air-wedge shearing interferometer and shadowgraph diagnostic provides measurements of electron density with a resolution of ∼40 μm. A ∼100-ns-long, monoenergetic electron bunch at 19.8 MeV and a current of 1.4 kA (8.5×1014 e−) heats 100-μm-thick aluminum (Al) foils in a 1-mm-spot to Te∼1 eV. A 5-ns-long, ∼60 mJ, frequency doubled Nd:YAG laser probes the dense Al plasma. Electron densities up to 1020cm−3 are resolved; the maximum resolvable density is limited by opacity, transmission, and spatial fringe density achievable with the detector. This diagnostic provides measurements of the total phase shift, transmission, and electron density. Several measurements at different time slices provide the ability to determine the velocity of the leading edge of the shadowgraph and compare it to the motion of different density shells. These measurements are also compared to radiation hydrodynamics simulations. A rough quantitative agreement is shown between the hydro simulations and the measurements; there are differences in the exact density distributions.

Coleman, J. E. (ORCID:0000000308873757)↗
DOE OSTI2022

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Al–W gradient density materials—Processing and dynamic ramp compression

Materials with high-density gradients are desired for controlling loading paths in dynamic compression, important for studying material properties in extreme conditions and inertial confinement fusion. The large density difference between Al and W makes them ideal choices for producing gradient density materials, but their extremely different melting temperatures make them challenging to fabricate simultaneously. We report a method for producing Al–W porosity-free materials with a fourfold increase in density (2.7–11 g/cm 3 ) across the composition range, from Al-rich to W-rich, without intermetallic phase formation. This was achieved by understanding the aluminum-dominated densification behavior and examining the influence of pressure and temperature on the densification of Al–W composites. Dynamic compression experiments conducted with the Al–W gradient density material produced shock ramp compressions as expected based on the designed composition, and the performed hydrodynamics simulations showed excellent agreement with experimental results. The results demonstrate that current activated pressure-assisted densification allows for the easy and rapid fabrication of gradient density materials with significant density gradients and tailored compositions, facilitating precise control of the loading paths. These materials have the potential to create customized pressure drives for advancing the fields of material science in extreme environments and dynamic compression.

Alloys↗
DOE OSTI2024

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗