Search NASA⌕ Search

SEARCH · Search NASA

Results for “water 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Reducing stomatal density by expression of a synthetic epidermal patterning factor increases leaf intrinsic water use efficiency and reduces plant water use in a C4 crop

Abstract Enhancing crop water use efficiency (WUE) is a key target trait for climatic resilience and expanding cultivation on marginal lands. Engineering lower stomatal density to reduce stomatal conductance (gs) has improved WUE in multiple C3 crop species. However, reducing gs in C3 species often reduces photosynthetic carbon gain. A different response is expected in C4 plants because they possess specialized anatomy and biochemistry which concentrates CO2 at the site of fixation. This modifies the relationship of photosynthesis (AN) with intracellular CO2 concentration (ci), such that photosynthesis is CO2 saturated and reductions in gs are unlikely to limit AN. To test this hypothesis, genetic strategies were investigated to reduce stomatal density in the C4 crop sorghum. Constitutive expression of a synthetic epidermal patterning factor (EPF) transgenic allele in sorghum led to reduced stomatal densities, reduced gs, reduced plant water use, and avoidance of stress during a period of water deprivation. In addition, moderate reduction in stomatal density did not increase stomatal limitation to AN. However, these positive outcomes were associated with negative pleiotropic effects on reproductive development and photosynthetic capacity. Avoiding pleiotropy by targeting expression of the transgene to specific tissues could provide a pathway to improved agronomic outcomes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Soil physical and chemical measurements for topsoils collected during NEON campaign in East River, CO (06/14/2018-06/28/2018)

The package is part of the DOE Watershed Function Science Focus Area (SFA) project and includes soil physical and chemical measurements from topsoils collected at the East River, Colorado, in conjunction with the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) survey conducted in June 2018. The soil measurements include soil bulk density, soil volumetric water content, soil microbial biomass C (Carbon), N (Nitrogen) and C:N (C to N ratio), soil DNA yield, soil total extractable organic C, soil total extractable N, soil extractable nitrate, soil extractable ammonium, soil dissolved inorganic N, soil dissolved organic N, soil pH, soil TOC400 (total organic carbon at 400°C), soil ROC (residual oxidizable carbon), soil TIC (total inorganic carbon), soil TOC (total organic carbon), soil TC (total carbon), soil N, soil OM (organic matter) loss on ignition. Additional associated site metadata can be found in the ESS-DIVE package 10.15485/1618130. The dataset includes (1) 2018_NEON_soil_physical_chemical_measurements.csv: soil physical and chemical measurements indexed by soil sample IGSNs; (2) samples.csv: sample metadata file used to register International Generic Sample Numbers (IGSNs); (3) flmd.csv: file level metadata file; and (4) dd.csv: data dictionary file. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS (Catchment Hydrology and ↗
DOE OSTI2024

Microbial abundances and carbon use under ambient temperature or experimental warming in a southern boreal peatland

Abstract Organic peat soils occupy relatively little of the global land surface area but store vast amounts of soil carbon in northern latitudes where climate is warming at a rapid pace. Warming may result in strong positive feedbacks of carbon loss and global climate change driven by microbial processes if warming alters the balance between primary productivity and decomposition. To elucidate effects of warming on the microbial communities mediating peat carbon dynamics, we explored the abundance of broad microbial groups and their source of carbon (i.e. old carbon versus more recently fixed photosynthate) using microbial lipid analysis (δ 13 C PLFA) of peat samples under ambient temperatures and before/after initiation of experimental peat warming (+ 2.25, + 4.5, + 6.75, and + 9 °C). This analysis occurred over a profile to 2 m depth in an undrained, ombrotrophic peat bog in northern Minnesota. We found that the total microbial biomass and individual indicator lipid abundances were stratified by depth and strongly correlated to temperature under ambient conditions. However, under experimental warming, statistically significant effects of temperature on the microbial community were sporadic and inconsistent. For example, 3 months after experimental warming the relative abundance of Gram-negative bacterial indicators across depth combined and > 50 cm depth and Gram-positive bacterial indicators at 20–50 cm depth showed significant positive relationships to temperature. At that same timepoint, however, the relative abundance of Actinobacterial indicators across depth showed a significant negative relationship to temperature. After 10 months of experimental warming, the relative abundance of fungal biomarkers was positively related to temperature in all depths combined, and the absolute abundance of anaerobic bacteria declined with increasing temperature in the 20–50 cm depth interval. The lack of observed response in the broader microbial community may suggest that at least initially, microbial community structure with peat depth in these peatlands is driven more by bulk density and soil water content than temperature. Alternatively, the lack of broad microbial community response may simply represent a lag period, with more change to come in the future. The long-term trajectory of microbial response to warming in this ecosystem then could either be direct, after this initial lag time, or indirect through other physical or biogeochemical changes in the peat profile. These initial results provide an important baseline against which to measure long-term microbial community and carbon-cycling responses to warming and elevated CO 2 .

Felice, Mark↗
DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The influence of submarine canyons-related processes on recent benthic foraminiferal distribution, Espírito Santo Basin, Southeastern Brazil

The complex topography of submarine canyons may result in different composition of benthic foraminifera assemblages. To understand how trophic, hydrological and sedimentological conditions in submarine canyons can influence the distribution of benthic foraminifera, and to use this information to corroborate paleoenvironmental interpretations for the Holocene, we investigated recent benthic foraminiferal assemblages (total fauna >63 μm) and sedimentological data in two canyons (CANWN and CAND) in the Espírito Santo Basin (ESB) between 18°20' and 21°20' S. Surface sediment samples (0–2 cm) were collected inside the canyons (150 to 1300 m water depth) and in shelf-slope adjacent transects (50 to 1300 m water depth). The density, taxonomic diversity, and assemblage composition of benthic foraminifera change with depth and location. The distinct ecological preferences of the most abundant taxa allowed us to recognize five benthic foraminiferal groups. Three groups (V, III, and I) are present in different bathymetric sectors; Group V: outer shelf (50 m), Group III: upper, and middle – lower slope (150 to 400 m), and Group I: middle – lower slope (1000 to 1300 m). Groups II and IV show no characteristic bathymetric distribution and are present only in CAND and in CANWN, respectively. Group V consists of Hanzawaia boueana, Peneroplis planatus, and Quinqueloculina lamarckiana; Group III is dominated by Globocassidulina rossensis and Trifarina spp.; Group I consists of Globocassidulina crassa, Bolivina lowmani, Gavelinopsis versiformis, Alabaminella weddellensis, and Epistominella exigua. The main species in Group II (CAND, 150, 1000 to 1300 m) are Trifarina angulosa, Globocassidulina subglobosa, and Discorbis vilardeboanus. Group IV (middle – lower CANWN, 1000 to 1300 m), consists mainly of agglutinated species Glomospira charoides, Rhabdammina abyssorum, and Psammosphaera fusca. Further, our data suggest that the quantity (and quality) of food supply, hydrodynamic conditions and sediment properties are the main drivers controlling the bathymetric distribution of benthic foraminiferal assemblages in both canyons. The middle – lower CANNW revealed unstable trophic conditions, related to terrigenous sediment input due to turbidity currents. In CAND, the foraminiferal assemblages ecology indicated sufficient organic matter supply that favors species establishment and diversity, indicating a more productive and less unstable environment than in CANWN.

58 GEOSCIENCES↗
DOE OSTI2021

A new generation of uranium coincidence fast neutron collars for assay of LWR fresh fuel assemblies

The active uranium neutron coincidence collar provides a means of non-destructively assaying the fissile linear density of Light Water Reactor fresh fuel assemblies containing low enriched uranium. These neutron collars can operate in two modes: a thermal and a fast mode. In fast mode, a neutron collar has an added cadmium (Cd) liner in the sample cavity of the detector to reduce the impact of the burnable poison (thermal neutron absorber) on the detector signal (doubles). The main advantage for operating in fast mode is a detected signal that is less dependent of the burnable neutron poison content and thus less dependent on facility operator declarations. The drawback is that operating in fast mode requires a longer measurement time (~hour vs tens of minutes for thermal mode) to achieve the statistically needed precision in the measurements. The trend in the modern reactor fuel assemblies is moving to higher burnup by using higher initial enrichment and, consequently, a higher number of burnable poison rods to compensate the initial neutron reactivity. The increase of the burnable poison loading has motivated the development of a new generation of high efficiency fast neutron collars to allow practical measurements in-field by nuclear inspectors. This paper describes the development and performance evaluation of a new generation of neutron collars, for both boiling water reactor (BWR) and pressurized water reactor (PWR) fuels, jointly developed between the US Department of Energy, through Los Alamos National Laboratory, and the Euratom Safeguards Directorate of the European Commission. In this work, we present here calibrations with reference fuel assemblies at Los Alamos National Laboratory as well as the results of in-field measurement campaigns in fuel fabrication plants with modern commercial fuel assemblies. The experimental results show that a typical PWR verification can be made in a total time of 30 min with an uncertainty in the measured mass of 2% at one standard deviation (1σ). A BWR verification can be made in 47 min with an uncertainty in the measured mass of 1.9% at 1σ, or a total time of 20 min with 1σ uncertainty in the measured mass of 2.5%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Visualizing degradation mechanisms in a gas-fed CO 2 reduction cell via operando X-ray tomography

We utilize operando X-ray computed tomography, coupled with real-time electrochemical analysis, to reveal the underlying failure mechanisms of membrane electrode assemblies (MEAs) for electrochemical CO 2 reduction (eCO 2 R). Through operando imaging, we can obtain unprecedented insights into the dynamic behavior of the MEA under different operating conditions, revealing critical changes in interface interactions, phase distribution, and structural integrity over time. Our findings identify phenomena giving rise to the transition from CO 2 R to the hydrogen evolution reaction (HER), as evidenced by shifts in cathode potential and CO 2 R selectivity. The formation of inhomogeneous precipitates at the gas diffusion electrode disrupts the CO 2 supply and reduces the active sites for eCO 2 R, resulting in a shift toward H2 production during low current density operation. Additionally, under high current density conditions, rapid water crossover up to the microporous layer/gas diffusion layer promotes the transition from CO 2 R to HER, further shifting cell potential toward anodic direction. Oscillating voltage conditions reveal the dissolution and regrowth of precipitates, providing direct visualization of the competing selectivity of CO 2 R and HER. This work offers new insight into the degradation mechanisms of MEAs, with implications for the design of more durable CO 2 R systems.

Lee, Sol A [California Institute of Technology (Ca↗
DOE OSTI2024

Chemical-specific Parameters Dataset

The chemical-specific parameters dataset is searchable for physicochemical information for multiple chemicals simultaneously. After selecting chemicals of interest and the desired parameters, the RAIS will generate a table containing the values, chosen according to an established hierarchy. Results can be downloaded in Excel format. Over 40 parameters are available, including melting point, boiling point, density, density, vapor pressure, water solubility, and Henry’s Law constants. Thirteen primary sources are used to populate the dataset of chemical-specific parameters. These values should be used in cancer risk and noncancer hazard assessments for the calculation of preliminary remediation goals (PRGs), hazard characterization, and transport modeling. Users can select up to 1000 chemicals per query. The dataset supports environmental risk assessments, regulatory decision-making, and environmental planning with tools for benchmarking against risk-based standards. This structured approach ensures a robust evaluation of environmental risks tailored to regulatory needs.

Dolislager, Fred [Oak Ridge National Laboratory (O↗
DOE OSTI2022

Micro-spectroscopic and freezing characterization of ice-nucleating particles collected in the marine boundary layer in the eastern North Atlantic

Abstract. Formation of atmospheric ice plays a crucial role in the microphysical evolution of mixed-phase and cirrus clouds and thus climate. How aerosol particles impact ice crystal formation by acting as ice-nucleating particles (INPs) is a subject of intense research activities. To improve understanding of atmospheric INPs, we examined daytime and nighttime particles collected during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign conducted in summer 2017. Collected particles, representative of a remote marine environment, were investigated for their propensity to serve as INPs in the immersion freezing (IMF) and deposition ice nucleation (DIN) modes. The particle population was characterized by chemical imaging techniques such as computer-controlled scanning electron microscopy with energy-dispersive X-ray analysis (CCSEM/EDX) and scanning transmission X-ray microscopy with near-edge X-ray absorption fine-structure spectroscopy (STXM/NEXAFS). Four major particle-type classes were identified where internally mixed inorganic–organic particles make up the majority of the analyzed particles. Following ice nucleation experiments, individual INPs were identified and characterized by SEM/EDX. The identified INP types belong to the major particle-type classes consisting of fresh sea salt with organics or processed sea salt containing dust and sulfur with organics. Ice nucleation experiments show IMF events at temperatures as low as 231 K, including the subsaturated regime. DIN events were observed at lower temperatures of 210 to 231 K. IMF and DIN observations were analyzed with regard to activated INP fraction, ice-nucleation active site (INAS) densities, and a water activity-based immersion freezing model (ABIFM) yielding heterogeneous ice nucleation rate coefficients. Observed IMF and DIN events of ice formation and corresponding derived freezing rates demonstrate that the marine boundary layer aerosol particles can serve as INPs under typical mixed-phase and cirrus cloud conditions. The derived IMF and DIN parameterizations allow for implementation in cloud and climate models to evaluate predictive effects of atmospheric ice crystal formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Organic carbon burial by river meandering partially offsets bank erosion carbon fluxes in a discontinuous permafrost floodplain

Abstract. Arctic river systems erode permafrost in their banks and mobilize particulate organic carbon (OC). Meandering rivers can entrain particulate OC from permafrost many meters below the depth of annual thaw, potentially enabling the production of greenhouse gases. However, the amount and fate of permafrost OC that is mobilized by river erosion is uncertain. To constrain OC fluxes due to riverbank erosion and deposition, we collected riverbank and floodplain sediment samples along the Koyukuk River, which meanders through discontinuous permafrost in the Yukon River watershed, Alaska, USA, with an average migration rate of 0.52 m yr−1. We measured sediment total OC (TOC) content, radiocarbon activity, water content, bulk density, grain size, and floodplain stratigraphy. Radiocarbon activity and TOC content were higher in samples dominated by silt as compared to sand, which we used to map OC content onto floodplain stratigraphy and estimate carbon fluxes due to river meandering. Results showed that the Koyukuk River erodes and re-deposits a substantial flux of OC each year due to its depth and high migration rate, generating a combined OC flux of a similar magnitude to the floodplain net ecological productivity. However, sediment being eroded from cutbanks and deposited as point bars had similar OC stocks (mean ± 1 SD of 125.3±13.1 kg OC m−2 in cutbanks versus 114.0±15.7 kg OC m−2 in point bars) whether or not the banks contained permafrost. We also observed radiocarbon-depleted biospheric OC in both cutbanks and permafrost-free point bars. These results indicate that a substantial fraction of aged biospheric OC that is liberated from floodplains by bank erosion is subsequently re-deposited in point bars rather than being oxidized. The process of aging, erosion, and re-deposition of floodplain organic material may be intrinsic to river–floodplain dynamics, regardless of permafrost content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparison of plug flow and multi-node stratified tank modeling approaches regarding computational efficiency and accuracy

Residential water heaters contain water stratified by temperature-driven density differences. This implies that a water tank can reach a state in which the top and bottom sections have different temperatures, unless mixing happens. A high degree of thermal stratification can improve the efficiency of some water heaters, by saving the amount of energy required for the heat-up process. Studies of stratification became popular in the 1970s and it remains an active research topic today. The research has led to the development of different models and techniques to better predict and define a stratified tanks behavior. By comparing these models and techniques used previously to describe thermal stratification, the phenomenon could be better understood, exploited, and used to increase efficiency and thermal energy capacity in modern water tanks. From the existing models, we found the one-dimensional standard plug-flow and a multi node model to be appropriate for analyzing the processes of the heat up and cool-down in a water tank. These two models are based on energy balances. This work involved comparing the accuracy and computational effort needed to implement these models. To assess accuracy, we compared both types of existing models to experimental data (also collected in this work) which included a heat up process using an external heat pump. This external process included a layering process that has an eddy diffusivity at five times the rate of thermal diffusion. For this project, we implemented the models in MATLAB, the multi-paradigm numerical computing environment. We quantified model accuracy using the root mean squared error between modeled data and experimental data for six measured tank temperatures. Comparing the accuracy and the computational time taken to run the simulation provides a method to contrast the performance of each model and a way to rate it. The multi node model was run using from 6 to 96 spatial nodes; the plug flow model was run using 1 to 0.001 º C temperature bin sizes. Additionally, timesteps were varied from 4 to 236 s. The results quantify the tradeoff between accuracy and computational time, providing guidance for simulations to intelligently select the best model type and simulation parameters. This research can be used to validate the pre-existing models and possibly improve the modern water tank.

Bulnes, Fernando Karg↗
DOE OSTI2020

Some Aspects of the Liquid Water Thermodynamic Behavior: From The Stable to the Deep Supercooled Regime

Liquid water is considered to be a peculiar example of glass forming materials because of the possibility of giving rise to amorphous phases with different densities and of the thermodynamic anomalies that characterize its supercooled liquid phase. In the present work, literature data on the density of bulk liquid water are analyzed in a wide temperature-pressure range, also including the glass phases. A careful data analysis, which was performed on different density isobars, made in terms of thermodynamic response functions, like the thermal expansion αP and the specific heat differences C P -C V , proves, exclusively from the experimental data, the thermodynamic consistence of the liquid-liquid transition hypothesis. The study confirms that supercooled bulk water is a mixture of two liquid “phases”, namely the high density (HDL) and the low density (LDL) liquids that characterize different regions of the water phase diagram. Furthermore, the C P -C V isobars behaviors clearly support the existence of both a liquid–liquid transition and of a liquid–liquid critical point.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Experimental Steps toward a Density Law for Chlorine-Crediting Criticality Models of Aqueous Plutonium solutions

Criticality accident prevention is an essential safety consideration for all operations involving fissionable and fissile materials. Predictive criticality calculations using advanced neutron transport codes, such as Monte Carlo N-Particle Transport (MCNP), are invaluable tools for designing and implementing operational limits. While indispensable, these tools are limited by the quality and accuracy of the inputs that the user provides to define the modeled system. Parameters such as atomic composition, shape, and material density must be accurately defined to obtain a meaningful result. In the case of material property details for a plutonium model, accurate material characterization data is sometimes sparse. Thus limitations in calculation accuracy can be reduced by improving our knowledge of material properties in the targeted fissile systems. One of the biggest remaining challenges to accurately defining fissile systems is the description of aqueous fissile solutions. Even simple properties, such as density, are not well-known for fissile solutions relevant to nuclear energy and security. Described here are initial efforts undertaken to improve criticality calculation inputs for fissile plutonium chloride solutions in water. This effort is focused on experimentally determining accurate solution characteristics for ternary plutonium chloride/hydrochloric acid/water systems, by measurement of water activity and solution density. The effect of inputting experimental densities for these solutions into MCNP criticality calculations is compared to the traditional approach of modeling an idealized (and fictitious) plutonium-water mixture. Expansion of these efforts to a working density law for aqueous plutonium chloride solutions is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Elucidating the temperature and density dependence of silver chloride hydration numbers in high-temperature water vapor: A first-principles molecular simulation study

Hydration numbers of metal complexes in low-density aqueous solutions are required for developing geochemical models for ore-forming metals and for designing supercritical desalination processes. In this study, we investigate the temperature and density dependence of the hydration numbers of silver chloride at 623 K, 673 K, and 713 K and densities of 10–100 kg/m 3 . Experimental estimates of the hydration number in the literature for AgCl at these conditions are inconclusive and possibly contradictory as to the temperature and density dependence. First-principles molecular simulation presents an attractive alternative to experimental measurements. Specifically, recent work shows that machine-learning-accelerated nested Monte Carlo simulations provide reliable estimates for the hydration numbers of CuCl at 623 K from 10 to 100 kg/m 3 . Using the same technique, we find a monotonic temperature dependence, with the hydration number decreasing slightly with increasing temperature. In addition, the simulation-predicted hydration numbers steadily increase with increasing density. These temperature and density trends are in agreement with certain experimental data sets. Therefore, this work demonstrates how first-principles Monte Carlo simulations assist in resolving discrepancies between experimental data sets. Our simulation results also correctly predict that the hydration number, and thus also the solubility, of AgCl is lower than that of CuCl under the same conditions. Furthermore, the bond length and angle formed between the water complex and AgCl differ from those for CuCl, consistent with the lower solubility of AgCl.

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

Blended fuel property analysis of butyl-exchanged polyoxymethylene ethers as renewable diesel blendstocks

Methyl-terminated polyoxymethylene ethers (MM-POMEs), having the formula CH 3 O-(CH 2 O)n-CH 3 (n = 3-6), are a class of oxygenates with desirable diesel-like fuel properties including high cetane number and low soot formation. However, their low energy density and high water-solubility present barriers to their adoption. Both concerns were recently addressed by our research group by synthesizing a mixture of POME structures having butyl end-groups and n = 1-6, termed B*POME1-6. B*POME1-6 maintained the advantageous properties of the parent MM-POMEs, and exhibited improved energy density and most notably, dramatically decreased water solubility. For evaluation against a set of criteria for a blended diesel blendstock, a 20 vol% blend of B*POME1-6 with a base diesel fuel was investigated here. Oxidation stability, cetane number, sooting tendency, lubricity and conductivity were improved in the B*POME1-6 blend compared with the base diesel, while also maintaining the flash point, cloud point, energy density, viscosity, and boiling point requirements. The B*POME1-6 product demonstrated a synergistic blending behavior at 10 vol% and a linear blending behavior at 20-30 vol% blends, in agreement with similar POME blends at comparable blend levels. Finally, common environmental and toxicity models performed on B*POME1-6 component molecules suggested they have a greater propensity to partition into the water compartment compared to a common diesel surrogate, but with a lower tendency to bioaccumulate.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Self-interaction correction in water–ion clusters

We study the importance of self-interaction errors in density functional approximations for various water–ion clusters. We have employed the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in conjunction with LSDA, PBE, and SCAN to describe binding energies of hydrogen-bonded water–ion clusters, i.e., water–hydronium, water– hydroxide, water–halide, as well as non-hydrogen-bonded water–alkali clusters. In the hydrogen-bonded water–ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules, because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water– ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher–level calculations. Here, the non-hydrogen-bonded water–alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies (≤ 50 meV/H 2 O) for hydrogen-bonded water–ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2023

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar

The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water at deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to vacuum, liquid water, and ice. In this work, we apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30,000 atoms. Future work will be directed towards the calculation of nucleation free energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces on ice nucleation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗