Search NASA⌕ Search

SEARCH · Search NASA

Results for “Supercooled water”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Observed Enhancement of Reflectivity and Electric Field in Long-Lived Florida Anvils

A study of two long-lived Florida anvils showed that reflectivity >20 dBZ increased in area, thickness and sometimes magnitude at mid-level well downstream of the convective cores. In these same regions electric fields maintained strengths >10 kV m1 for many tens of minutes and became quite uniform over tens of kilometers. Millimetric aggregates persisted at 9 to 10 km for extended times and distances. Aggregation of ice particles enhanced by strong electric fields might have contributed to reflectivity growth in the early anvil, but is unlikely to explain observations further out in the anvil. The enhanced reflectivity and existence of small, medium and large ice particles far out into the anvil suggest that an updraft was acting, perhaps in weak convective cells formed by instability generated from the evaporation and melting of falling ice particles. We conclude that charge separation must have occurred in these anvils, perhaps at the melting level but also at higher altitudes, in order to maintain fields >10 kV m 1 at 9 to 10 km for extended periods of time over large distances. We speculate that charge separation occurred as a result of ice-ice particle collisions (without supercooled water being present) via either a non-inductive or perhaps even an inductive mechanism, given the observed broad ice particle spectra, the strong pre-existing electric fields and the many tens of minutes available for particle interactions. The observations, particularly in the early anvil, show that the charge structure in these anvils was quite complex.

Dye, James E.↗
NASA NTRS2013

Understanding the Relationships Between Lightning, Cloud Microphysics, and Airborne Radar-derived Storm Structure During Hurricane Karl (2010)

This study explores relationships between lightning, cloud microphysics, and tropical cyclone (TC) storm structure in Hurricane Karl (16 September 2010) using data collected by the NASA DC-8 and Global Hawk (GH) aircraft during NASA's Genesis and Rapid Intensification Processes (GRIP) experiment. The research capitalizes on the unique opportunity provided by GRIP to synthesize multiple datasets from two aircraft and analyze the microphysical and kinematic properties of an electrified TC. Five coordinated flight legs through Karl by the DC-8 and GH are investigated, focusing on the inner-core region (within 50km of the storm center) where the lightning was concentrated and the aircraft were well coordinated. GRIP datasets are used to compare properties of electrified and nonelectrified inner-core regions that are related to the noninductive charging mechanism, which is widely accepted to explain the observed electric fields within thunderstorms. Three common characteristics of Karl's electrified regions are identified: 1) strong updrafts of 10-20ms21, 2) deep mixed-phase layers indicated by reflectivities.30 dBZ extending several kilometers above the freezing level, and 3) microphysical environments consisting of graupel, very small ice particles, and the inferred presence of supercooled water. These characteristics describe an environment favorable for in situ noninductive charging and, hence, TC electrification. The electrified regions in Karl's inner core are attributable to a microphysical environment that was conducive to electrification because of occasional, strong convective updrafts in the eyewall.

Cloud mircophysics↗
NASA NTRS2016

Altitude Effects on Thermal Ice Protection System Performance; a Study of an Alternative Approach

Research has been conducted to better understand the phenomena involved during operation of an aircraft's thermal ice protection system under running wet icing conditions. In such situations, supercooled water striking a thermally ice-protected surface does not fully evaporate but runs aft to a location where it freezes. The effects of altitude, in terms of air pressure and density, on the processes involved were of particular interest. Initial study results showed that the altitude effects on heat energy transfer were accurately modeled using existing methods, but water mass transport was not. Based upon those results, a new method to account for altitude effects on thermal ice protection system operation was proposed. The method employs a two-step process where heat energy and mass transport are sequentially matched, linked by matched surface temperatures. While not providing exact matching of heat and mass transport to reference conditions, the method produces a better simulation than other methods. Moreover, it does not rely on the application of empirical correction factors, but instead relies on the straightforward application of the primary physics involved. This report describes the method, shows results of testing the method, and discusses its limitations.

ice protection system performance↗
NASA NTRS2018

Design and Development of a Laboratory-Scale Ice Adhesion Testing Device

When an aircraft traverses through clouds containing supercooled water droplets, in-flight icing can occur that negatively affects vehicle performance by increasing weight and drag leading to loss of lift. Super-cooled water droplets present in clouds that impact vehicle surfaces can lead to inflight icing any time during the year.1 Most events occur at temperatures ranging from 0 to -20degC. Ice generated on the aircraft can vary between clear/glaze, rime, and mixed (Fig. 1) depending on air temperature (-5 to -20degC), liquid water content (0.3-0.6 g/m3), and droplet size (median volumetric diameter of 15-40 μm). Current strategies to remove ice are based on active technologies such as pneumatic boots, heated surfaces, and deicing agents (i.e., ethylene- and propylene-based glycols). The latter have potential environmental concerns. A passive approach to mitigate accreting ice that is actively being investigated are protective coatings. An ice mitigating coating could potentially be used as a stand-alone material, but more likely in combination with an active approach. In the latter scenario, potential reduction in power consumption by the active approach may be realized. To determine the ice adhesion strength of impact ice that is representative of the aircraft environment is not a trivial matter. Test methods utilizing slowly formed ice (i.e., freezer ice) do not accurately simulate this environment. Likewise, some testing methodologies involve sample relocation from the icing environment to the test chamber that can result in thermal shock to the sample, thus affecting the results. The Adverse Environment Rotor Test Stand (AERTS) located at Pennsylvania State University (PSU) has been demonstrated to simulate impact icing conditions within the icing envelope for the determination of ice adhesion shear strength (IASS) without removal/relocation of the sample.2 Due to the confidence in results obtained from AERTS, this instrument is in high demand and requires a significant amount of lead time and capital investment to obtain IASS results. As a solution for quickly and economically screening coatings in a controlled manner under impact icing conditions, a laboratory-scale ice adhesion test and dead blades were then removed from the rotor/blade assembly to obtain the final mass. The IASS of the live blade was determined from the difference in mass (before and after testing) of the live and dead blades, the ice shed area, and the rpm of the shed event. The same live blade sample was tested in triplicate at all three test temperatures. Surface roughness was determined using a Bruker Dektak XT Stylus Profilometer. Measurements were conducted using a 12.5 μm tip at a vertical range of 65.5 μm with an applied force of 3 mg. Data were collected over a 1.0 mm length at a resolution of 0.056 μm/point. Five single line scans at different locations were collected and processed using a two-point leveling subtraction. The resultant Ra (arithmetic roughness) and Rq (root mean square roughness) average values were calculated.

Smith, Joseph G., Jr.↗
NASA NTRS2019

The Effect of Stainless Steel 304 Surface Roughness on Ice Adhesion Shear Strength of Accreted Impact Ice

Aircraft in-flight icing is problematic due to the ad-verse effect on vehicle performance. It occurs when supercooled water droplets (SCWD) present in clouds, under the appropriate environmental conditions, col-lide with the aircraft surface resulting in accretion of ice (i.e., impact icing). Impact ice can range from clear/glaze to rime or a combination of the two (i.e., mixed) with the type determined by the air temperature (0 to -20°C), liquid water content (LWC, 0.3-0.6 g/cu.m), and droplet size [median volumetric diameter (MVD) of 15-40 µm] present during accretion.1 These impact icing events generally occur at temperatures ranging from 0 to -20°C. Below -20°C, ice crystals dominate the environment and typically do not adhere to the aircraft surface. A main difference between an impact icing occurrence and a slow growth icing (i.e., freezer ice) one is the speed of the icing event. Besides environmental conditions, ice adhesion strength (IAS) to a metallic substrate depends upon surface roughness. It is known that increasing surface roughness and decreasing temperature lead to in-creases in IAS.

Smith, Joseph G., Jr.↗
NASA NTRS2021

An Investigation of the Grain Size and Adhesion Strength of High-Speed Impact Ice

Ice accumulation on aircraft from supercooled water droplets (impact icing) can increase drag and decrease lift, leading to a decrease in fuel efficiency and a potential for unsafe situations. Anti-icing materials and coatings are of high interest to prevent aircraft icing. Because no anti-icing material suitable for in-flight situations has been realized, icing models aim to predict the accumulation and shape of the ice to understand its effects. In both cases, understanding the physical mechanism behind ice adhesion and shedding is necessary. A recent literature review has shown that the apparent ice adhesion strength to any specific material varies over several orders of magnitude across studies and test methods. One of the reasons for this high variance is likely to be the varying microstructure of the ice used for different studies. Although the microstructure of impact ice has been studied in the past, those conditions were not relevant to that of in-flight aircraft icing. The microstructure of the ice formed will determine both its bulk and surface properties, yet it has been overlooked in the ice adhesion community. Here, we use optical microscopy to study the non-uniform ice crystal grain size in impact ice accreted in NASA Glenn Research Center’s Icing Research Tunnel at various icing conditions. Similar to other environmental conditions studied in the literature, we show an increase in grain size at larger distances from the growth interface; however, the data shows substantially larger grains than what was expected from the prior literature. Ice accrued at lower temperatures has a smaller grain size at all distances from the growth interface. No change in grain size was noticed when varying only the wind speed. A correlational analysis shows that smaller grain sizes at the growth interface corresponded to an increase in the apparent ice adhesion strength when the adhesive failure was interfacial, but not when the failure was cohesive. These results may explain the overall trend in the literature of increasing apparent adhesion strength with decreasing temperature. Incorporating microstructure trends for various environmental conditions in icing models may increase the certainty in predicting ice adhesion strength and shedding, and will provide a better basis for comparison between adhesive tests on ice produced in different facilities and at different storage times.

• Ice Fabric↗
NASA NTRSDate not supplied

Plans for Ice Crystal Icing Tests Using a 3D Heated Test Article at the NASA Icing Research Tunnel

This presentation describes the plans for ice crystal icing and supercooled water icing tests to be conducted in Fall/Winter 2021 at the NASA Icing Research Tunnel. A test article, whose geometry is representative of an inter-compressor duct and strut region of a turbofan engine, has been designed and fabricated for the study of ice crystal icing physics tests. The Simulated Inter-compressor Duct Research Model (SIDRM) will have the ability to heat surfaces to simulate the warm surfaces of the inter-compressor duct within an engine. The test article is instrumented with heaters, heat flux gauges, thermocouples, and pressure taps to measure icing behavior. The aim of these tests is to generate ice shapes on the SIDRM test article under well-characterized conditions. The results of the testing with the SIDRM model will be used to provide icing physics and validation data for development of the engine icing simulation capabilities in GlennICE, the 3D ice prediction tool. This is intended to be a presentation only, and does not have an accompanying paper.

ice-crystal icing↗
NASA NTRS2021

Aerosol–Ice Formation Closure: A Southern Great Plains Field Campaign

Prediction of ice formation in clouds presents one of the grand challenges in the atmospheric sciences. Immersion freezing initiated by ice-nucleating particles (INPs) is the dominant pathway of primary ice crystal formation in mixed-phase clouds, where supercooled water droplets and ice crystals coexist, with important implications for the hydrological cycle and climate. However, derivation of INP number concentrations from an ambient aerosol population in cloud-resolving and climate models remains highly uncertain. We conducted an aerosol-ice formation closure pilot study using a field-observational approach to evaluate the predictive capability of immersion freezing INPs. The closure study relies on co-located measurements of the ambient size-resolved and single-particle composition and INP number concentrations. The acquired particle data serve as input in several immersion freezing parameterizations, that are employed in cloud-resolving and climate models, for prediction of INP number concentrations. We discuss in detail one closure case study in which a front passed through the measurement site, resulting in a change of ambient particle and INP populations. We achieved closure in some circumstances within uncertainties, but we emphasize the need for freezing parameterization of potentially missing INP types and evaluation of the choice of parameterization to be employed. Overall, this closure pilot study aims to assess the level of parameter details and measurement strategies needed to achieve aerosol-ice formation closure. The closure approach is designed to accurately guide immersion freezing schemes in models, and ultimately identify the leading causes for climate model bias in INP predictions.

ice formation↗
NASA NTRS2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed-Phase Arctic Stratus Cloud Layer

Mixed-phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long-lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice-nucleating particle (INP) reservoir dynamics governing immersion-mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol-cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time-independent (singular) number- and surface area-based descriptions and a time-dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT-based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

aerosol↗
NASA NTRS2023

CLASS Observation of Polarized Cloud Emission at Millimeter Wavelengths

The dynamic atmosphere imposes challenges to ground-based cosmic microwave background observation, especially for measurements on large angular scales. The hydrometeors in the atmosphere, mostly in the form of clouds, scatter the ambient thermal radiation and are known to be the main linearly polarized source in the atmosphere. This scattering-induced polarization is significantly enhanced for ice clouds due to the alignment of ice crystals under gravity, which are also the most common clouds seen at the millimeter-astronomy sites at high altitudes. This work presents a multifrequency study of cloud polarization observed by the Cosmology Large Angular Scale Surveyor experiment on Cerro Toco in the Atacama Desert of northern Chile, from 2016–2022, at the frequency bands centered around 40, 90, 150, and 220 GHz. Using a machine-learning-assisted cloud classifier, we made connections between the transient polarized emission found in all four frequencies with the clouds imaged by monitoring cameras at the observing site. The polarization angles of the cloud events are found to be mostly 90° from the local meridian, which is consistent with the presence of horizontally aligned ice crystals. The90 and 150 GHz polarization data are consistent with a power law with a spectral index of 3.90±0.06, while an excess/deficit of polarization amplitude is found at 40/220 GHz compared with a Rayleigh scattering spectrum. These results are consistent with Rayleigh-scattering-dominated cloud polarization, with possible effects from supercooled water absorption and/or Mie scattering from a population of large cloud particles that contribute to the220 GHz polarization.

Earth's clouds↗
NASA NTRS2002

Surface Crystallization of Cloud Droplets: Implications for Climate Change and Ozone Depletion

The process of supercooled liquid water crystallization into ice is still not well understood. Current experimental data on homogeneous freezing rates of ice nucleation in supercooled water droplets show considerable scatter. For example, at -33 C, the reported freezing nucleation rates vary by as much as 5 orders of magnitude, which is well outside the range of measurement uncertainties. Until now, experimental data on the freezing of supercooled water has been analyzed under the assumption that nucleation of ice took place in the interior volume of a water droplet. Here, the same data is reanalyzed assuming that the nucleation occurred "pseudoheterogeneously" at the air (or oil)-liquid water interface of the droplet. Our analysis suggest that the scatter in the nucleation data can be explained by two main factors. First, the current assumption that nucleation occurs solely inside the volume of a water droplet is incorrect. Second, because the nucleation process most likely occurs on the surface, the rates of nuclei formation could differ vastly when oil or air interfaces are involved. Our results suggest that ice freezing in clouds may initiate on droplet surfaces and such a process can allow for low amounts of liquid water (approx. 0.002 g per cubic meters) to remain supercooled down to -40 C as observed in the atmosphere.

Tabazadeh, A.↗
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 OSTI2023

Structural relaxation of water during rapid cooling from ambient temperatures

Experiments investigating the properties of deeply supercooled liquid water are needed to develop a comprehensive understanding of water’s anomalous properties. One approach involves transiently heating nanoscale water films into the supercooled region for several nanoseconds at a time and then interrogating the water films after they have quenched to cryogenic temperatures. To relate the results obtained with this approach to other experiments and simulations on supercooled water, it is important to understand how closely the quenched structure tracks the (metastable) equilibrium structure of water as a function of the transient heating temperature. A key step involves quantifying the extent to which water that is transiently heated to ambient temperatures [hyperquenched water (HQW)] subsequently relaxes toward the structure of low-density amorphous (LDA) ice as it cools. We analyzed the infrared reflection–absorption spectra of LDA, HQW, and crystalline ice films to determine their complex indices of refraction. With this information, we estimate that HQW retains ~50%–60% of a structural motif characteristic of water at high temperatures with the balance comprised of a low-temperature motif. This result, along with results from x-ray diffraction experiments on water and amorphous ices, allows one to quantify the fraction of the high-temperature motif at approximately zero pressure as a function of temperature from 150 to 350 K.

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

Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model

For the past 50 years, researchers have sought molecular models that can accurately reproduce water’s microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ~208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. Further, this work suggests that MB-pol may provide predictive capability for studies of liquid water’s physical properties across broad ranges of thermodynamic states, including the so-called water’s “no man’s land” which is difficult to probe experimentally.

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

Evaluating Cloud Properties at Scott Base: Comparing Ceilometer Observations With ERA5, JRA55, and MERRA2 Reanalyses Using an Instrument Simulator

This study compares CL51 ceilometer observations made at Scott Base, Antarctica, with statistics from the ERA5, JRA55, and MERRA2 reanalyses. To enhance the comparison we use a lidar instrument simulator to derive cloud statistics from the reanalyses which account for instrumental factors. The cloud occurrence in the three reanalyses is slightly overestimated above 3 km, but displays a larger underestimation below 3 km relative to observations. Unlike previous studies, we see no relationship between relative humidity and cloud occurrence biases, suggesting that the cloud biases do not result from the representation of moisture. We also show that the seasonal variation of cloud occurrence and cloud fraction, defined as the vertically integrated cloud occurrence, are small in both the observations and the reanalyses. We also examine the quality of the cloud representation for a set of weather states derived from ERA5 surface winds. The variability associated with grouping cloud occurrence based on weather state is much larger than the seasonal variation, highlighting weather state is a strong control of cloud occurrence. All the reanalyses continue to display underestimates below 3 km and overestimates above 3 km for each weather state. But the variability in ERA5 statistics matches the changes in the observations better than the other reanalyses. We also use a machine learning scheme to estimate the quantity of supercooled liquid water cloud from the ceilometer observations. Ceilometer low-level supercooled liquid water cloud occurrences are considerably larger than values derived from the reanalyses, further highlighting the poor representation of low-level clouds in the reanalyses.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

Using 3-D Numerical Weather Data in Piloted Simulations

This report describes the process of acquiring and using 3-D numerical model weather data sets in NASA Langley's Research Flight Deck (RFD). A set of software tools implement the process and can be used for other purposes as well. Given time and location information of a weather phenomenon of interest, the user can download associated numerical weather model data. These data are created by the National Oceanic and Atmospheric Administration (NOAA) High Resolution Rapid Refresh (HRRR) model, and are then processed using a set of Mathworks' Matlab(TradeMark) scripts to create the usable 3-D weather data sets. Each data set includes radar re ectivity, water vapor, component winds, temperature, supercooled liquid water, turbulence, pressure, altitude, land elevation, relative humidity, and water phases. An open-source data processing program, wgrib2, is available from NOAA online, and is used along with Matlab scripts. These scripts are described with sucient detail to make future modi cations. These software tools have been used to generate 3-D weather data for various RFD experiments.

Daniels, Taumi S.↗
DOE OSTI2024

Partition between supercooled liquid droplets and ice crystals in mixed-phase clouds based on airborne in situ observations

Abstract. The onset of ice nucleation in mixed-phase clouds determines the lifetime and microphysical properties of ice clouds. In this work, we develop a novel method that differentiates between various phases of mixed-phase clouds, such as clouds dominated by pure liquid or pure ice segments, compared with those having ice crystals surrounded by supercooled liquid water droplets or vice versa. Using this method, we examine the relationship between the macrophysical and microphysical properties of Southern Ocean mixed-phase clouds at −40 to 0 °C (e.g. stratiform and cumuliform clouds) based on the in situ aircraft-based observations during the US National Science Foundation Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) flight campaign. The results show that the exchange between supercooled liquid water and ice crystals from a macrophysical perspective, represented by the increasing spatial ratio of regions containing ice crystals relative to the total in-cloud region (defined as ice spatial ratio), is positively correlated with the phase exchange from a microphysical perspective, represented by the increasing ice water content (IWC), decreasing liquid water content (LWC), increasing ice mass fraction, and increasing ice particle number fraction (IPNF). The mass exchange between liquid and ice becomes more significant during phase 3 when pure ice cloud regions (ICRs) start to appear. Occurrence frequencies of cloud thermodynamic phases show a significant phase change from liquid to ice at a similar temperature (i.e. −17.5 °C) among three types of definitions of mixed-phase clouds based on ice spatial ratio, ice mass fraction, or IPNF. Aerosol indirect effects are quantified for different phases using number concentrations of aerosols greater than 100 or 500 nm (N>100 and N>500, respectively). N>500 shows stronger positive correlations with ice spatial ratios compared with N>100. This result indicates that larger aerosols potentially contain ice-nucleating particles (INPs), which facilitate the formation of ice crystals in mixed-phase clouds. The impact of N>500 is also more significant in phase 2 when ice crystals just start to appear in the mixed phase compared with phase 3 when pure ICRs have formed, possibly due to the competing aerosol indirect effects on primary and secondary ice production in phase 3. The thermodynamic and dynamic conditions are quantified for each phase. The results show stronger in-cloud turbulence and higher updraughts in phases 2 and 3 when liquid and ice coexist compared with pure liquid or ice (phases 1 and 4, respectively). The highest updraughts and turbulence are seen in phase 3 when supercooled liquid droplets are surrounded by ice crystals. These results indicate both updraughts and turbulence support the maintenance of supercooled liquid water amongst ice crystals. Overall, these results illustrate the varying effects of aerosols, thermodynamics, and dynamics through various stages of mixed-phase cloud evolution based on this new method that categorizes cloud phases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Liquid-Liquid Transition in Water from First Principles

A long-standing question in water research is the possibility that supercooled liquid water can undergo a liquid-liquid phase transition (LLT) into high- and low-density liquids. Here, we used several complementary molecular simulation techniques to evaluate the possibility of an LLT in an ab initio neural network model of water trained on density functional theory calculations with the SCAN exchange correlation functional. We conclusively show the existence of a first-order LLT and an associated critical point in the SCAN description of water, representing the first definitive computational evidence for an LLT in water from first principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗