Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrogen bond benchmarks”

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.

14 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2025

Spectroscopic Identification of Carbamate Formation and Synergistic Binding in Amide–CO$^{–}_{2}$ Complexes

Carbamate formation is an elementary step that governs nitrogencentered nucleophilic CO 2 capture across diverse environments, yet a direct, structure-specific understanding has been lacking. Here, we report the first gas-phase characterization of closed-shell carbamate formation by deprotonated amides, using cryogenic ion trap vibrational spectroscopy combined with quantum chemical calculations. Reactions of deprotonated benzamide and isophthalamide anions with CO 2 form carbamate species, with the amide nitrogen serving as the nucleophilic site. Diagnostic, strongly red-shifted antisymmetric CO 2 stretching vibrational bands, supported by a bonding analysis, establish chemisorption with substantial charge transfer. In the multiamide system, an intramolecular N−H···O hydrogen bond provides synergistic stabilization, correlating with larger red shifts and more exothermic binding. These structure-assigned benchmarks provide molecular-level insights into the binding motif, charge redistribution, and hydrogen bond-mediated stabilization of amide carbamates, which aid the characterization of carbamate formation in condensed phase.

Amides↗
DOE OSTI2026

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗
DOE OSTI2026

Size-dependent femtosecond proton transfer in protonated methanol clusters

Methanol clusters (CH 3 OH) n are prototypical hydrogen-bonded systems that bridge isolated molecules and bulk liquid methanol. Strong-field ionization of these clusters can trigger ultrafast charge migration, proton transfer, and isomerization; however, the onset timescales of these pathways remain largely unexplored despite extensive steady-state mass-spectrometric evidence for protonation and chemical rearrangement. Here we use femtosecond time-resolved strong-field ionization and disruptive probing to directly track the early-time dynamics of methanol clusters. Transient mass spectra yield size-dependent timescales for the formation of protonated clusters H + (CH 3 OH) n (n = 1–3) and an associated CHO˙ + radical cation. We find that both proton-transfer and subsequent stabilization become markedly faster with increasing cluster size. These measurements establish methanol clusters as a complementary benchmark to water for understanding ionization-initiated proton motion and chemical reorganization in hydrogen-bonded liquids.

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

Benchmarking of ENDF/B-VIII.1 Thermal Scattering Library for Hydrogen

In correctly characterizing the energy and momentum transfer at thermal and cold energies between neutrons and its interacting medium, thermal scattering libraries, which details energy states due to the intra- and inter-molecular bond effects for the medium materials, are applied in place of free-gas cross section libraries in a particle transport simulation code. They are essential to the neutron performance of a neutron facility like SNS, where thermalized neutrons from 20 K liquid hydrogen and ambient (~300 K) water are transported to the beamlines for neutron scattering experiments in studying materials. Recently, a new version of thermal scattering library for parahydrogen and orthohydrogen at 14-20 K was developed and to be released in ENDF/B-VIII.1. It is, therefore, important to benchmark its impacts on the prediction of moderator performance due to the updates in the thermal scattering library. In this study, the recent ENDF/B-VIII.1 thermal scattering library was compared to the current ENDF/B-VII.1 one in the neutron performance calculations for the decoupled and coupled hydrogen moderators at SNS under theorized and real working conditions. In addition, the predictions using both thermal scattering libraries were benchmarked to the measurements of moderator performance. The consistency between the libraries was observed mostly for parahydrogen and the difference in orthohydrogen at cold neutron energies was noted.

Lu, Wei [Oak Ridge National Laboratory (ORNL), Oak↗
DOE OSTI2026

Benchmarking quantum trial wavefunctions for phaseless auxiliary-field quantum Monte Carlo

The phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) method is a stochastic imaginary-time projection technique for computing ground-state properties of strongly correlated quantum systems, with accuracy that depends critically on the choice of trial wavefunction. Here, we investigate ph-AFQMC with trial states prepared using parameterized quantum circuits. In this work, we present a comprehensive benchmarking study of quantum trial wavefunctions spanning unitary coupled-cluster, Hamiltonian-informed, Jastrow-inspired, and adaptively constructed ansatze. The benchmarking evaluates accuracy, expressibility, and scalability of these ansatze within the QC-AFQMC framework. We test these ansatze on linear hydrogen chains under bond stretching and find that several ansatz families produce chemically accurate ph-AFQMC energies across the dissociation curve. We have performed simulations using the CUDA-Q quantum development platform on the GPU partition of the Perlmutter supercomputer. When comparing ansatze at similar numbers of variational parameters, we find that different ansatz families yield comparable ph-AFQMC results despite exhibiting substantially different variational energies, optimization costs, and circuit depths. Our results indicate that the variational energy of an ansatz is not always a reliable indicator of its quality for ph-AFQMC and reveal instances of over-parameterization. In the strongly correlated regime, trial wavefunctions obtained from adaptive ansatze, exemplified here by ADAPT-VQE with the UCCSD operator pool, can outperform their fixed-ansatz counterparts (UCCSD) in terms of projected energies while using substantially more compact circuits, providing a flexible route to optimize quantum resources within the ph-AFQMC framework.

Rofougaran, Rod [LBNL, Berkeley; Columbia U.; PNL,↗
DOE OSTI2025

Multi-Scale Modeling for Plasma-Enhanced Ammonia Decomposition over Carbides and Nitrides

Ammonia is a carbon-free hydrogen carrier, but its decomposition typically requires high temperatures over costly Ru-based catalysts due to the large barrier for N≡N bond formation. We develop a multiscale framework combining density functional theory, zero-dimensional plasma kinetics, and microkinetic modeling to elucidate how non-thermal plasma (NTP) enables low-temperature NH 3 decomposition over Co-based carbides and nitrides, benchmarked against Ru and Co. Under thermal conditions, all catalysts are limited by N≡N bond formation, with Co 3 C(001) most active owing to its negatively charged surface, strong N* binding, and low activation barriers of N≡N bond formation. Plasma-induced vibrational excitation of NH 3 and its reactive radicals promotes a radical-driven •NH 2 –N* coupling pathway that dominates on Co 3 C(001) and Co 3 N(001), shifting the rate-limiting step to NH 3 (v1) dissociation, increasing turnover frequencies by up to 6 orders of magnitude, and reducing the temperature needed to reach a turnover frequency of 5 s –1 from >680 °C (Ru and Co under thermal condition) to 267 °C (Co 3 C) and 415 °C (Co 3 N). These results identify Co-based carbides and nitrides as promising plasma-active catalysts for energy-efficient hydrogen production from ammonia.

ammonia decomposition↗
DOE OSTI2025

Quantitative and mechanistic insights into proton dynamics for fast energy storage

Proton conduction in hydrogen-bond-rich protic electrolytes enables fast mass and charge transport, crucial for electrochemical energy storage and power conversion. Such transport can give proton-based batteries exceptional rate capability and low-temperature operation beyond other working ions. Here we show that in phosphoric acid (H 3 PO 4 ) electrolytes, vehicular and structural proton transport coexist, and their contributions to conductivity can be quantitatively distinguished. We link structural diffusion directly to hydrogen-bond strength, enabling the precise tuning of proton migration. Guided by this, we reveal a double conductivity peak from regulated structural diffusion. The optimal electrolyte (5.8-M H 3 PO 4 ) achieves high overall (232.9 mS cm −1 ) and structural (164.9 mS cm−1) conductivity. A MoO 3 ‖CuFe-TBA battery with this electrolyte outperforms a deep-eutectic benchmark (8.3-M H 3 PO 4 ), delivering >17,474 W kg −1 at room temperature and retaining 15.1 Wh kg −1 at −75 °C. These findings provide a framework for designing advanced protic electrolytes across electrochemical systems.

Li, Ziyue [Fudan University, Shanghai (China, Peop↗
DOE OSTI2026

Exploring Local Vibrational Structure in Protein-Bound Chlorophyll a: Isotope-Enrichment Experiments and Electrostatic Analysis

Local chlorophyll (Chl) vibrations play an essential role in biological photosynthesis by facilitating rapid energy transfer between pigments. In research studies, they also provide a useful spectroscopic probe of the local protein environment that surrounds each pigment. However, measuring the complete vibrational absorption spectrum of a protein-bound Chl molecule is much more difficult than, for example, Chl in neat solvent due to overlap with protein vibrations that typically drown out Chl vibrational features. Resonance Raman and fluorescence spectroscopies provide a way around this problem for Franck–Condon active vibrations, but these often rely on cryogenic measurement conditions and fail to capture vibrational signatures from, e.g., ester-group vibrations that lack coupling to a convenient electronic transition. In the present contribution, we use 13 C-enrichment of the protein backbone to shift protein background signals to lower frequency, providing a largely clean spectral window in which to study local Chl a C═O stretch modes. A room-temperature absorption spectrum for Chl a in the water-soluble chlorophyll protein (WSCP) of Lepidium virginicum is thus extracted as a difference between protein-plus-pigment and protein-only vibrational spectra. Excellent agreement in the molecular fingerprint region with the vibrational spectrum of Chl a in organic solvents confirms that the resulting spectrum represents the response of the protein-bound Chl a molecule. Furthermore, the ester group resonance is observed to shift in response to the S53P mutation that eliminates a 17 3 ester group hydrogen bond from the protein environment. Finally, we analyze these experimental results using MD-based electrostatic analysis, finding that electric-field mapping at the C atom of the ester group provides a satisfactory explanation for the observed frequency shifts between organic solvent and protein environment. MD analysis further suggests that a red-shifted ester peak observed experimentally for the S53P mutant results from solvation of the 17 3 Chl ester group due increase water penetration into the Chl-binding pocket relative to wild-type WSCP. Here, we anticipate that these results will prove useful both for benchmarking future simulation work and as a reference for interpreting Chl vibrational spectra as a probe of pigment–protein interactions.

14 SOLAR ENERGY↗
DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

Cryogenic Vibrational Spectroscopy of the Deprotonated Dimer of Phosphoric Acid

Phosphate-containing molecules are ubiquitous in nature, where they play crucial roles in biochemical processes. Further, they are of technical importance, for example, in certain batteries and in fuel cells, where a unique property of phosphoric acid is exploited its exceptionally high proton conductivity. Proton transport in phosphoric acid is known to involve proton shuttling; however, the elementary steps involved are not clear. To elucidate the hydrogen bonding preferences of phosphoric acid, we investigate the dihydrogen phosphate anion as well as the deprotonated dimer of phosphoric acid (H 3 PO 4 ·H 2 PO 4 − ) in the gas phase using infrared action spectroscopy in helium nanodroplets and infrared D 2 -tagging photodissociation spectroscopy, and the experimental spectra are compared to theoretical ones. Theory finds for H 3 PO 4 ·H 2 PO 4 − two different structures that are predicted to be nearly isoenergetic. The comparison to the experimental spectra, however, allows for a clear assignment and structure identification. The resulting structure has an interesting binding motif, which might be of relevance to interactions of phosphoric acid in the condensed phase and which can serve as a benchmark for quantum chemical calculations.

Chemical structure↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2024

Promoted Ru/PrOx Catalysts for Mild Ammonia Synthesis

Ammonia synthesis is one of the most important chemical reactions. Due to thermodynamic restrictions and the reaction requirements of the current commercial iron catalysts, it is also one of the worst reactions for carbon dioxide emissions and energy usage. Ruthenium-based catalysts can substantially improve the environmental impact as they operate at lower pressures and temperatures. In this work, we provide a screening of more than 40 metals as possible promoter options based on a Ru/Pr2O3 catalyst. Cesium was the best alkali promoter and was held constant for the series of double-promoted catalysts. Ten formulations outperformed the Ru-Cs/PrOx benchmark, with barium being the best second promoter studied and the most cost-effective option. Designs of experiments were utilized to optimize both the pretreatment conditions and the promoter weight loadings of the doubly promoted catalyst. As a result, optimization led to a more than five-fold increase in activity compared to the unpromoted catalyst, therefore creating the possibility for low-ruthenium ammonia synthesis catalysts to be used at scale. Further, we have explored the roles of promoters using kinetic analysis, X-ray Photoelectron Spectroscopy (XPS), and in situ infrared spectroscopy. Here, we have shown that the role of barium is to act as a hydrogen scavenger and donor, which may permit new active sites for the catalyst, and have demonstrated that the associative reaction mechanism is likely used for the unpromoted Ru/PrOx catalyst with hydrogenation of the triple bond of the dinitrogen occurring before any dinitrogen bond breakage.

Chemistry↗
DOE OSTI2024

Techno-Economic Analysis of Green Hydrogen Energy Storage in a Cryogenic Flux Capacitor

Abstract The Cryogenic Flux Capacitor (CFC) is a cold, dense energy storage core that is being studied in the cryo-compressed, about 300 bar and 80K, region of gaseous hydrogen (GH2) storage and liquid hydrogen (LH2) region near the normal boiling point. Hydrogen storage is improved by physically bonding the molecules within the nanoscale pores of the aerogel composite blanket material. The process of bonding or debonding is governed by principles of physical adsorption (physisorption) and thermodynamics. The large surface area afforded by the nanoporous aerogel (∼1,000 m2/g) allows its storage performance to easily exceed capacities of high-pressure GH2 storage for an equivalent volume. With the integrated aerogel, subscale tests have shown that storage is increased by about 36% over a simple tank filled with GH2 at the same operating temperature and pressure. For LH2 conditions, the CFC is shown to operate at improved densities, but testing is ongoing. For the techno-economic analysis (TEA), the source of hydrogen is compared between onsite steam methane reforming (SMR) and onsite solar photovoltaic (PV) panels providing power to electrolyzers to produce green GH2. The TEA compares pure hydrogen produced at a small scale for a 25 MW power system and at a large scale in a 500 MW power system. The system allowed for hydrogen imports and exports at a set price with a tank sized for 10 hours of power production. The two power producing technologies are a combined cycle gas turbine (CCGT) and hydrogen fuel cells. The SMR system uses natural gas as an input and includes a carbon capture and storage (CCS) system. The levelized cost of electricity (LCOE), levelized cost of hydrogen (LCOH), and levelized cost of storage (LCOS) are developed based on the capital cost and operating cost of the systems. The results are shown for current costs using a 2021 benchmark and DOE projections for cost improvements by 2030. The TEA showed that onsite hydrogen generation from SMR has an LCOH of about 1.4 to 2 USD per kg over the life of the plant and the PV hydrogen production LCOH is about 5.2 to 5.5 USD per kg. The LCOS of conventional GH2 systems is estimated to be $210/MWh and cost of storage for LH2 systems is $205/MWh for fuel cell systems and $249/MWh for CCGT systems. CFC improved the LCOS of all these systems to $198/MWh, $191/MWh and $233/MWh respectively. The LCOE also improved with conventional systems between $171/MWh and $228/MWh improved by CFC to between $167/MWh and $212/MWh. Using projections for improvement in costs following DOE’s goals by 2030, green hydrogen improved to as low as $78/MWh LCOS and LCOE for conventional cases. CFC improved over conventional storage with the lowest LCOS being $62/MWh and the lowest LCOE being $73/MWh. These results correspond to an LCOH of $2/kg. Finally, the TEA shows how LCOE is improved for hydrogen conditioning and storage over conventional systems and caverns in the 10 to 50 hour range.

08 HYDROGEN↗