SEARCH · Search NASA

Results for “Supercooled phase transitions”

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

SourceResearch connectionFrom yearThrough year

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

69 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Undercooling minimization in ultrasound coupled DTA measurements of molten salts

Molten salts are ionic liquids that are used for the electrolytic pyroprocessing of metals and as heat transfer fluids in very high temperature processes. Recently, halide-based molten salt reactors (MSR) have gained momentum for high density and environmentally responsible electricity generation. The function of these reactors and their fuel cycle depend on a knowledge of the halide salt’s thermodynamic properties [1]. Therefore, high accuracy phase equilibria of MSR relevant base halide and actinide containing salts are needed. Modern thermal analysis of the phase transitions of halide salts is usually done during heating at relatively high scan rates with commercial devices. Normally such measurements are adequate for pure or pseudo-binary salts. However, as the number of components in the mixture increases, accurately resolving the liquidus becomes increasingly difficult. Reversing the scanning mode greatly increases the sensitivity of phase transition measurements but can decrease accuracy due to undercooling [2] — a common occurrence in molten halide systems [3]. This work presents the development of a differential thermal analysis (DTA) cell intended for use in radiological gloveboxes. Non-contact ultrasonic agitation of the halide salt is implemented to limit kinetic limitations on crystallization during cooling to minimize undercooling [4]. Measurements on halide salts are also presented to elucidate the effect of non-contact mixing on undercooling. Reference [1] S. Boyd and C. Taylor, “3 - Chemical fundamentals and applications of molten salts,” in Molten Salt Reactors and Thorium Energy, T. J. Dolan, Ed., Woodhead Publishing, 2017, pp. 29–91. doi: 10.1016/B978-0-08-101126-3.00003-8. [2] K. Nitsch, A. Cihlár, and M. Rodová, “Molten state and supercooling of lead halides,” J. Cryst. Growth, vol. 264, no. 1, pp. 492–498, Mar. 2004, doi: 10.1016/j.jcrysgro.2004.01.011. [3] L. Rycerz, “Practical remarks concerning phase diagrams determination on the basis of differential scanning calorimetry measurements,” J. Therm. Anal. Calorim., vol. 113, no. 1, pp. 231–238, Jul. 2013, doi: 10.1007/s10973-013-3097-0. [4] Md. H. Zahir, S. A. Mohamed, R. Saidur, and F. A. Al-Sulaiman, “Supercooling of phase-change materials and the techniques used to mitigate the phenomenon,” Appl. Energy, vol. 240, pp. 793–817, Apr. 2019, doi: 10.1016/j.apenergy.2019.02.045.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Phase transitions from the fifth dimension

We study the cosmological transition of 5D warped compactifications, from the high-temperature black-brane phase to the low-temperature Randall-Sundrum I phase. The transition proceeds via percolation of bubbles of IR-brane nucleating from the black-brane horizon. The violent bubble dynamics can be a powerful source of observable stochastic gravitational waves. While bubble nucleation is non-perturbative in 5D gravity, it is amenable to semiclassical treatment in terms of a “bounce” configuration interpolating between the two phases. We demonstrate how such a bounce configuration can be smooth enough to maintain 5D effective field theory control, and how a simple ansatz for it places a rigorous lower-bound on the transition rate in the thin-wall regime, and gives plausible estimates more generally. When applied to the Hierarchy Problem, the minimal Goldberger-Wise stabilization of the warped throat leads to a slow transition with significant supercooling. We demonstrate that a simple generalization of the Goldberger-Wise potential modifies the IR-brane dynamics so that the transition completes more promptly. Supercooling determines the dilution of any (dark) matter abundances generated before the transition, potentially at odds with data, while the prompter transition resolves such tensions. We discuss the impact of the different possibilities on the strength of the gravitational wave signals. Via AdS/CFT duality the warped transition gives a theoretically tractable holographic description of the 4D Composite Higgs (de)confinement transition. Our generalization of the Goldberger-Wise mechanism is dual to, and concretely models, our earlier proposal in which the composite dynamics is governed by separate UV and IR RG fixed points. The smooth 5D bounce configuration we introduce complements the 4D dilaton/radion dominance derivation presented in our earlier work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2021

Preconditioning of overcast-to-broken cloud transitions by riming in marine cold air outbreaks

Marine cold air outbreaks (CAOs) commonly form overcast cloud decks that transition into broken cloud fields downwind, dramatically altering the local radiation budget. In this study, we investigate the impact of frozen hydrometeors on these transitions. We focus on a CAO case in the NW Atlantic, the location of the multi-year flight campaign ACTIVATE (Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment). We use MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, version 2) reanalysis fields to drive large eddy simulations with mixed-phase two-moment microphysics in a Lagrangian framework. We find that transitions are triggered by substantial rain (rainwater paths >25 g m−2), and only simulations that allow for aerosol depletion result in sustained breakups, as observed. Using a range of diagnostic ice nucleating particle concentrations, Ninp, we find that increasing ice progressively accelerates transitions, thus abbreviating the overcast state. Ice particles affect the cloud-topped boundary layer evolution, primarily through riming-related processes prior to substantial rain, leading to (1) a reduction in cloud liquid water, (2) early consumption of cloud condensation nuclei, and (3) early and light precipitation cooling and moistening below cloud. We refer to these three effects collectively as “preconditioning by riming”. Greater boundary layer aerosol concentrations available as cloud condensation nuclei (CCN) delay the onset of substantial rain. However, cloud breakup and low CCN concentration final stages are found to be inevitable in this case, due, primarily, to liquid water path buildup. An ice-modulated cloud transition speed suggests the possibility of a negative cloud–climate feedback. To address prevailing uncertainties in the model representation of mixed-phase processes, the magnitude of ice formation and riming impacts and, thereby, the strength of an associated negative cloud–climate feedback process, requires further observational evaluation by targeting riming hot spots with in situ imaging probes that allow for both the characterization of ice particles and abundance of supercooled droplets.

Marine cold air outbreaks↗

NASA NTRS1978

Phase diagram of KHF2 and non-equilibrium effects

The equilibrium diagram for the KHF2-H2O system was constructed from cooling and heating curves for the compositions between 5 wt% and 40 wt% KHF2 and the results are shown. The phase diagrams shown is typical of that of a two component system with miscible liquid phases and whole solid phases consist of pure components. A eutectic point was found at approximately 15% KHF2 which remains completely liquid down to a temperature of -9.0 C. No hydrate formation was observed and no anomalous behavior such as the occurrence of solid transitions or metastable states was observed. The effect of rapid freezing on the equilibrium diagram did not appear, and cooling curves exhibited only one halt. Also, at rapid freezing rates, the supercooling of the solutions was smaller than those observed at the slow cooling rates. The existence of a eutectic composition and the slow rate of dissolution of the salt are used to interpret heat absorption behavior in practical applications of the KHF2-H2O system.

Hobson, M. C.↗

NASA NTRS1993

Thermal processing of omicron-terphenyl - A Raman study

Raman spectra from omicron-terphenyl sealed in glass capillary tubes were acquired over temperatures ranging from 233 to 343 K. Over this temperature range, omicron-terphenyl can exist in the solid-crystalline, molten, supercooled-liquid, and amorphous-glassy states. Representative structure-sensitive spectra are presented along with a qualitative description. The lattice spectral feature occurring at a Raman shift of 123/cm displays abrupt changes in profile and intensity as melting and glass transition occur. The intensity ratio of the C-C stretching 996/cm line and C-H stretching line at 1008/cm displays a change related to the transformation of the ordered-crystalline phase to the disordered-liquid or glassy state. At a Raman shift of 1162/cm the solid-crystalline state shows a splitting due to the effect of the crystal field. The possibility of using these structure-sensitive Raman characteristics for in situ diagnostics during materials processing is discussed.

Biswas, A.↗

NASA NTRS1993

Immiscible phase incorporation during directional solidification of hypermonotectics

Solidification processes in immiscible samples were investigated by directly observing the events taking place at the solid-liquid interface during directional solidification. Visualization of these events was made possible through the use of a transparent metal analog system and a temperature gradient stage assembly fitted to an optical microscope. The immiscible transparent analog system utilized was the succinonitrile-glycerol system. This system has been shown to exhibit the same morphological transitions as observed in metallic alloys of monotectic composition. Both monotectic and hypermonotectic composition samples were directionally solidified in order to gain an improved understanding of the manner in which the excess hypermonotectic liquid is incorporated into the solidifying structure. The processing conditions utilized prevented sedimentation of the excess hypermonotectic liquid by directionally solidifying the samples in very thin (13 microns), horizontally oriented cells. High thermal gradient to growth rate ratios (G/R) were used in an effort to prevent constitutional supercooling and the subsequent formation of L(sub 2) droplets in advance of the solidification front during the growth of fibrous composite structures. Results demonstrated that hypermonotectic composites could be produced in samples up to two weight percent off of the monotectic composition by using a G/R ratio greater than or equal to 4.6 x 10(exp 4) C(s)/mm(sup 2) to avoid constitutional supercooling. For hypermonotectic samples processed with G/R ratios below 4.6 x 10(exp 4) C(s)/mm(sup 2), constitutional supercooling occurred and resulted in slight interfacial instability. For these samples, two methods of incorporation of the hypermonotectic liquid were observed and are reported. The correlation between the phase spacing, lambda, and the growth rate, R, was examined and was found to obey a relationship generally associated with a diffusion controlled coupled growth process. For samples with compositions ranging from the monotectic composition up to 2 percent off of the monotectic composition, data indicated that the square of the phase spacing (lambda) varied linearly with the inverse of the growth rate (R).

Andrews, J. Barry↗

NASA NTRS1993

Homogeneous ice nucleation and supercooled liquid water in orographic wave clouds

This study investigates ice nucleation mechanisms in cold lenticular wave clouds, a cloud type characterized by quasi-steady-state air motions and microphysical properties. It is concluded that homogeneous ice nucleation is responsible for the ice production in these clouds at temperatures below about -33 C. The lack of ice nucleation observed above -33 C indicates a dearth of ice-forming nuclei, and hence heterogeneous ice nucleation, in these clouds. Aircraft measurements in the temperature range -31 to -41 C show the following complement of simultaneous and abrupt changes in cloud properties that indicate a transition from the liquid phase to ice: disappearance of liquid water; decrease in relative humidity from near water saturation to ice saturation; increase in mean particle size; change in particle concentration; and change in temperature due to the release of latent heat. A numerical model of cloud particle growth and homogeneous ice nucleation is used to aid in interpretation of our in situ measurements. The abrupt changes in observed cloud properties compare favorably, both qualitatively and quantitatively, with results from the homogeneous ice nucleation model. It is shown that the homogeneous ice nucleation rates from the measurements are consistent with the temperature-dependent rates employed by the model (within a factor of 100, corresponding to about 1 C in temperature) in the temperature range -35 deg to -38 C. Given the theoretical basis of the modeled rates, it may be reasonable to apply them throughout the -30 to -50 C temperature range considered by the theory.

Heymsfield, Andrew J.↗

NASA NTRSDate not supplied

Boundary-Layer Cloud Modeling Challenges on the North Slope of Alaska

The accurate modeling and prediction of cloud base heights is critical for energy balance calculations and aviation operations, alike. Low-level (i.e., boundary-layer) Arctic clouds can be difficult to model, making prediction of formation and dissipation challenging. Primarily mixed-phase, these clouds typically contain low quantities of supercooled liquid water and often slowly precipitate relatively small amounts of moderately and heavily rimed snow particles. While this appears to be the predominant cloudy state on the North Slope of Alaska (NSA), the delicate balance of microphysical, dynamical, radiative, surface coupling, and advective processes can rapidly shift to heavy snow (with various degrees of riming) or to a complete dissipation of the cloud layer without any precipitation, depending on the dominant processes. Here we strive to disentangle these various processes. First, we compare the predictive performances of four different numerical weather models in forecasting the presence and base-heights of low-level clouds: the High-Resolution Rapid Refresh - Alaska (HRRR-AK) model, the Polar Weather Research and Forecasting (Polar WRF) model, the Unified Model (UM), and the European Centre for Medium-range Weather Forecasting (ECMWF) model. Initial results comparing model output at two U.S. Department of Energy Atmospheric Radiation Measurement (AMT) NSA sites, during the fall season in 2019 and 2022, show that the UM slightly outperforms the HRRR-AK in terms of accurately forecasting the presence of a low-level cloud layer (89% of the time). All models have a significant bias of 300 to 800 meters in forecasting cloud base height (lower than is observed); however, the UM and ECMWF models have the lowest biases. Finally, a case study for a particularly challenging April 2017 thin-cloud event is presented, wherein we compare the performance of four different bulk microphysical parameterization schemes using a higher-resolution large eddy simulation (LES) model, the WRF-LES. Initial results show that the Thompson scheme was the only one able to reproduce and sustain a substantial supercooled liquid layer, but it was unable to reproduce the transition from a deep, liquid-rich cloud to a thin layer with moderately and heavily rimed precipitation. This is the first step in linking simulated LES-scale riming processes with those parameterized at a coarser mesoscale model scale. This has important implications for forecasting low-level clouds in an operational environment, given the efficiency of the riming process.

cloud base heights↗

DOE OSTI2025

Atomic Structure, Dynamics, Changes in Chemical Bonding and Semiconductor-Metal Transition in Sb 2 Se 3 : A Remarkable Material for Quantum Networks and Energy Applications

Antimony sesquiselenide has become an outstanding functional material for photovoltaics, energy storage and transformation, memory and photonic applications. Sb 2 Se 3 is one of the most successful emerging solar light absorbers and has also been identified as a highly promising ultralow-loss phase-change material (PCM) for next-generation coherent nanophotonic processors, photonic tensor cores, quantum and neuromorphic networks. Unlike benchmark telluride PCMs, Sb 2 Se 3 features a quasi-one-dimensional (1D) crystalline structure consisting of (Sb 4 Se 6 ) ∞ ribbons, lacks the typical PCM chemical bonding, and undergoes an extended semiconductor-metal transition above the melting point. Consequently, the origin of high optical contrast between crystalline (SET) and amorphous (RESET) logic states remains elusive and presents a significant challenge. Using high-energy X-ray diffraction and Raman spectroscopy over a wide temperature range, supported by first-principles simulations and complemented by thermal, optical and electrical measurements, as well as by 121 Sb-Mossbauer spectroscopy, the quasi-1D network of orthorhombic antimony sesquiselenide was found to undergo significant evolution in amorphous and supercooled Sb 2 Se 3 , leading to lower coordination, shorter interatomic distances and a higher p-electron density on antimony, indicating changes in chemical bonding. The observed novel Sb 2 Se 3 nanocrystalline polymorph, characterized by trigonal antimony coordination and more isolated Sb-Se ribbons, could help reduce multiple trapping defect states in the bandgap, which are typical of orthorhombic Sb 2 Se 3 , thereby enhancing the power-conversion efficiency of photovoltaic devices. Semimetallic and metallic liquid Sb 2 Se 3 exhibit a gradual transformation into a denser 2D and/or 3D network with higher antimony coordination. Localized electron states in the pseudogap are becoming extended, leading to an increase in electronic conductivity σ following the relationship σ ∝ N(E F ) 2 . Liquid Sb 2 Se 3 also appears to be strongly fragile, with a nonmonotonic change in viscosity and higher atomic mobility in the metallic liquid. Furthermore, these results explain extraordinary functionalities of Sb 2 Se 3 for photonic and energy applications.

antimony↗

NASA NTRS1995

Freezing Behavior of Stratospheric Sulfate Aerosols Inferred from Trajectory Studies

Based on the trajectory analysis presented in this paper, a new mechanism is described for the freezing of the stratospheric sulfate aerosols. Temperature histories based on 10-day back trajectories for six ER-2 flights during AASE-I (1989) and AAOE (1987) are presented. The mechanism requires, as an initial step, the cooling of a H2SO4/H2O aerosol to low temperatures. If a cooling cycle is then followed up by a warming to approximately 196-198 K, the aerosols may freeze due to the growth of the crystallizing embryos formed at the colder temperature. The HNO3 absorbed at colder temperatures may increase the nucleation rate of the crystalling embryos and therefore influence the crystallization of the supercooled aerosols upon warming. Of all the ER-2 flights described, only the polar stratospheric clouds (PSC), observed on the flights of January 24, and 25, 1989 are consistent with the thermodynamics of liquid ternary solutions of H2SO4/HNO3/H2O (type Ib PSCs). For those two days, back trajectories indicate that the air mass was exposed to sulfuric acid tetrahydrate (SAT) melting temperatures about 24 hours prior to being sampled by the ER-2. Temperature histories, recent laboratory measurements, and the properties of glassy solids suggest that stratospheric H2SO4 aerosols may undergo a phase transition to SAT upon warming at approximately 198 K after going through a cooling cycle to about 194 K or lower.

Tabazadeh, A.↗

DOE OSTI2024

Zero-Strain Metal-Insulator Transition by the Local Fluctuation of Cation Dimerization

The coupled electronic and structural transitions in metal-insulator transition (MIT) hinder ultrafast switching and ultimate endurance. Decoupling these transitions and achieving a zero-strain electronic MIT can overcome the fundamental limitations of MIT in solid materials. Here, this study demonstrates that iso-valent Ti dopants in supercooled VO 2 epitaxial films cause MIT with minimal hysteresis without changing unit-cell volume and crystal symmetry. The Ti dopants in the VO 2 lattice locally alter the configuration of V-V pairs, where the long-range ordering in V-V pairs is disrupted, and the nano-domains of V-V dimers are formed. Strikingly, these local V-V dimers persist even above the electronic transition temperature (T MI ), facilitating the zero-strain electronic MIT with nanoscale structural heterogeneity. The geometrically compatible interface between insulating and metallic phases drastically enhances switching speed and endurance during electrically and optically driven zero-strain MIT. In conclusion, this discovery offers a fresh perspective on the scientific understanding of MIT and the improved functionality in terms of device speed and reliability by decoupling electronic and structural transitions.

36 MATERIALS SCIENCE↗

NASA NTRS1995

Nucleation and growth of crystals under cirrus and polar stratospheric cloud conditions

Laboratory studies examine phase changes of hygroscopic substances which occur as aerosol in stratosphere and troposphere (sodium chloride, ammonium sulfate, ammonium bisulfate, nitric acid, sulfuric acid), under controlled conditions, in samples volume 1 to 10(exp -4) ml. Crystallization of salts from supersaturated solutions is examined by slowly evaporating a solution drop on a substrate, under controlled relative humidity, until self nucleation occurs; controlled nucleation of ice in a mm capillary U-tube gives a measured ice crystallization velocity at known supercooling. Two states of crystallization occur for regions where hydrates exist. It is inferred that all of the materials readily exist as supersaturated/supercooled solutions; the degree of metastability appears to be slightly enhanced by inclusion of aircraft produced soot. The crystallization velocity is taken as a measure of viscosity. Results suggest an approach to a glass transition at high molality, supersaturation and/or supercooling within the range of atmospheric interest. It is hypothesized that surface reactions occur more readily on solidified particles - either crystalline or glass, whereas volume reactions are more important on droplets with sufficiently low viscosity and volume diffusivity. Implications are examined for optical properties of such particles in the atmosphere. In a separate experiment, crystal growth was examined in a modified thermal vapor diffusion chamber over the range of cirrus temperature (-30 to -70 C) and under controlled supersaturation and air pressure. The crystals grew at a velocity of 1-2 microns/s, thickness 60-70 micron, in the form of thin column crystals. Design criteria are given for a system to investigate particle growth down to -100 C, (PSC temperatures) where nitric acid particles can be grown under similar control and in the form of hydrate crystals.

Hallett, John↗

DOE OSTI2024

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl 2 from T m + 85 K to T m + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl 2 . We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. Furthermore, this result may contribute to new understandings of the structural relaxation of other network liquids.

36 MATERIALS SCIENCE↗

NASA NTRS1999

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗

DOE OSTI2023

Melting curves of ice polymorphs in the vicinity of the liquid–liquid critical point

The possible existence of a liquid–liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164–168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid–liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid–liquid critical point lies between the melting curves of ice III and ice V. Herein we explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid–liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗