Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical molecular dynamics”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

SANE: strategic autonomous non-smooth exploration for multiple optima discovery in multi-modal and non-differentiable black-box functions

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and multimodal parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, material structure image spaces, and molecular embedding spaces. Often these systems are black-boxes and time-consuming to evaluate, which resulted in strong interest towards active learning methods such as Bayesian optimization (BO). However, these systems are often noisy which make the black box function severely multi-modal and non-differentiable, where a vanilla BO can get overly focused near a single or faux optimum, deviating from the broader goal of scientific discovery. To address these limitations, here we developed Strategic Autonomous Non-Smooth Exploration (SANE) to facilitate an intelligent Bayesian optimized navigation with a proposed cost-driven probabilistic acquisition function to find multiple global and local optimal regions, avoiding the tendency to becoming trapped in a single optimum. To distinguish between a true and false optimal region due to noisy experimental measurements, a human (domain) knowledge driven dynamic surrogate gate is integrated with SANE. We implemented the gate-SANE into pre-acquired piezoresponse spectroscopy data of a ferroelectric combinatorial library with high noise levels in specific regions, and piezoresponse force microscopy (PFM) hyperspectral data. SANE demonstrated better performance than classical BO to facilitate the exploration of multiple optimal regions and thereby prioritized learning with higher coverage of scientific values in autonomous experiments. Our work showcases the potential application of this method to real-world experiments, where such combined strategic and human intervening approaches can be critical to unlocking new discoveries in autonomous research.

Biswas, Arpan [University of Tennessee, Knoxville,↗
DOE OSTI2023

Reduced Order Model to Predict Dispersion of Flammable Refrigerant into a Space

As the HVAC&R industry mobilizes to deploy more low-GWP refrigerants, relevant standards are being continually reviewed and updated. Those include the general safety standards ISO 5149 and ASHRAE 15, and the equipment standards IEC and UL. The standards systematically set the allowable maximum amount of refrigerant that should be used in different equipment types and different applications. To do so, they rely on predictions of how a leaked refrigerant mass will disperse into a space. Dispersion characteristics, such as total flammable volume and its residence time, determine the risk associated with the presence of the flammable refrigerant. The standards have included provisions for the use of flammable refrigerants for approximately two decades. They relied on limited analytical analyses and test cases in their development. Dispersion of a refrigerant into a space is complex. Computational fluid dynamics (CFD) are the most accurate in predicting a given problem. However, CFD is computationally expensive and requires specialized expertise and resources and is not suitable for use by standards development working group as prediction tool. This paper presents the development of a reduced order model (ROM) that predicts the key dispersion characteristics relevant to the dispersion of a leaked refrigerant into a space for any combination of input variables. The inputs are the refrigerant release height, the total released refrigerant mass and its release flow rate, the refrigerant molecular weight, the ventilation flow rate, the floor area and height of the space, recirculation air flow rate, and the tightness of the space. The outputs are histograms of volume fraction of the room in prescribed concentration bins and the total mass of the refrigerant in each bin normalized by the total refrigerant charge at 13 prescribed simulation time stamps between 1 and 900 seconds. The ROM is constructed from a set of CFD simulations with carefully chosen combinations of input parameters. The selection if done using a multidimensional sparse grid which is a generalization of the classical tensor approach but offers additional flexibility and thus can be more carefully tuned towards a specific model. The tuning is done to improve the accuracy, measured in the difference between the output values of the ROM and the CFD model, while minimizing the computational cost, measured in number of CFD simulations which is orders of magnitude more expensive than the processing the training data.

Edwards, Dean↗
DOE OSTI2023

Quantum spin ice response to a magnetic field in the dipole-octupole pyrochlore Ce 2 Zr 2 O 7

The pyrochlore magnet Ce 2 Zr 2 O 7 has attracted much attention as a quantum spin ice candidate whose novelty derives in part from the dipolar-octupolar nature of the Ce 3+ pseudospin-1/2 degrees of freedom it possesses. We report heat capacity measurements on single crystal samples of Ce 2 Zr 2 O 7 down to T~0.1K in a magnetic field along the [$1,\bar{1},0$] direction. These measurements show that the broad hump in the zero-field heat capacity moves higher in temperature with increasing field strength and is split into two separate humps by the [$1,\bar{1},0$] magnetic field at ~2T. These separate features are due to the decomposition of the pyrochlore lattice into effectively decoupled chains for fields in this direction: One set of chains (α chains) is polarized by the field while the other (β chains) remains free. This situation is similar to that observed in the classical spin ices Ho 2 Ti 2 O 7 and Dy 2 Ti 2 O 7 , but with the twist that here the strong transverse exchange interactions produce substantial quantum effects. Our theoretical modeling suggests that the β chains are close to a critical state, with nearly-gapless excitations. We also report elastic and inelastic neutron scattering measurements on single crystal Ce 2 Zr 2 O 7 in [$1,\bar{1},0$] and [0,0,1] magnetic fields at temperatures down to T=0.03K. The elastic scattering behaves consistently with the formation of independent chains for a [$1,\bar{1},0$] field, while the [0,0,1] field produces a single field-induced elastic magnetic Bragg peak at (0,2,0) and equivalent wavevectors, indicating a polarized spin ice state for fields above ~3T. For both [$1,\bar{1},0$] and [0,0,1] magnetic fields, our inelastic neutron scattering results show an approximately dispersionless continuum of scattering that increases in both energy and intensity with increasing field strength. By modeling the complete set of experimental data using numerical linked cluster and semiclassical molecular dynamics calculations, we demonstrate the dominantly multipolar nature of the exchange interactions in Ce 2 Zr 2 O 7 and the smallness of the parameter θ, which controls the mixing between dipolar and octupolar degrees of freedom. Finally, these results support previous estimates of the microscopic exchange parameters and place strong constraints on the theoretical description of this prominent spin ice candidate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Singlet fission in core-linked terrylenediimide dimers

We have studied two regioisomeric terrylenediimide (TDI) dimers in which the 1-positions of two TDIs are linked via 1,3- or 1,4-phenylene spacers, mTDI 2 and pTDI, respectively. The nature and the dynamics of the multiexciton state are tuned by altering the through-bond electronic couplings in the ground and excited states and by changing the solvent environment. Our results show that controlling the electronic coupling between the two chromophores by an appropriate choice of linker can result in independent triplet state formation, even though the initial correlated triplet pair state is confined to a dimer. Moreover, even in polar solvents, if the electronic coupling is strong, the correlated triplet pair state is observed prior to symmetry-breaking charge separation. These results point out the close relationship between the singlet, correlated triplet pair, and charge transfer states in molecular dimers.

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

Application of machine learning in the determination of impact parameter in the 132 Sn+ 124 Sn system

Here, 132 Sn + 124 Sn collisions at a beam energy of 270 MeV/nucleon were performed at the Radioactive Isotope Beam Factory (RIBF) in RIKEN to investigate the nuclear equation of state. Reconstructing the impact parameter is one of the important tasks in the experiment as it relates to many observable. In this work, we employ three commonly used algorithms in machine learning, the artificial neural network (ANN), the convolutional neural network (CNN), and the light gradient boosting machine (LightGBM), to determine the impact parameter by analyzing either the charged particle spectra or several features simulated with events from the ultrarelativistic quantum molecular dynamics (UrQMD) model. To closely imitate experimental data and investigate the generalizability of the trained machine learning algorithms, incompressibility of nuclear equation of state and the in-medium nucleon-nucleon cross sections are varied in the UrQMD model to generate the training data. The mean absolute error Δb between the true and the predicted impact parameter is smaller than 0.45 fm if training and testing sets are sampled from the UrQMD model with the same parameter set. However, if training and testing sets are sampled with different parameter sets, Δb would increase to 0.8 fm. The generalizability of the trained machine learning algorithms suggests that these machine learning algorithms can be used reliably to reconstruct the impact parameter in experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Advanced Microscopy and Spectroscopy for Probing and Optimizing Electrode-Electrolyte Interphases in High Energy Lithium Batteries (Final Scientific/Technical Report)

The earliest developed LiCoO 2 layered oxide cathode material sparked the development of other layered cathode materials, dominating the positive electrode materials for lithium ion batteries. Within the practical operating conditions of today, the current generation layered oxide materials do not meet the future energy storage demands of 350 Wh kg -1 at cell level. Li-rich layered oxide (LRLO) materials have the potential to meet the high energy demands. Unlike the classical layered oxides, LRLO materials exhibit capacities that go beyond conventional topotactic mechanistic theoretical values. Despite its high capacities, this material has several challenges (voltage fading, structural instability, sluggish kinetics, cathode electrolyte interphase instability, etc.) that must be overcome in order to reach commercialization. In the past five years, our research group has made great progress on developing advanced characterization techniques coupled with atomic scale modeling to properly characterize the dynamic phenomena that govern the performance limitations of LRLO materials. Furthermore, our efforts have guided the material synthesis and surface modification to improve capacity retention. By carefully controlling the oxygen activities through the creation of uniform oxygen vacancies, we were able to avoid structural decomposition in LRLO. Our materials deliver a discharge capacity as high as 306 mAh g -1 with an initial coulombic efficiency of 90.6%. Furthermore, they do not exhibit obvious capacity decay with a reversible capacity over 300 mAh g -1 after 100 cycles at 0.1 C-rate. Our approach demonstrates the critical needs for advanced diagnosis and characterization. It is through the in-depth understanding of these high voltage cathode materials at atomistic and molecular level and their dynamic changes during the operation of batteries; we can successfully formulate strategies to further optimize this class of cathode materials, especially for the voltage instabilities. The diagnostic tools developed here can also be leveraged to study anode materials such as Li metal anode. The challenge of probing Li metal is the low dose tolerance (high beam sensitivity) of Li metal. We will have to apply cryogenic method (low-temperature) and low-dose electron microscopy and specialized camera to enable the characterization of the lithium metal.

25 ENERGY STORAGE↗
DOE OSTI2021

Uranium (VI) Adsorbate Structures on Portlandite [Ca(OH)2] Type Surfaces Determined by Computational Modelling and X-ray Absorption Spectroscopy

Portlandite [Ca(OH) 2 ] is a potentially dominant solid phase in the high pH fluids expected within the cementitious engineered barriers of Geological Disposal Facilities (GDF). This study combined X-ray Absorption Spectroscopy with computational modelling in order to provide atomic-scale data which improves our understanding of how a critically important radionuclide (U) will be adsorbed onto this phase under conditions relevant to a GDF environment. Such data are fundamental for predicting radionuclide mass transfer. Surface coordination chemistry and speciation of uranium with portlandite [Ca(OH) 2 ] under alkaline groundwater conditions (ca. pH 12) were determined by both in situ and ex situ grazing incidence extended X-ray absorption fine structure analysis (EXAFS) and by computational modelling at the atomic level. Free energies of sorption of aqueous uranyl hydroxides, [UO 2 (OH) n ] 2–n (n = 0–5) with the (001), (100) and (203) or (101) surfaces of portlandite are predicted from the potential of mean force using classical molecular umbrella sampling simulation methods and the structural interactions are further explored using fully periodic density functional theory computations. Although uranyl is predicted to only weakly adsorb to the (001) and (100) clean surfaces, there should be significantly stronger interactions with the (203/101) surface or at hydroxyl vacancies, both prevalent under groundwater conditions. The uranyl surface complex is typically found to include four equatorially coordinated hydroxyl ligands, forming an inner-sphere sorbate by direct interaction of a uranyl oxygen with surface calcium ions in both the (001) and (203/101) cases. In contrast, on the (100) surface, uranyl is sorbed with its axis more parallel to the surface plane. The EXAFS data are largely consistent with a surface structural layer or film similar to calcium uranate, but also show distinct uranyl characteristics, with the uranyl ion exhibiting the classic dioxygenyl oxygens at 1.8 Å and between four and five equatorial oxygen atoms at distances between 2.28 and 2.35 Å from the central U absorber. These experimental data are wholly consistent with the adsorbate configuration predicted by the computational models. These findings suggest that, under the strongly alkaline conditions of a cementitious backfill engineered barrier, there would be significant uptake of uranyl by portlandite to inhibit the mobility of U(VI) from the near field of a geological disposal facility.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Para-particle oscillator simulations on a trapped-ion quantum computer

Deformed oscillators allow for a generalization of the standard fermions and bosons, namely, for the description of para-particles. Such particles remain hypothetical and unobserved in nature; yet, they can model physical phenomena, such as topological phases of matter. Here, we report the digital quantum simulation of para-particle oscillators by mapping para-particle states to the state of a qubit register, which allows us to identify the para-particle oscillator Hamiltonian as an XY model and further digitize the system onto a universal set of gates. In both instances, the gate depth grows polynomially with the number of qubits used. To establish the validity of our results, we experimentally simulate the dynamics of para-fermions and para-bosons, demonstrating full control of para-particle oscillators on a quantum computer. Furthermore, we compare the overall performance of the digital simulation of dynamics of the driven para-Fermi oscillator to a recent analog quantum simulation result.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The Role of Cation Coordination in the Electrical and Optical Properties of Amorphous Transparent Conducting Oxides

Amorphous oxide semiconductor materials have demonstrated numerous advantages without compromise of electrical properties as compared to their crystalline counterparts, yet understanding of the fundamental principles allowing this has remained elusive. To study the origins of enhanced optoelectronic properties, we apply high-throughput, combinatorial sputtering, structural and spectral mapping, and computationally intensive ab initio molecular dynamics simulations with density functional theory to a ternary, post-transition metal oxide system, namely, zinc tin oxide. The deposited thin films exhibit a high figure of merit, achieving carrier densities in the range of 1019 to 1020 cm–3 and carrier mobilities up to 35 cm2/Vs. These results highlight the role of local distortions and cation coordination in determining the microscopic origins of carrier generation and transport. In particular, we identify the strong likelihood of Sn undercoordination in both Zn-poor and Zn-rich phases leading to the high carrier concentrations observed. This not only diverges from the still widespread historical indictment of oxygen vacancies controlling carrier population in crystalline oxides but also provides a comprehensive framework to describe the unique structure–property relationships using specific structural and electronic descriptors in disordered phase materials.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗
DOE OSTI2020

Selectivity of nitrate and chloride ions in microporous carbons: the role of anisotropic hydration and applied potentials

Understanding ion transport in porous carbons is critical for a wide range of technologies, including supercapacitors and capacitive deionization for water desalination, yet many details remain poorly understood. For instance, an atomistic understanding of how ion selectivity is influenced by the molecular shape of ions, morphology of the micropores and applied voltages is largely lacking. In this work, we combined molecular dynamics simulations with enhanced sampling methods to elucidate the mechanism of nitrate and chloride selectivity in subnanometer graphene slit-pores. We demonstrate that nitrate is preferentially adsorbed over chloride in the slit-like micropores. This preferential adsorption was found to stem from the weaker hydration energy and unique anisotropy of the ion solvation of nitrate. Beside the effects of ion dehydration, we found that applied potential plays an important role in determining the ion selectivity, leading to a lower selectivity of nitrate over chloride at a high applied potential. We conclude that the measured ion selectivity results from a complex interplay between voltage, confinement, and specific ion effects-including ion shape and local hydration structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicted thermophysical properties of UN, PuN and (U,Pu)N

Molecular dynamics and density functional theory simulations are used to predict the lattice and electronic contributions of thermophysical properties for UN, PuN, and mixed (U,Pu)N systems. The properties predicted include the lattice parameter, linear thermal expansion, enthalpy, and specific heat capacity, as a function of temperature. The simulation predictions for high temperature specific heat capacity are compared against experimental measurements to understand the behavior, and why differences in the experimental measurements are observed. The influence of adding U vacancies, N interstitials, and Pu to UN is also examined. For this, a new PuN potential parameter set is developed and used with the Kocevski UN potential, enabling the dynamics of mixed (U,Pu)N systems to be studied. How defects impact the thermophysical properties is important for understanding fuel behavior under different reactor conditions, and these mechanistic predictions can be used to support fuel performance codes where data is scarce.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Comparative analysis of characteristics of a supersonic cw chemical HF laser on molecular fluorine and nitrogen trifluoride

We have experimentally studied the characteristics of a supersonic cw chemical HF laser with a flat nozzle array based on a nozzle – nozzle reagent mixing scheme operating with the use of F{sub 2} – D{sub 2} – He and NF{sub 3} – D{sub 2} – He fuel compositions in an atomic fluorine generator. The active medium flow field pattern, its gas-dynamic characteristics, and the laser radiation spectrum are considered. The comparison of these characteristics has made it possible to establish a strong impact on the flow field (lasing zone length) of the chemical composition of the combustion products from the atomic fluorine generator. The gas-dynamic parameters of the active medium formed by using a NF{sub 3} – D{sub 2} – He fuel have turned out to be less acceptable in terms of the gas flow pressure recovery than in the case of F{sub 2} – D{sub 2} – He fuel, while the radiation spectrum was more preferable from the viewpoint of passing through the atmosphere. (paper)

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular Modes of Attosecond Charge Migration

Here, first-principles calculations are employed to elucidate the modes of attosecond charge migration (CM) in halogenated hydrocarbon chains. We use constrained density functional theory (DFT) to emulate the creation of a localized hole on the halogen and follow the subsequent dynamics via time-dependent DFT. We find low-frequency CM modes (~ 1 eV) that propagate across the molecule and study their dependence on length, bond order, and halogenation. We observe that the CM speed (~ 4Å/ fs) is largely independent of molecule length, but is lower for triple-bonded versus double-bonded molecules. Additionally, as the halogen mass increases, the hole travels in a more particlelike manner as it moves across the molecule. These heuristics will be useful in identifying molecules and optimal CM detection methods for future experiments, especially for halogenated hydrocarbons which are promising targets for ionization-triggered CM.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Crystallization of binary nanocrystal superlattices and the relevance of short-range attraction

The synthesis of binary nanocrystal superlattices (BNSLs) enables the targeted integration of orthogonal physical properties, such as photoluminescence and magnetism, into a single superstructure, unlocking a vast design space for multifunctional materials. However, the formation mechanism of BNSLs remains poorly understood, restricting the prediction of the structure and properties of superlattices. Here we use a combination of in situ scattering and molecular simulation to elucidate the self-assembly of two common BNSLs (AlB 2 and NaZn 13 ) through emulsion templating. Our self-assembly experiments reveal that no intermediate structures precede the formation of the final binary phases, indicating that their formation proceeds through classical nucleation. Using simulations, we find that, despite the formation of AlB 2 and NaZn 13 typically being attributed to entropy, their self-assembly is most consistent with the nanocrystals possessing short-range interparticle attraction, which we find can accelerate nucleation kinetics in BNSLs. Here, we also find homogeneous, classical nucleation in simulations, corroborating our experiments. These results establish a robust correspondence between experiment and theory, paving the way towards prediction of BNSLs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Biomimetic Light Harvesting Complexes Based on Self-Assembled Dye-DNA Nanostructures (Final Technical Report)

A fundamental limitation in the development of artificial molecular light harvesting structures is that generally if one places pigment molecules close enough to achieve very high energy transfer efficiency and absorbance cross section, the dyes interact in such a way that self-quenching occurs. This limits both the optical cross section and the overall efficiency of light capture in such systems. Nature overcomes this problem and achieves efficient energy transfer without self-quenching by controlling the environment, vibration characteristics, orientation and dynamics of the pigments involved in light harvesting. Natural light harvesting systems have evolved highly intricate chromophore architectures of protein-bound organic molecules that absorb and funnel solar energy with high efficiency and fast transfer rates using elegant excitonic circuitries. Nature achieves this precise control over the properties of pigment complexes by using specific proteins as templates to guide the formation of very well-defined pigment aggregates, groups of pigment molecules that are in close enough contact to undergo excitonic interactions resulting in new electronic properties that were not present in the monomer. In synthetic systems, pigment aggregates are well known, and J-aggregates and H-aggregates have been long studied for their unique properties. In Nature, the chlorosomes of green bacteria form what appear to be classical J-aggregates and the light harvesting complexes of purple nonsulfur bacteria form ring structures that have many of the properties of J-aggregates such as a red-shifted absorbance spectrum and excited states that are delocalized over multiple pigment molecules. These photosynthetic systems leverage the spectral and enhanced transport properties of J-like aggregates to achieve fast exciton transport within the antenna and thus high light-harvesting efficiency. The close interaction of the pigments results in substantial electronic coupling that lifts the degeneracy of excitonic states and therefore modifies the absorption and emission spectra of the monomer. In the case of J-aggregation, lower energy excitonic states acquire significant oscillator strength, resulting in a characteristic bathochromic shift of the absorption and fluorescence and often superradiant delocalized exciton states.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Spectrally filtered ps–ns emission dynamics of atmospheric-pressure nanosecond pulsed plasmas

The temporal dynamics of atmospheric-pressure nanosecond pulsed plasma discharges in a pin-to-pin electrode configuration are studied using streak-camera line imaging of the interelectrode gap with a time resolution as short as ~25 ps. Discharge emission initiates homogeneously throughout the interelectrode gap with no detectable streamer propagation and then temporally decays in two distinct phases. Plasma emission bands attributed to various electronic transitions are tracked for single discharges in air and N 2 . Spectral filtering of the excited molecular states reveals that the N 2 (C–B) and N 2 (B–A) emission bands evolve in distinct early and late phases, respectively, with a time separation of ~15–20 ns. Furthermore, significant differences in the temporal dynamics of plasma discharges in air and N2 are observed. High levels of excited-state atomic oxygen and NO appear after the initial decay of the N 2 (C) state and coincide primarily with the latter phases of plasma evolution in air environments. From temporal traces of discharge emission, the formation and relaxation timescales of the electronically excited states of N 2 are quantified in pure N2 and air environments with sub-nanosecond resolution. The streak-OES (optical emission spectroscopy) technique enables quantitative time-resolved studies of key chemical species for model validation in ultra-short-pulsed plasmas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Isotope effects in supercooled H 2 O and D 2 O and a corresponding-states-like rescaling of the temperature and pressure

Water shows anomalous properties that are enhanced upon supercooling. The unusual behavior is observed in both H 2 O and D 2 O, however, with different temperature dependences for the two isotopes. It is often noted that comparing the properties of the isotopes at two different temperatures (i.e., a temperature shift) approximately accounts for many of the observations—with a temperature shift of 7.2 K in the temperature of maximum density being the most well-known example. However, the physical justification for such a shift is unclear. Motivated by recent work demonstrating a “corresponding-states-like” rescaling for water properties in three classical water models that all exhibit a liquid–liquid transition and critical point, the applicability of this approach for reconciling the differences in the temperature- and pressure-dependent thermodynamic properties of H 2 O and D 2 O is investigated here. Utilizing previously published data and equations-of-state for H 2 O and D 2 O, we show that the available data and models for these isotopes are consistent with such a low temperature correspondence. These observations provide support for the hypothesis that a liquid–liquid critical point, which is predicted to occur at low temperatures and high pressures, is the origin of many of water’s anomalies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗