Search NASA⌕ Search

SEARCH · Search NASA

Results for “minimum”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Advancing the Performance of Lithium-Rich Oxides in Cooperation with Inherent Complexities: Engineering Domain Structure

The nanodomain structure of lithium- and manganese-rich, composite cathode materials has been systematically altered through synthesis conditions while keeping larger-scale morphological differences to a minimum. Clear changes in electrochemical performance across the samples studied are observed, especially with respect to the anomalous impedance at low states-of-charge. Atomic-scale modeling, coupled to electrochemical measurements and physical characterization, reveals the local consequences of the high-voltage activation charge process as a function of specific domain structures. Furthermore, the results explain the influence of different postactivation structures on the insertion of Li-ions and the associated impedance during discharge. This work adds to our series of studies on lithium- and manganese-rich oxides, demonstrating that the inherent performance of LMRs can be greatly enhanced through rational, highly controlled synthetic strategies based on an understanding of synthesis–structure–property relationships across length scales.

25 ENERGY STORAGE↗
DOE OSTI2025

Deepening the LUMO: Brominated Naphthalene Diimide Electron Transport Layers for Low-Hysteresis Perovskite Solar Cells

Precise energy level alignment at the interfaces between the charge transport layers, active layer, and electrodes plays a key role in maximizing photovoltaic performance and operational stability in perovskite solar cells (PSCs). Organic electron transport layers (ETLs) have received little attention compared to their inorganic counterparts but offer the unique advantage of facile functionalization for fine-tuning of electronic properties. Here, we report the design, synthesis, and characterization of two benzyl-phosphonic acid (BnPA)-functionalized naphthalene diimide (NDI) derivatives, NDI-(BnPA) 2 and Br 2 -NDI-(BnPA) 2 as organic ETLs for PSCs in an n-i-p device configuration. Bromination of the NDI core at the 4,9-positions deepens the lowest unoccupied molecular orbital (LUMO) in Br 2 -NDI-(BnPA) 2 by 0.29 eV, enabling improved energy alignment with the conduction band minimum of the perovskite. Devices incorporating Br 2 -NDI-(BnPA) 2 demonstrated enhanced short-circuit current density (JSC), reduced hysteresis, and a maximum power conversion efficiency of 13.67%, compared to 13.20% for the unsubstituted NDI-(BnPA) 2 analog. The deeper LUMO of Br 2 -NDI-(BnPA) 2 is hypothesized to facilitate more efficient electron extraction, suppress interfacial charge accumulation, and reduce field-driven ion migration, collectively contributing to the observed reduction in hysteresis. These results highlight the effectiveness of combining molecular-level LUMO tuning with robust interfacial anchoring to advance the performance and durability of organic ETLs in PSCs.

Deposition↗
DOE OSTI2024

Numerical Simulations in Support of a Long-Term Test of Gas Production From Hydrate Accumulations on the Alaska North Slope: Water Production and Associated Design and Management Issues

Here, we investigated numerical simulation strategies for a long-term test of depressurization-induced gas production from the B1 Sand of Unit B at the Hydrate-01 Stratigraphic Test Well. The main objective of this study was to estimate fluid production rates (with emphasis on water production) under a variety of conditions and production scenarios and contribute new insights to the design and management of the field test. In the first part of the study, we investigated the system response to a three-step depressurization process using two limiting sets of flow properties─the expected maximum and minimum intrinsic and effective permeabilities─for the very heterogeneous reservoir. In the second part, we investigated the effect of the production interval length and placement within the formation relative to the boundaries of the hydrate-bearing unit. The best performing well configuration was used in the third part of the study, which used the most representative subsurface flow properties to investigate the effect of the depressurization strategy on the production performance. The best overall performance (largest gas production with modest water production and a strong response at the observation wells) was obtained with a 10 m-long well situated 3 m below the top of the formation and a three-step depressurization scheme at 15-day intervals to a terminal bottomhole pressure of 2.8 MPa. The overall production performance was enhanced by a faster rate of depressurization. Estimated water production rates in all cases were limited and easily manageable. None of the tested well configurations or depressurization strategies significantly reduced water production without also severely reducing gas production. In all the investigated cases, 95% of the long-term fraction of produced water was replenished by inflows from the boundaries and could not be reduced. These substantial water inflows are an unavoidable feature of HU-B and cannot be easily mitigated by a hydraulic control.

02 PETROLEUM↗
DOE OSTI2025

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Process Feasibility Analysis of Waste Biomass Valorization to Biochar and Bio-Oil via Slow and Fast Pyrolysis

The United States has abundant biomass and waste feedstock to support the nation's energy addition and affordability targets. Pyrolysis, a thermochemical conversion process, decomposes lignocellulosic feedstocks into liquid, solid, and gaseous fuels that can contribute to the domestic production of biofuels, biopower, and bioproducts. Growing private sector interest in this technology is a key motivation for this comprehensive techno-economic process modeling analysis of a respective biorefinery that includes feedstock preprocessing, slow and fast pyrolysis, and product separation to bio-oil, biochar, and syngas hydrocarbons. Results show that biochar from slow pyrolysis could achieve minimum selling prices (MSPs) of $\$$188-$\$$260/t, competitive with reported market values, while bio-oil from fast pyrolysis is estimated to yield MSPs of $\$$6.49-$\$$9.68/GGE, approximately twice conventional fuel benchmarks. Sensitivity analysis identifies feedstock cost, product yield, and scale as primary cost drivers, while scenarios involving biochar carbon credits and high value applications may substantially improve economics. Overall, these results suggest that continued innovation in feedstock logistics, process integration, and market development will be critical to achieving economically viable and scalable bioproducts.

09 BIOMASS FUELS↗
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

Influence of CO 2 -Regenerative Film Properties in Enhancing C 2+ Products Selectivity While Mitigating CO 2 Crossover

Zero-gap anion-exchange membrane electrode assembly (AEMEA) electrolyzers operating in alkaline media face challenges such as CO 2 crossover and salting-out. The bipolar membrane electrode assembly (BPMEA) electrolyzer, using DI water as the electrolyte, addresses both CO 2 crossover and salting-out issues. As cations, which are crucial in stabilizing the CO 2 R intermediates, are absent in the electrolyte, cations embedded in the AEM play an important role in dictating the CO 2 R selectivity and activity in BPMEA systems. So far, no systematic study has been conducted on the influence of cations embedded in AEMs on CO 2 selectivity and activity in BPMEA systems. Moreover, BPMEA systems impose an additional challenge: low stability due to delamination of the bipolar membrane caused by CO 2 regeneration at the membrane-membrane interface. To enhance the stability of the electrolyzer, a simple yet highly reproducible strategy for coating a porous CO 2 regenerative film on smooth Nafion 117 is demonstrated in this work, along with a systematic study of four different commercially available AEMs, composed of different cations and cation densities, for use in combination with Nafion 117 and copper catalyst at the cathode. We found that the PiperION membrane delivers selectivity and activity comparable to those of alkali-metal cations, owing to the enhanced local electric field resulting from the combined effects of a high positive charge on the N atom of the piperidinium cation and the high ion-exchange capacity. Further, we studied the influence of the thickness of the PiperION porous layer over Nafion 117 on CO 2 R selectivity and found that 70 μm is the minimum thickness to achieve maximum C 2+ products selectivity, reduced the CO 2 crossover to 5% from 25% at 4 SCCM and 150 mA cm –2 , and was stable for operation beyond 100 h. The selectivity of this system, compared with AEMEA, and stability outperformed both AEMEA and BPMEA. This study helps design more effective BPMs to inhibit CO 2 crossover while enabling stable and selective electrochemical CO 2 reduction to C 2+ hydrocarbons.

Cations↗
DOE OSTI2026

Superstructure Optimization of Waste Plastic Pyrolysis, Integrating Thermal, Catalytic, and Plasma Technologies with Machine Learning

Global plastic waste generation exceeds 430 million tonnes per year, yet fewer than 9% are recycled in the United States. Pyrolysis offers a chemical recycling route at scale, but existing techno-economic and life cycle assessments fix product yields to single pure polymers, producing economic and environmental outputs that break down when the feed composition changes. Here, we present a superstructure optimization framework that addresses this by embedding a composition-aware random forest yield predictor, trained on 566 pyrolysis experiments, within a full-scale process simulation. Product distributions update automatically as feed allocation shifts across four reactor chemistries: conventional thermal, catalytic (HZSM-5), thermal oxo-degradation, and nonequilibrium CO2 plasma. The optimal superstructure achieves minimum selling prices of −0.56 to −0.76/kg feed and global warming potentials of −0.276 to −0.322 kg CO2-eq/kg feed across four commodity price scenarios, confirming profitable, carbon-negative operation without tipping fees. Carbon abatement costs of $\$$0.46 to $\$$1.25/kg CO2-eq are competitive with direct air capture. Sensitivity analysis shows that the catalytic-plasma split fraction is the single largest driver of both economic and climate performance, while hydrocracking allocation in the wax upgrading stage is emission-neutral across the full variable range. Mixed plastic waste streams, evaluated as composition-variable feedstocks rather than pure resins, are profitable and carbon-negative across realistic market conditions. These results give a quantitative basis for reactor selection, circular economy investment, and policy design targeting chemical recycling on a large scale.

Life cycle assessment↗
DOE OSTI2025

Sulfate Promotes Compact CaCO 3 Formation and Protects Portland Cement from Supercritical CO 2 Attack

Supercritical (sc) CO 2 in geologic carbon sequestration (GCS) can chemically and mechanically deteriorate wellbore cement, raising concerns for long-term operations. In contrast to the conventional view of “sulfate attack” on cement, we found that adding 0.15 M sulfate to the acidic brine can significantly reduce the impact of scCO 2 attack on Portland cement, resulting in stronger cement than that found in a sulfate-free system. Scanning electron microscopy revealed a decreased total attack depth in reacted cement in the presence of sulfate. With a newly defined minimum porosity term in reactive transport modeling, our model suggests that sulfate caused CaCO 3 to fill more nanopore spaces in the cement. Small angle X-ray scattering experiments also showed that sulfate can decrease the pore sizes of the carbonate layer. The results suggest that the interactions between sulfate and cement can generate a less porous CaCO 3 layer, which better resists acidic brine. Using this mechanism as a proof-of-concept, we tested the incorporation of sodium sulfate into Portland cement and synthesized new cement composites that show stronger resistance against scCO 2 attacks. Finally, these newly discovered interfacial interactions between CaCO 3 and sulfate provide new insights into engineering mechanically strong and green materials for safer GCS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Techno-Economic Analysis and Life Cycle Assessment of Alternative Fuels for Locomotives in the U.S. Freight Rail Sector

Freight rail is more energy-efficient than truck transport over long-haul distances, offering a low-energy and emissions-intensive option for transporting freight. This study evaluates techno-economic analysis and life cycle assessment of seven alternative unblended fuels for freight locomotive engines─biodiesel, renewable diesel (RD), bio-oils, methanol, dimethyl ether (DME), ethanol, and ammonia─across 16 fuel pathways utilizing soybean, corn, woody biomass, renewable hydrogen, and waste sources, e.g., sludge, manure, and industrial CO 2 , and compares these to conventional diesel. The minimum fuel selling price (MFSP) ranged from $\$2.05$ to $\$8.27$ per diesel gallon equivalent (2020 US dollars), with biocrude and RDs produced from hydrothermal liquefaction (HTL) of sludge having the lowest MFSPs due to coproduct credits and avoided waste treatment cost. Life cycle GHG emissions ranged from −41 to 53 g of CO 2 e/MJ. RD from waste via HTL achieves negative emissions by diverting sludge/manure from GHG-intensive conventional management. Few pathways such as biocrude, methanol, and DME require additional control for SO X emissions in the refinery, while ethanol, FT-diesel, and bio-oil require additional control for particulate matter emissions. Bio-oil and RD from sludge have lower marginal abatement cost or MAC (–$\$38$/tonne CO 2 lowest) while methanol and ammonia with renewable hydrogen have higher MAC ($\$490$/tonne CO 2 maximum).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Cost and Carbon Intensity Implications of Coprocessing Sustainable Aviation Fuel at Petroleum Refineries

Sustainable aviation fuel (SAF) will play a critical role in decarbonizing the aviation industry. Among SAF production pathways, alcohol-to-jet (ATJ) stands out for its scalability, supported by abundant feedstock availability and a well-established bioethanol industry. However, significant reductions in SAF carbon intensity (CI) require the use of future feedstocks (e.g., cellulosic) whose adoption is hindered by high capital costs for feedstock processing and ethanol upgrading. Here, we evaluate the financial viability and environmental implications of integrating an ATJ SAF biorefinery within a petroleum refinery, utilizing miscanthus and switchgrass as example feedstocks. Three scenarios are evaluated: standalone (benchmark), colocated, and repurposing (coprocessing SAF within the petroleum refinery). Results show repurposing reduces baseline capital costs by 36% and SAF minimum selling price (MSP) by 12% to 8.14 USD·gal −1 ; the superior performance of repurposing is consistent across both feedstocks. Integration has a limited effect on SAF CI, which remains stable across scenarios, whereas using cellulosic feedstocks reduces CI by over 70% relative to corn, with baseline values of 17.01 g CO 2 e· MJ −1 for miscanthus and 12.23 g CO 2 e·MJ −1 for switchgrass. Global sensitivity analysis reveals MSP declines with greater coprocessing levels.

09 BIOMASS FUELS↗
DOE OSTI2026

A Kinetic Model-Driven Techno-Economic Analysis of Plastic Pyrolysis: Linking Process Dynamics to Economic Viability

This study employs a kinetic model integrated into Aspen Plus to predict pyrolysis product distribution under various conditions. A techno-economic assessment calculated the minimum selling price (MSP) of pyrolysis oil under different operating conditions for the baseline capacity of 100 kta, and across eight processing capacities ranging from 30 to 150 kta. The lowest MSP under the baseline capacity is estimated at $\$$420/ton, which is 33% lower than the 2023 average US crude oil price ($\$$74.6/bbl, equivalent to $\$$634/ton based on the density of pyrolysis oil). Under Monte Carlo simulation, accounting for variability in key economic and technical parameters, the mean MSP is estimated at $\$$1137/ton. The economic viability depends on feedstock price remaining below $\$$320/ton, defining the break-even feedstock price threshold. Sensitivity analysis further identifies capital investment and transportation cost as key economic drivers. Capacities beyond 90 kta show limited economies of scale benefits. Reducing product storage time cuts capital costs by 7% but raises operational risk. Uncertainty analysis suggests the economic feasibility of pyrolysis oil is unlikely to compete with crude oil without policy incentives.

petrochemicals↗
DOE OSTI2025

Fermi’s Golden Rule Rate Expression for Transitions Due to Nonadiabatic Derivative Couplings in the Adiabatic Basis

Starting from a general molecular Hamiltonian expressed in the basis of adiabatic electronic and nuclear position states, where a compact and complete expression for the nonadiabatic derivative coupling (NDC) Hamiltonian term is obtained, we provide a general analysis of the Fermi’s golden rule (FGR) rate expression for nonadiabatic transitions between adiabatic states. We then consider a quasi-adiabatic approximation that uses crude adiabatic states and NDC couplings, both evaluated at the minimum potential energy configuration of the initial adiabatic state, for the definition of the zeroth and first-order terms of the Hamiltonian. Although the application of this approximation is rather limited, it allows deriving a general FGR rate expression without further approximation while accounting for non-Condon contribution to the FGR rate arising from momentum operators of NDC terms and its coupling with vibronic displacements. For a generic and widely used model where all nuclear degrees of freedom and environmental effects are represented as linearly coupled harmonic oscillators, we derive a closed-form FGR rate expression that requires only Fourier transform. The resulting rate expression includes quadratic contributions of NDC terms and their couplings to Franck–Condon modes, which require evaluation of two additional bath spectral densities in addition to the conventional one that appears in a typical FGR rate theory based on the Condon approximation. Model calculations for the case where nuclear vibrations consist of both a sharp high-frequency mode and an Ohmic bath spectral density illustrate new features and implications of the rate expression. We then apply our theoretical expression to the nonradiative decay from the first excited singlet state of azulene, which illustrates the utility and implications of our theoretical results.

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

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2024

Comparative Study of Neutral and Cationic Sn 2 H 2 : Toward Laboratory Detection of the Cation

Group 14 M 2 H 2 isomers (M = Si, Ge, Sn, and Pb) have attracted interest due to their radically differing electronic structures from acetylene. To better understand the Sn–H interactions of the neutral and cationic Sn 2 H 2 structures, we present the most rigorous study of these systems to date. CCSD(T)/cc-pwCVTZ harmonic frequencies are presented as the first predictions for the neutral and cationic species to date. CCSDT(Q)/CBS relative energies are reported using the focal point approach, confirming the butterfly isomer as the global minimum on the potential energy surface for both the neutral and cationic species. In all, there exist 7 minima and 15 transition states. NBO analysis is also performed to elucidate the changes in bond order going from neutral to cation across all isomers of Sn 2 H 2 . Our results provide insights into the important Sn–H interaction and provide guidance for future work that may detect Sn 2 H$^{+}_{2}$ in the laboratory for the first time.

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

Toward a Machine Learning Approach to Interpreting X-ray Spectra of Trace Impurities by Converting XANES to EXAFS

The fact that the photoabsorption spectrum of a material contains information about the atomic structure, commonly understood in terms of multiple scattering theory, is the basis of the popular extended X-ray absorption spectroscopy (EXAFS) technique. How much of the same structural information is present in other complementary spectroscopic signals is not obvious. Here we use a machine learning approach to demonstrate that within theoretical models that accurately predict the EXAFS signal, the extended near-edge region does indeed contain the EXAFS-accessible structural information. We do this by exhibiting deep operator neural networks (DeepONets) that have learned the relationship between the extended and near edge portions of the X-ray absorption spectrum to predict the former from the latter. We find that we can accurately predict the EXAFS spectrum between 6 and 14 Å –1 from the first 6 Å –1 (≈100 eV) of the absorption spectrum of Cu 2 + substitutional defects in the Fe 3+ mineral hematite (α-Fe 2 O 3 ). This surprising finding implies that theoretical analyses of X-ray absorption spectra could be implemented that extract the same conclusions as high-quality EXAFS studies from spectra collected over a much smaller range of photon energies. This relaxes a host of experimental limitations related to the X-ray source and measurement sample, including collection time, minimum dopant concentration, source brilliance, and energy range. We describe the theoretical data sets and DeepONet construction and show that the resulting DeepONets produce EXAFS that recovers linear combination fits to experimental data with accuracy approaching the original ab initio calculations. We discuss the implications of our findings for minor constituent characterization and for understanding the information content of spectroscopic data more broadly, including how this approach might be applied to measured experimental spectra. In conclusion, to encourage similar efforts, the simulated X-ray spectra, machine learning, and fitting code are publicly available.

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

Low-Energy Isomers of the Magic Number H + (H 2 O) 21 Cluster

Electronic structure calculations are used to characterize low-energy isomers of H + (H 2 O) 21 . Eleven different classes of isomers, based on the (H 2 O) 20 pentagonal dodecahedron with the excess proton localized on the surface (as a hydrated hydronium ion) and the “extra” water molecule located in the interior of the cluster, are characterized. In 10 of these classes, the internal water molecule is engaged in six 5-membered rings, but in the remaining class, which is predicted to start at only 0.6 kcal/mol above the global minimum, the internal water is engaged in a 4- membered ring, an additional 6-membered ring, and four 5- membered rings. In addition, isomers with two 4-membered rings and two 6-membered rings on the cluster surface are predicted to start at only ∼1.3 kcal/mol above the lowest-energy dodecahedralbased structure.

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

A Comparison of Electronic Structure Methods for Predicting the Hydrogenation Energies of Candidate Molecules for Hydrogen Storage

The development of novel energy materials and fuels is required to expand current available energy sources. Aiming to reach this goal, there is growing interest in using molecular hydrogen as an energy carrier due to its abundance and high energy density. Liquid organic hydrogen carriers (LOHCs) are a promising route to the large-scale storage and transport of hydrogen for use in the energy economy. The search for thermodynamically viable LOHC molecules for real world use has led to a set of constraints on the dehydrogenation enthalpy and the minimum gravimetric hydrogen capacity. These constraints allow one to formulate the search for an ideal LOHC candidate molecule as an optimization problem well suited to the strengths of machine learning and artificial intelligence computational approaches. A critical barrier to a large-scale, high-throughput screening of LOHC candidate molecules is the lack of reliable training data. Computational electronic structure methods including density functional theory, coupled cluster approximations, and diffusion Monte Carlo can be used to provide training data where experimental data are either unreliable or do not exist. In this work, we use these methods to calculate the dehydrogenation energies and enthalpies of candidate LOHC molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗