Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density Functionals”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2024

Using ground state and excited state density functional theory to decipher 3d dopant defects in GaN

Abstract Using ground state density functional theory (DFT) and implementing an occupation-constrained DFT (occ-DFT) for self-consistent excited state calculations, we decipher the electronic structure of the Mn dopant and other 3 d defects in GaN across the band gap. Our analysis, validated with broad agreement with defect levels (ground-state calculations) and photoluminescence data (excited-state calculations), mandates reinterpretation and reassignment of 3 d defect data in GaN. The Mn Ga defect is determined to span stable charge states from (1−) in n -type GaN through (2+) in p -type GaN. The Mn(2+) is predicted to be a d 2 ground state spin triplet defect with a singlet excited state, isoelectronic with the defect associated with the 1.19 eV photoluminescence in n -type GaN. The combined analysis of defect levels and excited states invites reassessment of all d 2 -capable dopants in GaN. We demonstrate that the 1.19 eV defect, a candidate defect for optically controlled quantum applications, cannot be the Cr(1+) assumed in literature and instead must be the V(0). The combined ground-state/excited-state DFT analysis is shown to be able to chemically fingerprint defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1994

A Comparison of Density Functional Theory with Ab initio Approaches for Systems Involving First Transition Row Metals

Density functional theory (DFT) is found to give a better description of the geometries and vibrational frequencies of FeL and FeL(sup +) systems than second order Moller Plesset perturbation theory (MP2). Namely, the DFT correctly predicts the shift in the CO vibrational frequency between free CO and the Sigma(sup -) state of FeCO and yields a good result for the Fe-C distance in the quartet states of FeCH4(+) 4 These are properties where the MP2 results are unsatisfactory. Thus DFT appears to be an excellent approach for optimizing the geometries and computing the zero-point energies of systems containing first transition row atoms. Because the DFT approach is biased in favor of the 3d(exp 7) occupation, whereas the more traditional approaches are biased in favor of the 3d(exp 6) occupation, differences are found in the relative ordering of states. It is shown that if the dissociation is computed to the most appropriate atomic asymptote and corrected to the ground state asymptote using the experimental separations, the DFT results are in good agreement with high levels of theory. The energetics at the DFT level are much superior to the MP2 and in most cases in good agreement with high levels of theory.

Ricca, Alessandra↗
DOE OSTI2024

QRCODE: Massively parallelized real-time time-dependent density functional theory for periodic systems

We present a new software module, QRCODE (Quantum Research for Calculating Optically Driven Excitations), for massively parallelized real-time time-dependent density functional theory (RT-TDDFT) calculations of periodic systems in the open-source Qbox software package. Our approach utilizes a custom implementation of a fast Fourier transformation scheme that significantly reduces inter-node message passing interface (MPI) communication of the major computational kernel and shows impressive scaling up to 16,344 CPU cores. In addition to improving computational performance, QRCODE contains a suite of various time propagators for accurate RT-TDDFT calculations. As benchmark applications of QRCODE, we calculate the current density and optical absorption spectra of hexagonal boron nitride (h-BN) and photo-driven reaction dynamics of the ozone-oxygen reaction. We also calculate the second and higher harmonic generation of monolayer and multi-layer boron nitride structures as examples of large material systems. Our optimized implementation of RT-TDDFT in QRCODE enables large-scale calculations of real-time electron dynamics of chemical and material systems with enhanced computational performance and impressive scaling across several thousand CPU cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory↗
DOE OSTI2025

Global optimization of harmonic oscillator basis in covariant density functional theory

The present investigation focuses on the improvement of the accuracy of the description of binding energies within moderately sized fermionic basis. Using the solutions corresponding to infinite fermionic basis it was shown that in the case of meson exchange (ME) covariant energy density functionals (CEDFs) the global accuracy of the description of binding energies in the finite $N_F$ = 16 - 20 bases can be drastically (by a factor ranging from ~3 up to ~9 dependent on the functional and $N_F$) improved by a global optimization of oscillator frequency of the basis. This is a consequence of the unique feature of the ME functionals in which with increasing fermionic basis size fermionic and mesonic energies approach the exact (infinite basis) solution from above and below, respectively. As a consequence, an optimal oscillator frequency $\hbar\omega_0$ of the basis can be defined which provides an accurate reproduction of exact total binding energies by the ones calculated in truncated basis. This leads to a very high accuracy of the calculations in moderately sized $N_F=20$ basis when mass dependent oscillator frequency is used: global rms differences $\delta B_{rms}$ between the binding energies calculated in infinite and truncated bases are only 0.025 MeV and 0.031 MeV for the NL5(Z) and DD-MEZ functionals, respectively. Optimized values of the oscillator frequency $\hbar\omega_0$ are provided for three major classes of CEDFs, i.e. for density dependent meson exchange functionals, nonlinear meson exchange ones and point coupling functionals.

Binding energy & masses↗
DOE OSTI2025

Time-Dependent Density Functional Theory Description of 238 U⁡(n,f), 240,242 Pu⁢(n,f), and 237 Np(n,f) Reactions

In nuclei with an odd nucleon number the nonvanishing spin number density is the source of a pseudomagnetic field, which favors the splitting of the nucleon Cooper pairs. Such a pseudomagnetic field is generated always in the dynamics of any nucleus, but its effects on Cooper pairs are significantly enhanced in the dynamic evolution of nuclei with an odd number of nucleons. We present for the first time a microscopic study of the induced fission of the odd neutron compound nuclei 239 U, 241,243 Pu, and the odd proton, odd neutron compound nucleus 238 Np, performed within the time-dependent density functional theory extended to superfluid fermion systems, without any simplifying assumptions, with controlled numerical approximations, and for a very large number of initial conditions. Because of the presence of the unpaired odd nucleon(s), the time-reversal symmetry of the fission compound nucleus is spontaneously broken, an aspect routinely neglected in the most advanced microscopic approaches of the past. The emerging fission fragment properties are quite similar to the properties of fission fragments of neighboring even-even nuclei. The time from saddle-to-scission is often significantly longer in odd-odd or odd-mass nuclei than for even-even nuclei, since systems with unpaired nucleons are easier to excite and the potential energy surfaces of these nuclei have more structure, often resembling a very complicated obstacle course, rather than a more direct evolution of the nuclear shape from the top of the outer fission barrier to the scission configuration. The Pauli blocking approximation, often invoked in the literature, expected to inhibit the fission of nuclei with unpaired nucleons, is surprisingly strongly violated during the fission dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Laser pulse induced second- and third-harmonic generation of gold nanorods with real-time time-dependent density functional tight binding (RT-TDDFTB) method

In this paper, we investigate second- and third-harmonic generation processes in Au nanorod systems using the real-time time-dependent density functional tight binding method. Our study focuses on the computation of nonlinear signals based on the time dependent dipole response induced by linearly polarized laser pulses interacting with nanoparticles. We systematically explore the influence of various laser parameters, including pump intensity, duration, frequency, and polarization directions, on harmonic generation. We demonstrate all the results using Au nanorod dimer systems arranged in end-to-end configurations, and disrupting the spatial symmetry of regular single nanorod systems is crucial for second-harmonic generation processes. Furthermore, we study the impact of nanorod lengths, which lead to variable plasmon energies, on harmonic generation, and estimates of polarizabilities and hyper-polarizabilities are provided.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗
NASA NTRS1995

Probability Density Functions of Observed Rainfall in Montana

The question of whether a rain rate probability density function (PDF) can vary uniformly between precipitation events is examined. Image analysis on large samples of radar echoes is possible because of advances in technology. The data provided by such an analysis easily allow development of radar reflectivity factors (and by extension rain rate) distribution. Finding a PDF becomes a matter of finding a function that describes the curve approximating the resulting distributions. Ideally, one PDF would exist for all cases; or many PDF's that have the same functional form with only systematic variations in parameters (such as size or shape) exist. Satisfying either of theses cases will, validate the theoretical basis of the Area Time Integral (ATI). Using the method of moments and Elderton's curve selection criteria, the Pearson Type 1 equation was identified as a potential fit for 89 percent of the observed distributions. Further analysis indicates that the Type 1 curve does approximate the shape of the distributions but quantitatively does not produce a great fit. Using the method of moments and Elderton's curve selection criteria, the Pearson Type 1 equation was identified as a potential fit for 89% of the observed distributions. Further analysis indicates that the Type 1 curve does approximate the shape of the distributions but quantitatively does not produce a great fit.

Larsen, Scott D.↗
DOE OSTI2025

Probability Density Function Control of Frequency Fluctuations in Renewable-Rich Power Systems

The stochastic nature of renewable energy sources (RESs) necessitates treating power system frequency response as a random process with a nonstationary probability density function (PDF). Here, based upon the stochastic distribution control theory originated by the second author, this paper proposes a novel stochastic controller to improve the frequency PDF in power grids when integrating a large amount of RESs, thereby minimizing the effects of uncertainties and enhancing overall system stability. The key idea is to manipulate the controllable power generation resources so that the frequency PDF is make to follow a target PDF by using the stochastic distribution control theory originated by the second author. The proposed method can easily be plugged into existing automatic generation controls for multi-area transmission grids. The proposed method is validated via a modified Kundar's two area system and 240-bus Western Electricity Coordinating Council systems. The simulation results show that the proposed control shapes the frequency PDF narrower and sharper, leading to a notable improvement toward minimizing the effects of randomness and uncertainty during grid operation.

frequency↗
DOE OSTI2024

Tunneling magnetoresistance in MgO tunnel junctions with Fe-based leads in empirically corrected density functional theory

The minority-spin Fe/MgO interface states are at the Fermi level in density functional theory (DFT), but experimental evidence and GW calculations place them slightly higher in energy. This small shift can strongly influence tunneling magnetoresistance (TMR) in junctions with a thin MgO barrier and its dependence on the concentration of Co in the electrodes. Here, in this study, an empirical potential correction to DFT is introduced to shift the interface states up to match the tunnel spectroscopy data. With this shift, TMR in Fe/MgO/Fe junctions exceeds 800% and 3000% at 3 and 4 monolayers (ML) of MgO, respectively. We further consider the effect of alloying of the Fe electrodes with up to 30% Co or 10% V, treating them in the coherent potential approximation (CPA). Alloying with Co broadens the interface states and brings a large incoherent minority-spin spectral weight to the Fermi level. Alloying with V brings the minority-spin resonant states close to the Fermi level. However, in both cases the minority-spin spectral weight at the Fermi level resides primarily at the periphery of the Brillouin zone, which is favorable for spin filtering. Using convolutions of k II -resolved barrier densities of states calculated in CPA, it is found that TMR is strongly reduced by alloying with Co or V but still remains above 500% at 4 ML of MgO up to 30% of Co or 5% V. At 5 ML, the TMR increases above 1000% in all systems considered. However, while TMR declines sharply with increasing bias up to 0.2 eV in the tunnel junctions with pure Fe leads, it remains almost constant up to 0.5 eV if leads are alloyed with Co.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Halogen Thermochemistry Assessed with Density Functional Theory: Systematic Errors, Swift Corrections and Effects on Electrochemistry

Abstract Despite its sizable errors, density functional theory (DFT) is extensively used to evaluate thermochemical properties of gases, liquids and their interfaces with solids. As numerous halogen‐containing compounds appear as reactants, products and/or electrolytes in electrochemical reactions, and ionic effects are currently an active area of research, it is important to evaluate the accuracy of DFT for halogen thermochemistry. Herein, we assess the formation energies of interhalogens, hydrogen halides, diatomic and atomic halogens and their ions using six widespread functionals at the GGA, meta‐GGA and hybrid levels. We observe that DFT errors with respect to experiments are correlated with the electronegativity of the species and there are systematic trends across functionals, such that swift corrections were devised. Specifically, the average of the mean absolute errors for the six functionals decreased from 0.19 eV before the corrections to 0.08 eV after them. Besides, the overall maximum absolute error (MAX) decreased from 0.76 to 0.44 eV and the average of the MAXs decreased from 0.51 to 0.24 eV. Finally, we illustrate the qualitative and quantitative impact of gas‐phase errors on the predictions of surface Pourbaix diagrams.

Chemistry↗
DOE OSTI2025

Synergistic Density Functional Theory and Molecular Dynamics Approach to Elucidate PNIPAM–Water Interaction Mechanisms

This study employs Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations to investigate interactions between water molecules and Poly(N-isopropylacrylamide) (PNIPAM). DFT reveals preferential water binding sites, with enhanced binding energy observed in the linker zone. Quantum Theory of Atoms in Molecules (QTAIM) and electron localization function (ELF) analyses highlight the roles of hydrogen bonding and steric hindrance. MD simulations unveil temperature-dependent hydration dynamics, with structural transitions marked by changes in the radius of gyration (Rg) and the radial distribution function (RDF), aligning with DFT findings. Our work goes beyond prior studies by combining a DFT, QTAIM and MD simulations approach across different PNIPAM monomer-to-30mer structures. It introduces a systematic quantification of pseudo-saturation thresholds and explores water clustering dynamics with structural specificity, which have not been previously reported in the literature. These novel insights establish a more complete molecular-level picture of PNIPAM hydration behavior and temperature responsiveness, emphasizing the importance of amide hydrogen and carbonyl oxygen sites in hydrogen bonding, which weakens above the lower critical solution temperature (LCST), resulting in increased hydrophobicity and paving the way for understanding water sorption mechanisms, offering guidance for future applications such as dehumidification and atmospheric water harvesting.

08 HYDROGEN↗
NASA NTRS1995

Using Density Functional Theory (DFT) for the Calculation of Atomization Energies

The calculation of atomization energies using density functional theory (DFT), using the B3LYP hybrid functional, is reported. The sensitivity of the atomization energy to basis set is studied and compared with the coupled cluster singles and doubles approach with a perturbational estimate of the triples (CCSD(T)). Merging the B3LYP results with the G2(MP2) approach is also considered. It is found that replacing the geometry optimization and calculation of the zero-point energy by the analogous quantities computed using the B3LYP approach reduces the maximum error in the G2(MP2) approach. In addition to the 55 G2 atomization energies, some results for transition metal containing systems will also be presented.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2025

Recent progress in global optimizations of covariant energy density functionals

The recent progress on global optimizations of covariant energy den­sity functionals (CEDFs) and global calculations of binding energies within the covariant density functional theory (CDFT) has been analyzed and reviewed. Recently developed anchor-based optimization approach of Ref. [1] allows global optimizations of CEDFs at a reasonable numerical cost. Moreover, it permits such optimizations in a very large fermionic basis with a proper extrap­olation to an infinite one. This allows to accurately estimate global calculation errors due to use of truncated fermionic basis and neglect of some contributions to binding energies (such as total electron binding energy).

Afanasjev, A. V. [Mississippi State University, MS↗