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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Organic Photovoltaic Efficiency Predictor: Data-Driven Models for Non-Fullerene Acceptor Organic Solar Cells

In the design of organic solar cells, there has been a need for materials with high power conversion efficiencies. Scharber’s model is commonly used to predict efficiency; however, it exhibits poor performance with new non-fullerene acceptor (NFA) devices, since it was designed for fullerene-based devices. In this work, an empirical model is proposed that can be a more accurate alternative for NFA organic solar cells. Additionally, many screening studies use computationally expensive methods. A model based on using semiempirical simplified time-dependent density functional theory (sTD-DFT) as an alternative method can accelerate the calculations and yield a similar accuracy. The models presented in this paper, termed organic photovoltaic efficiency predictor (OPEP) models, have shown significantly lower errors than previous models, with OPEP/B3LYP yielding errors of 1.53% and OPEP/sTD- DFT of 1.55%. As a result, the proposed computational models can be used for the fast and accurate screening of new high-efficiency NFAs/donor pairs.

14 SOLAR ENERGY↗
DOE OSTI2020

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Connection between coherent phonons and electron-phonon coupling in Sb (111)

We report time- and angle-resolved photoemission spectroscopy (trARPES) measurements on the Sb(111) surface. We observe band- and momentum-dependent binding-energy oscillations in the bulk and surface bands driven by $A_{1g}$ and $E_{g}$ coherent phonons. While the bulk band shows simultaneous $A_{1g}$ and $E_{g}$ oscillations, the surface bands show either $A_{1g}$ or $E_{g}$ oscillations. The observed behavior is reproduced by frozen-phonon calculations based on density-functional theory. Furthermore, this evidences the connection between electron-phonon coupling and coherent binding energy dynamics tied to lattice vibration and confirms that band-, momentum-, and mode-dependent electron-phonon coupling can indeed be probed by trARPES in the low fluence limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Simulated field-modulated x-ray absorption in titania

Here, we present a method to compute the x-ray absorption near-edge structure (XANES) spectra of solid-state transition metal oxides using real-time time-dependent density functional theory, including spin–orbit coupling effects. This was performed on bulk-mimicking anatase titania (TiO 2 ) clusters, which allows for the use of hybrid functionals and atom-centered all electron basis sets. Furthermore, this method was employed to calculate the shifts in the XANES spectra of the Ti L-edge in the presence of applied electric fields to understand how external fields can modify the electronic structure, and how this can be probed using x-ray absorption spectroscopy. Specifically, the onset of t 2g peaks in the Ti L-edge was observed to red shift and the e g peaks were observed to blue shift with increasing fields, attributed to changes in the hybridization of the conduction band (3d) orbitals.

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

Effect of Fluoroethylene Carbonate Additives on the Initial Formation of the Solid Electrolyte Interphase on an Oxygen-Functionalized Graphitic Anode in Lithium-Ion Batteries

The formation of a solid electrolyte interphase (SEI) at the electrode/electrolyte interface substantially affects the stability and lifetime of lithium-ion batteries (LIBs). One of the methods to improve the lifetime of LIBs is by the inclusion of additive molecules to stabilize the SEI. To understand the effect of additive molecules on the initial stage of SEI formation, in this work we compare the decomposition and oligomerization reactions of a fluoroethylene carbonate (FEC) additive on a range of oxygen-functionalized graphitic anodes to those of an ethylene carbonate (EC) organic electrolyte. A series of density functional theory (DFT) calculations augmented by ab initio molecular dynamics (AIMD) simulations reveal that EC decomposition on an oxygen-functionalized graphitic ($11\bar{2}0$) edge facet through a nucleophilic attack on an ethylene carbon site (C E ) of an EC molecule (S2 mechanism) is spontaneous during the initial charging process of LIBs. However, decomposition of EC through a nucleophilic attack on a carbonyl carbon (C C ) site (S1 mechanism) results in alkoxide species regeneration that is responsible for continual oligomerization along the graphitic surface. In contrast, FEC prefers to decompose through an S1 pathway, which does not promote alkoxide regeneration. Including FEC as an additive is thus able to suppress alkoxide regeneration and results in a smaller and thinner SEI layer that is more flexible toward lithium intercalation during the charging/discharging process. In addition, we find that the presence of different oxygen functional groups at the surface of graphite dictates the oligomerization products and the LiF formation mechanism in the SEI.

25 ENERGY STORAGE↗
DOE OSTI2022

Evidence of Two-Source King Plot Nonlinearity in Spectroscopic Search for New Boson

Optical precision spectroscopy of isotope shifts can be used to test for new forces beyond the standard model, and to determine basic properties of atomic nuclei. We measure isotope shifts on the highly forbidden 2 S 1/2 → 2 F 7/2 octupole transition of trapped 168,170, 172, 174, 176 Yb ions. When combined with previous measurements in Yb + and very recent measurements in Yb, the data reveal a King plot nonlinearity of up to 240σ. Furthermore, the trends exhibited by experimental data are explained by nuclear density functional theory calculations with the Fayans functional. We also find, with 4.3σ confidence, that there is a second distinct source of nonlinearity, and discuss its possible origin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M –1 cm –1 ) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Detonation-induced transformation of graphite to hexagonal diamond

In this work, we explore the structural evolution of highly oriented pyrolytic graphite (HOPG) under detonation-induced shock conditions using in situ synchrotron x-ray diffraction in the ns timescale. We observe the formation of hexagonal diamond (lonsdaleite) at pressures above 50 GPa, in qualitative agreement with recent gas gun experiments. First-principles density functional calculations reveal that under uniaxial compression, the energy barrier for the transition toward hexagonal diamond is lower than that for cubic diamond. Finally, no indication of cubic diamond formation was observed up to >70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Signatures of van Hove singularities in the anisotropic in-plane optical conductivity of the topological semimetal Nb 3 ⁢SiTe 6

In this work, we present a temperature-dependent infrared spectroscopy study on the layered topological semimetal Nb 3 ⁢SiTe 6 combined with density-functional theory calculations of the electronic band structure and optical conductivity. Our results reveal an anisotropic behavior of the in-plane (ac-plane) optical conductivity, with three pronounced excitations located at around 0.15, 0.28, and 0.41 eV for the polarization of the incident radiation along the c axis. These excitations are well reproduced in the theoretical spectra. Based on the ab initio results, the excitations around 0.15 and 0.28 eV are interpreted as fingerprints of van Hove singularities in the electronic band structure and compared to the findings for other topological semimetals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Effects of carbon concentration on the local atomic structure of amorphous GST

Ge-Sb-Te (GST) alloys are leading phase-change materials for data storage due to the fast phase transition between amorphous and crystalline states. Ongoing research aims at improving the stability of the amorphous phase to improve retention. This can be accomplished by the introduction of carbon as a dopant to Ge 2 Sb 2 Te 5 , which is known to alter the short- and mid-range structure of the amorphous phase and form covalently bonded C clusters, both of which hinder crystallization. The relative importance of these processes as a function of C concentration is not known. Here we used molecular dynamics simulation based on density functional theory to study how carbon doping affects the atomic structure of GST-C. Carbon doping results in an increase in tetrahedral coordination, especially of Ge atoms, and this is known to stabilize the amorphous phase. We observe an unexpected, non-monotonous trend in the number of tetrahedral bonded Ge with the amount of carbon doping. Our simulations show an increase in the number of tetrahedral bonded Ge up to 5 at.% C, after which the number saturates and begins to decrease above 14 at.% C. The carbon atoms aggregate into clusters, mostly in the form of chains and graphene flakes, leaving less carbon to disrupt the GST matrix at higher carbon concentrations. Different degrees of carbon clustering can explain divergent experimental results for recrystallization temperature for carbon doped GST.

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

Interfacial charge transfer and interaction in the MXene/TiO 2 heterostructures

Hybrid materials of MXenes [two-dimensional (2D) carbides and nitrides] and transition-metal oxides have shown great promise in electrical energy storage (EES) and 2D heterostructures have been proposed as the next-generation electrode materials to expand the limits of current technology. Here we use first principles density functional theory to investigate the interfacial structure, energetics, and electronic properties of the heterostructures of MXenes (Ti n+1 C n T 2 ; T = terminal groups) and anatase TiO 2 . We find that the greatest work-function differences are between OH-terminated MXene (1.6 eV) and anatase TiO 2 (101) (6.4 eV), resulting in the largest interfacial electron transfer (~0.9e/nm 2 across the interface) from MXene to the TiO 2 layer. This interface also has the strongest adhesion and is further strengthened by hydrogen bond formation. For O–, F–, or mixed O–/F– terminated Ti n+1 C n MXenes, electron transfer is minimal and interfacial adhesion is weak for their heterostructures with TiO 2 . The strong dependence of the interfacial properties of the MXene/TiO 2 heterostructures on the surface chemistry of the MXenes will be useful to tune the heterostructures for EES applications.

2-dimensinal systems↗
DOE OSTI2022

Unoccupied electronic structure of actinide dioxides

Ligand field density functional theory calculations of the dioxides of thorium, uranium, and plutonium have been combined with high-energy-resolution fluorescence detection (HERFD) in x-ray absorption spectroscopy and inverse photoelectron spectroscopy (IPES) measurements to provide powerful insight into the underlying composition of the unoccupied 5f electronic structure in these 5f localized systems. Fine structure in the 5f 5/2 transitions in HERFD can be directly correlated with the fine structure in the leading edge of the IPES. The shapes, intensities, and systematics in HERFD and IPES are explained in a consistent and rigorous fashion in terms of the j-specific 5f electronic structure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A -type antiferromagnetic order in the Zintl-phase insulator EuZn 2 P 2

Zintl phases, containing strongly covalently bonded frameworks with separate ionically bonded ions, have emerged as a critical materials family in which to couple magnetism and strong spin-orbit coupling to drive diverse topological phases of matter. Here we report the single-crystal synthesis, magnetic, thermodynamic, transport, and theoretical properties of the Zintl compound EuZn 2 P 2 that crystallizes in the anti-La 2 O 3 (CaAl 2 Si 2 ) P-3m1 structure, containing triangular layers of Eu 2+ ions. In-plane resistivity measurements reveal insulating behavior with an estimated activation energy of E g = 0.11eV. Specific heat and magnetization measurements indicate antiferromagnetic ordering at T N = 23K. Curie-Weiss analysis of in-plane and out of plane magnetic susceptibility from T = 150 to 300 K yields p eff = 8.61 for μ 0 H⊥c and p eff = 7.74 for μ 0 H//c, close to the expected values for the 4f 7 J = S = 7/2 Eu 2+ ion and indicative of weak anisotropy. Below T N , a significant anisotropy of χ ⊥ /χ // ≈ 2.3 develops, consistent with A-type magnetic order as observed in isostructural analogs and as predicted by the density functional theory calculations reported herein. The positive Weiss temperatures of θ W =19.2K for μ 0 H⊥c and θ W =41.9K for μ 0 H//c show a similar anisotropy and suggest competing ferromagnetic and antiferromagnetic interactions. Comparing Eu magnetic ordering temperatures across trigonal EuM 2 X 2 (M= divalent metal, X= pnictide) shows that EuZn 2 P 2 exhibits the highest ordering temperature, with variations in T N correlating with changes in expected dipolar interaction strengths within and between layers and independent of the magnitude of electrical conductivity. These results provide experimental validation of the crystochemical intuition that the cation Eu 2+ layers and the anionic (M 2 X 2 ) 2– framework can be treated as electronically distinct subunits, enabling further predictive materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Charting C–C coupling pathways in electrochemical CO 2 reduction on Cu(111) using embedded correlated wavefunction theory

The electrochemical CO 2 reduction reaction (CO 2 RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C 2+ ) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO 2 RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C 2+ product formation via CO 2 RR electrocatalysts.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Bandgap Engineering of Ga 2 O 3 by MOCVD Through Alloying with Indium

Ga 2 O 3 and In 2 O 3 are vital semiconductors with current and future electronic device applications. Here, we study the alloying of In 2 O 3 and Ga 2 O 3 (IGO) and the associated changes in structure, morphology, band gap, and electrical transport properties. Undoped films of IGO were deposited on sapphire substrates with varying indium (In) percentage from zero to 100% by metal-organic chemical vapor deposition (MOCVD). Some films were annealed in H 2 to induce electrical conductivity. The measurements showed the optical band gap decreased by adding In; this was confirmed by density functional (DFT) calculations, which revealed that the nature of the valence band maximum and conduction band minimum strongly relate to the chemistry and that the band gap drops by adding In. The as-grown films were highly resistive except for pure In 2 O 3 , which possesses p-type conductivity, likely arising from In vacancy-related acceptor states. N-type conductivity was induced in all films after H-anneal. DFT calculations revealed that the presence of In decreases the electron effective mass, which is consistent with the electrical transport measurements that showed higher electron mobility for higher In percentage. The work revealed the successful band gap engineering of IGO and the modification of its band structure while maintaining high-quality films by MOCVD.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic and vibrational properties of bulk Cr 2 Ge 2 Te 6 from first-principles calculations

Here we use density functional theory based ab initio calculations to investigate the structural, vibrational, magnetic, and electronic properties of the layered ferromagnet Cr 2 Ge 2 Te 6 (CGT) that has attracted attention for potential spintronic applications. We optimize the structure for a fixed c/a ratio. Our results are in excellent agreement with experimental data on structure, phonons, and electronic properties. The use of van der Waals interactions and relativistic spin-orbit coupling yields accurate lattice constants and interlayer distances. Simultaneously, the adopted theoretical methods lead to a rigorous description of the vibrational normal modes, as well as the valence bands, that are in excellent agreement with Raman spectroscopy and angle-resolved photoemission spectroscopy spectra, respectively. The magnetic moment is slightly overestimated, and the magnetic anisotropy has the correct sign but is greater in magnitude than that in experiment. We find that inclusion of the on-site Coulomb repulsion on Cr−3d orbitals worsens agreement with experiment, in contrast to previous studies. Our results provide a path toward ab initio analysis of magnetic heterostructures based on CGT.

36 MATERIALS SCIENCE↗