Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional approximations”

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

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

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

Stabilizing lattice oxygen redox in layered sodium transition metal oxide through spin singlet state

Abstract Reversible lattice oxygen redox reactions offer the potential to enhance energy density and lower battery cathode costs. However, their widespread adoption faces obstacles like substantial voltage hysteresis and poor stability. The current research addresses these challenges by achieving a non-hysteresis, long-term stable oxygen redox reaction in the P3-type Na 2/3 Cu 1/3 Mn 2/3 O 2 . Here we show this is accomplished by forming spin singlet states during charge and discharge. Detailed analysis, including in-situ X-ray diffraction, shows highly reversible structural changes during cycling. In addition, local CuO 6 Jahn-Teller distortions persist throughout, with dynamic Cu-O bond length variations. In-situ hard X-ray absorption and ex-situ soft X-ray absorption study, along with density function theory calculations, reveal two distinct charge compensation mechanisms at approximately 3.66 V and 3.99 V plateaus. Notably, we observe a Zhang-Rice-like singlet state during 3.99 V charging, offering an alternative charge compensation mechanism to stabilize the active oxygen redox reaction.

25 ENERGY STORAGE↗
DOE OSTI2021

Computing molecular excited states on a D-Wave quantum annealer

Abstract The possibility of using quantum computers for electronic structure calculations has opened up a promising avenue for computational chemistry. Towards this direction, numerous algorithmic advances have been made in the last five years. The potential of quantum annealers, which are the prototypes of adiabatic quantum computers, is yet to be fully explored. In this work, we demonstrate the use of a D-Wave quantum annealer for the calculation of excited electronic states of molecular systems. These simulations play an important role in a number of areas, such as photovoltaics, semiconductor technology and nanoscience. The excited states are treated using two methods, time-dependent Hartree–Fock (TDHF) and time-dependent density-functional theory (TDDFT), both within a commonly used Tamm–Dancoff approximation (TDA). The resulting TDA eigenvalue equations are solved on a D-Wave quantum annealer using the Quantum Annealer Eigensolver (QAE), developed previously. The method is shown to reproduce a typical basis set convergence on the example $$\hbox {H}_2$$ H 2 molecule and is also applied to several other molecular species. Characteristic properties such as transition dipole moments and oscillator strengths are computed as well. Three potential energy profiles for excited states are computed for $$\hbox {NH}_3$$ NH 3 as a function of the molecular geometry. Similar to previous studies, the accuracy of the method is dependent on the accuracy of the intermediate meta-heuristic software called qbsolv.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

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

Determination of single-crystal elastic moduli of Li RE F 4 ( RE =Y, Gd, and Tb) by resonant ultrasound spectroscopy

The tetragonal fluoro-scheelite Li RE F 4 compounds (RE = rare earth) have been shown to exhibit a variety of useful optical and magnetic properties. While LiYF 4 has been widely studied, many of the fundamental thermodynamic properties of other members of this family remain unknown. Here, we report the complete elastic tensors (C ij ) of single-crystalline LiYF 4 , LiTbF 4 , and LiGdF 4 using resonant ultrasound spectroscopy and density functional theory (DFT). Here we compare the results for LiYF 4 with prior experimental results using time-of-flight ultrasound methods. This is the first report, however, of the experimental elastic tensors of LiTbF 4 and LiGdF 4 . The present results point to a softening of the elastic moduli of the Li RE F 4 system when Y is replaced by the larger ionic radius of Tb or Gd. Furthermore, we find that just 0.3% doping with Nd on the Y site also leads to a slight softening of the moduli. The variation of the elastic moduli as a function of temperature up to 216°C was also measured. A nearly linear softening of all seven independent elastic moduli was observed with increasing temperature. Phonon dispersions and phonon density of states obtained by DFT support the experimental finding of a significantly higher sound velocity due to lighter Y atoms in LiYF 4 , as compared to heavier LiTbF 4 and LiGdF 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

The electronic structure, crystal fields, and magnetic anisotropy in RECo 5 magnets

The current progress in describing rare-earth-based magnets' electronic structure and magnetic properties is discussed. We use several currently popular electronic structure methods to show the typical values of critical parameters that define the physics of RECo 5 (RE = rare earth atom) materials. The magnetic moments and magnetic anisotropy of 4f atoms are obtained using several approaches, including anisotropic 4f-charge density-constrained DFT and DFT+HI methods. We also suggest the introduction of "penalty" functional for obtaining correct variational total energy in the traditional Hund's rule-constrained DFT-based techniques. The applicability and future extensions are discussed. The proposed combination of methods is potentially suitable for high-throughput computational searches of new rare-earth-containing magnetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗
DOE OSTI2020

Gamow-Teller Strength in Ca 48 and Ni 78 with the Charge-Exchange Subtracted Second Random-Phase Approximation

Here, we develop a fully self-consistent subtracted second random-phase approximation for charge-exchange processes with Skyrme energy-density functionals. As a first application, we study Gamow-Teller excitations in the doubly magic nucleus 48 Ca, the lightest double-β emitter that could be used in an experiment, and in 78 Ni, the single-beta-decay rate of which is known. The amount of Gamow-Teller strength below 20 or 30 MeV is considerably smaller than in other energy-density-functional calculations and agrees better with experiment in 48 Ca, as does the beta-decay rate in 78 Ni. These important results, obtained without ad hoc quenching factors, are due to the presence of two-particle–two-hole configurations. Their density progressively increases with excitation energy, leading to a long high-energy tail in the spectrum, a fact that may have implications for the computation of nuclear matrix elements for neutrinoless double-β decay in the same framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Protein Molecular Dynamics Simulations with Approximate QM: What Can We Learn?

Classical force fields are essential for computer simulations of proteins and are typically parameterized to reproduce secondary and tertiary structure of isolated proteins. However, while protein-protein interactions are ubiquitous in nature, they are not considered in parameterization efforts and are far less understood than isolated proteins. A better characterization of intermolecular interactions is widely recognized as a key to revolutionizing drug and therapeutic developments with high-throughput computational screening. Urgently needed is a critical assessment of the performance of modern protein force fields against first-principles electronic structure methods and experiments. In a daring step toward this goal, we here describe a comparison of peptide folding dynamics as predicted by a molecular mechanics force field on the one hand and by an approximate electronic structure quantum mechanical (QM) method based on density-functional tight-binding (DFTB) on the other. We further compare the dynamics from straightforward DFTB simulations with a near-linear scaling version of DFTB for massively parallel computation based on the fragment molecular orbital (FMO-DFTB) method. We illustrate differences between the phenomenology of the folding dynamics from these three methods for a small model peptide, as well as charge polarization and dynamic fluctuations, point out possible correlations and implications for force field developers, and discuss the lessons learned that might become applicable to future predictive high-throughput computer screening for personalized neoantigen cancer therapy.

Irle, Stephan↗
DOE OSTI2023

Interaction energy and isosteric heat of adsorption between hydrogen and magnesium diboride

Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).

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

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗
DOE OSTI2024

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.

algorithmic development↗
DOE OSTI2020

First Principle Investigation of Electronic, Transport, and Bulk Properties of Zinc-Blende Magnesium Sulfide

We have studied electronic, structural, and transport properties of zinc-blende magnesium sulfide (zb-MgS). We employed a local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO) method. Our computational method is able to reach the ground state of a material, as dictated by the second theorem of density functional theory (DFT). Consequently, our findings have the physical content of DFT and agree with available, corresponding experimental ones. The calculated band gap of zb-MgS, a direct gap equal to 4.43 eV, obtained at the experimental lattice constant of 5.620 Å, completely agrees with the experimental band gap of 4.45 ± 0.2 eV. We also report total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated pDOS also agree with the experiment for the description of the states at the top and the bottom of the valence and conduction bands, respectively.

42 ENGINEERING↗
DOE OSTI2024

Density functional theory of material design: fundamentals and applications—II

Abstract This is the second and the final part of the review on density functional theory (DFT), referred to as DFT-II. In the first review, DFT-I, we have discussed wavefunction-based methods, their complexity, and basics of density functional theory. In DFT-II, we focus on fundamentals of DFT and their implications for the betterment of the theory. We start our presentation with the exact DFT results followed by the concept of exchange-correlation (xc) or Fermi-Coulomb hole and its relationship with xc energy functional. We also provide the exact conditions for the xc-hole, xc-energy and xc-potential along with their physical interpretation. Next, we describe the extension of DFT for non-integer number of electrons, the piecewise linearity of total energy and discontinuity of chemical potential at integer particle numbers, and derivative discontinuity of the xc potential, which has consequences on fundamental gap of solids. After that, we present how one obtains more accurate xc energy functionals by going beyond the LDA. We discuss the gradient expansion approximation (GEA), generalized gradient approximation (GGA), and hybrid functional approaches to designing better xc energy functionals that give accurate total energies. However, these functionals fail to predict properties like the ionization potential and the band gap. Thus, we next describe different methods of modelling these potentials and results of their application for calculation of the band gaps of different solids to highlight accuracy of different xc potentials. Finally, we conclude with a glimpse on orbital-free density functional theory and the machine learning approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Why Are 5-Thioglycopyranosyl Donors More Axially Selective than their Glycopyranosyl Counterparts? A Low and Variable Temperature NMR Spectroscopy and Computational Study

5-Thioglycopyranosyl donors differ in reactivity and selectivity from simple glycopyranosyl donors. An extensive study has been conducted on the nature and stability of the reactive intermediates generated on the activation of per-O-acetyl and per-O-methyl 5-thioglucopyranosyl donors and the corresponding glucopyranosyl donors. Variable temperature nuclear magnetic resonance (NMR) studies with per-O-methylated or per-O-acetyl glycosyl sulfoxides and trichloroacetimidates on activation with trifluoromethanesulfonic anhydride or trimethylsilyl triflate are reported. These show that following initial adduct formation with the promoter conversion of the 5-thioglucopyranosyl donors to the 5-thioglucopyranosyl triflates requires higher temperatures than that of the glucopyranosyl donors to the glucopyranosyl triflates. It is demonstrated that neighboring group participation is a less important phenomenon for the peracetylated thioglucosyl donors than for the peracetylated glucosyl donors. A simple thiocarbenium ion was generated by protonation of 2,3-dihydro-4H-thiopyran at low temperature and characterized by NMR spectroscopy. However, the corresponding 5-thioglucopyranosyl thenium ions were not observed in any of the NMR studies of the 5-thiopyranosyl donors: the electron-withdrawing C–O bonds around the thiopyranoside core discourage thiocarbenium ion formation, just as they discourage oxocarbenium ion formation. Density functional theory (DFT) calculations reveal the tetrahydrothiopyranyl thiocarbenium ion to be approximately 2.5 kcal/mol lower in energy than the corresponding tetrahydropyranyl oxocarbenium ion relative to the corresponding covalent triflates. However, the computations reveal a 5.8 kcal/mol activation barrier for conversion of the tetrahydrothiopyranyl triflate to the thiocarbenium ion, while formation of the oxocarbenium ion–triflate ion pair from tetrahydropyranyl triflate requires only 2.6 kcal·mol –1 . Overall, the greater axial selectivity of 5-thioglycopyranosyl donors compared to analogous glycopyranosyl donors derives from (i) the lower kinetic reactivity necessitating higher reaction temperatures, (ii) the greater stability of the thiocarbenium ion over the oxocarbenium ion facilitating equilibration under thermodynamic conditions, (iii) the greater magnitude of the anomeric effect in the 5-thiosugars, and (iv) decreased neighboring group participation in the per-esterified 5-thiosugars.

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

Re-examining the silicon self-interstitial charge states and defect levels: A density functional theory and bounds analysis study

The self-interstitial atom (SIA) is one of two fundamental point defects in bulk Si. Isolated Si SIAs are extremely difficult to observe experimentally. Even at very low temperatures, they anneal before typical experiments can be performed. Given the challenges associated with experimental characterization, accurate theoretical calculations provide valuable information necessary to elucidate the properties of these defects. Previous studies have applied Kohn–Sham density functional theory (DFT) to the Si SIA, using either the local density approximation or the generalized gradient approximation to the exchange-correlation (XC) energy. The consensus of these studies indicates that a Si SIA may exist in five charge states ranging from -2 to +2 with the defect structure being dependent on the charge state. This study aims to re-examine the existence of these charge states in light of recently derived “approximate bounds” on the defect levels obtained from finite-size supercell calculations and new DFT calculations using both semi-local and hybrid XC approximations. We conclude that only the neutral and +2 charge states are directly supported by DFT as localized charge states of the Si SIA. Within the current accuracy of DFT, our results indicate that the +1 charge state likely consists of an electron in a conduction-band-like state that is coulombically bound to a +2 SIA. Furthermore, the -1 and -2 charge states likely consist of a neutral SIA with one and two additional electrons in the conduction band, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Comparison of GW band structure to semiempirical approach for an FeSe monolayer

We present the G 0 W 0 band structure, core levels, and deformation potential of monolayer FeSe in the paramagnetic phase based on a starting mean field of the Kohn-Sham density functional theory (DFT) with the Perdew, Burke, and Ernzerhof functional. We find the GW correction increases the bandwidth of the states forming the M pocket near the Fermi energy, while leaving the Γ pocket roughly unchanged. We then compare the G 0 W 0 quasiparticle band energies with the band structure from a simple empirical + A approach, which was recently proposed to capture the renormalization of the electron-phonon interaction going beyond DFT in FeSe, when used as a starting point in density functional perturbation theory. We show that this empirical correction succeeds in approximating the GW nonlocal and dynamical self-energy in monolayer FeSe and reproduces the GW band structure near the Fermi surface, the core energy levels, and the deformation potential (electron-phonon coupling).

36 MATERIALS SCIENCE↗