Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Measuring Intermolecular Excited State Geometry for Favorable Singlet Fission in Tetracene

Singlet fission (SF) is the process of converting an excited singlet to a pair of excited triplets. Harvesting two charges from a single photon has the potential to increase photovoltaic device efficiencies. Acenes, such as tetracene and pentacene, are model molecules for studying SF. Despite SF being an endoergic process for tetracene and exoergic for pentacene, both acenes exhibit near unity SF quantum efficiencies, raising questions about how tetracene can overcome the energy barrier. Here, we use recently developed instrumentation to measure inelastic neutron scattering (INS) while optically exciting the model molecules using two different excitation energies. The spectroscopic results reveal intermolecular structural relaxation due to the presence of a triplet excited state. Here, the structural dynamics of the combined excited state molecule and surrounding tetracene molecules are further studied using time-dependent density functional theory (TD-DFT), which shows that the singlet and triplet levels shift due to the excited state geometry, reducing the uphill energy barrier for SF to within kT.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogenation of ethylene over molybdenum–sulfur complexes supported on UiO-66

Development of supported single-site catalysts using small metal sulfide complexes could significantly help in the development of cost-effective catalytic materials to drive selective hydrogenation and hydrogenolysis. The goal of this study is to contribute to the development of metal sulfide catalysts by calculating the thermodynamics of a catalytic cyclic involving a metal organic framework functionalized by insertion of metal sulfide. Anchored metal sulfide complexes can potentially be designed with ligands with distinctly different electronic and catalytic properties for specific catalytic applications. Here we examine the hydrogenation of ethylene as a model. We use density functional theory to investigate molybdenum–sulfur complexes as active catalysts anchored on the metal–organic framework UiO-66 as a stable support. Our calculations show that the anchored complexes with more than two sulfur ligands are unfavorable for ethylene adsorption, so we study complexes with one or two sulfur ligands. Hydrogenation of the unsaturated carbon double bond requires the transfer of two hydrogen atoms, which can occur via heterolytic activation of hydrogen to form a Mo-hydride and a protonated sulfur – either by hydride transfer followed by proton transfer or via proton transfer followed by hydride transfer, and we find that both mechanisms proceed via two-state reactivity involving two spin states along the reaction path. Of the two catalysts studied in gas the phase, the MoS single-sulfur–ligand complex with lower oxidation states produces thermodynamically more favorable intermediates along the pathway for the first hydrogen transfer for both the hydride-first mechanism and the proton-first mechanism. As a result, the quantum mechanical calculations provide experimentally inaccessible partial atomic charges and geometries of the various intermediates encountered along the steps of the reaction mechanisms.

Kermani, Maryam Mansoori [University of Minnesota,↗
DOE OSTI2026

Neutron skins probed in proton knockout from neutron-rich nuclei

Proton-induced quasifree knockout reactions provide a powerful probe of nuclear single-particle structure and reaction dynamics in both stable and neutron-rich nuclei. Here, in this work, we develop a unified theoretical framework for the calculation of inclusive (𝑝, 2⁢𝑝) and sequential (𝑝, 3⁢𝑝) reaction cross sections and fragment momentum distributions at intermediate and relativistic energies. The approach is based on a probabilistic extension of Glauber multiple-scattering theory combined with microscopic nuclear densities obtained from Hartree-Fock-Bogoliubov calculations using Skyrme energy-density functionals. We focus in particular on the sensitivity of total cross sections and longitudinal momentum dispersions to neutron-skin thickness along isotopic chains. Our results indicate that both (𝑝, 2⁢𝑝) and (𝑝, 3⁢𝑝) reactions exhibit a systematic decrease of cross section and momentum width with increasing neutron excess, reflecting enhanced attenuation and surface bias induced by neutron skins. The effect is significantly stronger for two-proton removal, suggesting that (𝑝, 3⁢𝑝) reactions may offer enhanced sensitivity to isovector nuclear structure. These findings establish proton-induced knockout reactions as complementary hadronic probes of neutron skins and the density dependence of the nuclear symmetry energy.

direct reactions↗
DOE OSTI2023

DFT (and DMFT) for DFT users

The interesting properties of f -electron systems comes from a simultaneous existence of free and discrete level electrons. The challenge we have is to describe free electrons and discrete level electrons equally well in a unified picture. The ELF index can be used to find regions in real space where discrete level physics needs to be taken into account. We have identified the HO gas as a model system that can be used to gain insight about this kind of physics. We will use the HO gas model system for creating a functional suitable for these systems via the subsystem functional scheme. With a predictive capability for the rare earths and lanthindes, we can identify new materials with new/improved properties. This is crucial for, e.g., developing material and technology substitutes for addressing the challenges associated with the use of rare earth elements and other critical materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reining in Radium for Nuclear Medicine: Extra-Large Chelator Development for an Extra-Large Ion

Targeted α therapy (TAT) of soft-tissue cancers using the α particle-emitting radionuclide 223 Ra holds great potential because of its favorable nuclear properties, adequate availability, and established clinical use for treating metastatic prostate cancer of the bone. Despite these advantages, the use of 223 Ra has been largely overshadowed by other α emitters due to its challenging chelation chemistry. A key criterion that needs to be met for a radionuclide to be used in TAT is its stable attachment to a targeting vector via a bifunctional chelator. The low charge density of Ra 2+ arising from its large ionic radius weakens its electrostatic binding interactions with chelators, leading to insufficient complex stability in vivo. In this study, we synthesized and evaluated macropa-XL as a novel chelator for 223 Ra. It bears a large 21-crown-7 macrocyclic core and two picolinate pendent groups, which we hypothesized would effectively saturate the large coordination sphere of the Ra 2+ ion. The structural chemistry of macropa-XL was first established with the nonradioactive Ba 2+ ion using X-ray diffraction and X-ray absorption spectroscopy, which revealed the formation of an 11-coordinate complex in a rare anti pendent-arm configuration. Subsequently, the stability constant of the [Ra(macropa-XL)] complex was determined via competitive cation exchange with 223 Ra and 224 Ra radiotracers and compared with that of macropa, the current state-of-the-art chelator for Ra 2+ . A moderate log K ML value of 8.12 was measured for [Ra(macropa-XL)], which is approximately 1.5 log K units lower than the stability constant of [Ra(macropa)]. This relative decrease in Ra 2+ complex stability for macropa-XL versus macropa was further probed using density functional theory calculations. Additionally, macropa-XL was radiolabeled with 223 Ra, and the kinetic stability of the resulting complex was evaluated in human serum. Furthermore, although macropa-XL could effectively bind 223 Ra under mild conditions, the complex appeared to be unstable to transchelation. Collectively, this study sheds additional light on the chelation chemistry of the exotic Ra 2+ ion and contributes to the small, but growing, number of chelator development efforts for 223 Ra-based TAT.

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

Chemical manipulation of hydrogen induced high p-type and n-type conductivity in Ga 2 O 3

Advancement of optoelectronic and high-power devices is tied to the development of wide band gap materials with excellent transport properties. However, bipolar doping (n-type and p-type doping) and realizing high carrier density while maintaining good mobility have been big challenges in wide band gap materials. Here P-type and n-type conductivity was introduced in β-Ga 2 O 3 , an ultra-wide band gap oxide, by controlling hydrogen incorporation in the lattice without further doping. Hydrogen induced a 9-order of magnitude increase of n-type conductivity with donor ionization energy of 20 meV and resistivity of 10 -4 Ω.cm. The conductivity was switched to p-type with acceptor ionization energy of 42 meV by altering hydrogen incorporation in the lattice. Density functional theory calculations were used to examine hydrogen location in the Ga 2 O 3 lattice and identified a new donor type as the source of this remarkable n-type conductivity. Positron annihilation spectroscopy measurements confirm this finding and the interpretation of the experimental results. This work illustrates a new approach that allows a tunable and reversible way of modifying the conductivity of semiconductors and it is expected to have profound implications on semiconductor field. At the same time, it demonstrates for the first time p-type and remarkable n-type conductivity in Ga 2 O 3 which should usher in the development of Ga 2 O 3 devices and advance optoelectronics and high-power devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mean force kinetic theory applied to self-diffusion in supercritical Lennard-Jones fluids

A new kinetic-theory-based calculation of the self-diffusion coefficient for dense supercritical Lennard-Jones fluids is presented. The mean force kinetic theory, which was recently developed for transport in dense plasmas, is applied for the calculation of diffusion in dense neutral fluids. The calculation only requires the pair distribution function, a quantity that is readily calculable from equilibrium statistical mechanics for many systems, including the Lennard-Jones fluid. The self-diffusion coefficients are compared with calculations from molecular dynamics simulations, and good agreement at high density is demonstrated, even in the vicinity of the solid–fluid coexistence line. A comparison of different kinetic models with molecular dynamics simulations demonstrates that the transport coefficients have important contributions due to particle interaction via a potential of mean force and local correlations, which increase the collision rate. The new calculations compare well to those from free-volume theory and overcome a limitation of this theory that prevents its use in systems that interact via long range monotonic potentials. It is expected that this approach will also apply to other systems, including neutral-plasma and neutral-electrolyte mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

HfNb_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-niobium (Hf-Nb). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Developments in Physics-Based Modeling and Machine Learning for Environmental Effects in Nuclear Materials [Slides]

This presentation is part of an annual review for the new Advanced Materials and Manufacturing Technologies program under the Office of Nuclear Energy within the Department of Energy. This presentation covers advances in physics-based modeling and machine learning for environmental effects in nuclear materials. It discusses efforts in combining density functional theory with machine learning to predict material properties, studying irradiation defect recombination, and improving MOOSE-based tools for crystal plasticity modeling and the Stochastic Tools Module.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Accelerated discovery of a large family of quaternary chalcogenides with very low lattice thermal conductivity

The development of efficient thermal energy management devices such as thermoelectrics and barrier coatings often relies on compounds having low lattice thermal conductivity (κ l ). Here, we present the computational discovery of a large family of 628 thermodynamically stable quaternary chalcogenides, AMM'Q 3 (A = alkali/alkaline earth/post-transition metals; M/M' = transition metals, lanthanides; Q = chalcogens) using high-throughput density functional theory (DFT) calculations. We validate the presence of low κ l in these materials by calculating κ l of several predicted stable compounds using the Peierls–Boltzmann transport equation. Our analysis reveals that the low κ l originates from the presence of either a strong lattice anharmonicity that enhances the phonon-scatterings or rattler cations that lead to multiple scattering channels in their crystal structures. Our thermoelectric calculations indicate that some of the predicted semiconductors may possess high energy conversion efficiency with their figure-of-merits exceeding 1 near 600 K. Our predictions suggest experimental research opportunities in the synthesis and characterization of these stable, low κ l compounds.

36 MATERIALS SCIENCE↗
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

Charge Density Wave Order and Electronic Phase Transitions in a Dilute d ‐Band Semiconductor

Abstract As one of the most fundamental physical phenomena, charge density wave (CDW) order predominantly occurs in metallic systems such as quasi‐1D metals, doped cuprates, and transition metal dichalcogenides, where it is well understood in terms of Fermi surface nesting and electron–phonon coupling mechanisms. On the other hand, CDW phenomena in semiconducting systems, particularly at the low carrier concentration limit, are less common and feature intricate characteristics, which often necessitate the exploration of novel mechanisms, such as electron–hole coupling or Mott physics, to explain. In this study, an approach combining electrical transport, synchrotron X‐ray diffraction, and density‐functional theory calculations is used to investigate CDW order and a series of hysteretic phase transitions in a dilute d ‐band semiconductor, BaTiS 3 . These experimental and theoretical findings suggest that the observed CDW order and phase transitions in BaTiS 3 may be attributed to both electron–phonon coupling and non‐negligible electron–electron interactions in the system. This work highlights BaTiS 3 as a unique platform to explore CDW physics and novel electronic phases in the dilute filling limit and opens new opportunities for developing novel electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
DOE OSTI2020

Computational study of TATB

We have performed a computational study of TATB. The study is composed of two parts, one where we perform static T = 0 K calculations with several different DFT functionals, to investigate structural properties and obtain a cold curve. Even though the functionals used in this study give poor results at ambient pressure, they perform much better at compression, when the dominance of the van der Waals’ forces in the binding is replaced by more normal interactions. In particular AM05 gives good results, mainly because it doesn’t include any “faulty” van der Waals’. Indeed the van der Waals’ corrected functional (PBE+D2) results shows the worst performance on structural properties. For the cold curve the PBE+D2 functional gives the best results compared to experiments (room temperature) but since the structural properties were bad with this functional we cannot for sure say that this is in any way a validation. It is shown that scaling the AM05 results with a common factor reproduces the experimental data. We leave the question on why to a follow up study but speculate that maybe the experiment and our calculations do not represent the same system set-up. From the structural properties we identify a site in the TATB lattice that contains a high concentration of oxygen and hydrogen. The consequences of this finding, if any, are left to a future project. In the second part of the study we perform extensive high-quality DFT-MD (molecular dynamics with DFT forces) calculations to obtain points on the Hugoniot. This part of the study did not get finished within this project and a more extensive summary of this part will be given in a follow up report at a future date. However, we show that obtaining calculated data is essential for equations of state development since the existing experimental Hugoniot data barely deviates from room temperature isotherm data and thus give no information of the high temperature dependency of materials properties. Our conclusion so far in the project is that being able to calculate materials properties is essential for high explosives since experiments are unable to give information in the high temperature part of phase space. While blind use of DFT cannot (yet) give unambiguous answers, we can still explore qualitatively and use such information to guide us in EOS development. However, the main conclusion is that we need better functionals in DFT to easily extract needed properties from calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Axially-deformed solution of the Skyrme-Hartree-Fock-Bogoliubov equations using the transformed harmonic oscillator basis (IV) HFBTHO (v4.0): A new version of the program

We describe the new version 4.0 of the code HFBTHO that solves the nuclear Hartree-Fock-Bogoliubov problem by using the deformed harmonic oscillator basis in cylindrical coordinates. In the new version, we have implemented the restoration of rotational, particle number, and reflection symmetry for even-even nuclei. The restoration of rotational symmetry does not require using bases closed under rotation. Furthermore, we added the SeaLL1 functional and improved the calculation of the Coulomb potential. Finally, we refactored the code to facilitate maintenance and future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Zero-Field NMR and Millitesla-SLIC Spectra for >200 Molecules from Density Functional Theory and Spin Dynamics

NMR is usually performed at magnetic fields of 1 T and above to obtain sufficient sensitivity and spectral dispersion to identify chemicals based on chemical shifts and J couplings. At lower fields, the advent of hyperpolarization technologies and sensitive detectors can address sensitivity concerns. However, it remains disputed whether spectral signatures at zero and ultra-low fields are sufficient for chemical identification. Here, we report an all–electron DFT-based batch calculation of J-coupling constants, which are used to generate J coupling NMR spectra at zero field and 6.5 mT for over 200 small molecules. In the developed computational tool chain, we first used the all-electron FHI-aims code to calculate the molecular J couplings and chemical shifts. We then fed the calculated NMR parameters into the NMR simulation package SPINACH to simulate both heteronuclear J coupling spectra at zero-field, and homonuclear J coupling spectra as spin-lock induced crossing (SLIC) spectra at ultra-low field (6.5 mT). The resulting spectra demonstrate that zero and ultra-low field NMR spectra can represent unique identifiers of chemical structure for small molecules.

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

Effect of H 2 O on the ethylene glycol/alkali dismantling of bagasse for high-value conversion

Using a high-boiling alcohol system to dismantle main components of biomass is a feasible technology. Reducing dismantle operating costs and improving dismantle efficiency are essential for promoting the green, economical, and sustainable development of biomass refining. Therefore, based on the low cost and chemical properties of H 2 O at high temperature, the effects of different H 2 O dosages in NaOH-catalyzed ethylene glycol (HBAA) system on the dismantling efficiency of bagasse, surface lignin coverage, recovered-lignin activity and enzymatic hydrolysis efficiency were investigated. Compared with the HBAA dismantling system without H 2 O, the HBAA system with 60% w/v H 2 O can obviously increase the removal rates of lignin and hemicellulose, while recovering up to 99% of cellulose and significantly declining surface lignin coverage, thus enhancing the enzymatic hydrolysis efficiency. Additionally, the results of density functional theory calculations and 2D HSQC NMR analysis prove that the synergy between H 2 O and ethylene glycol can promote the esterification reaction occurrence at the α-C carbon cation in β-O-4 structure of lignin, thereby protecting the β-O-4 aromatic ether bond. Simultaneously, when the H 2 O dosage increase from 0% to 60%, the enzymatic yield increases from 84.51% to 93.74% with an enzyme load of 10 FPU/g. Based on experimental results, this study conducted a techno-economic analysis of bagasse dismantling for ethanol and co-production of lignin, achieving a minimum ethanol selling price of $\$$1.07 per kg. Here, in this study, a green and economical solution for dismantling the main components of bagasse is developed, which is important for the high-value conversion of bagasse.

Ethylene glycol↗
DOE OSTI2021

Discovery of Energy Storage Molecular Materials Using Quantum Chemistry-Guided Multiobjective Bayesian Optimization

Redox flow batteries (RFBs) are a promising technology for stationary energy storage applications due to their flexible design, scalability, and low cost. In RFBs, energy is carried in flowable redoxactive materials (redoxmers) which are stored externally and pumped to the cell during operation. Further improvements in the energy density of RFBs necessitates redoxmer designs with wider redox potential windows and higher solubility. Additionally, designing redoxmers with a fluorescence-enabled self-reporting functionality allows monitoring of the state of health of RFBs. To accelerate the discovery of redoxmers with desired properties, state-of-the- art machine learning ( ML) methods, such as multiobjective Bayesian optimization (MBO), are useful. Here, we first employed density functional theory calculations to generate a database of reduction potentials, solvation free energies, and absorption wavelengths for 1400 redoxmer molecules based on a 2,1,3-benzothiadiazole (BzNSN) core structure. From the computed properties, we identified 22 Pareto-optimal molecules that represent best trade-off among all of the desired properties. We further utilized these data to develop and benchmark an MBO approach to identify candidates quickly and efficiently with multiple targeted properties. With MBO, optimal candidates from the 1400-molecule data set can be identified at least 15 times more efficiently compared to the brute force or random selection approach. Importantly, we utilized this approach for discovering promising redoxmers from an unseen database of 1 million BzNSN-based molecules, where we discovered 16 new Pareto-optimal molecules with significant improvements in properties over the initial 1400 molecules. Furthermore, we anticipate that this active learning technique is general and can be utilized for the discovery of any class of functional materials that satisfies multiple desired property criteria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗