Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Comparing the Expense and Accuracy of Methods to Simulate Atomic Vibrations in Rubrene

Atomic vibrations can inform about materials properties from hole transport in organic semiconductors to correlated disorder in metal–organic frameworks. Currently, there are several methods for predicting these vibrations using simulations, but the accuracy–efficiency tradeoffs have not been examined in depth. In this work, rubrene is used as a model system to predict atomic vibrational properties using six different simulation methods: density functional theory, density functional tight binding, density functional tight binding with a Chebyshev polynomial-based correction, a trained machine learning model, a pretrained machine learning model called ANI-1, and a classical forcefield model. The accuracy of each method is evaluated by comparison to the experimental inelastic neutron scattering spectrum. All methods discussed here show some accuracy across a wide energy region, though the Chebyshev-corrected tight-binding method showed the optimal combination of high accuracy with low expense. We then offer broad simulation guidelines to yield efficient, accurate results for inelastic neutron scattering spectrum prediction.

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

Building materials genome from ground‐state configuration to engineering advance

Individual phases are commonly considered as the building blocks of materials. However, the accurate theoretical prediction of properties of individual phases remains elusive. The top-down approach by decoding genomic building blocks of individual phases from experimental observations is nonunique. The density functional theory (DFT), as a state-of-the-art solution of quantum mechanics, prescribes the existence of a ground-state configuration at 0 K for a given system. It is self-evident that the ground-state configuration alone is insufficient to describe a phase at finite temperatures as symmetry-breaking non-ground-state configurations are excited statistically at temperatures above 0 K. Our multiscale entropy approach (recently terms as Zentropy theory) postulates that the entropy of a phase is composed of the sum of the entropy of each configuration weighted by its probability plus the configurational entropy among all configurations. Consequently, the partition function of each configuration in statistical mechanics needs to be evaluated by its free energy rather than total energy. The combination of the ground-state and symmetry-breaking non-ground-state configurations represents the building blocks of materials and can be used to quantitatively predict free energy of individual phases with the free energy of each configuration predicted from DFT as well as all properties derived from free energy of individual phases.

CALPHAD↗
DOE OSTI2024

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

Initial Fermi orbital descriptors for FLOSIC calculations: The quick-FOD method

Fermi orbital descriptors (FODs) play a key role in Fermi-Löwdin orbital self-interaction correction (FLOSIC) calculations used to remove self-interaction from approximate density functionals. Optimal FODs are obtained by minimizing the self-interaction-corrected total energy, and, in this process, identifying initial sets of FODs becomes crucial. Here we propose, implement and test a novel method for automatically initializing FODs, quick-FOD, based on the minimization of an empirical energy expression that involves a Coulomb-like FOD-electron density attraction, an FOD-FOD short-range repulsion, and an exchange-like FOD-density repulsion. Quick-FOD successfully reproduces FOD arrangements in qualitatively good agreement with Lewis theory and with full-fledged FLOSIC calculations.

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

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2024

A reactive molecular dynamics model for uranium/hydrogen containing systems

Uranium-based materials are valuable assets in the energy, medical, and military industries. However, understanding their sensitivity to hydrogen embrittlement is particularly challenging due to the toxicity of uranium and the computationally expensive nature of quantum-based methods generally required to study such processes. In this regard, we have developed a Chebyshev Interaction Model for Efficient Simulation (ChIMES) that can be employed to compute energies and forces of U and UH3 bulk structures with vacancies and hydrogen interstitials with accuracy similar to that of Density Functional Theory (DFT) while yielding linear scaling and orders of magnitude improvement in computational efficiency. Here, we show that the bulk structural parameters, uranium and hydrogen vacancy formation energies, and diffusion barriers predicted by the ChIMES potential are in strong agreement with the reference DFT data. We then use ChIMES to conduct molecular dynamics simulations of the temperature-dependent diffusion of a hydrogen interstitial and determine the corresponding diffusion activation energy. Our model has particular significance in studies of actinides and other high-Z materials, where there is a strong need for computationally efficient methods to bridge length and time scales between experiments and quantum theory.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electron Binding Energy Spectra of Al n Pt – Clusters—A Combined Experimental and Computational Study

Results of size-selected electron photo-detachment experiments and density functional theory calculations on anionic Al n Pt – , n = 1–7, clusters are presented and analyzed. The measured and calculated spectra of electron binding energies are, overall, in excellent accord with each other. The analysis reveals the general importance of accounting for the multiplicity of structural forms of a given-size cluster that can contribute to its measured spectrum, especially when the clusters are fluxional and/or the conditions of the experiment allow for structural transitions. Here, we show that for the systems studied here, the size-specific peculiarities of the measured spectra can be understood in terms of the combined contributions of corresponding different accessible stable equilibrium conformations, bona-fide transition-state configurations, and electronic-crossing structures that may play the role of effective barriers in electronically nonadiabatic dynamics.

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

Understanding the Existence of a Na 2 Dimer in a High-Spin State

The recent observation of a high-spin Na 2 dimer formed on the surface of liquid helium nanodroplets raises some fundamental questions, as the ground state of Na 2 is known to have zero spin. Is it protected against spontaneous dissociation? What is its binding energy and interatomic distance? Is it stable at a higher temperature? Using calculations based on density functional theory (with and without long-range interaction) and coupled cluster methods, CCSD(T), we show that the bonding in the high-spin Na 2 dimer is governed by van der Waals interaction with binding energy (bond length) varying between −0.030 eV (5.108 Å) and −0.192 eV (4.231 Å), depending on the computational method used. Thus, the experimental method used by Kresin and coworkers can be very useful to study larger metastable high-spin clusters such as Li 4 which was predicted in 1985 to have a tetrahedral structure carrying a magnetic moment of 2 μ B , while its ground state is planar and nonmagnetic.

Basis sets↗
DOE OSTI2024

Nonprecious Single Atom Catalyst for Methane Pyrolysis

The development of a suitable catalytic system for methane pyrolysis reactions requires a detailed investigation of the activation energy of C-H bonds on catalysts, as well as their stability against sintering and coke formation. In this work, both single-metal Ni atoms and small clusters of Ni atoms deposited on titanium nitride (TiN) plasmonic nanoparticles were characterized for the C-H bond activation of a methane pyrolysis reaction using ab initio spin-polarized density functional theory (DFT) calculations. The present work shows the complete reaction pathway, including energy barriers for C-H bond activation and dehydrogenated fragments, during the methane pyrolysis reaction on catalytic systems. Interestingly, the C-H bond activation barriers were low for both Ni single-atom and Ni-clusters, showing the energy barriers of ~1.10 eV and ~0.88 eV, respectively. Additionally, single-atom Ni-TiN showed weaker binding to adsorbates, and a net endothermic reaction pathway indicated that the single-atom Ni-TiN was expected to resist coke formation on its surface. However, these Ni single-atom catalysts can sinter, aggregate into a small cluster, and form a coke layer from the highly exothermic reaction pathway that the cluster takes despite the facile reaction pathway.

08 HYDROGEN↗
DOE OSTI2024

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Finally, other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

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

Defect levels and self-compensation in iodine-doped CdTe single crystals

This study systematically documents defect levels in n-type iodine-doped cadmium telluride (CdTe:I) crystals as measured by thermoelectric effect spectroscopy, Hall-effect measurements, and photoluminescence and calculated by density functional theory. The primary donor, I Te , was identified with an activation energy of ∼0.05 eV (measured)/∼0.13 eV (calculated). Deep acceptor states, V Cd and I Te –V Cd complex, exhibited activation energies of ∼0.1 eV (measured)/∼0.19 eV (calculated) and ∼0.12 eV (measured)/∼0.24 eV (calculated), respectively. Self-compensation and Fermi level pinning were observed at ∼0.7 eV (measured)/∼1.08 eV (calculated), contributing to high resistivity in as-grown samples. Post-growth Cd annealing effectively removes compensation centers, releases iodine donors, and significantly reduces resistivity, resulting in an ideal n-type solar cell material.

Cd annealing↗
DOE OSTI2024

Structure, Bonding, and Vibrational Dynamics of a Triamine High Energy Density Material under Pressure

High energy density materials have complex intermolecular interactions which influence their stability and performance. We used a combination of synchrotron X-ray diffraction, synchrotron infrared spectroscopy, and Raman vibrational spectroscopy, supplemented by density functional theory calculations, to probe pressure-induced changes in structure and intermolecular interactions of 1H,4'H-[3,3'-bis(1,2,4-triazole)]-4',5,5'-triamine as a model high energy density material up to 40 GPa. We find that compression of the triamine is accompanied by increased intermolecular interactions that give rise to an interesting evolution of the structure and bonding with pressure. Analysis of the equation of state determined from the X-ray diffraction indicates a change in compression mechanism near 19 GPa consistent with changes in vibrational spectra that provide evidence for a structural rearrangement associated with changes in hydrogen bonding near that pressure. As a result, the overall compressional behavior calculated theoretically agrees with that observed experimentally though differences are found that indicate the need for improved treatment of the intermolecular interactions including hydrogen bonding under pressure.

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

Projector‐Based Quantum Embedding Study of Iron Complexes

Projection‐based embedding theory (PBET) is used to calculate and assess the challenging spin‐crossover energies for a selection of small Fe‐containing systems by embedding the metal center into the frozen potential of the ligands. MP2, CCSD, and CCSD(T) are embedded in potentials from the SCAN and r 2 SCAN functionals and compared with the canonical values for the constituent methods and previously reported reference values. Considering the PBET calculations as a correction for the underlying DFT, the embedding calculations are able to provided improvement for most cases. In some cases, the PBET methods are able to compensate for limitations in the wave function methods and produce results similar to more rigorous calculations from the literature. For the systems with spin‐crossover energies near zero, the current methodology fails to provide consistent improvement. In conclusion, the isolated recalculation of the electronic structure around the metal center when embedded into a DFT treatment of the ligand field shows promise as a pragmatic and lower cost treatment compared to the canonical treatment of the whole system of the difficult class of spin‐crossover complexes.

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

Characterizing the surface compositions of supported bimetallic PtSn clusters: Effects of cluster-support interactions and surface adsorbates

PtSn bimetallic clusters on TiO2(110) and highly oriented pyrolytic graphite (HOPG) surfaces have been characterized by scanning tunneling microscopy, low energy ion scattering (LEIS), Xray photoelectron spectroscopy, and temperature programmed desorption (TPD); density functional theory (DFT) calculations have also been performed to better understand adsorption of CO and D2 on the PtSn surfaces. On TiO2 at coverages of 2 ML of Pt and 2 ML of Sn, exclusively bimetallic clusters are formed for both orders of deposition because clusters of the first metal completely cover the surface such that all atoms of the second metal are incorporated into the existing clusters. In contrast, on HOPG, the high mobility and weak cluster-support interactions on HOPG result in much larger 2 ML monometallic clusters (~30 Å high) that do not completely cover the surface, and deposition of the second metal produces larger clusters as well as smaller ones. Despite the difference in cluster morphologies for the different orders of deposition and supports, the LEIS experiments demonstrate that in all cases, the PtSn clusters are rich in Sn at the surface, as expected based on the lower surface free energy for Sn compared to Pt. Furthermore, the +0.2 eV shift in the Sn(3d5/2) binding energy observed on all surfaces in the presence of Pt is consistent with PtSn alloy formation. Deposition of 2 ML of Sn on TiO2 produces two-dimensional clusters with oxidation of Sn and reduction of titania at the clustersupport interface, but addition of Pt to the Sn clusters causes Sn to diffuse away from this interface, leaving Sn in the metallic state. TPD experiments on 2 ML Pt/TiO2 with increasing coverages of Sn show that the number of adsorption sites for D2 sharply decreases to nearly zero at 0.5 ML, while CO adsorption decreases to zero only at much higher Sn coverages of 2 ML. DFT studies for Sn modified Pt surfaces and bulk structures demonstrate that for CO adsorption at low Sn coverages (<0.25 ML), the strong Pt-CO interactions induce diffusion of Pt to the cluster surface and the formation of a bulk Pt3Sn alloy, whereas D2 adsorption does not lead to interactions with the Pt surface that are strong enough to induce alloy formation. A single Sn adatom prevents D2 adsorption on four neighboring Pt atoms via site-blocking and the donation of electron density to Pt.

Li, Fangliang↗
DOE OSTI2021

Supramolecular Self-Assembled Multi-Electron-Acceptor Organic Molecule as High-Performance Cathode Material for Li-Ion Batteries

Organic electrode materials possess many advantages such as low toxicity, sustainability, and chemical/structural tunability toward high energy density. However, to compete with inorganic-based compounds, crucial aspects such as redox potential, capacity, cycling stability, and electronic conductivity need to be improved. Herein, a comprehensive strategy on the molecular design of small organic electron-acceptor-molecule—hexaazatrianthranylene (HATA) embedded quinone (HATAQ) is reported. By introducing conjugated quinone moieties into the electron-deficient hexaazatriphenylene-derivative core, HATAQ with highly extended π-conjugation can yield extra-high capacity for lithium storage, delivering a capacity of 426 mAh g -1 at 200 mA g -1 (0.4C). At an extremely high rate of 10 A g -1 (19C), a reversible capacity of 209 mAh g -1 corresponding to nearly 85% retention is obtained after 1000 cycles. A unique network of unconventional lock-and-key hydrogen bonds in the solid-state facilitates favorable supramolecular 2D layered arrangement, enhancing cycling stability. To the best of the authors’ knowledge, the capacity and rate capability of HATAQ are found to be the best ever reported for organic small-molecule-based cathodes. Finally, these results together with density functional theory studies provide proof-of-concept that the design strategy is promising for the development of organic electrodes with exceptionally high energy density, rate capability, and cycling stability.

25 ENERGY STORAGE↗
DOE OSTI2024

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Benchmark calculations of pure neutron matter with realistic nucleon-nucleon interactions

Here, we report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner–Bethe–Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne v 18 interaction and two of its simplified versions, and four of the new Norfolk Δ -full chiral effective field theory potentials. Primarily because of the advancements in the auxiliary-field diffusion Monte Carlo, we observe good agreement among the three many-body techniques up to nuclear saturation density—the maximum difference in the energy per particle is within 1.5 MeV for all the potentials we consider. At higher densities, the divergences become more important, and are mainly driven by the Fermi hypernetted chain/single-operator calculations. We also study the connection between nucleon-nucleon scattering data and the energy per particle of pure neutron matter. Our results suggest that fitting to higher-energy nucleon-nucleon scattering helps reduce the spread of energies among the models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Development of interatomic potential and effect of ordering on defect properties in CrMnV

Developing materials that can withstand extreme environments, such as high radiation doses and elevated temperatures, is crucial for next-generation particle accelerators, including the 2.4 MW Long-Baseline Neutrino Facility. High-Entropy Alloys have emerged as promising candidates for beam window materials due to their superior mechanical strength, corrosion resistance, and radiation tolerance. In this study, we focus on the Cr–Mn–V alloy system, developing and employing machine-learning interatomic potentials (MLIPs) to investigate the formation of an ordered phase and its influence on defect properties. Using hybrid Monte Carlo-Molecular Dynamics simulations, we observe the formation of a B2-ordered phase at lower temperatures, consistent with Density Functional Theory (DFT) predictions. Ordered structures display a bimodal distribution of migration energies and reduced mean square displacement values, indicating suppressed vacancy diffusion. Our results also show that the migration energy barrier varies based on the atomic species, with Mn and V exhibiting the highest and lowest average barriers, respectively. These findings suggest that atomic ordering inhibits defect mobility, potentially enhancing the radiation resistance of CrMnV alloys. The validated MLIP provides a reliable framework for simulations that are faster than traditional DFT while maintaining the accuracy required to study defect and ordering properties.

36 MATERIALS SCIENCE↗