Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 19 records

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↗

Benchmarking Guanidinium Organosulfonate Hydrogen-Bonded Frameworks for Structure Determination of Encapsulated Guests

Single crystal X-ray diffraction (SCXRD) is arguably the most definitive method for molecular structure determination, but it is often challenged by compounds that are liquids or oils at room temperature or do not form crystals adequate for analysis. Our laboratory previously reported a simple, cost-effective, single-step crystallization method based on guanidinium organosulfonate (GS) hydrogen bonded frameworks for structure determination of a wide range of encapsulated guest molecules, including assignment of the absolute configuration of chiral centers. Herein, we expand on those results and report a head-to-head comparison of the GS method with adamantoid “molecular chaperones”, which have been reported to be useful hosts for structure determination. Inclusion compounds limited to only two GS hosts are characterized by low R1 values and Flack parameters, infrequent disorder of the host and guest, and manageable disorder when it does exist. The structures of some target molecules that were not included or resolved using the adamantoid chaperones were successfully included and resolved by the GS hosts, and vice versa. Of the 32 guests attempted by the GS method, 31 inclusion compounds afforded successful guest structure solutions, a 97% success rate. The GS hosts and adamantoid chaperones are complementary with respect to guest inclusion, arguing that both should be employed in the arsenal of methods for structure determination. Furthermore, the low cost of organosulfonate host components promises an accessible route to molecular structure determination for a wide range of users.

36 MATERIALS SCIENCE↗

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↗

On the temperature dependence of H-U{sub iso} in the riding hydrogen model

The temperature dependence of hydrogen U{sub iso} and parent U{sub eq} in the riding hydrogen model is investigated by neutron diffraction, aspherical-atom refinements and QM/MM and MO/MO cluster calculations. Fixed values of 1.2 or 1.5 appear to be underestimated, especially at temperatures below 100 K. The temperature dependence of H-U{sub iso} in N-acetyl-l-4-hydroxyproline monohydrate is investigated. Imposing a constant temperature-independent multiplier of 1.2 or 1.5 for the riding hydrogen model is found to be inaccurate, and severely underestimates H-U{sub iso} below 100 K. Neutron diffraction data at temperatures of 9, 150, 200 and 250 K provide benchmark results for this study. X-ray diffraction data to high resolution, collected at temperatures of 9, 30, 50, 75, 100, 150, 200 and 250 K (synchrotron and home source), reproduce neutron results only when evaluated by aspherical-atom refinement models, since these take into account bonding and lone-pair electron density; both invariom and Hirshfeld-atom refinement models enable a more precise determination of the magnitude of H-atom displacements than independent-atom model refinements. Experimental efforts are complemented by computing displacement parameters following the TLS+ONIOM approach. A satisfactory agreement between all approaches is found.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Molecular ageing: Free radical initiated epimerization of thymopentin – A case study

The epimerization of amino acid residues increases with age in living organisms. In the present study, the structural consequences and thermodynamic functions of the epimerization of thymopentin (TP-5), the active site of the thymic hormone thymopoietin, were studied using molecular dynamics and density functional theory methods. The results show that free radical-initiated D-amino acid formation is energetically favoured (−130 kJmol{sup −1}) for each residue and induces significant changes to the peptide structure. In comparison to the wild-type (each residue in the L-configuration), the radius of gyration of the D-Asp{sup 3} epimer of the peptide decreased by 0.5 Å, and disrupted the intramolecular hydrogen bonding of the native peptide. Beyond establishing important structural, energetic and thermodynamic benchmarks and reference data for the structure of TP-5, these results disseminate the understanding of molecular ageing, the epimerization of amino acid residues.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches

Determining the influence of the solvent on electrochemical reaction energetics is a central challenge in our understanding of electrochemical interfaces. To date, it is unclear how well existing methods predict solvation energies at solid/liquid interfaces, since they cannot be assessed experimentally. Ab initio molecular dynamics (AIMD) simulations present a physically highly accurate, but also a very costly approach. In this work, we employ extensive AIMD simulations to benchmark solvation at charge-neutral metal/water interfaces against commonly applied continuum solvent models. We consider a variety of adsorbates including *CO, *CHO, *COH, *OCCHO, *OH, and *OOH on Cu, Au, and Pt facets solvated by water. The surfaces and adsorbates considered are relevant, among other reactions, to electrochemical CO 2 reduction and the oxygen redox reactions. We determine directional hydrogen bonds and steric water competition to be critical for a correct description of solvation at the metal/water interfaces. As a consequence, we find that the most frequently applied continuum solvation methods, which do not yet capture these properties, do not presently provide more accurate energetics over simulations in vacuum. We find most of the computed benchmark solvation energies to linearly scale with hydrogen bonding or competitive water adsorption, which strongly differ across surfaces. Furthermore, we determine solvation energies of adsorbates to be non-transferable between metal surfaces, in contrast to standard practice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Appointing silver and bronze standards for noncovalent interactions: A comparison of spin-component-scaled (SCS), explicitly correlated (F12), and specialized wavefunction approaches

A systematic examination of noncovalent interactions as modeled by wavefunction theory is presented in comparison to gold-standard quality benchmarks available for 345 interaction energies of 49 bimolecular complexes. Quantum chemical techniques examined include spin-component-scaling (SCS) variations on second-order perturbation theory (MP2) [SCS, SCS(N), SCS(MI)] and coupled cluster singles and doubles (CCSD) [SCS, SCS(MI)]; also, method combinations designed to improve dispersion contacts [DW-MP2, MP2C, MP2.5, DW-CCSD(T)-F12]; where available, explicitly correlated (F12) counterparts are also considered. Dunning basis sets augmented by diffuse functions are employed for all accessible ζ-levels; truncations of the diffuse space are also considered. After examination of both accuracy and performance for 394 model chemistries, SCS(MI)-MP2/cc-pVQZ can be recommended for general use, having good accuracy at low cost and no ill-effects such as imbalance between hydrogen-bonding and dispersion-dominated systems or non-parallelity across dissociation curves. Moreover, when benchmarking accuracy is desirable but gold-standard computations are unaffordable, this work recommends silver-standard [DW-CCSD(T**)-F12/aug-cc-pVDZ] and bronze-standard [MP2C-F12/aug-cc-pVDZ] model chemistries, which support accuracies of 0.05 and 0.16 kcal/mol and efficiencies of 97.3 and 5.5 h for adenine·thymine, respectively. Choice comparisons of wavefunction results with the best symmetry-adapted perturbation theory [T. M. Parker, L. A. Burns, R. M. Parrish, A. G. Ryno, and C. D. Sherrill, J. Chem. Phys. 140, 094106 (2014)] and density functional theory [L. A. Burns, Á. Vázquez-Mayagoitia, B. G. Sumpter, and C. D. Sherrill, J. Chem. Phys. 134, 084107 (2011)] methods previously studied for these databases are provided for readers' guidance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On single-crystal neutron-diffraction in DACs: quantitative structure refinement of light elements on SNAP and TOPAZ

Quantitative single crystal neutron-diffraction in diamond anvil cells has so far been limited by the neutron flux available at the various neutron sources. As a result, highly precise measurements of the exact position of light elements have not been possible preventing, for example, structural studies of hydrogen and hydrogen bonds under pressure. Here we report experiments carried out on SNAP at the Spallation Neutron Source (ORNL, TN, USA) to explore the possibility and current limits of such studies. Furthermore, we benchmarked the obtained data quality with reference experiments carried out on TOPAZ, a dedicated single-crystal instrument.We show that measuring single-crystal diffraction intensities on SNAP is possible to such a precision that we are able to resolve the hydrogen bonds in potassium dideuterium phosphate (DKDP) as well as in ice VI.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Benchmarking the Performance of the ReaxFF Reactive Force Field on Hydrogen Combustion Systems

A thorough understanding of the kinetics and dynamics of combusting mixtures is of considerable interest, especially in regimes beyond the reach of current experimental validation. The ReaxFF reactive force field method has provided a way to simulate large-scale systems of hydrogen combustion via a parametrized potential that can simulate bond breaking. This modeling approach has been applied to hydrogen combustion, as well as myriad other reactive chemical systems. In this work, we benchmark the performance of several common parametrizations of this potential against higher-level quantum mechanical (QM) approaches. We demonstrate instances where these parametrizations of the ReaxFF potential fail both quantitatively and qualitatively to describe reactive events relevant for hydrogen combustion systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Order and correlation contributions to the entropy of hydrophobic solvation

The entropy of hydrophobic solvation has been explained as the result of ordered solvation structures, of hydrogen bonds, of the small size of the water molecule, of dispersion forces, and of solvent density fluctuations. We report a new approach to the calculation of the entropy of hydrophobic solvation, along with tests of and comparisons to several other methods. The methods are assessed in the light of the available thermodynamic and spectroscopic information on the effects of temperature on hydrophobic solvation. Five model hydrophobes in SPC/E water give benchmark solvation entropies via Widom’s test-particle insertion method, and other methods and models are tested against these particle-insertion results. Entropies associated with distributions of tetrahedral order, of electric field, and of solvent dipole orientations are examined. We find these contributions are small compared to the benchmark particle-insertion entropy. Competitive with or better than other theories in accuracy, but with no free parameters, is the new estimate of the entropy contributed by correlations between dipole moments. Dipole correlations account for most of the hydrophobic solvation entropy for all models studied and capture the distinctive temperature dependence seen in thermodynamic and spectroscopic experiments. Entropies based on pair and many-body correlations in number density approach the correct magnitudes but fail to describe temperature and size dependences, respectively. Hydrogen-bond definitions and free energies that best reproduce entropies from simulations are reported, but it is difficult to choose one hydrogen bond model that fits a variety of experiments. The use of information theory, scaled-particle theory, and related methods is discussed briefly. Our results provide a test of the Frank-Evans hypothesis that the negative solvation entropy is due to structured water near the solute, complement the spectroscopic detection of that solvation structure by identifying the structural feature responsible for the entropy change, and point to a possible explanation for the observed dependence on length scale. Our key results are that the hydrophobic effect, i.e. the signature, temperature-dependent, solvation entropy of nonpolar molecules in water, is largely due to a dispersion force arising from correlations between rotating permanent dipole moments, that the strength of this force depends on the Kirkwood g-factor, and that the strength of this force may be obtained exactly without simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum cluster equilibrium model of N-methylformamide–water binary mixtures

The established quantum cluster equilibrium (QCE) approach is refined and applied to N-methylformamide (NMF) and its aqueous solution. The QCE method is split into two iterative cycles: one which converges to the liquid phase solution of the QCE equations and another which yields the gas phase. By comparing Gibbs energies, the thermodynamically stable phase at a given temperature and pressure is then chosen. The new methodology avoids metastable solutions and allows a different treatment of the mean-field interactions within the gas and liquid phases. These changes are of crucial importance for the treatment of binary mixtures. For the first time in a QCE study, the cis-trans-isomerism of a species (NMF) is explicitly considered. Cluster geometries and frequencies are calculated using density functional theory (DFT) and complementary coupled cluster single point energies are used to benchmark the DFT results. Independent of the selected quantum-chemical method, a large set of clusters is required for an accurate thermodynamic description of the binary mixture. The liquid phase of neat NMF is found to be dominated by the cyclic trans-NMF pentamer, which can be interpreted as a linear trimer that is stabilized by explicit solvation of two further NMF molecules. This cluster reflects the known hydrogen bond network preferences of neat NMF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab Initio Molecular Dynamics Spectra for Characterization of Hydrated Al 2 O 3 Supported MoO x

The structure of supported MoO x /Al 2 O 3 catalysts is investigated using ab initio molecular dynamics (AIMD) and density-functional theory (DFT). Phase diagrams were computed to understand the hydroxylation coverage as a function of the temperature, H 2 O partial pressure, and MoO x loading. We relate the shifts in experimental Raman vibrational frequencies under hydrated conditions to hydrogen bonding interactions and the MoO x anchoring location. We showcase that the use of AIMD as a benchmark for DFT can provide insight into how, under certain hydrated conditions, DFT excels as an inexpensive method for computing vibrational frequencies, while, under dehydrated conditions, it is susceptible to the largest errors. Additionally, to facilitate the analysis of the hydroxyl region of experimental infrared spectra, we compute the power spectra of individual hydroxyls and see a strong relationship between OH stretching frequencies and the surrounding coordination environment, which are also impacted by anchored MoO x . Here, we compare computational and experimental IR and Raman spectra of catalysts synthesized herein under the same conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

An accurate and linear-scaling method for calculating charge-transfer excitation energies and diabatic couplings

Quantum-mechanical methods that are both computationally fast and accurate are not yet available for electronic excitations having charge transfer character. In this work, we present a significant step forward towards this goal for those charge transfer excitations that take place between non-covalently bound molecules. In particular, we present a method that scales linearly with the number of non-covalently bound molecules in the system and is based on a two-pronged approach: The molecular electronic structure of broken-symmetry charge-localized states is obtained with the frozen density embedding formulation of subsystem density-functional theory; subsequently, in a post-SCF calculation, the full-electron Hamiltonian and overlap matrix elements among the charge-localized states are evaluated with an algorithm which takes full advantage of the subsystem DFT density partitioning technique. The method is benchmarked against coupled-cluster calculations and achieves chemical accuracy for the systems considered for intermolecular separations ranging from hydrogen-bond distances to tens of Angstroms. Numerical examples are provided for molecular clusters comprised of up to 56 non-covalently bound molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Crosslinked Ionic Organic Framework for Efficient Iodine and Iodide Remediation in Water

Iodine is widely used as an antimicrobial reagent for water disinfection in the wilderness and outer space, but residual iodine and iodide need to be removed for health reasons. Currently, it is challenging to remove low concentrations of iodine and iodide in water (~5 ppm). Furthermore, the remediation of iodine and iodide across a broad temperature range (up to 90 °C) has not previously been investigated. In this work, we report a nitrate dimer-directed synthesis of a single-crystalline ionic hydrogen-bonded crosslinked organic framework (HCOF-7). HCOF-7 removes iodine and iodide species in water efficiently through halogen bonding and anion exchange, reducing the total iodine concentration to 0.22 ppm at room temperature. Packed HCOF-7 columns were employed for iodine/iodide breakthrough experiments between 23 and 90 °C, and large breakthrough volumes were recorded (≥18.3 L/g). Finally, the high iodine/iodide removal benchmarks recorded under practical conditions make HCOF-7 a promising adsorbent for water treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗