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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2022

High Thermoelectric Performance in Chalcopyrite Cu 1-x Ag x GaTe 2 –ZnTe: Nontrivial Band Structure and Dynamic Doping Effect

The understanding of thermoelectric properties of ternary I–III–VI 2 type (I = Cu, Ag; III = Ga, In; and VI = Te) chalcopyrites is less well developed. Although their thermal transport properties are relatively well studied, the relationship between the electronic band structure and charge transport properties of chalcopyrites has been rarely discussed. In this study, we reveal the unusual electronic band structure and the dynamic doping effect that could underpin the promising thermoelectric properties of Cu 1–x Ag x GaTe 2 compounds. Density functional theory (DFT) calculations and electronic transport measurements suggest that the Cu 1–x Ag x GaTe 2 compounds possess an unusual non-parabolic band structure, which is important for obtaining a high Seebeck coefficient. Moreover, a mid-gap impurity level was also observed in Cu 1–x Ag x GaTe 2 , which leads to a strong temperature-dependent carrier concentration and is able to regulate the carrier density at the optimized value for a wide temperature region and thus is beneficial to obtaining the high power factor and high average ZT of Cu 1–x Ag x GaTe 2 compounds. We also demonstrate a great improvement in the thermoelectric performance of Cu 1–x Ag x GaTe 2 by introducing Cu vacancies and ZnTe alloying. The Cu vacancies are effective in increasing the hole density and the electrical conductivity, while ZnTe alloying reduces the thermal conductivity. As a result, a maximum ZT of 1.43 at 850 K and a record-high average ZT of 0.81 for the Cu 0.68 Ag 0.3 GaTe 2 –0.5%ZnTe compound are achieved.

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

Feasibility study of Mg storage in a bilayer silicene anode via application of an external electric field

With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K

25 ENERGY STORAGE↗
DOE OSTI2022

Revisiting Theoretical Tools and Approaches for the Valorization of Recalcitrant Lignocellulosic Biomass to Value-Added Chemicals

Biorefinery processes for converting lignocellulosic biomass to fuels and chemicals proceed via an integrated series of steps. Biomass is first pretreated and deconstructed using chemical catalysts and/or enzymes to liberate sugar monomers and lignin fragments. Deconstruction is followed by a conversion step in which engineered host organisms assimilate the released sugar monomers and lignin fragments, and produce value-added fuels and chemicals. Over the past couple of decades, a significant amount of work has been done to develop innovative biomass deconstruction and conversion processes that efficiently solubilize biomass, separate lignin from the biomass, maximize yields of bioavailable sugars and lignin fragments and convert the majority of these carbon sources into fuels, commodity chemicals, and materials. Herein, we advocate that advanced in silico approaches provide a theoretical framework for developing efficient processes for lignocellulosic biomass valorization and maximizing yields of sugars and lignin fragments during deconstruction and fuel and chemical titers during conversion. This manuscript surveys the latest developments in lignocellulosic biomass valorization with special attention given to highlighting computational approaches used in process optimization for lignocellulose pretreatment; enzyme engineering for enhanced saccharification and delignification; and prediction of the genome modification necessary for desired pathway fine-tuning to upgrade products from biomass deconstruction into value-added products. Physics-based modeling approaches such as density functional theory calculations and molecular dynamics simulations have been most impactful in studies aimed at exploring the molecular level details of solvent-biomass interactions, reaction mechanisms occurring in biomass-solvent systems, and the catalytic mechanisms and engineering of enzymes involved in biomass degradation. More recently, with ever increasing amounts of data from, for example, advanced mutli-omics experiments, machine learning approaches have begun to make important contributions in synthetic biology and optimization of metabolic pathways for production of biofuels and chemicals.

09 BIOMASS FUELS↗
DOE OSTI2021

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

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

A family of ionic supersalts with covalent-like directionality and unconventional multiferroicity

Ionic crystals composed of elemental ions such as NaCl are non-polar due to directionless ionic bonding interactions. Here, we show that these can develop polarity by changing their building blocks from elemental ions to superalkali and superhalogen cluster-ions, which mimic the chemistry of alkali and halogen atoms, respectively. Due to the non-spherical geometries of these cluster ions, corresponding supersalts form anisotropic polar structures with ionic bonding, yet covalent-like directionality, akin to sp 3 hybridized systems. Using density functional theory and extensive structure searches, we predict a series of stable ferroelectric/ferroelastic supersalts, PnH 4 MX 4 (Pn = N, P; M = B, Al, Fe; X = Cl, Br) composed of superalkali PnH 4 and superhalogen MX 4 ions. Unlike traditional ferroelectric/ferroelastic materials, the cluster-ion based supersalts possess ultra-low switching barrier and can endure large ion displacements and reversible strain. In particular, PH 4 FeBr 4 exhibits triferroic coupling of ferroelectricity, ferroelasticity, and antiferromagnetism with controllable spin directions via either ferroelastic or 90-degree ferroelectric switching.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pathways to fluorescence via restriction of intramolecular motion in substituted tetraphenylethylenes

The design of materials with enhanced luminescence properties is a fast-developing field due to the potential applicability of these materials as light-emitting diodes or for bioimaging. A transparent way to enhance the emission properties of interesting molecular candidates is blocking competing and unproductive non-radiative relaxation pathways by the restriction of intramolecular motions. Rationalized functionalization is an important possibility to achieve such restrictions. Using time-dependent density functional theory (TD-DFT) based on the ωB97XD functional and the semiempirical tight-binding method including long-range corrections (TD-LC-DFTB), this work investigates the effect of functionalization of the paradigmatic tetraphenylethylene (TPE) on achieving restricted access to conical intersections (RACI). Photodynamical surface hopping simulations have been performed on a larger set of compounds including TPE and ten functionalized TPE compounds. Functionalization has been achieved by means of electron-withdrawing groups, bulky groups which block the relaxation channels via steric hindrance and groups capable of forming strong hydrogen bonds, which restrict the motion via the formation of hydrogen bond channels. Most of the investigated functionalized TPE candidates show ultrafast deactivation to the ground state due to their still existing structural flexibility, but two examples, one containing –CN and –CF 3 groups and a second characterized by a network of hydrogen bonds, have been identified as interesting candidates for creating efficient luminescence properties in solution.

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

Density Functional Theory Study of the Adsorption and Dissociation of Copper(I) Acetamidinates on Ni(110): The Effect of the Substrate

In order to assess the role of a metal substrates on the thermal chemistry of adsorbed acetamidinate metalorganic compounds, we have studied the surface chemistry of copper(I)–N,N'-dimethylacetamidinate on Ni(110) using density functional theory and contrasted it with similar calculations we previously carried out on Cu(110). At low coverages, it was found that, in its most stable configuration, the molecular adsorption of copper(I)–N,N'-dimethylacetamidinate dimers occurs with the Cu atoms occupying surface hollow sites. In this work, the ligands reorient away from those metal centers, and the N atoms develop new direct bonds with surface Ni atoms. However, this configuration is not stable and decomposes by losing both ligands to the surface. In the final state, the two ligands bind via their N atoms to Ni sites one lattice space away from the sites where the Cu atoms remain, with their molecular planes perpendicular to that of the surface. This is in contrast with what happens in the case of adsorption on Cu(110), where the Cu atoms from the metalorganic complex still occupy hollow sites but where only one of the ligands breaks away and binds directly to the surface; the other remains on top of the two Cu ions. In terms of the energetics of adsorption and decomposition, the reactions on Ni(110) are much more exothermic than on Cu(110). Further analysis of the distribution of charge within the adsorbates shows a minor reduction of the Cu atoms of the dimer upon interaction with the surface; full reduction to metallic copper is complete only when both ligands have fully migrated to their new Ni surface sites.

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

Synergistic effects of mixing and strain in high entropy spinel oxides for oxygen evolution reaction

Developing stable and efficient electrocatalysts is vital for boosting oxygen evolution reaction (OER) rates in sustainable hydrogen production. High-entropy oxides (HEOs) consist of five or more metal cations, providing opportunities to tune their catalytic properties toward high OER efficiency. This work combines theoretical and experimental studies to scrutinize the OER activity and stability for spinel-type HEOs. Density functional theory confirms that randomly mixed metal sites show thermodynamic stability, with intermediate adsorption energies displaying wider distributions due to mixing-induced equatorial strain in active metal-oxygen bonds. The rapid sol-flame method is employed to synthesize HEO, comprising five 3d-transition metal cations, which exhibits superior OER activity and durability under alkaline conditions, outperforming lower-entropy oxides, even with partial surface oxidations. The study highlights that the enhanced activity of HEO is primarily attributed to the mixing of multiple elements, leading to strain effects near the active site, as well as surface composition and coverage.

08 HYDROGEN↗