Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 OSTI2021

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

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

Energy Barriers for H 2 Uptake and Release by Defects in Boron Nitride

Density functional theory simulations (DFT) have been used in combination with the climbing image nudged elastic band method (CI-NEB) to determine the energy barriers required to dissociate H 2 at four types of h-BN defects and evaluate the potential for reacted H species to hop to adjacent available sites. The defects were selected based on previous work in which ab initio thermodynamics was used to identify specific defects with optimal thermodynamics for reversible H 2 storage, that is, the defects would dissociate molecular H 2 while not binding H species too strongly so that hydrogen release is favorable. The four defects investigated are the N and B monovacancies, which only offer either B or N defective sites respectively, as well as larger cluster vacancy defects named 3V(1B2N) and hexagonal 6V(3B3N), which present a mixture of B and N defect sites available for H binding. While homogeneous H 2 dissociation by B and N monovacancies have energy barriers larger than 0.4 eV, we found that mixed BN terminated defects have energy barrier lower than 0.4 eV for H 2 dissociation. In particular, H 2 dissociation by the hexagonal 6V(3B3N) defect, exclusively made of Frustrated Lewis Pairs (FLPs), is almost barrierless. The energy barrier for subsequent H hopping to available adjacent B or N sites of the defect were also calculated. It was found that the hopping of H species generally involves high energy barriers ranging from 0.26 eV to 2.57 eV. In conclusion, this work completes previous thermodynamic theoretical investigations and provides a more comprehensive picture of H 2 reaction with defective h-BN for potential H 2 storage and recovery.

Sassi, Michel J. [Pacific Northwest National Labor↗
DOE OSTI2024

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗
DOE OSTI2023

Active Site Engineering via Optimizing the Heterogeneous Support Structure for Single-Atom Catalysis

Supported single-atom catalysts show a large range of activities and selectivities that depend on the local environment of the catalytic sites. Here a theory-based optimization strategy is presented that is based on a density functional theory determination of the transition states and intermediates for a low-dimensional coordinate representation of the heterogeneity of the active sites. The approach is applied to a vanadium catalyst on an amorphous SiO 2 support that involves a large kinetic network described using a full chemistry model. Without assuming a priori scaling relations or mechanism reduction, the optimal state of heterogeneity is found to lie at atomic configurations where the activation energies for two distinct key chemical processes are equal. It is found a posteriori that the behavior of the system is consistent with linear free energy scaling relations in the randomness parameters. The energetic span theory proves quite useful in reducing the full chemistry model to a small number of key reactions. The use of a nonlinear optimization algorithm in combination with energetic span theory provides significant simplification in treating disordered systems.

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

Structural anisotropy in Sb thin films

Sb thin films have attracted wide interest due to their tunable band structure, topological phases, high electron mobility, and thermoelectric properties. We successfully grow epitaxial Sb thin films on a closely lattice-matched GaSb(001) surface by molecular beam epitaxy. We find a novel anisotropic directional dependence on their structural, morphological, and electronic properties. The origin of the anisotropic features is elucidated using first-principles density functional theory (DFT) calculations. The growth regime of crystalline and amorphous Sb thin films was determined by mapping the surface reconstruction phase diagram of the GaSb(001) surface under Sb 2 flux, with confirmation of structural characterizations. Crystalline Sb thin films show a rhombohedral crystal structure along the rhombohedral (211) surface orientation parallel to the cubic (001) surface orientation of the GaSb substrate. At this coherent interface, Sb atoms are aligned with the GaSb lattice along the [1̄10] crystallographic direction but are not aligned well along the [110] crystallographic direction, which results in anisotropic features in reflection of high-energy electron diffraction patterns, misfit dislocation formation, surface morphology, and transport properties. Our DFT calculations show that the preferential orientation of the rhombohedral Sb (211) plane may originate from the GaSb surface, where Sb atoms align with the Ga and Sb atoms on the reconstructed surface. The formation energy calculations confirm the stability of the experimentally observed structures. Our results provide optimal film growth conditions for further studies of novel properties of Bi 1-x Sb x thin films with similar lattice parameters and an identical crystal structure, as well as functional heterostructures of them with III–V semiconductor layers along the (001) surface orientation, supported by a theoretical understanding of the anisotropic film orientation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shape Matters: Understanding the Effect of Electrode Geometry on Cell Resistance and Chemo-Mechanical Stress

Rechargeable batteries that incorporate shaped three-dimensional electrodes have been shown to have increased power and energy densities when compared to a conventional geometry, i.e. a planar cathode and anode that sandwich an electrolyte. Electrodes can be shaped to enable a higher active material loading, while keeping ion transport distances small. However, the relationship between electrical and mechanical performance of shaped electrodes remains poorly understood. Many electrode designs have been explored, where the electrodes are individually shaped or intertwined, and advances in manufacturing and shape/topology optimization have made such designs a reality. Here, we explore sinusoidal half cells and interdigitated full cells. First, we use a simple electrostatics model to understand the cell resistance as a function of shape. Here we focus on low-temperature conditions, where the electrolyte conductivity decreases relative to that of the electrode; here, LiPF 6 EC:DMC electrolyte and MnO 2 electrode are considered. Next, we use a chemo-mechanics model to examine the stress that arises due to intercalation-driven volume expansion. We show that shaped electrodes provide a significant reduction in resistance in low-temperature conditions, however, they exhibit unfavorable stress concentrations. Overall, we find that the fully interdigitated electrodes may provide the best balance with respect to this resistance-stress trade-off.

25 ENERGY STORAGE↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

The Community Land Model Version 5: Description of New Features, Benchmarking, and Impact of Forcing Uncertainty

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

58 GEOSCIENCES↗
DOE OSTI2023

Non-reactive facet specific adsorption as a route to remediation of chlorinated organic contaminants

The present work quantifies metal-contaminant interactions between palladium substrates and three salient chlorinated organic contaminants, namely trichloroethylene 1,3,5-trichlorobenzene (TCB), and 3,3′,4,4′-tetrachlorobiphenyl (PCB77). Given that Pd is one of the conventional catalytically active materials known for contaminant removal, maximizing catalytic efficiency through optimal adsorption dynamics reduces the cost of remediation of contaminants that are persistent water pollutants chronically affecting public health. Adsorption efficiency analyses from all-atom molecular dynamics (MD) simulations advance the understanding of reaction mechanisms available from density functional theory (DFT) calculations to an extractable feature scale that can fit the parametric design of supported metal catalytic systems and feed into high throughput catalyst selection. Data on residence time, site-specific adsorption, binding energies, packing geometries, orientation profiles, and the effect of adsorbate size show the anomalous behaviour of organic contaminant adsorption on the undercoordinated {110} surface as compared to the {111} and {100} surfaces. The intermolecular interaction within contaminants from molecular dynamics simulation exhibits refreshing results than ordinary single molecule density functional theory calculation. Since complete adsorption and dechlorination is an essential step for chlorinated organic contaminant remediation pathways, the presented profiles provide essential information for designing efficient remediation systems through facet-controlled palladium nanoparticles.

Guo, Hao↗
DOE OSTI2022

Critical role of acceptor dopants in designing highly stable and compatible proton-conducting electrolytes for reversible solid oxide cells

Proton-conducting electrolytes are receiving increasing attention due to their high ionic conductivity at intermediate temperatures, enabling the operation of solid oxide cells with high energy efficiency at low cost. However, the effect of B-site dopants on the properties of doped barium hafnate-cerate electrolyte materials, especially in single cells under operating conditions, has not been systematically studied. Here we report our findings in the development of a series of proton-conducting electrolytes with a general formula of BaHf 0.1 Ce 0.7 R 0.2 O 3–δ (BHCR172, R = Yb, Er, Y, Gd, Sm). Here, the results reveal that electrical conductivity, ionic transference number, chemical stability against steam and CO 2 , and compatibility with NiO during sintering are all closely correlated with the dopant size. In particular, the reaction with NiO is found to strongly affect the properties of the electrolytes and hence cell performance. Among all tested compositions, BaHf 0.1 Ce 0.7 Yb 0.2 O 3–δ (BHCYb172) shows excellent chemical stability and minimal reactivity towards NiO, as predicted from density functional theory (DFT)-based calculations and confirmed by experimental results. In addition, proton-conducting reversible solid oxide cells (P-ReSOCs) based on the optimized electrolyte composition, BHCYb172, demonstrate exceptional performance and stability, achieving a remarkable peak power density of 1.74 W cm –2 (O 2 as the oxidant) at 600 °C in the fuel cell mode and a high current density of 2.0 A cm –2 at 1.3 V and 600 °C in the steam electrolysis mode while maintaining excellent durability for over 1000 h.

25 ENERGY STORAGE↗
DOE OSTI2025

Processability and Material Behavior of NiTi Shape Memory Alloys Using Wire Laser-Directed Energy Deposition (WL-DED)

Utilizing additive manufacturing (AM) techniques with shape memory alloys (SMAs) like NiTi shows great promise for fabricating highly flexible and functionally superior 3D metallic structures. Compared to methods relying on powder feedstocks, wire-based additive manufacturing processes provide a viable alternative, addressing challenges such as chemical composition instability, material availability, higher feedstock costs, and limitations on part size while simplifying process development. This study presented a novel approach by thoroughly assessing the printability of Ni-rich Ni55.94Ti (Wt. %) SMA using the wire laser-directed energy deposition (WL-DED) technique, addressing the existing knowledge gap regarding the laser wire-feed metal additive manufacturing of NiTi alloys. For the first time, the impact of processing parameters—specifically laser power (400–1000 W) and transverse speed (300–900 mm/min)—on single-track fabrication using NiTi wires in the WL-DED process was examined. An optimal range of process parameters was determined to achieve high-quality prints with minimal defects, such as wire dripping, stubbing, and overfilling. Building upon these findings, we printed five distinct cubes, demonstrating the feasibility of producing nearly porosity-free specimens. Notably, this study investigated the effect of energy density on the printed part density, impurity pick-up, transformation temperature, and hardness of the manufactured NiTi cubes. The results from the cube study demonstrated that varying energy densities (46.66–70 J/mm3) significantly affected the quality of the deposits. Lower to intermediate energy densities achieved high relative densities (>99%) and favorable phase transformation temperatures. In contrast, higher energy densities led to instability in melt pool shape, increased porosity, and discrepancies in phase transformation temperatures. These findings highlighted the critical role of precise parameter control in achieving functional NiTi parts and offer valuable insights for advancing AM techniques in fabricating larger high-quality NiTi components. Additionally, our research highlighted important considerations for civil engineering applications, particularly in the development of seismic dampers for energy dissipation in structures, offering a promising solution for enhancing structural performance and energy management in critical infrastructure.

Dabbaghi, Hediyeh↗
DOE OSTI2023

Crystal structure of meglumine diatrizoate, (C 7 H 18 NO 5 )(C 11 H 8 I 3 N 2 O 4 )

The crystal structure of meglumine diatrizoate has been solved and refined using synchrotron X-ray powder diffraction data and optimized using density functional theory techniques. Meglumine diatrizoate crystallizes in space group P2 1 (#4) with a = 10.74697(4), b = 6.49364(2), c = 18.52774(7) Å, β = 90.2263(3), V = 1292.985(5) Å 3 , and Z = 2. Two different crystal structures, which yielded essentially identical refinement residuals and positions of the non-H atoms, were obtained. The differences were in the H atom positions and the hydrogen bonding. One structure was 123.0 kJ/mol/cell lower in energy than the other and was adopted for the final description. The crystal structure consists of alternating double layers of cations and anions along the c-axis. The hydrogen bonds link the cations and anions into a three-dimensional framework. Each of the hydrogen atoms on the ammonium nitrogen of the cation acts as a donor in a strong N–H∙∙∙O hydrogen bond. One of these is to a hydroxyl group of another cation, and the other is to the carboxylate group of the anion. Each of the amide nitrogen atoms of the anion forms a strong N–H∙∙∙O intermolecular hydrogen bond, one to a carbonyl and the other to a carboxylate group. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigating magnetic van der Waals materials using data-driven approaches

In this work, we investigate magnetic monolayers of the form A i A ii B 4 X 8 based on the well-known intrinsic topological magnetic van der Waals (vdW) material MnBi 2 Te 4 (MBT) using first-principles calculations and machine learning techniques. We select an initial subset of structures to calculate the thermodynamic properties, electronic properties, such as the band gap, and magnetic properties, such as the magnetic moment and magnetic order using density functional theory (DFT). Data analytics approaches are used to gain insight into the microscopic origin of materials’ properties. The dependence of materials’ properties on chemical composition is also explored. For example, we find that the formation energy and magnetic moment depend largely on A and B sites whereas the band gap depends on all three sites. Finally, we employ machine learning tools to accelerate the search for novel vdW magnets in the MBT family with optimized properties. Finally, this study creates avenues for rapidly predicting novel materials with desirable properties that could enable applications in spintronics, optoelectronics, and quantum computing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water

The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Achieving performance portability in Gaussian basis set density functional theory on accelerator based architectures in NWChemEx

The numerical integration of the exchange–correlation (XC) potential is one of the primary computational bottlenecks in Gaussian basis set Kohn–Sham density functional theory (KS-DFT). To achieve optimal performance and accuracy, care must be taken in this numerical integration to preserve local sparsity as to allow for near linear weak scaling with system size. This leads to an integration scheme with several performance critical kernels which must be hand optimized for each architecture of interest. As the set of available accelerator hardware goes more diverse, a key challenge for developers of KS-DFT software is to maintain performance portability across a wide range of computational architectures. In this article, we examine a modular software design pattern which decouples the implementation details of performance critical kernels from the expression of high-level algorithmic workflows in a device-agnostic language such as C++; thus allowing for developers to target existing and emerging accelerator hardware within a single code base. We consider the efficacy of such a design pattern in the numerical integration of the XC potential by demonstrating its ability to achieve performance portability across a set of accelerator architectures which are representative of those on current and future U.S. Department of Energy Leadership Computing Facilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Facilitating Hydrogen Dissociation over Dilute Nanoporous Ti–Cu Catalysts

The dissociation of H 2 is an essential elementary step in many industrial chemical transformations, typically requiring precious metals. Here, we report a hierarchical nanoporous Cu catalyst doped with small amounts of Ti (npTiCu) that increases the rate of H 2 –D 2 exchange by approximately one order of magnitude compared to the undoped nanoporous Cu (npCu) catalyst. The promotional effect of Ti was measured via steady-state H 2 –D 2 exchange reaction experiments under atmospheric pressure flow conditions in the temperature range of 300–573 K. Pretreatment with flowing H 2 is required for stable catalytic performance, and two temperatures, 523 and 673 K, were investigated. The experimentally determined H 2 –D 2 exchange rate is 5–7 times greater for npTiCu vs the undoped Cu material under optimized pretreatment and reaction temperatures. The H 2 pretreatment leads to full reduction of Cu oxide and partial reduction of surface Ti oxide species present in the as-prepared catalyst as demonstrated using in situ ambient pressure X-ray photoelectron spectroscopy and X-ray absorption spectroscopy. The apparent activation energies and pre-exponential factors measured for H 2 –D 2 exchange are substantially different for Ti-doped vs undoped npCu catalysts. Density functional theory calculations suggest that isolated, metallic Ti atoms on the surface of the Cu host can act as the active surface sites for hydrogen recombination. Furthermore, the increase in the rate of exchange above that of pure Cu is caused primarily by a shift in the rate-determining step from dissociative adsorption on Cu to H/D atom recombination on Ti-doped Cu, with the corresponding decrease in activation entropy that it produces.

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

A combined theoretical and experimental study of small anthracene–water clusters

Water-cluster interactions with polycyclic aromatic hydrocarbons (PAHs) are of paramount interest in many chemical and biological processes. Here we report a study of anthracene monomers and dimers with water (up to four)-cluster systems utilizing molecular beam vacuum-UV photoionization mass spectrometry and density functional calculations. Structural loss in photoionization efficiency curves when adding water indicates that various isomers are generated, while theory indicates only a slight shift in energy in photoionization states of different isomers. Calculations reveal that the energetic tendency of water is to remain clustered and not to disperse around the PAH. Theoretically, we observe water confinement exclusively in the case of four water clusters and only when the anthracenes are in a cross configuration due to optimal OH∙∙∙π interactions, indicating dependence on the size and structure of the PAH. Furthermore theory sheds light on the structural changes that occur in water upon ionization of anthracene, due to the optimal interactions of the resulting hole and water hydrogen atoms.

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

A charge-density-based general cation insertion algorithm for generating new Li-ion cathode materials

Future lithium (Li) energy storage technologies, in particular solid-state configurations with a Li metal anode, opens up the possibility of using cathode materials that do not necessarily contain Li in its as-made state. To accelerate the discovery and design of such materials, we develop a general, chemically, and structurally agnostic methodology for identifying the optimal Li sites in any crystalline material. For a given crystal structure, we attempt multiple Li insertions at symmetrically in-equivalent positions by analyzing the electronic charge density obtained from first-principles density functional theory. In this report, we demonstrate the effectiveness of this procedure in successfully identifying the positions of the Li ion in well-known cathode materials using only the empty host (charged) material as guidance. Furthermore, applying the algorithm to over 2000 candidate cathode empty host materials we obtain statistics of Li site preferences to guide future developments of novel Li-ion cathode materials, particularly for solid-state applications.

25 ENERGY STORAGE↗