SEARCH · Search NASA

Results for “Liquid structure simulation”

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.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Experimentally Driven Automated Machine-Learned Interatomic Potential for a Refractory Oxide

Understanding the structure and properties of refractory oxides is critical for high temperature applications. In this work, a combined experimental and simulation approach uses an automated closed loop via an active learner, which is initialized by x-ray and neutron diffraction measurements, and sequentially improves a machine-learning model until the experimentally predetermined phase space is covered. Furthermore, a multiphase potential is generated for a canonical example of the archetypal refractory oxide, HfO 2 , by drawing a minimum number of training configurations from room temperature to the liquid state at similar to 2900 degrees C. The method significantly reduces model development time and human effort.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Femtosecond x-ray diffraction reveals a liquid–liquid phase transition in phase-change materials

In phase-change memory devices, a material is cycled between glassy and crystalline states. The highly temperature-dependent kinetics of its crystallization process enables application in memory technology, but the transition has not been resolved on an atomic scale. Using femtosecond x-ray diffraction and ab initio computer simulations, we determined the time-dependent pair-correlation function of phase-change materials throughout the melt-quenching and crystallization process. We found a liquid–liquid phase transition in the phase-change materials Ag 4 In 3 Sb 67 Te 26 and Ge 15 Sb 85 at 660 and 610 kelvin, respectively. The transition is predominantly caused by the onset of Peierls distortions, the amplitude of which correlates with an increase of the apparent activation energy of diffusivity. This reveals a relationship between atomic structure and kinetics, enabling a systematic optimization of the memory-switching kinetics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Studying the SBND Cathode Structure

The Short-Baseline Near Detector (SBND) is a Liquid Argon Time Projection Chamber (LArTPC) located at Fermilab. It utilizes electrons ionized by passing particles inside the detector to create 3D reconstructions of interactions. This is achieved by drifting ionized electrons towards the instrumented wire planes. Structurally, the cathode can be mapped using cosmic ray muons by using information of their trajectories and their timing. This study shows that the cathode has rich structure. Though preliminary, this work shows evidence of non-negligible effects and warrants further investigation of these structures which could lead to improved detector simulations and more reliable event reconstructions.

Cataldo, Aedan [Notre Dame U.]↗

DOE OSTI2022

Liquid Core Detection and Strand Condition Monitoring in a Continuous Caster Using Optical Fiber

Real-time monitoring of the liquid core position during the continuous casting of steel has been demonstrated using low-cost distributed optical-fiber-based strain sensors. These sensors were installed on the containment roll support structures in the segments of a production continuous caster to detect the position of the solid–liquid interface and monitor the strand condition during the continuous casting. Distributed Fiber Bragg Grating sensors (FBGs) were used in this work to monitor strain at six roll positions in the caster. The sensor performance was first validated by comparing optical strain measurements with conventional strain gauge measurements in the lab. Next, optical strain measurements were performed on an isolated caster segment in a segment maintenance facility using hydraulic jacks to simulate the presence of a liquid core under the roll. Finally, the sensors were evaluated during caster operation. The sensors successfully detected the load increase associated with the presence of a liquid core under each instrumented roll location. Incidents of bulging and roll eccentricity were also detected using frequency analysis of the optical strain signal. The liquid core position measurements were compared using predictions from computer models (digital twins) in use at the production site. The measurements were in good agreement with the model predictions, with a few exceptions. Under certain transient caster operating conditions, such as spraying practice changes and SEN exchanges, the model predictions deviated slightly from the liquid core position determined from strain measurements.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Identifying synoptic controls on boundary layer thermodynamic and cloud properties in a regional forecast model

Although most of our understanding of boundary layer cloudiness is based on idealized, subtropical, barotropic marine environments, boundary layer clouds exist across a range of conditions. In this study, we use the Naval Research Laboratory's Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) and an automated cold-front-relative analysis framework to explore the boundary layer structure associated with low clouds across a transect through the cold front of a midlatitude synoptic cyclone. The model credibly captures boundary layer structure in line with conceptual models of midlatitude cyclones and ground-based observations at Graciosa Island in the Azores. The warm sector is conditionally unstable, with clouds that are too shallow and with too little liquid water, compared to cloud property retrievals from satellite and surface-based instruments. The cold-frontal region exhibits convection associated with weak stability and ascent. Northwest of the cold front, the boundary layer is well mixed, deeper, and capped by a strong inversion maintained by large-scale subsidence. Simulated clouds in frontal and post-frontal regions are mostly too thick, with too much liquid water and too little cloud base drizzle. The post-frontal clouds are associated with grid-scale updrafts, which appear to be the model's attempt to represent mesoscale organization of cellular convection typically observed in the cold sector of midlatitude cyclones. The deep, well-mixed post-frontal boundary layers and cloudiness are maintained by strong surface fluxes, as in cold air outbreaks. Our analysis framework serves as a unique approach to model verification, and our results offer insights into differences in boundary layer cloud behavior between subtropical and synoptic cold-sector regimes.

Eissner, Jordan M. [Univ. of Kansas, Lawrence, KS ↗

NASA NTRS2013

Flight Performance Feasibility Studies for the Max Launch Abort System

In 2007, the NASA Engineering and Safety Center (NESC) initiated the Max Launch Abort System Project to explore crew escape system concepts designed to be fully encapsulated within an aerodynamic fairing and smoothly integrated onto a launch vehicle. One objective of this design was to develop a more compact launch escape vehicle that eliminated the need for an escape tower, as was used in the Mercury and Apollo escape systems and what is planned for the Orion Multi-Purpose Crew Vehicle (MPCV). The benefits for the launch vehicle of eliminating a tower from the escape vehicle design include lower structural weights, reduced bending moments during atmospheric flight, and a decrease in induced aero-acoustic loads. This paper discusses the development of encapsulated, towerless launch escape vehicle concepts, especially as it pertains to the flight performance and systems analysis trade studies conducted to establish mission feasibility and assess system-level performance. Two different towerless escape vehicle designs are discussed in depth: one with allpropulsive control using liquid attitude control thrusters, and a second employing deployable aft swept grid fins to provide passive stability during coast. Simulation results are presented for a range of nominal and off-nominal escape conditions.

Tarabini, Paul V.↗

DOE OSTI2022

Short- to Intermediate-Range Structure, Transport, and Thermophysical Properties of LiF–NaF–ZrF4 Molten Salts

LiF–NaF–ZrF 4 multicomponent molten salts are identified as promising candidates for coolant salts in molten salt reactors and advanced high-temperature reactors. This study focused on low-melting point salt compositions of interest: 38LiF–51NaF–11ZrF 4 , 42LiF–29NaF–29ZrF 4 , and 26LiF–37NaF–37ZrF 4 . Ab-initio molecular dynamics (AIMD) calculations were performed and compared with available experimental data to assess the ability of rigid ion models (RIM) to reproduce short- to intermediate-range structure, transport, and thermophysical properties of the LiF–NaF–ZrF 4 salt mixtures. It is found that as ZrF 4 mol% increases, the average cation–anion coordination number (CN) of monovalent cations (Li + , Na + ) obtained from RIM calculations decreases, while multivalent Zr 4+ CN varied from 15% to 19% in comparison to corresponding AIMD values. In addition, RIM is found to predict the existence of 7, 8, and 9 coordinated fluorozirconate complexes, while AIMD and the available experimental data showed an occurrence of 6, 7, and 8 coordinated complexes in the melt. The intermediate-range structure analysis revealed that while the RIM parameters are able to reproduce a local structure for lower ZrF 4 mol% salts such as in 38LiF–51NaF–11ZrF 4 , an extensive fluorozirconate network formation is observed in RIM simulations for higher ZrF 4 mol% compositions. The network generated by RIM parameters is found to be mainly connected by “corner-sharing” fluorozirconate complexes as opposed to both “edge-sharing” and “corner-sharing” connectively portrayed by AIMD. It is found that a close agreement between AIMD and the RIM salt structure for the 11-mol% ZrF 4 salt resulted in good agreement in the calculated Zr diffusivities and the viscosity values. However, due to the inaccurate short- to intermediate-range structure prediction by RIM for higher ZrF 4 mol% compositions, thermophysical properties such as densities and heat capacity differ by up to 26% and 27%, respectively, upon comparison with AIMD and experimental values. Also, the network-dominated properties such as diffusion coefficients and viscosities differed by up to two and three orders of magnitude, respectively. This study signifies the importance of accurate salt structure generation for an accurate prediction of transport and thermophysical properties of multicomponent molten salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

Experimental and Theoretical Investigation of the Ion Conduction Mechanism of Tris(adiponitrile)perchloratosodium, a Self-Binding, Melt-Castable Crystalline Sodium Electrolyte

Sodium perchlorate (NaClO4) crystallizes with adiponitrile (ADN) as a 1:3 solvate to produce (ADN) 3 NaClO 4 , a solid electrolyte for sodium ion conduction. The solid possesses high thermal stability (up to 150 °C) and the ability to be melt-cast (T m = 81 °C). The pressed solid has a high ionic conductivity of 2.2 × 10 –4 S cm –1 at room temperature with a low activation barrier for ion conduction of 22 kJ mol –1 . The high conductivity is the result of low-affinity ion-conduction channels in the bulk based on the X-ray crystal structure, and by low grain-boundary resistance and possibly a grain-boundary percolating network due to a fluidlike nanoliquid layer between the grains, observable by scanning electron microscopy and differential scanning calorimetry. When the liquid nanolayer is rinsed away or removed by excessive drying, the bulk room temperature ionic conductivity is 4 × 10 –5 S cm –1 , activation energy for ionic conduction for an organic solid is 37 kJ mol –1 , and the sodium ion transference number is 0.71. Scanning electron microscopy and classical molecular dynamics simulations suggest that these cocrystals form a fluid layer of ADN at the surface, which facilitates the Na+ ion migration between the grains. Finally, density functional theory calculations are consistent with the possibility of ion conduction via a solvent–anion coordinated transition state through vacancy defects in the three symmetry-equivalent ion channels along separate directions, suggesting the possibility of ionic conductivity in three dimensions.

25 ENERGY STORAGE↗

DOE OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Large-Eddy Simulations of Marine Boundary Layer Clouds Associated with Cold-Air Outbreaks during the ACTIVATE Campaign. Part I: Case Setup and Sensitivities to Large-Scale Forcings

Large-eddy simulation (LES) is able to capture key boundary layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LESs. We derive these quantities from measurements and reanalysis obtained for two cold-air outbreak (CAO) events during Phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE) in February–March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from ERA5 data agree reasonably well with those derived from ACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time-evolving heat fluxes, wind, and advective tendencies profiles from ERA5 data to drive the LES. We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LESs. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path simulated in LES agrees reasonably well with the ACTIVATE measurements. This study paves the path to investigate aerosol–cloud–meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Lattice-scale variations in viscosity are correlated with solution structure at mineral-water interfaces

At solid-liquid interfaces, the viscosity increases markedly from the bulk due to the collective interactions of ions and water molecules, influencing phenomena relevant to nanofluidics, colloidal dynamics, and electrochemistry. Here, in this study, we investigated dissipative forces at the boehmite-water interface using 3D atomic force microscopy. We observed an increase in interfacial solution viscosity, η, by 10-100-fold as the nanoprobe approached the surface in normal direction, with up to four oscillatory features showing average peaks of η/η bulk = 44-71. Moreover, the viscosity showed sub-nanometer variations within 0.5 nm from the interface, templated by the underlying crystal lattice and correlated with the interfacial solution structure. Beyond a near-wall region of approximately 1.2 nm, the dissipative response was comparable to that in bulk solution. Molecular dynamics simulations, along with statistical mechanical analyses, provided details on hydrodynamic structures near the interface. Specifically, the lattice-dependent dissipative responses are correlated with extensive hydrogen bonding by interfacial water molecules, which increased friction, particularly along the [001] direction. These results demonstrate how solution viscosity at mineral-water interfaces is anisotropic and correlated with the local solution structure, providing insights into the dynamics of nanocrystal attachment.

Viscosity↗

NASA NTRS1981

Design and fabrication of brazed Rene 41 honeycomb sandwich structural panels for advanced space transportation systems

The design and fabrication of two large brazed Rene 41 honeycomb panels, the establishment of a test plan, the design and fabrication of a test fixture to subject the panels to cyclic thermal gradients and mechanical loads equivalent to those imposed on an advanced space transportation vehicle during its boost and entry trajectories are discussed. The panels will be supported at four points, creating three spans. The outer spans are 45.7 cm (18 in.) and the center span 76.2 cm (30 in). Specimen width is 30.5 cm (12 in.). The panels were primarily designed by boost conditions simulated by subjecting the panels to liquid nitrogen, 77K (-320 F) on one side and 455K (360 F) on the other side and by mechanically imposing loads representing vehicle fuel pressure loads. Entry conditions were simulated by radiant heating to 1034K (1400 F). The test program subjected the panels to 500 boost thermal conditions. Results are presented.

Hepler, A. K.↗

NASA NTRS2010

MCNP Simulations of Measurement of Insulation Compaction in the Cryogenic Rocket Fuel Tanks at Kennedy Space Center by Fast/Thermal Neutron Techniques

MCNP simulations have been run to evaluate the feasibility of using a combination of fast and thermal neutrons as a nondestructive method to measure of the compaction of the perlite insulation in the liquid hydrogen and oxygen cryogenic storage tanks at John F. Kennedy Space Center (KSC). Perlite is a feldspathic volcanic rock made up of the major elements Si, AI, Na, K and 0 along with some water. When heated it expands from four to twenty times its original volume which makes it very useful for thermal insulation. The cryogenic tanks at Kennedy Space Center are spherical with outer diameters of 69-70 feet and lined with a layer of expanded perlite with thicknesses on the order of 120 cm. There is evidence that some of the perlite has compacted over time since the tanks were built 1965, affecting the thermal properties and possibly also the structural integrity of the tanks. With commercially available portable neutron generators it is possible to produce simultaneously fluxes of neutrons in two energy ranges: fast (14 Me V) and thermal (25 me V). The two energy ranges produce complementary information. Fast neutrons produce gamma rays by inelastic scattering, which is sensitive to Fe and O. Thermal neutrons produce gamma rays by prompt gamma neutron activation (PGNA) and this is sensitive to Si, Al, Na, K and H. The compaction of the perlite can be measured by the change in gamma ray signal strength which is proportional to the atomic number densities of the constituent elements. The MCNP simulations were made to determine the magnitude of this change. The tank wall was approximated by a I-dimensional slab geometry with an 11/16" outer carbon steel wall, an inner stainless wall and 120 cm thick perlite zone. Runs were made for cases with expanded perlite, compacted perlite or with various void fractions. Runs were also made to simulate the effect of adding a moderator. Tallies were made for decay-time analysis from t=0 to 10 ms; total detected gamma-rays; detected gamma-rays from thermal neutron reactions d. detected gamma-rays from non-thermal neutron reactions and total detected gamma-rays as a function of depth into the annulus volume. These indicated a number of possible independent metrics of perlite compaction. For example the count rate for perlite elements increased from 3600 to 8500 cps for an increase in perlite density from 6 lbs/lcf to 16.5 lbs/cf. Thus the MCNP simulations have confirmed the feasibility of using neutron methods to map the compaction of perlite in the walls of the cryogenic tanks.

Livingston, R. A.↗

DOE OSTI2020

Dehydration entropy drives liquid-liquid phase separation by molecular crowding

Complex coacervation driven liquid-liquid phase separation (LLPS) of biopolymers has been attracting attention as a novel phase in living cells. Studies of LLPS in this context are typically of proteins harboring chemical and structural complexity, leaving unclear which properties are fundamental to complex coacervation versus protein-specific. This study focuses on the role of polyethylene glycol (PEG)—a widely used molecular crowder—in LLPS. Significantly, entropy-driven LLPS is recapitulated with charged polymers lacking hydrophobicity and sequence complexity, and its propensity dramatically enhanced by PEG. Experimental and field-theoretic simulation results are consistent with PEG driving LLPS by dehydration of polymers, and show that PEG exerts its effect without partitioning into the dense coacervate phase. Furthermore, it is then up to biology to impose additional variations of functional significance to the LLPS of biological systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Dynamic phases, stratification, laning, and pattern formation for driven bidisperse disk systems in the presence of quenched disorder

Using numerical simulations, we examine the dynamics of driven two-dimensional bidisperse disks flowing over quenched disorder. The system exhibits a series of distinct dynamical phases as a function of applied driving force and packing fraction, including a phase-separated state as well as a smectic state with liquid-like or polycrystalline features. At low driving forces, we find a clogged phase with an isotropic density distribution, while at intermediate driving forces the disks separate into bands of high and low densities with either liquid-like or polycrystalline structure in the high-density bands. In addition to the density phase separation, we find that in some cases there is a fractionation of the disk species, particularly when the disk size ratio is large. The species phase-separated regimes form a variety of patterns such as large disks separated by chains of smaller disks. Here, our results show that the formation of laning states can be enhanced by tuning the ratio of disk radius of the two species, due to the clumping of small disks in the interstitial regions between the large disks. This system could be experimentally realized using sterically interacting colloidal particles suspended in a viscous fluid driven over random pinning arrays or granular matter suspended in fluid moving over a random landscape.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Structural insights into hydrolytic defluorination of difluoroacetate by microbial fluoroacetate dehalogenases

Fluorine forms the strongest single bond to carbon with the highest bond dissociation energy among natural products. However, fluoroacetate dehalogenases (FADs) have been shown to hydrolyze this bond in fluoroacetate under mild reaction conditions. Furthermore, two recent studies demonstrated that the FAD RPA1163 fromRhodopseudomonas palustriscan also accept bulkier substrates. In this study, we explored the substrate promiscuity of microbial FADs and their ability to defluorinate polyfluorinated organic acids. Enzymatic screening of eight purified dehalogenases with reported fluoroacetate defluorination activity revealed significant hydrolytic activity against difluoroacetate in three proteins. Product analysis using liquid chromatography‐mass spectrometry identified glyoxylic acid as the final product of enzymatic DFA defluorination. The crystal structures of DAR3835 fromDechloromonas aromaticaand NOS0089 fromNostocsp. were determined in the apo‐state along with the DAR3835 H274N glycolyl intermediate. Structure‐based site‐directed mutagenesis of DAR3835 demonstrated a key role for the catalytic triad and other active site residues in the defluorination of both fluoroacetate and difluoroacetate. Computational analysis of the dimer structures of DAR3835, NOS0089, and RPA1163 indicated the presence of one substrate access tunnel in each protomer. Moreover, protein‐ligand docking simulations suggested similar catalytic mechanisms for the defluorination of both fluoroacetate and difluoroacetate, with difluoroacetate being defluorinated via two consecutive defluorination reactions producing glyoxylate as the final product. Thus, our findings provide molecular insights into substrate promiscuity and catalytic mechanism of FADs, which are promising biocatalysts for applications in synthetic chemistry and bioremediation of fluorochemicals.

Biochemistry & Molecular Biology↗

NASA NTRSDate not supplied

Spacesuit Sizing by Multivariate Simulation: Case Study for Liquid Cooling and Ventilation Garment

A technique was developed to assess the fit of spacesuit garment liquid cooling and ventilation garment (LCVG). A fit scoring scheme was built from the spacesuit sizing requirements and guidelines, and the accommodated proportion of the wearer population was estimated from statistically synthesized pseudo population data. The method evaluated different scenarios of LCVG sizing systems including 1) baseline 7-size system, 2) 7-sizes with target body dimensions iteratively adjusted for maximum fit, and 3) 9-sizes with iterative adjustments. Overall, the baseline sizing system showed 81% accommodation rate, while the hypothetical new systems in-tegrating iterative adjustments for 7- and 9-sizes predictively indicated that accommodation could improve up to 92 and 97%. The proposed method provided a structured metric to quantify fit performance of a garment sizing system, which can be further applied to different spacesuit component, PPE and garment designs.

Han Kim↗

DOE OSTI2019

High-Level VSCF/VCI Calculations Decode the Vibrational Spectrum of the Aqueous Proton

We report that the hydrated excess proton is a common species in aqueous chemistry which complexes with water in a variety of structures. The infrared spectrum of the aqueous proton is particularly sensitive to this array of structures, which manifests as continuous IR absorption from 1000-3000 cm -1 known as the “proton continuum”. Because of the extreme breadth of the continuum and strong anharmonicity of the involved vibrational modes, this spectrum has eluded straightforward interpretation and simulation. Using protonated water hexamer clusters from reactive molecular dynamics trajectories, and focusing on their central H + (H 2 O) 2 structures’ spectral contribution, we reproduce the linear IR spectrum of the aqueous proton with a high-level local- monomer quantum method and highly accurate many-body potential energy surface. The accuracy of this approach is first verified in the vibrational spectra of the two isomers of the protonated water hexamer in the gas phase. We then apply this approach to 800 H + (H 2 O) 6 clusters, also written as [H + (H 2 O) 2 ](H 2 O) 4 , drawn from MS-EVB simulations of the bulk liquid to calculate the infrared spectrum of the aqueous proton complex. Incorporation of anharmonic effects to the vibrational potential and quantum mechanical treatment of the proton produce better agreement to the infrared spectrum compared to the double-harmonic approximation. We assess the correlation of the pro- ton stretching mode with different atomistic coordinates, finding the best correlation with $\langle$ROH$\rangle$, the expectation value of the proton-oxygen distance ROH. Finally, we also decompose the IR spectrum based on normal mode vibrations and $\langle$ROH$\rangle$ to provide insight on how different frequency regions in the continuum report on different configurations, vibrational modes, and mode couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗