SEARCH · Engineering Papers

Results for “Enthalpy”

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 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Thermodynamics of cesium lead halide (CsPbX 3 , x= I, Br, Cl) perovskites

Inorganic cesium lead halide (CsPbX 3 , X = I, Br, Cl) perovskites have sparked broad interest for photovoltaic applications, which, however, may be plagued by chemical and thermodynamic instabilities. This work reports calorimetric studies on the formation enthalpies of CsPbX 3 perovskites, confirming decreasing thermodynamic stability in the order CsPbCl 3 , CsPbBr 3 and CsPbI 3 . Cesium lead iodide (CsPbI 3 ), one of the most studied inorganic lead halide perovskites, is shown to be only barely stable from a thermodynamic perspective. This order of increasing instability is consistent with the variation of Goldschmidt tolerance factor. The perovskites are high temperature entropy-stabilized phases with respect to low temperature polymorphs, with the thermodynamics of phase transitions reflecting the above thermodynamic trends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Effect of pressure on TATB and LX-17 thermal decomposition

Thermal decomposition of TATB (1,3,5-triamino-2,4,6-trinitrobenzene) and its formulation LX-17 is studied at pressures from 0.1 to 7 MPa for both isothermal heating at 340 K and ramped heating at 1 to 6 K min -1 . Conditions that eliminate self-heating are thoroughly explored to avoid experimental artifacts. The increase in pressure accelerates the rate of decomposition by only about 10 %, but it substantially increases the enthalpy of the reaction, presumably because of longer volatile product residence times in the heated zone. The narrowness of the decomposition profile and the acceleratory phase during isothermal pyrolysis are consistent with a multistep autocatalytic mechanism, and the lack of a significant pressure effect suggests that the autocatalytic species have low volatility. The apparent activation energy is about 210 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Using a Chromatographic Pseudophase Model To Elucidate the Mechanism of Olefin Separation by Silver(I) Ions in Ionic Liquids

Silver(I) ions undergo selective olefin complexation and have been utilized in various olefin/paraffin separation techniques such as argentation chromatography and facilitated transport membranes. Ionic liquids (ILs) are solvents known for their low vapor pressure, high thermal stability, low melting points, and ability to promote a favorable solvation environment for silver(I) ion–olefin interactions. To develop highly selective separation systems, a fundamental understanding of analyte partitioning to the stationary phase and the thermodynamic driving forces behind solvation is required. In this study, a chromatographic model treating silver(I) ions as a pseudophase is constructed and employed for the first time to investigate the olefin separation mechanism in silver(I) salt/IL mixtures. Stationary phases containing varying amounts of noncoordinated silver(I) salt ([Ag + ][NTf 2 – ]) dissolved in the 1-decyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C 10 MIM + ][NTf 2 – ]) IL are utilized to determine the partition coefficients of various analytes including alkanes, alkenes, alkynes, aromatics, aldehyde, esters, and ketones. As ligand coordination to silver(I) ions is known to lower its olefin complexation capability, this study also examines two different types of coordinated silver(I) ion pseudophases, namely, monocoordinated silver(I) salt ([Ag + (1-decyl-2-methylimidazole, DMIM)][NTf 2 – ]) and dicoordinated silver(I) salt ([Ag + (1-methylimidazole, MIM)(DMIM)][NTf 2 – ]). The extent of olefin partitioning to the coordinated silver(I) ion pseudophases over the carrier gas and IL decreased by up to two orders of magnitude. Values for enthalpy, entropy, and free energy of solvation were determined for the three silver(I) ion-containing systems. Olefin retention was observed to be enthalpically dominated, while ligand coordination to the silver(I) ion pseudophase resulted in variations for both enthalpic and entropic contributions to the free energy of solvation. Furthermore, the developed model can be used to study chemical changes that occur in silver(I) ions over time as well as identify optimal silver(I) salt/IL mixtures that yield high olefin selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Voltammetric Measurement of Rates and Energetics for Surface Methoxylation of Si(100) in Methanol with Dissolved Electron Acceptors Using Si Ultramicroelectrodes

The steady-state voltammetric responses of n-type Si(100) semiconductor ultramicroelectrodes (SUMEs) immersed in air- and water-free methanolic electrolytes have been measured. The response characteristics of these SUMEs in the absence of illumination were modeled and understood through a framework that describes the distribution of the applied potential across the semiconductor/electrolyte contact using four discrete regions: the semiconductor space charge, surface, Helmholtz, and diffuse layers. The latter region was described by the full Gouy–Chapman model. This framework afforded insight on how relevant parameters such as the semiconductor band edge potentials, the reorganization energies for charge transfer, the standard potential of redox species in solution, the density and energy of surface state populations, and the presence of an insulating (tunneling) layer individually and collectively dictate the observable current–potential responses. With this information, the methoxylation of Si surfaces was evaluated by analysis of the change in voltammetric responses during the course of prolonged immersion in methanol. Here, the electrochemical data were consistent with a surface methoxylation mechanism that depended on the standard potential of redox species dissolved in solution. Estimates of the enthalpies of adsorption as well as the potential-dependent rate constant for surface methoxylation were obtained. Collectively, these measurements supported the contention that the rates of Si surface reactions can be systematically tuned by exposure to dissolved outer-sphere electron acceptors. Moreover, the data represent the quantitative utility of voltammetry with SUMEs for the measurement of semiconductor/liquid contacts.

14 SOLAR ENERGY↗

DOE OSTI2024

Assessing Gas-Phase Ion Reactivity of 50 Elements with NO and the Direct Application for 239 Pu in Complex Matrices Using ICP-MS/MS

Understanding the reactivity of metal cations with various reaction gases in ICP-MS/MS is important to determine the best gas to use for a given analyte/interference pair. In this study, nitric oxide (NO) was investigated as the reaction gas following previous experimental designs. The reactions with 50 elements were investigated to examine periodic trends in reactivity, validate theoretical modeling of reaction enthalpies as a method to screen reactant gases, and provide a baseline for potential in-line gas separation methods. ICP-MS/MS studies involving actinides are typically limited to Th, U and Pu, with analyses of Np and Am rarely reported in the literature. To date, only two previous methods have investigated the use of NO in ICP-MS/MS analyses. To showcase the utility of NO, a method was developed to measure 239 Pu in the presence of environmental matrix constituent and other actinides, like what could be expected from post-detonation debris, with no chemical separation prior to analysis. 239 Pu + was reacted to form 239 Pu 16 O + , eliminating interferences derived from the sample matrix by measuring 239 Pu + intensity at m/z = 255 ( 239 Pu 16 O + ). To validate NO for 238 U 1 H + interference removal in environmental matrices, standard reference materials were diluted to 1 mg per gram of solution and spiked to 0.05 pg/g of 239 Pu and 1 µg/g 238 U (Pu/U = 5 x 10 -8 ). Measured 239 Pu concentrations were within 6% of the spiked value. In conclusion, these results demonstrate that reliable 239 Pu measurements can be made at levels relevant to nuclear forensics without the need for extensive chemical matrix separation prior to analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Volatility Basis Set Distributions and Viscosity of Organic Aerosol Mixtures: Insights from Chemical Characterization Using Temperature-Programmed Desorption–Direct Analysis in Real-Time High-Resolution Mass Spectrometry

Quantitative assessment of gas-particle partitioning of individual components within complex atmospheric organic aerosols (OA) mixtures is critical for predicting and comprehending the formation and evolution of OA particles in the atmosphere. This investigation leverages previously documented data obtained through a temperature programmed desorption - direct analysis in real time – high resolution mass spectrometry (TPD-DART-HRMS) platform. This methodology facilitates the bottom-up construction of volatility basis set (VBS) distributions for constituents found in three biogenic secondary organic aerosol (SOA) mixtures produced through the ozonolysis of a-pinene, limonene, and ocimene. The apparent enthalpies (ΔH*, kJ mol -1 ) and saturated vapor mass concentrations (C T *, µg∙m -3 ) of individual SOA components, determined as a function of temperature (T, K), facilitated an assessment of changes in VBS distributions and gas-particle partitioning with respect to T and atmospheric total organic mass loadings (tOM, µg∙m -3 ). Further, the VBS distributions reveal distinct differences in volatilities among monomers, dimers, and trimers, enabling their categorization into separate volatility bins. At the ambient temperature of T = 298 K, only monomers efficiently partition between gas and particle phases across a broad range of atmospherically relevant total organic mass loadings (tOM) values of 1–100 µg∙m -3 . Partitioning of dimers and trimers becomes notable only at T > 360 K and T > 420 K, respectively. The viscosity of SOA mixtures is assessed using a bottom-up calculation approach, incorporating the input of elemental formulas, ΔH*, C T *, and particle-phase mass fractions of the SOA components. Through this approach, we are able to accurately estimate the variations in SOA viscosity that result from the evaporation of its components. These variations are, in turn, influenced by atmospherically relevant changes in tOM and T. Comparison of the calculated SOA viscosity and diffusivity values with literature reported experimental results shows close agreement, thereby validating the employed calculation approach. These findings underscore the significant potential for TPD-DART-HRMS measurements in enabling the untargeted analysis of organic molecules within OA mixtures. This approach facilitates quantitative assessment of their gas-particle partitioning and allows for the estimation of their viscosity and condensed-phase diffusion, thereby contributing valuable insights to atmospheric models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

High-Pressure Synthesis of Bulk Cobalt Cementite, Co 3 C

Transition metal carbides find use in a wide range of advanced high-resilience applications including high-strength steels, heat shields, and deep-earth drills. However, carbides of the mid-to-late transition metals remain difficult to isolate and characterize on account of their metastability, which precludes the preparation of high-quality bulk single crystal samples using traditional solid-state methods. Herein, we report a combined computational and experimental survey of the cobalt-carbon binary system under high pressures and demonstrate that pressure offers a route toward the bulk synthesis of the metastable cementite-type cobalt carbide, Co 3 C, which under ambient conditions can only be prepared in low-dimensional thin film or nanoparticle forms. First-principles calculations reveal two competitive low-energy stoichiometric phases under ambient pressures - Pnnm-Co 2 C (Fe 2 C-type) and Pnma-Co 3 C (Fe 3 Ctype) - consistent with the known low-dimensional phases that have been studied for their promising magnetic properties. However, the calculated formation enthalpy of Pnma-Co 3 C decreases steadily with the applied pressure, while that of Pnnm-Co 2 C increases. Here, we pursue these results using high-pressure laser-heated synthesis methods coupled with in situ X-ray diffraction and observe the formation of Pnma-Co 3 C above 4.8 GPa. We determine the experimental bulk modulus of Co 3 C to be K o = 237 GPa (K p = 4.0). First-principles calculations of the phonon modes in Co 3 C reveal dynamical instabilities at ambient pressure that are absent under compression. These results offer a promising new route for the synthesis of rare-earth-free magnets.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Generic Behavior of Ultrastability and Anisotropic Molecular Packing in Codeposited Organic Semiconductor Glass Mixtures

Vapor-deposited glass mixtures of organic semiconductors commonly serve as active layers in organic electronic devices, whose lifetime and performance are strongly influenced by the stability and structure of these mixed glasses. Here, we study the stability and anisotropic molecular packing of six co-deposited organic semiconductor glass mixtures with 50:50 weight ratio, by differential scanning calorimetry and spectroscopic ellipsometry. We also find that all six binary systems exhibit high kinetic stability and significantly reduced enthalpy relative to the corresponding liquid-cooled glassy mixtures (ultrastable behavior), even for systems where the glass transition temperatures of the components differ by more than 90 K. Furthermore, we demonstrate that the birefringence of a co-deposited glass mixture, a measure of its anisotropic packing, can be predicted from the birefringence of glasses of the two pure components. These results for stability and structure are expected to be applicable to other co-deposited organic semiconductor glass mixtures, so long as the two components mix well in the glass and individually can form ultrastable glasses. Therefore, our findings are significant for designing novel electronic devices with enhanced device lifetime and increased operational efficiency.

36 MATERIALS SCIENCE↗

DOE OSTI2026

A Review on Direct Air Capture of Carbon Dioxide: Sorbent Materials, Process Engineering, Industrial Scale-Up, and Future Perspectives

The relentless accumulation of anthropogenic greenhouse gases has driven atmospheric carbon dioxide concentrations to approximately 426 ppm, necessitating the aggressive deployment of negative-emission technologies to achieve net zero by 2050. Direct air capture (DAC) offers a scalable, location-independent approach to atmospheric carbon removal; however, it is fundamentally constrained by the significant thermodynamic barriers associated with capturing CO 2 from ultradilute ambient conditions, requiring minimum thermodynamic energy inputs substantially higher than those for postcombustion point sources. This comprehensive review critically examines the technological landscape of DAC, focusing on the interdependent triad of sorbent material design, contactor engineering, and regeneration thermodynamics. We evaluate the fundamental boundaries of adsorption, emphasizing that an optimal adsorption enthalpy and isosteric heat of adsorption must balance the high CO 2 uptake capacity with the energetic penalties of sorbent regeneration. A systematic, comparative analysis of state-of-the-art sorbents is presented, encompassing mesoporous silicas, zeolites, carbon-based materials (CBMs), metal−organic frameworks (MOFs), porous organic polymers (POPs), and polymeric membranes. Special attention is devoted to surface functionalization strategies, particularly amine grafting and impregnation, which transition capture mechanisms from physisorption to chemisorption to enhance selectivity under ambient moisture and low partial pressures. Furthermore, we assess the operational merits of various reactor configurations, including gas−solid, gas−liquid, and membrane contactors, alongside regeneration cycles such as temperature, vacuum, pressure, and moisture swing adsorption. Finally, the review bridges fundamental materials science with industrial application by chronicling the scale-up milestones of pioneering entities and providing a strategic roadmap for advancing DAC technology readiness levels toward global deployment.

Adsorption↗

DOE OSTI2026

Optimal Design and Techno-Economic Analysis of 3D-Printed, Intensified Packings for Absorbers and Strippers in Solvent-Based CO 2 Capture

A potential technology for the CO 2 absorption process is utilizing intensified structured packing with embedded cooling/heating channels for continuous heat exchange, which can overcome limitations of discrete methods, such as discrete intercooling and centralized reboilers, to aid in reducing energy consumption and decreasing costs. This work investigates the modeling of intensified packing (IP) for the stripper tower, extending on previous work for the absorber, which distributes heat internally within the column, improving the thermodynamics for the solvent regeneration process. The model includes submodels for steam turbine extraction to produce steam at various qualities as well as a surrogate model for calculating steam enthalpy. A cost model for a plant-scale absorption capture process was developed, allowing for the design of the plant to be optimized, subject to minimizing capture cost using two different power plant flue gas sources. In this optimization, the placement of IP in both towers is optimized to balance the trade-off between enhanced heat transfer and reduced mass transfer volume. For natural gas combined cycle flue gas, the standard process configuration had a minimum cost of $\$$65.40/tonne CO 2 , and considering IP, the minimum capture cost is reduced to $\$$62.73/tonne, with utilization in the stripper column, which reduces yearly costs by up to $\$$2.67 MM/yr. Cooling the absorber through IP, or intercoolers, was only found to be beneficial at higher capture rates, with IP in both towers having a cost of capture of $\$$68.08/tonne at 99.9% capture, a reduction of $\$$12.64/tonne when using only intercoolers at the same capture rate. When capturing from pulverized-coal power plants, the minimum cost of capture when using IP in both towers is $\$$44.18/tonne (at 97% capture), while the standard configuration with and without intercoolers was $\$$45.69 and $\$$47.22 per tonne, respectively. This results in a reduction in yearly costs of $\$$16.98 MM/yr from the base-case configuration. At this higher CO 2 concentration, cooling in the absorber from the IP becomes extremely beneficial, reducing energy consumption by up to 6%.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2023

Comparison of Heterogeneous Langmuirian Models for Mixed-Gas Adsorption Equilibria

Mixed-gas adsorption remains a very active area of research, where newly proposed adsorption equilibrium models need to be compared to standards in the field. The generalized Langmuir (gL) isotherm, proposed by Hamid et al. to address adsorption of different size molecules on heterogeneous adsorbents without requiring evaluation of the Gibbs isotherm, has thus far been tested mainly against the adsorbed solution theory (AST). The gL isotherm was shown to be thermodynamically consistent and gave either better or similar representation than AST models. Here, the gL isotherm is contrasted with the Dual-Process Langmuir isotherm popularized by Ritter and his co-workers and the Loading Ratio Correlation popularized for simple mixtures, both of which are considered the “standard” models for process simulation. Furthermore, the comparison is based on their theoretical roots, interpretations, and representation of the equilibrium loading and the isosteric enthalpy of adsorption of unary, binary, and ternary mixtures of various gaseous adsorbates on different adsorbents over a range of experimental pressures and temperatures.

42 ENGINEERING↗

DOE OSTI2021

The Role of Water and Hydroxyl Groups in the Structures of Stetindite and Coffinite, MSiO 4 (M = Ce, U)

Orthosilicates adopt the zircon structure types ( I 4 1 /amd ), consisting of isolated SiO 4 tetrahedra joined by A-site metal cations, such as Ce and U. They are of significant interest in the fields of geochemistry, mineralogy, nuclear waste form development, and material science. Stetindite (CeSiO 4 ) and coffinite (USiO 4 ) can be formed under hydrothermal conditions despite both being thermodynamically metastable. Water has been hypothesized to play a significant role in stabilizing and forming these orthosilicate phases, though little experimental evidence exists. To understand the effects of hydration or hydroxylation on these orthosilicates, in situ high-temperature synchrotron and laboratory-based X-ray diffraction was conducted from 25 to ~850 °C. Stetindite maintains its I 4 1 /amd symmetry with increasing temperature but exhibits a discontinuous expansion along the a- axis during heating, presumably due to the removal of water confined in the [001] channels, which shrink against thermal expansion along the a -axis. Furthermore, additional in situ high-temperature Raman and Fourier transform infrared spectroscopy also confirmed the presence of the confined water. Coffinite was also found to expand nonlinearly up to 600 °C and then thermally decompose into a mixture of UO 2 and SiO 2 . A combination of dehydration and dehydroxylation is proposed for explaining the thermal behavior of coffinite synthesized hydrothermally. Additionally, we investigated high-temperature structures of two coffinite-thorite solid solutions, uranothorite (U x Th 1– x SiO 4 ), which displayed complex variations in composition during heating that was attributed to the negative enthalpy of mixing. Lastly, for the first time, the coefficients of thermal expansion of CeSiO 4 , USiO 4 , U 0.46 Th 0.54 SiO 4 , and U 0.9 Th 0.1 SiO 4 were determined to be α V = 14.49 × 10 –6 , 14.29 × 10 –6 , 17.21 × 10 –6 , and 17.23 × 10 –6 °C –1 , respectively.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Contrasting Trivalent Lanthanide and Actinide Complexation by Polyoxometalates via Solution-State NMR

Deciphering the solution chemistry and speciation of actinides is inherently difficult due to radioactivity, rarity, and cost constraints, especially for transplutonium elements. In this context, the development of new chelating platforms for actinides and associated spectroscopic techniques is particularly important. In this study, we investigate a relatively overlooked class of chelators for actinide binding, namely, polyoxometalates (POMs). We provide the first NMR measurements on americium–POM and curium–POM complexes, using one-dimensional (1D) 31 P NMR, variable-temperature NMR, and spin-lattice relaxation time (T 1 ) experiments. The proposed POM–NMR approach allows for the study of trivalent f-elements even when only microgram amounts are available and in phosphate-containing solutions where f-elements are typically insoluble. The solution-state speciation of trivalent americium, curium, plus multiple lanthanide ions (La 3+ , Nd 3+ , Sm 3+ , Eu 3+ , Yb 3+ , and Lu 3+ ), in the presence of the model POM ligand PW 11 O 39 7– was elucidated and revealed the concurrent formation of two stable complexes, [M III (PW 11 O 39 )(H 2 O) x ] 4– and [M III (PW 11 O 39 ) 2 ] 11– . Interconversion reaction constants, reaction enthalpies, and reaction entropies were derived from the NMR data. The NMR results also provide experimental evidence of the weakly paramagnetic nature of the Am 3+ and Cm 3+ ions in solution. Furthermore, the study reveals a previously unnoticed periodicity break along the f-element series with the reversal of T 1 relaxation times of the 1:1 and 1:2 complexes and the preferential formation of the long T 1 species for the early lanthanides versus the short T 1 species for the late lanthanides, americium, and curium. Furthermore, given the broad variety of POM ligands that exist, with many of them containing NMR-active nuclei, the combined POM–NMR approach reported here opens a new avenue to investigate difficult-to-study elements such as heavy actinides and other radionuclides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Metal Atom (Dis)Order and Superconductivity in YCaH n ( n = 8–20) High-Pressure Superhydrides

High-pressure superhydrides have attracted much attention due to their high superconducting critical temperatures (T c s). Our density functional theory (DFT) calculations, focusing on YCaH n (n = 8–20) compositions, found a number of nearly isoenthalpic YCaH 8 phases, differing only in the arrangement of the metal atoms, suggesting the potential stability of metal alloy superhydrides. The computed T c s of the considered YCaH 8 phases were higher than those of the isostructural I 4/ mmm MH 4 parent compounds. DFT enthalpies suggested that YCaH 12 could also be disordered; however, the T c s of the ordered variants spanned a wide range from 105 to 253 K at 200 GPa, showing that alloying could either mildly enhance or drastically reduce T c from that of the $Im\bar{3}m$ MH 6 parents. Finally, for YCaH 18 and YCaH 20 , only a single dynamically stable ordered superhydride was found, which we attribute to the differences in the structures of the most stable MH 9 and MH 10 binary hydride parents.

anions↗

DOE OSTI2020

Dynamics of Cation-Induced Conformational Changes in Nanometer-Sized Uranyl Peroxide Clusters

Conformational changes of the pyrophosphate (Pp)-functionalized uranyl peroxide nanocluster [(UO 2 ) 24 (O 2 ) 24 (P 2 O 7 ) 12 ] 48– ({U 24 Pp 12 }), dissolved as a Li/Na salt, can be induced by the titration of alkali cations into solution. The most symmetric conformer of the molecule has idealized octahedral ( O h ) molecular symmetry. One-dimensional 31 P NMR experiments provide direct evidence that both K + and Rb + ions trigger an O h -to- D 4 h conformational change within {U 24 Pp 12 }. Variable-temperature 31 P NMR experiments conducted on partially titrated {U 24 Pp 12 } systems show an effect on the rates; increased activation enthalpy and entropy for the D 4 h -to- O h transition is observed in the presence of Rb + compared to K + . Two-dimensional, exchange spectroscopy 31 P NMR revealed that magnetization transfer links chemically unique Pp bridges that are present in the D 4 h conformation and that this magnetization transfer occurs via a conformational rearrangement mechanism as the bridges interconvert between two symmetries. The interconversion is triggered by the departure and reentry of K (or Rb) cations out of and into the cavity of the cluster. This rearrangement allows Pp bridges to interconvert without the need to break bonds. Cs ions exhibit unique interactions with {U 24 Pp 12 } clusters and cause only minor changes in the solution 31 P NMR signatures, suggesting that O h symmetry is conserved. Single-crystal X-ray diffraction measurements reveal that the mixed Li/Na/Cs salt adopts D 2 h molecular symmetry, implying that while solvated, this cluster is in equilibrium with a more symmetric form. Thus, these results highlight the unusually flexible nature of the actinide-based {U 24 Pp 12 } and its sensitivity to countercations in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to stop unphysical rearrangements during geometry relaxation. A simple fixed-atom or “coordinate-lock” approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

QC-GN 2 oMS 2 : a Graph Neural Net for High Resolution Mass Spectra Prediction

Predicting the mass spectrum of a molecular ion is often accomplished via three generalized approaches: rules-based methods for bond breaking, deep learning, or quantum chemical (QC) modeling. Rules-based approaches are often limited by the conditions for different chemical subspaces and perform poorly under chemical regimes with few defined rules. QC modeling is theoretically robust but requires significant amounts of computational time to produce a spectrum for a given target. Among deep learning techniques, graph neural networks (GNNs) have performed better than previous work with fingerprint-based neural networks in mass spectra prediction. To explore this technique further, we investigate the effects of including quantum chemically derived information as edge features in the GNN to increase predictive accuracy. The models we investigated include categorical bond order, bond force constants derived from extended tight-binding (xTB) quantum chemistry, and acyclic bond dissociation energies. Throughout this work, we evaluated these models against a control GNN with no edge features in the input graphs. Bond dissociation enthalpies yielded the best improvement with a cosine similarity score of 0.462 relative to the baseline model (0.437). In this work we also apply dynamic graph attention which improves performance on benchmark problems and supports the inclusion of edge features. Between implementations, we investigate the nature of the molecular embedding for spectra prediction and discuss the recognition of fragment topographies in distinct chemistries for further development in tandem mass spectrometry prediction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗