Search NASA⌕ Search

SEARCH · Search NASA

Results for “Isobaric”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Beam-spin asymmetry Σ for Σ – hyperon photoproduction off the neutron

We report a new measurement of the beam-spin asymmetry $\boldsymbol{\Sigma}$ for the $\vec{\gamma} n \rightarrow K^+\Sigma^-$ reaction using quasi-free neutrons in a liquid-deuterium target. The new dataset includes data at previously unmeasured photon energy and angular ranges, thereby providing new constraints on partial wave analyses used to extract properties of the excited nucleon states. The experimental data were obtained using the CEBAF Large Acceptance Spectrometer (CLAS), housed in Hall B of the Thomas Jefferson National Accelerator Facility (JLab). The CLAS detector measured reaction products from a liquid-deuterium target produced by an energy-tagged, linearly polarised photon beam with energies in the range 1.1 to 2.3 GeV. Predictions from an isobar model indicate strong sensitivity to $N(1720)3/2^+$, $\Delta(1900)1/2^-$, and $N(1895)1/2^-$, with the latter being a state not considered in previous photoproduction analyses. When our data are incorporated in the fits of partial-wave analyses, one observes significant changes in $\gamma$-$n$ couplings of the resonances which have small branching ratios to the $\pi N$ channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

57 Zn β-delayed proton emission establishes the 56 Ni rp-process waiting point bypass

We measured the 57 Zn β-delayed proton (βp) and γemission at the National Superconducting Cyclotron Laboratory. We find a 57 Zn half-life of 43.6 ± 0.2ms, βp branching ratio of (84.7 ±1.4)%, and identify four transitions corresponding to the exotic β-γ-pdecay mode, the second such identification in the fp-shell. The p/γratio was used to correct for isospin mixing while determining the 57 Zn mass via the isobaric multiplet mass equation. Previously, it was uncertain as to whether the rp-process flow could bypass the textbook waiting point 56 Ni for astrophysical conditions relevant to Type-I X-ray bursts. Our results definitively establish the existence of the 56 Ni bypass, with 14-17% of the rp-process flow taking this route.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A pathway to unveiling neutrinoless ββ decay nuclear matrix elements via γγ decay

We investigate the experimental feasibility of detecting second-order double-magnetic dipole (γγ-M1M1) decays from double isobaric analog states (DIAS), which have recently been found to be strongly correlated with the nuclear matrix elements of neutrinoless ββ decay. Using the nuclear shell model, we compute theoretical branching ratios for γγ-M1M1 decays and compare them with other competing processes, such as single-γ decay and proton emission, which represent the dominant decay channels. We also estimate the potential competition from internal conversion and internal pair creation, which can influence the decay dynamics. Additionally, we propose an experimental strategy based on using LaBr3 scintillators to identify γγ-M1M1 transitions from the DIAS amidst the background of the competing processes. Our approach emphasizes the challenges of isolating the rare γγ-M1M1 decay and suggests ways to enhance the experimental detection sensitivity. Our simulations suggest that it may be possible to access experimentally γγ-M1M1 decays from DIAS, shedding light on the neutrinoless ββ decay nuclear matrix elements.

Romeo, Beatriz↗
DOE OSTI2022

Elucidating the challenges in extracting ultra-slow flame speeds in a closed vessel - A CH 2 F 2 microgravity case study using optical and pressure-rise data

Refrigerants with a low global warming potential (GWP) possess mild flammability. Hence, a fundamental understanding of their combustion characteristics is required to assess their fire-hazardous potential. The laminar burning velocity is one fundamental property to describe these under safety evaluation aspects. However, measuring laminar burning velocities for slow-burning components, such as low-GWP alternatives, is challenging. This is due to influencing effects such as buoyant flame deformation, stretch, radiation, and confinement, especially at the flammability limits. Here, in the present study, we investigated near-limit flames of the representative refrigerant difluoromethane (CH 2 F 2 ) with nitrogen-enriched oxidizer mixtures to elucidate the potentials and limitations at ultra-slow flame speeds of two widely used flame speed measurement methods in the closed vessel configuration: the optical flame speed measurement in the early quasi-isobaric regime and the flame speed determination from the pressure-rise during the later phase of the isochoric combustion process. Experiments were performed under terrestrial- and microgravity, conducted in a specially designed setup suitable for drop tower applications, using both methods in parallel. Recommendations for the extraction of flame speed data are derived from the non-buoyant reference investigations under microgravity and transferred to ultra-slow flames at terrestrial gravity. Near-limit flames under microgravity were found to be significantly affected by stretch effects for the optical method so that flame speed extrapolations to unstretched conditions are prone to errors. Stretch also affects the pressure method’s lower evaluation limit, and errors can be estimated by combining flame speed data from the pressure method and Markstein lengths from optical experiments.

42 ENGINEERING↗
DOE OSTI2023

Plasma-assisted deflagration to detonation transition in a microchannel with fast-frame imaging and hybrid fs/ps coherent anti-Stokes Raman scattering measurements

Our study examines kinetic enhancement by nanosecond dielectric barrier discharge (ns-DBD) plasma on fuellean dimethyl ether (DME), oxygen (O 2 ), and argon (Ar) premixtures during deflagration to detonation transition (DDT) experiments in a microchannel. Non-equilibrium plasma produces active species and radicals and creates fast and slow heating of a mixture to promote ignition due to electronic and vibrational excitation. Experiments have been conducted to examine the influence of the plasma discharge on the premixture and on the resultant deflagration to detonation transition (DDT) onset time and distance through the use of high speed imaging and one-dimensional, two-beam, femtosecond/picosecond, coherent anti-Stokes Raman scattering (CARS). A highspeed camera is used to trace the time histories of flame front position and velocity and to identify the dynamics and onset of DDT. The results show that plasma discharge can nonlinearly affect the onset time and distance of DDT. It is shown that a small number of plasma discharge pulses prior to ignition result in reduced DDT onset time and distance by 60% and 40%, respectively, when compared to the results without pre-excitation by ns discharges. The results also show that an increase of plasma discharge pulses results in an extended DDT onset time and distance of 224% and 94%, respectively. Time history of the deflagration wave speed of DME and the analysis of ignition timescale under the choking condition of the deflagration front suggest low temperature ignition may play a role for DME near the isobaric choking condition of the burned gas and the DDT. Plasmainduced conversion of the reactive mixture was assessed via the O 2 to CO 2 ratio as measured through fs/ps CARS during the DBD discharges. The present experiments demonstrate the ability of non-equilibrium plasma to alter the chemistry of DME/O 2 /Ar premixtures in order to control DDT for applications in advanced propulsion engines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gas-phase ion-molecule interactions in a collision reaction cell with triple quadrupole-inductively coupled plasma mass spectrometry: Investigations with N 2 O as the reaction gas

Nitrous oxide (N 2 O) was used as a reaction gas to investigate the gas phase ion-molecule interactions using the Agilent 8900 QQQ-ICP-MS. A multi-element standard containing 45 elements with masses ranging from 9 to 208 u was measured in the presence and absence of N 2 O. The main product ion species observed were oxides and nitrides. Comparison of the N 2 O reaction results with similar measurements conducted with O 2 revealed that N 2 O was more effective at forming oxides in general: the elements Cd and Pb were shown to produce oxides with N 2 O where the reaction did not occur with O 2 . Nitrous oxide was also shown to produce a significant amount of nitride species in a few cases. The general reactivity was shown to be consistent with density functional theory (DFT)-predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels at least 1% of the unreacted ion. Our results show that reaction enthalpy is a reasonable predictor of reactivity with N2O on the timescales of the interactions in non-thermal ICP-MS/MS systems. Our work demonstrates the utility of two relatively new platforms (commercial elemental ICP-MS/MS and EMSL Arrows interface to the NWChem program suite), which allows for the study of a large number of elements within a short period. While DFT with the basis sets utilized here is not the most accurate computational method, it is also not computationally expensive and is shown to be suitable for predicting gas phase reactivity in the QQQ-ICP-MS for the majority of ions studied. Here, the ease and rapidity of data collection and DFT calculations has the potential to be very impactful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

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

How equation of state selection impacts accuracy near the critical point: Forced convection supercritical CO 2 flow over a cylinder

Increasing application of supercritical fluids leads to increased use of numerical simulations to guide designs. However, the impact of Equations of State (EoS) selection on computational cost and accuracy is not well understood. Laminar flow around a heated cylinder is used as a test case. A robust matrix of 99-simulations is used to analyze the flow of carbon dioxide just above the critical point at the isobaric condition of 8 MPa and temperatures between 305 and 390 K. Four EoS were studied: Ideal Gas Law, cubic EoS’s Peng-Robinson and Soave-Redlich-Kwong, and the Span-Wagner EoS. Results were identical for isothermal cases, but for heated or cooled cases, the high-fidelity Span-Wagner EoS simulations presented drastically different wake structures. Compared to the high-fidelity SW EoS, cEoS’s underestimated Nu values by up to 38% and cEoS’s underestimated CD values by 18–42%. In conclusion, the Span-Wagner EoS required only 17% more computation time then the cubic EoS simulations, making the Span-Wagner EoS recommended for engineering design.

42 ENGINEERING↗
DOE OSTI2020

Automated Coupling of Nanodroplet Sample Preparation with Liquid Chromatography–Mass Spectrometry for High-Throughput Single-Cell Proteomics

Single-cell proteomics can provide critical biological insight into the cellular heterogeneity that is masked by bulk-scale analysis. Here, we have developed a nanoPOTS (nanodroplet processing in one pot for trace samples) platform and demonstrated its broad applicability for single-cell proteomics. However, because of nanoliter-scale sample volumes, the nanoPOTS platform is not compatible with automated LC-MS systems, which significantly limits sample throughput and robustness. To address this challenge, we have developed a nanoPOTS autosampler allowing fully automated sample injection from nanowells to LC-MS systems. We also developed a sample drying, extraction, and loading workflow to enable reproducible and reliable sample injection. The sequential analysis of 20 samples containing 10 ng tryptic peptides demonstrated high reproducibility with correlation coefficients of >0.995 between any two samples. The nanoPOTS autosampler can provide analysis throughput of 9.6, 16, and 24 single cells per day using 120, 60, and 30 min LC gradients, respectively. As a demonstration for single-cell proteomics, the autosampler was first applied to profiling protein expression in single MCF10A cells using a label-free approach. At a throughput of 24 single cells per day, an average of 256 proteins was identified from each cell and the number was increased to 731 when the Match Between Runs algorithm of MaxQuant was used. Using a multiplexed isobaric labeling approach (TMT-11plex), ~77 single cells could be analyzed per day. We analyzed 152 cells from three acute myeloid leukemia cell lines, resulting in a total of 2558 identified proteins with 1465 proteins quantifiable (70% valid values) across the 152 cells. These data showed quantitative single-cell proteomics can cluster cells to distinct groups and reveal functionally distinct differences.

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

Simultaneous Isotopic Analysis of U, Pu, and Am in Spent Nuclear Fuel by Resonance Ionization Mass Spectrometry

Solid samples of spent nuclear fuel were analyzed for actinide isotopic composition by resonance ionization mass spectrometry. Isotopes of U, Pu, and Am were simultaneously quantified using a new method that removes and/or resolves the isobaric interferences at 238 U/ 238 Pu and 241 Pu/ 241 Am without sample preparation other than cutting and mounting small (~10 μm) samples. Trends in burnup and neutron capture product distributions were correlated with the sampling positions inside the reactor. Overall, the results show the skin effect, in which the core and near-edge regions of a fuel pellet exhibit strong differences in actinide concentrations and isotope distributions due to differences in the neutron energy spectra between the pellet rim and the core. While no elemental concentration measurements were made, the ability to measure the 238 Pu/ 239 Pu ratio in the presence of a 7400× excess of 238 U enabled an estimate of the enhancement in Pu concentration due to the skin effect at the rim of the pellet.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Evaluating Linear Ion Trap for MS3-Based Multiplexed Single-Cell Proteomics

There is a growing demand to develop high-throughput and high-sensitivity mass spectrometry methods for single-cell proteomics. The commonly used isobaric labeling-based multiplexed single-cell proteomics approach suffers from distorted protein quantification due to co-isolated interfering ions during MS/MS fragmentation, also known as ratio compression. We reasoned that the use of MS3-based quantification could mitigate ratio compression and provide better quantification. However, previous studies indicated reduced proteome coverages in the MS3 method, likely due to long duty cycle time and ion losses during multilevel ion selection and fragmentation. Here, in this paper, we described an improved MS acquisition method for MS3-based single-cell proteomics by employing a linear ion trap to measure reporter ions. We demonstrated that linear ion trap can increase the proteome coverages for single-cell-level peptides with even higher gain obtained via the MS3 method. The optimized real-time search MS3 method was further applied to study the immune activation of single macrophages. Among a total of 126 single cells studied, over 1200 and 1000 proteins were quantifiable when at least 50 and 75% nonmissing data were required, respectively. Our evaluation also revealed several limitations of the low-resolution ion trap detector for multiplexed single-cell proteomics and suggested experimental solutions to minimize their impacts on single-cell analysis.

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

A Dual-Gated Structures for Lossless Ion Manipulations-Ion Mobility Orbitrap Mass Spectrometry Platform for Combined Ultra-High-Resolution Molecular Analysis

High-resolution ion mobility spectrometry-mass spectrometry (HR-IMS-MS) instruments have enormously advanced the ability to characterize complex biological mixtures. Unfortunately, HR-IMS and HR-MS measurements are typically performed independently due to mismatches in analysis time scales. Here we overcome this limitation by using a dual-gated ion injection approach to couple an 11-meter path length structures for lossless ion manipulations (SLIM) module to a Q-Exactive Plus Orbitrap MS. The dual-gate setup was implemented by placing one ion gate before the SLIM module and a second ion gate after. The dual-gated ion injection approach allowed the new SLIM-Orbitrap platform to simultaneously perform an 11-meter SLIM separation, Orbitrap mass analysis using the highest selectable mass resolution setting (up to 140k), and high-energy collision induced dissociation (HCD) in ~25 minutes over an $m/z$ range of ~1500 amu. The SLIM-Orbitrap was initially characterized using a mixture of standard phosphazene cations and demonstrated an average SLIM CCS resolving power (Rp CCS ) of ~218 and SLIM peak capacity of ~156 while simultaneously obtaining high mass resolutions. SLIM-Orbitrap analysis with fragmentation was then performed on mixtures of standard peptides and two reverse peptides (SDGRG 1+ , GRGDS 1+ , Rp CCS = 305) to demonstrate the utility of combined HR-IMS-MS/MS measurements for peptide identification. Our new HR-IMS-MS/MS capability was further demonstrated by analyzing a complex lipid mixture and showcasing SLIM separations on isobaric lipids. In conclusion, this new SLIM-Orbitrap platform demonstrates a critical new capability for proteomics and lipidomics applications, and the high-resolution multimodal data obtainable with this system establishes the foundation for reference-free identification of unknown ion structures.

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

Structure–Transport Relationships of Water–Organic Solvent Co-transport in Carbon Molecular Sieve (CMS) Membranes

We explore the effects of the carbon molecular sieve (CMS) microstructure on the separation performance and transport mechanism of water–organic mixtures. Specifically, we utilize PIM-1 dense films and integrally skinned asymmetric hollow fiber membranes as polymer precursors for the CMS materials. The PIM-1 membranes were pyrolyzed under several different pyrolysis atmospheres (argon, carbon dioxide, and diluted hydrogen gas) and at multiple pyrolysis temperatures. Detailed gas physisorption measurements reveal that membranes pyrolyzed under 4% H 2 and CO 2 had broadened ultramicropore distributions (pore diameter <7 Å) compared to Ar pyrolysis, and pyrolysis under CO 2 increased ultramicropore volume and broadened micropore distributions at increased pyrolysis temperatures. Gravimetric water and p-xylene sorption and diffusion measurements reveal that the PIM-1-derived CMS materials are more hydrophilic than other CMS materials that have been previously studied, which leads to sorption-diffusion estimations showing water-selective permeation. Water permeation in the vapor phase, pervaporation, and liquid-phase hydraulic permeation reveal that the isobaric permeation modes (vapor permeation and pervaporation) are reasonably well predicted by the sorption-diffusion model, whereas the hydraulic permeation mode is significantly underpredicted (>250×). Conversely, the permeation of p-xylene is well predicted by the sorption-diffusion model in all cases. The collection of pore size analysis, vapor sorption and diffusion, and permeation in different modalities creates a picture of a combined transport mechanism in which water–under high transmembrane pressures–permeates via a Poiseuille-style mechanism, whereas p-xylene solutes in the mixture permeate via sorption-diffusion.

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

In Situ High-Temperature Raman Spectroscopy of UCl 3 : A Combined Experimental and Theoretical Study

Uranium trichloride (UCl 3 ) has received growing interest for its use in uranium-fueled molten salt reactors and in the pyrochemical processing of used fuel. In this paper, we report for the first time the experimentally determined Raman spectra of UCl 3 , at both ambient condition and in situ high temperatures up to 871 K. The frequencies of five of the Raman-active vibrational modes (v i ) of UCl 3 exhibit a negative temperature derivative ((∂ν i /∂T) P ) with increasing temperature. Here, this red-shift behavior is likely due to the elongation of U–Cl bonds. The average isobaric mode Grüneisen parameter (γ iP = 0.91 ± 0.02) of UCl 3 was determined through use of the coefficient of thermal expansion published in Vogel et al. (2021) and the (∂ νi /∂T) P values determined in this study. These results are in general agreement with those calculated here by density functional theory (DFT+U). Finally, a comparison of the ambient band positions of UCl 3 to those of isostructural lanthanide (La–Eu) and actinide chlorides (Am–Cf) has been made.

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

Molecular Simulations Probing the Adsorption and Diffusion of Ammonia, Nitrogen, Hydrogen, and Their Mixtures in Bulk MFI Zeolite and MFI Nanosheets at High Temperature and Pressure

Recent advances in the synthesis of MFI zeolite nanosheets have led to highly selective membranes that are promising candidates for small-scale ammonia separation from ammonia/nitrogen/hydrogen mixtures in distributed green ammonia production plants. Using force-field-based molecular simulations, we evaluate the performance of bulk all-silica MFI zeolite and a 3 nm thick all-silica MFI nanosheet for ammonia/nitrogen/hydrogen separations at moderate temperature and pressure conditions (T = 373 K, p = 5 bar) as well as conditions relevant for a membrane-based reactor–separator process (T = 523 or 623 K, p = 80 bar). Isobaric–isothermal Gibbs ensemble Monte Carlo (GEMC) simulations were carried out to understand the selective adsorption behavior for these mixtures. Molecular dynamics simulations in the canonical ensemble were performed to examine the transport behavior of ammonia in a mixture with nitrogen or hydrogen at loadings obtained from the GEMC simulations. Our results show that both bulk MFI and the nanosheet with explicit surface silanols are highly selective toward ammonia adsorption, but the adsorption selectivity decreases by factors of 4 and 10 for ammonia/nitrogen and ammonia/hydrogen mixtures as the temperature is increased from 373 to 623 K. Conversely, the diffusion selectivities toward ammonia are more favorable at process-relevant temperatures and for the MFI nanosheet. Here, at T = 523 K and p = 80 bar, our simulations predict overall separation factors of 3.8 ± 0.2 and 4.5 ± 0.3 for ammonia/nitrogen and ammonia/hydrogen mixtures, respectively, in the bulk MFI zeolite, and separation factors of 8.6 ± 0.4 and 12.5 ± 0.8, respectively, for the MFI nanosheet.

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

Python Library for Monte Carlo Simulations with Ab Initio and Machine-Learned Interatomic Potentials

There is a growing need in the simulation community for software that provides a transparent, reproducible, usable, and extensible (TRUE) Monte Carlo (MC) simulation framework employing energies from ab initio methods and machine-learning interatomic potentials (MLIPs). We introduce a Python library (ASE-MC) that adds Monte Carlo functionality to the Atomic Simulation Environment (ASE) package. Now, we can combine the powerful tools used to build systems and perform ab initio and MLIP in ASE with MC simulation algorithms to sample the configurational space with a concise Python script. After presenting the design philosophy, we demonstrate the flexibility of our approach using selected examples. These example simulations include liquid water described with a message-passing MLIP in the canonical and isothermal–isobaric ensembles, sampling the characteristic dihedral angle of biphenyl and comparing an MLIP to first-principles calculations, and a grand canonical Monte Carlo simulation of ammonia adsorption on Pt(111). These examples showcase the main features of the software, which include flexibility in the choice of ab initio or MLIP engine, ab initio or MLIP grand canonical MC with cavity bias insertions and deletions, the ability to add custom MC moves to the move set, and how users can condense complex MC workflows into a single Python script. Finally, this library serves as a framework for reproducible Monte Carlo simulations, facilitating easy reproduction of the work and application to new systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Phase Change-Driven Negative Activation Energies in Pd/Carbon-Based/Organic Getter Hydrogenation Reactions

In this study, the hydrogenation of 1,4-diphenylbutadiyne (DPB) blended with carbon-supported Pd (DPB–Pd/C) in the form of pellets was investigated by isothermal–isobaric experiments at 1333 Pa of H 2 and in the temperature range of 291–315 K. The extracted kinetics were then used in conjunction with a complementary constant rate of H 2 input experimentation to model the performance of a DPB-catalysis/support system as a function of temperature and H 2 partial pressure. First-principles density functional theory (DFT) calculations were also performed to shed light on the molecular level energetics of DPB and its intermediate states. A seemingly puzzling formation of alternate positive activation energy barrier (higher reaction rate with higher temperature) and negative activation energy barrier (higher reaction rate with lower temperature) zones during the hydrogenation process was discovered. However, this observed phenomenon can be logically explained in terms of the associated phase changes and H 2 transport in the material. This work provides a good illustration of a rarely encountered chemical process with a negative activation energy barrier.

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

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

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

Accurate Modeling of Bromide and Iodide Hydration with Data-Driven Many-Body Potentials

Ion–water interactions play a central role in determining the properties of aqueous systems in a wide range of environments. However, a quantitative understanding of how the hydration properties of ions evolve from small aqueous clusters to bulk solutions and interfaces remains elusive. Here, we introduce the second generation of data-driven many-body energy (MBnrg) potential energy functions (PEFs) representing bromide–water and iodide–water interactions. The MB-nrg PEFs use permutationally invariant polynomials to reproduce two-body and three-body energies calculated at the coupled cluster level of theory, and implicitly represent all higher-body energies using classical many-body polarization. Further, a systematic analysis of the hydration structure of small Br - (H 2 O) n and I - (H 2 O) n clusters demonstrates that the MBnrg PEFs predict interaction energies in quantitative agreement with “gold standard” coupled cluster reference values. Importantly, when used in molecular dynamics simulations carried out in the isothermal-isobaric ensemble for single bromide and iodide ions in liquid water, the MB-nrg PEFs predict extended X-ray absorption fine structure (EXAFS) spectra that accurately reproduce the experimental spectra, which thus allows for characterizing the hydration structure of the two ions with high level of confidence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗