SEARCH · Engineering Papers
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.
Correction for ‘Computational approaches for XANES, VtC-XES, and RIXS using linear-response time-dependent density functional theory based methods’ by Daniel R. Nascimento et al. , Phys. Chem. Chem. Phys. , 2022, 24 , 14680–14691, https://doi.org/10.1039/D2CP01132H.
The lack of particle-hole attraction and orbital relaxation within time-dependent density functional theory (TDDFT) lead to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n - 1-electron system as the reference, building orbital relaxation and a proper hole into the initial density. Here, our approach is an exact generalization of the static-exchange approximation that ameliorates the particle-hole interaction error associated with the adiabatic approximation and reduces errors in TDDFT XAS by orders of magnitude. With a statistical performance of just 0.5 eV root-mean-square error and the same computational scaling as TDDFT under the core-valence separation approximation, we anticipate that this approach will be of great utility in XAS calculations of large systems.
Fundamental knowledge on the energetics at the interface between a water layer and a metal catalyst is essential so as to understand the roles that water can play in the synthesis, activation and regeneration of noble metal-based catalysts. Here, we report the direct measurement of the enthalpy of water adsorption (Δh ads ) on activated carbon (C) and activated C-supported Ru nanoparticles, which are promising catalyst as applied to the hydrogenation/hydrodeoxygenation (HDO) of oxygenates (phenolics, aldehydes, etc.). Specifically, the near-zero coverage enthalpy of water adsorption on a C-supported Ru catalyst is -75.3 ± 0.4 kJ/(mol water), suggesting favorable water–metal binding. This is much more exothermic than that on C, which has an enthalpy of adsorption of -50.3 ± 1.3 kJ/(mol water). Despite the favorable initial binding, the magnitudes of enthalpies of water condensation on C and Ru-C indicate that overall, their surfaces are both hydrophobic. Moreover, the experimentally-measured near-zero coverage water adsorption enthalpy at the Ru sites is in very good agreement with our density functional theory based calculations. At low coverages, we obtain a water binding energy of -61.7 kJ/(mol water), which increases to -78.1 kJ/(mol water) at saturation. Complementary results are also obtained from a thermal analysis, which employed a thermogravimetric analysis–differential scanning calorimetry–mass spectrometry (TG-DSC-MS), a spectroscopic investigation using ex situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and a morphological evaluation with transmission electron microscopy (TEM). We point out that in carbon-supported metal catalysts, such as Ru-C, a strong hydration at near-zero coverage and relative weak water-surface interactions occurs upon saturation. Such heterogeneity is essential and crucial for catalytic hydrogenation/HDO reactions that involve balanced interactions among the water-rich reactant mixture and nonpolar organic products.
Here, using the density functional theory-coupled Monte Carlo approach, we explored the chemical short-range order (SRO) and element segregation in equimolar CrMnFeCoNi alloy. We found that state-of-the-art approximation of random element distribution is only applicable at > 1100 K close to the melting temperature, while the Cr-Cr repulsion driving the system stabilization and accompanying the formation of cubic Cr sublattice, and mild Ni-Ni attraction are the most prominent pair interactions at lower temperatures. Chemical potential and vacancy formation energy calculations indicate that Cr is most sensitive to the local chemical environment, making Cr atoms most stabilized when the preferred SRO is introduced. While the vacancy formation is predicted equally probable among five constituting elements in the random solid solution, Cr and Ni atoms show the lowest vacancy formation energies in the structure with SRO. Furthermore, distinct element segregation was predicted in the vicinity of planar defects, including symmetric tilt grain boundary and stacking fault, which we correlated to the site- and chemistry-dependent atomic volume and bond lengths. It suggests that the local mechanical strain and bond energy induce the SRO development and element segregation: Namely, the system takes advantage of segregation of Ni atoms having large atomic volume or Cr-Cr pairs having elongated bond lengths to fill in the excess volume at defects that relaxes the mechanical strain field and optimizes bond energy distribution. The correlation between the SRO and properties of CrMnFeCoNi alloy needs further investigations, which is expected to greatly help understand and control the properties of high-entropy alloys.
Understanding structure–property relationships in atomically precise metal nanoclusters is vital in finding selective and tunable catalysts. Here, in this study, density functional theory (DFT) was used to benchmark seven exchange correlation functionals at different basis sets for 17 atomically precise nanoclusters against experimentally determined geometries, band gaps, and optical gaps. The set contains both monometallic and bimetallic clusters that possess at least two types of 3d transition metals (specifically, Cu, Ni, Fe, or Co). The benchmark highlights that PBE0 is a good functional to use regardless of the basis set, and Minnesota functionals do well with respect to specific metals. Further, while long-range corrected functionals overestimate band and optical gaps, they model absorption features better than the other considered functionals. The study additionally looks at the photoinduced hydrogen evolution reaction (HER) and the CO 2 reduction mechanism on nanoclusters reported from the literature.
The delta phase of Pu is stabilized by Ga doping, but the mechanism of this stabilization remains an open question. Density functional theory calculations focused on how Ga doping affects the phonons sheds some light on the phonons' contribution to the stabilization. Here, the calculated phonon modes of Ga-doped delta phase Pu fall into two distinct types: localized, high frequency Ga-dominated phonon modes, and Pu-dominated modes at lower frequencies. Increasing the Ga concentration has an effect on the Pu-dominated phonon modes opposite to that of compression: higher-frequency modes soften, and lower-frequency modes stiffen. The latter provides an indication that the stabilization mechanism is not due to a thermodynamic contribution from the phonons. Furthermore, the stiffened phonon modes include candidate modes that describe possible pathways into low-temperature phases, suggesting that doping with Ga could impede such pathways.
We report Dynamical Mean Field Theory (DMFT) is a successful method to compute the electronic structure of strongly correlated materials, especially when it is combined with density functional theory (DFT). Here, we present an open-source computational package (and a library) combining DMFT with various DFT codes interfaced through the Wannier90 package. The correlated subspace is expanded as a linear combination of Wannier functions introduced in the DMFT approach as local orbitals. In particular, we provide a library mode for computing the DMFT density matrix. This library can be linked and then internally called from any DFT package, assuming that a set of localized orbitals can be generated in the correlated subspace. The existence of this library allows developers of other DFT codes to interface with our package and achieve the charge-self-consistency within DFT+DMFT loops. To test and check our implementation, we computed the density of states and the band structure of well-known solid-state correlated materials, namely LaNiO, SrVO, and NiO. The obtained results are compared to those obtained from other DFT+DMFT implementations.
Several mechanical and thermophysical properties are required as a function of non-stoichiometry for the successful implementation of YH 2-x for nuclear reactor moderator applications. Density functional theory calculations, in combination with neutron diffraction experiments, were used to study the structural and mechanical properties of YH 2-x . Point defect analysis indicated H occupation primarily at the tetrahedral site within an fcc Y sub-lattice, confirming a fluorite YH 2 structure. The small positive formation energy for H vacancies under Y-rich conditions predicted that hypo-stoichiometry is accommodated by Y+YH 2 at ambient conditions and by H vacancies in the YH 2-x single phase that is relevant to high temperatures. Neutron diffraction studies were used to confirm both the occupation of H on tetrahedral sites and the near-stoichiometric composition of the hydride phase in the two-phase Y+YH 2 region of the phase diagram that dominates at room temperature. Energy minimized special-quasirandom-structures of H vacancies were used to calculate lattice parameters, elastic constants, and several other properties as a function of composition for the single phase YH 2-x . The lattice parameter of YH 2-x decreased by only 0.004 Å with increasing H/Y for the range 1.31 ≤ H / Y ≤ 2.0 indicating a negligible effect on lattice parameters due to vacancy formation. In the two-phase region, however, calculations predicted the density of z Y + ( 1 - z ) Y H 2 - x to increase with decreasing H/Y at lower temperatures due to the increased fraction of high-density Y metal. Additionally, for the high temperature single phase, decreasing H/Y reduced the density as a consequence of the lattice expansion associated with vacancy formation. All elastic constants and moduli increased with increasing hydrogen content in single-phase YH 2-x .
Herein, we examine the effect of removing self-interaction error (SIE) on the calculation of molecular polarizabilities in the local spin density (LSDA) and generalized gradient approximations (GGA). To this end, we utilize a database of 132 molecules taken from a recent benchmark study [Hait and Head-Gordon, Phys. Chem. Chem. Phys. , 2018, 20 , 19800] to assess the influence of SIE on polarizabilities by comparing results with accurate reference data. Our results confirm that the general overestimation of molecular polarizabilities by these density functional approximations can be attributed to SIE. However, removing SIE using the Perdew–Zunger self-interaction-correction (PZ-SIC) method, implemented using the Fermi–Löwdin Orbital SIC approach, leads to an underestimation of molecular polarizabilities, showing that PZ-SIC overcorrects when combined with LSDA or GGA. Application of a recently proposed locally scaled SIC [Zope, et al., J. Chem. Phys. , 2019, 151 , 214108] is found to provide more accurate polarizabilities. We attribute this to the ability of the local scaling scheme to selectively correct for SIE in the regions of space where the correction is needed most.
Abstract Barium titanate (BTO) is a widely researched ferroelectric useful for energy storage. While BTO’s surface chemistry is commonly studied using density functional theory, little has been published on the TiO2 surface. Here, we determined that BTO’s surface response can be decoupled from the ferroelectric response by using a pre-optimized ferroelectric slab and allowing only the top three atomic z-layers to respond to ligand binding. Multiple favorable binding modes were identified for hydrogen, hydroxyl, water, and tert-butyl phosphonic acid on BTO’s TiO 2 surface. Of these ligands, tBuPA dominates surface binding with binding energies as low as − 2.61 eV for its nine configurations. Graphical abstract
We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (~1.2 × 10 -3 E h / a 0 ) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.
The funding from this award supported the development of new theoretical methods that both increase the accuracy of quantum simulations beyond that of density functional theory, as well as reduce the cost of such simulations, for example, through machine learning and algorithmic improvements. The improvements lead to, in some cases, orders of magnitude increases in simulation speed as well as orders of magnitude increase in accuracy, compared to before this work was started. This translates to both greater certainty in making predictions about chemical phenomena that are currently studied computationally, as well as laying the ground-work for the study of new areas of chemistry that have not traditionally been simulated, because of the new time-scales and accuracies that can be reached.
In order to provide a comprehensive theoretical description of MgSiO 3 at extreme conditions, we combine results from path integral Monte Carlo and density functional molecular dynamics simulations and generate a consistent equation of state for this material. We consider a wide range of temperature and density conditions from 10 4 to 10 8 K and from 0.321 to 64.2 g/cm -3 (0.1- to 20-fold the ambient density). We study how the L and K shell electrons are ionized with increasing temperature and pressure. We derive the shock Hugoniot curve and compare with experimental results. Our Hugoniot curve is in good agreement with the experiments, and we predict a broad compression maximum that is dominated by the K shell ionization of all three nuclei while the peak compression ratio of 4.70 is obtained when the Si and Mg nuclei are ionized. To conclude, we analyze the heat capacity and structural properties of the liquid.
Here, we present the theory, implementation, and benchmarking of a real-time time-dependent density functional theory (RT-TDDFT) module within the RMG code, designed to simulate the electronic response of molecular systems to external perturbations. Our method offers insights into nonequilibrium dynamics and excited states across a diverse range of systems, from small organic molecules to large metallic nanoparticles. Benchmarking results demonstrate excellent agreement with established TDDFT implementations and showcase the superior stability of our time integration algorithm, enabling long-term simulations with minimal energy drift. The scalability and efficiency of RMG on massively parallel architectures allow for simulations of complex systems, such as plasmonic nanoparticles with thousands of atoms. Future extensions, including nuclear and spin dynamics, will broaden the applicability of this RT-TDDFT implementation, providing a powerful toolset for studies of photoactive materials, nanoscale devices, and other systems where real-time electronic dynamics is essential.
Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.
The incorporation of uranium, plutonium and technetium in the negative thermal expansion (NTE) α-Zr (WO 4 ) 2 has been investigated within the framework of density functional theory (DFT). It is found that the vacancy formation energies of the charged vacancies are overall larger than that of its counterpart neutral Frenkel defects and Schottky defects. DFT calculations suggest that U and Pu substitutions for the Zr site are preferred in α-Zr (WO 4 ) 2 . In case of Tc substitution, both Tc(IV) for the Zr site and Tc(VII) for the W site are considered under oxygen-poor and oxygen-rich conditions, while Tc(VII) substitution can be improved significantly by including Y 2 O 3 (charge compensation).
Ni coarsening in the fuel electrode of solid oxide cells (SOCs) is known to be significantly faster under a humid atmosphere. In this talk, density-functional theory and kinetic Monte Carlo methods are used to explore the hypothesis that the surface diffusion of NiH on Ni may promote Ni coarsening in the SOC operated in electrolysis cell mode. Defining the surface diffusivity as the product of the fractional surface coverage and single-molecule surface diffusivity, the surface diffusivity of NiH on Ni (111) is found to be large enough under a significant overpotential to support the above hypothesis. However, as NiH could dissociate on Ni (111), more work is needed to show that NiH may promote Ni coarsening.
In this work, the initial kinetics of water dissociation on two facets of δ-plutonium, δ-Pu(111) and δ-Pu(100), are explored through density functional theory in order to understand how water dissociation occurs on these facets. We explored the dissociation of water via the formation of hydroxyls, atomic hydrogen, and atomic oxygen species on each facet. We calculate low energetic barriers for water to split to adsorbed hydrogen and hydroxyl species at 0.19 eV for δ-Pu(111) and 0.07 eV for δ-Pu(100). The hydroxyl has a barrier of 0.64 and 0.37 eV to cleave the hydrogen–oxygen bond on δ-Pu(111) and δ-Pu(100), respectively. Due to the highly exergonic adsorption free energy of atomic oxygen of −2.10 eV, the metallic surfaces are found to be fully covered in oxygen, even with the inclusion of oxygen lateral interactions. When combined with the reaction thermodynamics, this free energy of adsorption forms a molecular oxygen desorption barrier greater than 9.51 eV (918 kJ/mol). These results, combined with simulated temperature programed desorption spectra, indicate that oxygen formed via water dissociation on δ-Pu (111) and (100) facets induces an irreversible poison that prevents further reaction of water directly on metallic plutonium surfaces, most likely due to the strong hybridization of the Pu and O valence states. Therefore, these results imply that another mechanism is responsible for the continuously experimentally measured water dissociation, which is most likely due to the Pu oxide.