Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory simulations”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling bicarbonate formation in an alkaline solution with multi-level quantum mechanics/molecular dynamics simulations

Understanding carbonate speciation and how it may be modulated is essential for the advancement of carbon dioxide (CO 2 ) capture and storage technologies, which often rely on the transformation of CO 2 into carbonate, e.g. via the formation of carbonate minerals. To date, few atomic-level, quantum-mechanics-based simulations have been carried out to characterize how carbonic acid (H 2 CO 3 ) and bicarbonate ($HCO^{-}_{3}$) form in aqueous solution, and how pH affects this process. Recently, Martirez and Carter utilized rare-event sampling density functional theory molecular dynamics simulations in combination with multi-level embedded correlated wavefunction theory, thus accounting for both solvent dynamics and electron correlation accurately, to elucidate the mechanism of H 2 CO 3 formation in neutral solution (J. Am. Chem. Soc., 145, 12561, 2023). Here, we perform a complementary simulation using the same method to map out the energetics of $HCO^{-}_{3}$ formation from dissolved CO 2 in basic solution. We find that, as in H 2 CO 3 formation, including water dynamics is important to obtain an accurate prediction of the energetics for the aforementioned reaction. Furthermore, only with MD did we identify the correct pathway for the reaction, in which water – not hydroxide – acts as the initial nucleophile and only at the transition state does it lose a proton.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Best practices for first-principles simulations of epitaxial inorganic interfaces

Abstract At an interface between two materials physical properties and functionalities may be achieved, which would not exist in either material alone. Epitaxial inorganic interfaces are at the heart of semiconductor, spintronic, and quantum devices. First principles simulations based on density functional theory (DFT) can help elucidate the electronic and magnetic properties of interfaces and relate them to the structure and composition at the atomistic scale. Furthermore, DFT simulations can predict the structure and properties of candidate interfaces and guide experimental efforts in promising directions. However, DFT simulations of interfaces can be technically elaborate and computationally expensive. To help researchers embarking on such simulations, this review covers best practices for first principles simulations of epitaxial inorganic interfaces, including DFT methods, interface model construction, interface structure prediction, and analysis and visualization tools.

Physics↗
DOE OSTI2020

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗
DOE OSTI2021

Energetics of polymeric carbon monoxide

The transformation of carbon monoxide (CO) from a molecular liquid to a polymeric solid under isothermal compression at room temperature is investigated using first principles theory. We report structural and thermodynamic properties from ambient density up to 2.45 g/cc obtained using density functional theory molecular dynamics simulations, including hybrid exchange corrections. The theoretical results are compared with newly obtained polymeric CO samples, synthesized in a large volume press. The explosive performance of polymeric CO is predicted and discussed. Under most favorable assumptions, it is found to be comparable to trinitrotoluene.

25 ENERGY STORAGE↗
DOE OSTI2024

Effects of carbon concentration on the local atomic structure of amorphous GST

Ge-Sb-Te (GST) alloys are leading phase-change materials for data storage due to the fast phase transition between amorphous and crystalline states. Ongoing research aims at improving the stability of the amorphous phase to improve retention. This can be accomplished by the introduction of carbon as a dopant to Ge 2 Sb 2 Te 5 , which is known to alter the short- and mid-range structure of the amorphous phase and form covalently bonded C clusters, both of which hinder crystallization. The relative importance of these processes as a function of C concentration is not known. Here we used molecular dynamics simulation based on density functional theory to study how carbon doping affects the atomic structure of GST-C. Carbon doping results in an increase in tetrahedral coordination, especially of Ge atoms, and this is known to stabilize the amorphous phase. We observe an unexpected, non-monotonous trend in the number of tetrahedral bonded Ge with the amount of carbon doping. Our simulations show an increase in the number of tetrahedral bonded Ge up to 5 at.% C, after which the number saturates and begins to decrease above 14 at.% C. The carbon atoms aggregate into clusters, mostly in the form of chains and graphene flakes, leaving less carbon to disrupt the GST matrix at higher carbon concentrations. Different degrees of carbon clustering can explain divergent experimental results for recrystallization temperature for carbon doped GST.

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

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗
DOE OSTI2026

Insights into the Reactivity of Brookite TiO 2 Nanorods in Liquid Water from Ab Initio Molecular Dynamics Simulations

Brookite TiO 2 , a rare natural polymorph of TiO 2 , has been reported to be an excellent photocatalyst for the production of hydrogen from water and aqueous alcohol solutions, especially when it is reduced and synthesized in the form of nanorods. Here, we investigate the reactivity of stoichiometric and reduced brookite nanorods in liquid water using ab initio molecular dynamics and hybrid density functional theory calculations. Our simulations show a much higher water dissociation fraction on reduced nanorods than on stoichiometric ones, with an accumulation of the resulting bridging hydroxyls (O br H) and terminal hydroxyls (Ti –OH) on different facets of the nanorod. ObrH groups accumulate preferentially on low-energy (210) facets, where they are stabilized by adjacent reduced Ti (Ti 3+ ) sites, while Ti –OH groups prefer to form at the four-fold coordinated Ti atoms on high-energy (010) facets. This hydroxylation pattern also favors the spatial localization of excited holes on the (010) facets. This coupling between water-induced surface chemistry and charge separation underpins the enhanced photocatalytic activity of brookite nanorods, providing useful information for the design of more efficient TiO2-based nanostructures for solar-driven hydrogen evolution.

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

Atomistic simulations of the thermal conductivity of liquids

We present a method based on sinusoidal approach to equilibrium molecular dynamics (SAEMD) to compute the thermal conductivity of liquids Similar to nonequilibrium molecular dynamics, and unlike equilibrium simulations based on the Green-Kubo formalism, the method only requires the calculation of forces and total energies. The evaluation of heat fluxes and energy densities is not necessary, thus offering the promise of efficiently implementing first principles simulations based on density functional theory or deep molecular dynamics. Our approach is a generalization of SAEMD for solids, where the thermal conductivity is computed in the steady state, instead of a transient regime, thus properly taking into account diffusive terms in the heat equation. We present results for liquid water at ambient conditions and under pressure and discuss simulation requirements to obtain converged values of the thermal conductivity as a function of size and simulation time.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Kinetics of low-temperature methane activation on IrO 2 (1 1 0): Role of local surface hydroxide species

The ability of the IrO 2 (1 1 0) surface to promote CH 4 activation at low temperatures (~150 K) suggests possibilities for developing IrO 2 -based catalysts to selectively convert light alkanes to high-value chemicals. In this study, we present experimental results and microkinetic simulations, based on density functional theory (DFT) calculations, of temperature programmed reaction spectra (TPRS) obtained for CH 4 σ-complexes adsorbed on IrO 2 (1 1 0), and focus on clarifying how surface OH groups modify the branching between CH 4 desorption, activation and subsequent reactions during TPRS. Our DFT results predict that surface OH groups strongly destabilize CH 4 σ-complexes on IrO 2 (1 1 0) in addition to deactivating surface O-atoms that are needed to achieve CH 4 activation at low temperature. We demonstrate that a microkinetic model that incorporates the influence of surface OH groups on the CH 4 binding and reactivity reproduces experimental TPRS results which show that adsorbed CH 4 σ-complexes preferentially dissociate at low CH 4 coverage, but that an increasing fraction of the adsorbed CH 4 desorbs at low temperature (~120 K) with increasing initial CH 4 or H coverage. The simulations reveal that low-temperature CH4 desorption during TPRS arises primarily from CH 4 σ-complexes that are kinetically-trapped between adjacent surface OH groups, and are thus destabilized and unable to access reactive O-atoms. In addition, when the effect of adjacent OH groups is incorporated we find that the PBE-D3 functional incorporating dispersion provides CH 4 binding energies that agree more closely with that obtained from the TPRS experiments. Overall, our results demonstrate that the local effect of adjacent OH groups must be incorporated into any microkinetic models to properly capture the selectivity between extensive and partial oxidation in alkane conversion on IrO 2 (1 1 0) under reaction conditions.

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

Aggregation of Size-Selected Oxide Clusters Deposited onto Au(111)

Kinetic Monte Carlo (kMC) simulations along with density functional theory (DFT) calculations were used to investigate the aggregation of size-selected Nb 3 O y (y = 5, 6, 7) clusters deposited onto the Au(111) surface. Recent STM experiments showed that the cluster binding sites and sizes of the cluster assemblies on the Nb 3 O y /Au(111) surfaces strongly depend on the stoichiometry of the clusters, i.e., the oxygen-to-niobium ratio. To better understand the origins of these differences, kMC simulations of the nucleation and growth of cluster assemblies were performed using energy barriers for diffusion and intercluster interactions estimated from DFT calculations of cluster binding and dimerization energies, respectively. Comparisons of the kMC simulations with STM images of the as-deposited Nb 3 O y /Au(111) surfaces at RT and after high temperature annealing were used to further optimize the energetics and gauge the importance of nearest neighbor interactions. The kMC simulations demonstrate that the assembly of Nb 3 O y clusters on Au(111) are largely controlled by the magnitude of the barriers for diffusion and interparticle-bond formation, while changes at higher temperatures are sensitive to the binding energies between nearest neighbors. Simulations for the Nb 3 O 5 and Nb 3 O 6 clusters, which exhibit smaller cluster assembly sizes in STM, required larger diffusion barriers as well as different barriers for interparticle binding, which reflected differences in DFT calculated dimerization energies. Furthermore, the results demonstrate the effectiveness of combined DFT and kMC calculations for understanding how the stoichiometry affects the aggregation of small oxide clusters on a metal surface.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/$\textit{t}$ algorithm, in particular the short iteration performance. Further, a proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li 0.5 La 0.5 TiO 3 . This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li 0.5 La 0.5 TiO 3 are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li 0.5 La 0.5 TiO 3 system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

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

Development and Applications of an eReaxFF Force Field for Graphitic Anodes of Lithium-Ion Batteries

Graphene is one of the most promising materials for lithium-ion battery anodes due to its superior electronic conductivity, high surface area for lithium intercalation, fast ionic diffusivity and enhanced specific capacity. A reliable description of many battery processes requires an explicit description of electrochemical interactions involving electrons. A detailed atomistic modeling of electronic conduction and non-zero voltage simulations of graphitic materials require the inclusion of an explicit electronic degree of freedom. To enable large length- and time-scale simulations of electron conduction in graphitic anodes, we developed an eReaxFF force field concept describing graphitic materials with an explicit electron. The newly developed force field, verified against quantum chemistry-based data describing, amongst others, electron affinities and equation of states, reproduces the qualitative behavior of electron conductivity in pristine and imperfect graphitic materials at different applied temperatures and voltages. In addition, excess electron localization near a defect site estimated from eReaxFF simulations agree quite well with the corresponding density functional theory calculations. Here, our eReaxFF simulations show the initiation of lithium-metal-plating driven by electron transfer from the graphene surface to the exposed lithium ions demonstrating the method’s potential for studying lithium-graphene interactions with explicit electrons and explain many unresolved electrode and electrode-electrolyte interface processes.

25 ENERGY STORAGE↗
DOE OSTI2023

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗
DOE OSTI2020

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗
DOE OSTI2026

Thermally Stable Co@C3N4 Single-Atom Catalysts for CO Oxidation: Atomic-Level Insights into Structure and Activity

Single Co atoms supported on C3N4 (Co@C3N4) have demonstrated high activity and selectivity in photocatalysis. However, the investigation of structure–function relationships and reaction mechanisms under photocatalytic conditions is very challenging due to the complex conditions of light absorption, charge transfer, and catalysis. In this study, we employed thermal CO oxidation as a prototypical probe reaction to benchmark the intrinsic catalytic performance and track the active-site evolution of Co@C3N4. Single Co atoms were identified and shown to be the catalytically active sites for CO oxidation based on control experiments and isotope-labeling experiments. The Co sites remained atomically dispersed before, during, and after the reaction with temperatures up to 400 °C, as established by in situ X-ray absorption fine structure (XAFS) combined with density functional theory (DFT), FDMNES simulations, and dynamic-time-warping (DTW)-assisted X-ray absorption near edge structure (XANES) matching. Together with theoretical calculations, the integrated analysis reveals a stable coordination environment under reaction conditions, which correlates with sustained activity, establishing Co@C3N4 single-atom catalysts as thermally stable CO oxidation catalysts. Beyond these findings, the current study provides a workflow for unambiguously assigning active sites in Co@C3N4 for thermal CO oxidation. This workflow will aid the understanding of their behavior in photocatalysis in the future, where light-driven dynamics obscure direct structure–function links. Notably, this study provides fundamental insights for the rational design of robust single-atom catalysts and a foundation for the broader application of Co@C3N4 catalysts in oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗