Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Time-dependent density functional theory applied to average atom opacity

We focus on studying the opacity of iron, chromium, and nickel plasmas at conditions relevant to experiments carried out at Sandia National Laboratories. In this work, we calculate the photoabsorption cross sections and subsequent opacity for plasmas using linear-response time-dependent density functional theory (TD-DFT). Our results indicate that the physics of channel mixing accounted for in linear-response TD-DFT leads to an increase in the opacity in the bound-free quasicontinuum, where the Sandia experiments indicate that models underpredict iron opacity. However, the increase seen in our calculations is only in the range of 5%–10%. Further, we do not see any change in this trend for chromium and nickel. This behavior indicates that channel mixing effects do not explain the trends in opacity observed in the Sandia experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Time-dependent density functional theory of charge, energy, and spin transport, and quantum continuum mechanics (Final Report)

The overarching goal of this project, which unfolded over a period of seventeen years, including two no-cost extensions in the last two years, was to extend and deepen our ability to use time-dependent density functional theory and related theoretical methods for the description of the transport of charge, spin, or energy in quantum electronic systems at the nanoscale. The output of our efforts is contained in 52 papers, listed below, which we have published since 2006. There have also been many presentations at international conferences, workshops, colloquia and seminars. This body of work has been carried out with the essential input of 8 postdoctoral students (who were supported by the grant) and several external senior collaborators and their postdocs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography

Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

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

Strain-Induced Shifts in Defective Graphite Phonon Modes Predicted by Density Functional Theory

Carbon fiber composites have gained attention as a structural material because of their high strength-to-weight ratio, and understanding the effect of defects on reactivity and mechanical properties is important for the longevity and safety of the composite. Although it is known that strain causes the underlying graphitic vibrational modes to redshift, it is not clear how strain may alter reactivity and defect-induced vibrational changes. To investigate the strain-induced phonon changes of defective carbon fiber composites, density functional theory calculations of graphite are used, including intercalated hydrogen and fluorine defects. By comparing changes in the bond lengths, formation energies, and phonon density of states for uniaxially and biaxially strained graphite, strain was found to generally make defect formation more favorable and the specific behavior changes are dependent on the strain direction and defect identity. Specifically, intercalated fluorine phonons are more sensitive to strain than hydrogen intercalation phonons, and strain applied along the zigzag direction alters the calculated properties more than strain along the armchair direction. Finally, these results highlight the importance of understanding the microstructural effect of deviations from the ideal material because small changes in strain or defect type can significantly alter the behavior of the carbon fiber composite core.

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

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C 2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B 2 , C 2 , O 2 , and F 2 in the local spin density approximation (LSDA), the Perdew–Burke–Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C 2 . Here, this work shows that spin symmetry breaking in singlet C 2 , which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

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

Real-Space Constrained Density Functional Theory Investigation of Site-Specific, Interfacial Charge Recombination Dynamics Across the Au Nanoparticle/TiO 2 Heterojunction

Au nanoparticle (NP)/TiO 2 heterojunction is a representative system to study interfacial charge transfer in photocatalysis and photovoltaics, where suppressing recombination from TiO 2 to Au can enhance hot carrier extraction. We apply real-space constrained density functional theory (CDFT) with Marcus theory to quantify charge recombination time scales across Au/TiO 2 . This approach enables direct control and visualization of charge-separated states, aligning with site-specific probes like time-resolved X-ray photoelectron spectroscopy (trXPS). We find that the charge-separated state features a bipolaron, with recombination dominated by TiO 2 LUMO to Au HOMO transitions, primarily at interfacial Au sites. Marcus rate predictions are benchmarked with surface hopping methods, quantifying differences in time scales and computational efficiency. Lastly, we examine how the Au cluster size affects the free energy change (ΔG) and reorganization energy (λ), explaining trends in closed-shell systems and highlighting challenges for open-shell extrapolations. Overall, CDFT + Marcus theory provides efficient, mechanistically transparent interfacial charge transfer modeling, and we clearly defined its applicability and limitation.

Glenna, Drew M. [Univ. of Idaho, Idaho Falls, ID (↗
DOE OSTI2020

Ni Nanoparticles on CeO 2 (111): Energetics, Electron Transfer, and Structure by Ni Adsorption Calorimetry, Spectroscopies, and Density Functional Theory

The morphology, interfacial bonding energetics and charge transfer of Ni clusters and nanoparticles on slightly-reduced CeO 2-x (111) surfaces at 100 to 300 K have been studied using single crystal adsorption calorimetry (SCAC), low-energy ion scattering spectroscopy (LEIS), Xray photoelectron spectroscopy (XPS), low energy electron diffraction (LEED) and density functional theory (DFT). The initial heat of adsorption of Ni vapor decreased with the extent of pre-reduction (x) of the CeO 2-x (111), showing that stoichiometric ceria adsorbs Ni more strongly than oxygen vacancies. On CeO 1.95 (111) at 300 K, the heat dropped quickly with coverage in the first 0.1 ML, attributed to nucleation of Ni clusters on stoichiometric steps, followed by the Ni particles spreading onto less favorable terrace sites. At 100 K, the clusters nucleate on terraces due to slower diffusion. Adsorbed Ni monomers are in the +2 oxidation state, and they bind by ~45 kJ/mol more strongly to step sites than terraces. The measured heat of adsorption versus average particle size on terraces is favorably compared to DFT calculations. The Ce 3d XPS lineshape showed an increase in Ce 3+ /Ce 4+ ratio with Ni coverage, providing the number of electrons donated to the ceria per Ni atom. The charge transferred per Ni is initially large but strongly decreases with increasing cluster size for both experiments and DFT, and shows large differences between clusters at steps versus terraces. This charge is localized on the interfacial Ni and Ce atoms in their atomic layers closest to the interface. This knowledge is crucial to understanding the nature of the active sites on the surface of Ni-CeO 2 catalysts for which metal-oxide interactions play a very important role in the activation of O-H and C-H bonds. In conclusion, the changes in these interactions with Ni particle size (metal loading) and the extent of reduction of the ceria help to explain how previously reported catalytic activity and selectivity change with these same structural details.

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

Density Functional Theory Study of Controllable Optical Absorptions and Magneto-Optical Properties of Magnetic CrI 3 Nanoribbons: Implications for Compact 2D Magnetic Devices

A chromium triiodide (CrI 3 ) monolayer has an interesting ferromagnetic ground-state. In this work, we calculate band structures and magnetic moments of tensile-strained and bent zigzag CrI 3 nanoribbons with density functional theory. The edge iodine atoms form flat low-lying conduction bands and couple with chromium atoms ferromagnetically, while the non-edge iodine atoms weakly couple antiferromagnetically. Narrow CrI 3 nanoribbons have two locally stable magnetic moment orientations, namely out-of-plane and in-plane (along the nanoribbon periodic direction) configurations. This enables four magnetization states in CrI 3 nanoribbons, including two out-of-plane ones (up and down) and two in-plane ones (forward and backward along the nanoribbon periodical direction), increasing the operating controllability. Based on the one-dimensional Ising spin chain model, the spin correlation length of the narrow CrI 3 nanoribbon is estimated as about 10 Å at its estimated Curie temperature of 27 Kelvin, which is lower than the measured 45 Kelvin of the monolayer CrI 3 . The optical absorption and magneto-optical properties of CrI 3 nanoribbons are investigated with many-body perturbation GW-BSE (Bethe-Salpeter equation), including magnetic dichroism, Faraday and magneto-optical Kerr effects. Here, the low-energy dark excitons are mainly from transitions between electrons and holes with unlike spins and are non-Frenkel-like, while the bright excitons have mixed spin configurations. The intrinsic lifetime of excitons can be over one nanosecond, suitable for quantum information processes. Tensile strains and bending manifestly modulate the absorption spectra and magneto-optical properties of CrI 3 nanoribbons within a technologically important photon energy-range of ~1.0-2.0 eV. The CrI 3 nanoribbons can be used in 1D or 2D magnetic storage nanodevices, tunable magnetic optoelectronics, and spin-based quantum information controls.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of single atom Platinum (Pt) doping and facet dependent on the electronic structure and light absorption of Lanthanum Titanium Oxide (La 2 Ti 2 O 7 ): A Density Functional Theory study

Charge generation and separation are two key processes for semiconductor photocatalysis. Here, we use Pt as a single atom catalyst to systematically examine the facet-dependent electronic band structure and light absorption of the layered perovskite-type wide-gap semiconductor, lanthanum titanium oxide (La 2 Ti 2 O 7 , LTO) by means of density functional theory simulations. In this study, it is found that single Pt atom doping of different LTO surfaces (here, (100), (101) and (001)) can not only create states in the bandgap that would promote the formation of recombination centers, but also shift the optical absorption edge to the visible region. Interestingly, the Pt doping forms a heterojunction, with the valence band maximum consisting of Pt states and the conduction band minimum consisting of LTO states, respectively. To study the facet dependent surface activity, acetic acid (CH 3 COOH) was used as a model molecule to investigate the adsorption and charge transfer on the (101), (100) and (001) Pt-LTO surface facets. The results show that the (101) facet could enable stronger adsorption of CH 3 COOH by promoting more electron transfer during the interfacial interaction. Our theoretical findings aim to promote the design and optimization of the single atom catalysts (SACs) for photocatalytic applications and other broad catalysis systems.

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

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

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

Excitation energy of fission fragments within nuclear time-dependent density functional theory

The number and properties of the neutrons and photons emitted in nuclear fission are directly related to the excitation energy of the fission fragments when they are formed at scission. Though not observable experimentally because of the extremely short timescales, the excitation energy of fission fragments can be predicted by microscopic theory based on time-dependent density functional theory (TDDFT). Initial results on the value of the total kinetic energy of fission reactions were very promising, but could not probe all possible fragmentations. Here, in this work, we perform large-scale TDDFT calculations in 240 Pu enabled by the development of a new TDDFT solver. We obtain TDDFT trajectories covering nearly all possible fragmentations. We find that the total kinetic energy is close to experimental values only for the most likely fission while it is severely underestimated at both small and large asymmetries. This conclusion seems rather independent of the parametrization of the energy functional, both in its particle-hole and particle-particle channels.

Physics - Nuclear physics and radiation physics↗
DOE OSTI2021

A probability density function model describing height estimation uncertainty due to image pixel intensity noise in digital fringe projection measurements

Digital fringe projection is a surface-profiling technique used for highly accurate non-contact measurements. As with any measurement technique, a variety of sources degrade to the measurement accuracy of the method. Here, this paper presents an analytically-derived probability density function that explicitly models the surface height measurement error due to inevitable phase measurement error, and it includes the specific case of pixel noise inducing the phase measurement error that ultimately leads to the height estimation error. The accuracy of the model was validated through Monte-Carlo simulations of resultant height distributions subject to arbitrarily correlated pixel intensity noise and experimental digital fringe projection measurements where the pixel-by-pixel height uncertainty estimations were compared to the predictions of the derived model.

42 ENGINEERING↗
DOE OSTI2025

Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Adsorption↗
DOE OSTI2021

van der Waals corrected density functionals for cylindrical surfaces: Ammonia and nitrogen dioxide adsorbed on a single-walled carbon nanotube

In this work, we extend the damped Zaremba-Kohn model (dZK) for long-range dispersion interaction between a molecule and a planar surface to molecules adsorbed on a curved cylindrical surface, and employ this extended model as an additive correction to the semilocal density functionals PBE (Perdew-Burke-Ernzerhof) and SCAN (strongly constrained and appropriately normed). The resulting PBE+vdW (van der Waals)-dZK and SCAN+vdW-dZK are applied to two systems, NH3 and NO2 molecules adsorbed on a single-wall carbon nanotube (CNT), for calculations of binding energies and equilibrium distances. For comparison, the results from vdW nonlocal functionals, such as SCAN+rVV10 and PBE+rVV10, are also presented. The binding energies from PBE+rVV10 (Vydrov and Van Voorhis), SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK are about 70–115 meV for the system of CNT + NH 3 and 300–500 meV for the system of CNT + NO 2 . The results from PBE+vdW-dZK and SCAN+vdW-dZK are closer to each other than those from PBE+rVV10 and SCAN+rVV10 are. The relatively closer results from PBE+vdW-dZK and SCAN+vdW-dZK indicate the consistency of our developed vdW–dZK model for cylindrical surfaces. All methods, including PBE, SCAN, PBE+rVV10, SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK, give approximately the same binding energy differences between two adsorption configurations (types I and II) for the two systems. This implies that the two adsorption sites have approximately the same adsorption stability. The exponent of the vdW interaction power law from our vdW-dZK model for the two systems is about 0 at short distance, largely due to the damping factor, and tends slowly to –4 to –4.5 at distances D about 20–50 Å. At even larger distances, the vdW power-law exponent approaches –5. This feature is very similar to the one calculated with random-phase approximation and renormalization group approaches, supporting the applicability of our methods. Our developed vdW-dZK method provides a highly efficient and reliable method for large systems with cylindrical surfaces, such as vdW interactions with nanotubes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Density Functional Theory Investigation of the NiO@Graphene Composite as a Urea Oxidation Catalyst in the Alkaline Electrolyte

Developing efficient and low-cost urea oxidation reaction (UOR) catalysts is a promising but still challenging task for environment and energy conversion technologies such as wastewater remediation and urea electrolysis. In this work, NiO nanoparticles that incorporated graphene as the NiO@Graphene composite were constructed to study the UOR process in terms of density functional theory. The single-atom model, which differed from the previous heterojunction model, was employed for the adsorption/desorption of urea and CO 2 in the alkaline media. As demonstrated from the calculated results, NiO@Graphene prefers to adsorb the hydroxyl group than urea in the initial stage due to the stronger adsorption energy of the hydroxyl group. After NiOOH@Graphene was formed in the alkaline electrolyte, it presents excellent desorption energy of CO 2 in the rate-determining step. Electronic density difference and the d band center diagram further confirmed that the Ni(III) species is the most favorable site for urea oxidation while facilitating charge transfer between urea and NiO@Graphene. Moreover, graphene provides a large surface for the incorporation of NiO nanoparticles, enhancing the electron transfer between NiOOH and graphene and promoting the mass transport in the alkaline electrolyte. Notably, this work provides theoretical guidance for the electrochemical urea oxidation work.

25 ENERGY STORAGE↗
DOE OSTI2020

Toward Engineering of Solution Microenvironments for the CO 2 Reduction Reaction: Unraveling pH and Voltage Effects from a Combined Density-Functional–Continuum Theory

Engineering the electrolyte microenvironment represents an attractive route to tuning the selectivity of electrocatalytic reactions beyond catalyst composition and morphology. However, harnessing the full potential of this approach requires understanding the interplay between voltage, electrolyte composition, and adsorbate binding within the electrical double layer, which is absent from the usual theoretical approaches. In this study, we apply a recently developed density functional theory (DFT)–continuum approach based on the effective screening medium method and reference interaction site model (ESM–RISM) to explore electrolyte effects with an enhanced description of the electrochemical interface. Furthermore, applying this method to the binding of CO adsorbates in potassium-containing electrolytes on copper, a problem of direct relevance to CO 2 electroreduction to value-added products, we show that the interdependence of voltage and pH leads to an unexpected change in adsorption site preference on Cu(001) terraces. Our findings highlight the often-overlooked importance of the electrical double-layer structure for predicting catalyst operation.

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

Critical Assessment of the Thermodynamics of Vacancy Formation in Fe 2 O 3 Using Hybrid Density Functional Theory

Fe 2 O 3 hematite is a technologically important material with applications in energy storage as well as being the key phase formed in the rust of iron-based materials. Despite this central importance, there is much that is still unknown regarding the properties of defects in Fe 2 O 3 and consequently oxide growth. In this work, using screened hybrid density functional theory (HSE06), we consider the thermodynamics of vacancies in Fe 2 O 3 , considering the effects of ionic and electronic chemical potentials on both iron and oxygen vacancy formation. We find that, in the oxygen-rich limit, iron vacancies are easier to form, though the difference in formation energy between the two vacancies is only about 1 eV. In contrast, in the Fe-rich limit, oxygen vacancies have an extremely low formation energy, only 0.07 eV at mid-gap, and would spontaneously form as the Fermi level is reduced, while Fe vacancies require over 5 eV to form. Consistent with experiment, this indicates that Fe 2 O 3 is relatively easily reduced but not oxidized. However, the theoretical picture is very different when considering other exchange–correlation functionals (GGA + U or SCAN), emphasizing the critical role of the exchange–correlation functional in describing this system accurately.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗