Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles density functional theory”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Polar Magnetic and Insulating Double Corundum Oxide: Mn2MnSbO6 with Ordered Mn(II) and Mn(III) Ions

A new magnetic insulator Mn2MnSbO6 with a polar crystal structure and an ordered Mn2+ and Mn3+ arrangement was synthesized under a high pressure of 7.5 GPa and 1300 degrees C. The crystal structure of Mn2MnSbO6, investigated by synchrotron powder X-ray diffraction, was found to be isomorphous with that of Ni3TeO6-type, space group R3. The non-centrosymmetric structure was confirmed by the second-harmonic generation measurements. The X-ray absorption near-edge spectroscopy measurement confirmed the nominal oxidation states of Mn22+Mn3+SbO6. Magnetic measurements indicate that Mn2MnSbO6 orders antiferromagnetically below 44 K and undergoes a field-induced spin-flop transition at 5 K. First-principles calculations indicate an antiferromagnetic ground state with up/up/up/down/down/down (uuuddd) spin configuration of the six crystallographically unique Mn ions in the c-axis doubled magnetic structure. The density functional theory calculations also substantiate the experimentally observed charge ordering of the Mn2+/Mn3+ ions and the insulating behavior due to a bandgap of 0.52 eV. To the best of our knowledge, this is the first double corundum oxide containing Jahn-Teller active Mn3+ ions.

Feng, Hai L.↗
DOE OSTI2021

Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model

Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multiscale simulations of electron and ion dynamics in self-irradiated silicon

The interaction of energetic ions with the electronic and ionic system of target materials is an interesting but challenging multiscale problem, and understanding of the early stages after impact of heavy, initially charged ions is particularly poor. At the same time, energy deposition during these early stages determines later formation of damage cascades. In this work, we address the multiscale character by combining real-time time-dependent density functional theory for electron dynamics with molecular dynamics simulations of damage cascades. Our first-principles simulations prove that core electrons affect electronic stopping and have an unexpected influence on the charge state of the projectile. We show that this effect is absent for light projectiles, but dominates the stopping physics for heavy projectiles. By parametrizing an inelastic energy loss friction term in the molecular dynamics simulations using our first-principles results, we also show a qualitative influence of electronic stopping physics on radiation-damage cascades.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

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

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

First-principles ionized-impurity scattering and charge transport in doped materials

Scattering of carriers with ionized impurities governs charge transport in doped semiconductors. However, electron interactions with ionized impurities cannot be fully described with quantitative first-principles calculations, so their understanding relies primarily on simplified models. Here we show an ab initio approach to compute the interactions between electrons and ionized impurities or other charged defects. It includes the short- and long-range electron-defect (e-d) interactions on equal footing and allows for efficient interpolation of the e-d matrix elements. Here we combine the e-d and electron-phonon interactions in the Boltzmann transport equation to compute the carrier mobilities in doped silicon over a wide range of temperature and doping concentrations, seamlessly spanning the defect- and phonon-limited transport regimes. The individual contributions of the defect- and phonon-scattering mechanisms to the carrier relaxation times and mean-free paths are analyzed. Our method provides a powerful tool to study electronic interactions in doped materials. It broadens the scope of first-principles transport calculations, enabling studies of a wide range of doped semiconductors and oxides with application to electronics, energy, and quantum technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles study of L -shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Catalyst Design in Nitrate Removal

Based on the volcano plot developed by Dr. Goldsmith group (Report linked in submission), we utilized DFT (density functional theory) calculations to search for bimetallic materials in the application of catalysts in aqueous nitrate removal. The calculations are conducted via the high-throughput automated workflow package developed by our group (Github linked in submission) using VASP commercial first-principles calculation software.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

A first-principles study of hydrogen surface coverage on δ-Pu (100), (111), and (110) surfaces

Hydriding corrosion of plutonium leads to surface cracking, pitting, and ultimately structural failure. Laboratory experiments demonstrate that hydriding begins on the surface or near the subsurface of plutonium. However, there has not yet been a systematic evaluation of hydrogen surface coverage on plutonium. In this work, we compute the surface energies of the low facet surfaces of face-centered cubic δ-Pu. The adsorption free energies of expected hydrogen structures at low and high coverage are presented along with the likely progression for filling sites as the H 2 partial pressure increases. Implications for near-equilibrium pressure hydride nucleation and non-equilibrium millibar pressure hydriding are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
DOE OSTI2021

Carbon complexes in highly C-doped GaN

We investigate the properties of heavily C-doped GaN grown by hydride vapor phase epitaxy using both optical experiments and hybrid density functional theory calculations. Previous work has established that carbon acceptors (C N ) give rise to a yellow luminescence band near 2.2 eV along with a blue luminescence band near 2.9 eV. Photoluminescence measurements show the yellow band shifting as a function of carbon concentration, suggesting a change in the behavior of carbon species as carbon content increases. With hybrid density functional theory we calculate the electrical and optical behavior of carbon centers containing multiple carbon impurities, which may arise in heavily doped material. We compare the behavior of these complexes to the isolated centers, and find that the dicarbon donor-acceptor (C Ga –C N ) complex is a candidate to explain the shift in the yellow luminescence peak. Tricarbon complexes have high formation energies and modest binding energies, and also give rise to optical transitions that are inconsistent with the observed spectra. We also identify the split dicarbon interstitial on the gallium site as a low-energy species with a large binding energy that may act to compensate carbon acceptors. Here, local vibrational modes are calculated for carbon impurity centers, and we compare these results to recent experiments. Dicarbon and tricarbon complexes involving C Ga and C N exhibit modes that are only slightly higher than the isolated species, while carbon interstitials and related complexes give rise to vibrational modes much higher than C Ga and C N .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dispersion interactions in proposed covalent superhydride superconductors

Recent developments in high temperature superconducting materials under high pressure have made numerical evaluation of the superconducting transition temperature (T c ) of predicted materials critically important as a means to identification. Existing methods of calculating T c often do not agree with each other or experiment, often due to the large number of complex factors that contribute to this property; among them is the neglect of dispersion interactions in commonly used density functionals. We evaluate the effect of including dispersion interactions on the predicted superconducting properties of two examples of the covalent superhydride class of very high-T c superconducting materials. In both cases, dispersion is found to have sizable effects, increasing the electron-phonon coupling as compared to the reference case of elemental niobium. A detailed investigation traces the origin of this effect in a 270 GPa $\mathrm {R3m CSH_7}$ structure to structural distortions driven by long-range electron-phonon interactions rather than novel bonding networks.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Electronic and total energy properties of ternary and quaternary semiconductor compounds, alloys, and superlattices: Theoretical study of Cu/graphite bonding

The goals of the research were to provide a fundamental science basis for why the bonding of Cu to graphite is weak, to critically evaluate the previous analysis of the wetting studies with particular regard to the values used for the surface energies of Cu and graphite, and to make recommendations for future experiments or other studies which could advance the understanding and solution of this technological problem. First principles electronic structure calculations were used to study the problem. These are based on density functional theory in the local density approximation and the use of the linear muffin-tin orbital band structure method. Calculations were performed for graphite monolayers, single crystal graphite with the hexagonal AB stacking, bulk Cu, Cu(111) surface, and Cu/graphite superlattices. The study is limited to the basal plane of graphite because this is the graphite plane exposed to Cu and graphite surface energies and combined with the measured contact angles to evaluate the experimental adhesion energy.

Lambrecht, Walter R. L.↗
DOE OSTI2021

Elucidating the local atomic and electronic structure of amorphous oxidized superconducting niobium films

Qubits made from superconducting materials are a mature platform for quantum information science application, such as quantum computing. However, material-based losses are now a limiting factor in reaching the coherence times needed for applications. In particular, knowledge of the atomistic structure and properties of the circuit materials is needed to identify, understand, and mitigate material-based decoherence channels. Here, we characterize the atomic structure of the native oxide film formed on Nb resonators by comparing fluctuation electron microscopy experiments to density functional theory calculations, finding that an amorphous layer is consistent with an Nb 2 O 5 stoichiometry. Comparing x-ray absorption measurements at the Oxygen K edge with first-principles calculations, we find evidence of d-type magnetic impurities in our sample, known to cause impedance in proximal superconductors. This work identifies the structural and chemical composition of the oxide layer grown on Nb superconductors and shows that soft x-ray absorption can fingerprint magnetic impurities in these superconducting systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Version 2.0.0 - M-SPARC: Matlab-Simulation Package for Ab-initio Real-space Calculations

M-SPARC is a Matlab code for performing ab initio Kohn–Sham Density Functional Theory simulations. Version 2.0.0 of the software further extends its capability to include relativistic effects, dispersion interactions, and advanced semilocal/nonlocal exchange–correlation functionals. These features significantly increase the fidelity of first principles calculations that can be performed using M-SPARC.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Control of polymorphism during epitaxial growth of hyperferroelectric candidate LiZnSb on GaSb (111)B

A major challenge for ferroelectric devices is the depolarization field, which competes with and often destroys long-range polar order in the limit of ultrathin films. Recent theoretical predictions suggest a new class of materials, termed hyperferroelectics, that should be robust against the depolarization field and enable ferroelectricity down to the monolayer limit. Here, the authors demonstrate the epitaxial growth of hexagonal LiZnSb, one of the hyperferroelectric candidate materials, by molecular-beam epitaxy on GaSb (111)B substrates. Due to the high volatility of all three atomic species, they find that LiZnSb can be grown in an adsorption-controlled window, using an excess zinc flux. Within this window, the desired polar hexagonal phase is stabilized with respect to a competing cubic polymorph, as revealed by x-ray diffraction and transmission electron microscopy measurements. First-principles calculations show that for moderate amounts of epitaxial strain and moderate concentrations of Li vacancies, the cubic LiZnSb phase is lower in formation energy than the hexagonal phase, but only by a few millielectronvolts per formula unit. Therefore, they suggest that kinetics plays a role in stabilizing the desired hexagonal phase at low temperatures. Their results provide a path toward experimentally demonstrating ferroelectricity and hyperferroelectricity in a new class of ternary intermetallic compounds.

42 ENGINEERING↗
DOE OSTI2022

Stabilization of Dinuclear Rhodium and Iridium Clusters on Layered Titanate and Niobate Supports

Atomically dispersed organometallic clusters can provide well-defined nuclearity of active sites for both fundamental studies as well as new regimes of activity and selectivity in chemical transformations. More recently, dinuclear clusters adsorbed onto solid surfaces have shown novel catalytic properties resulting from the synergistic effect of two metal centers to anchor different reactant species. Difficulty in synthesizing, stabilizing, and characterizing isolated atoms and clusters without agglomeration challenges allocating catalytic performance to atomic structure. Here, we explore the stability of dinuclear rhodium and iridium clusters adsorbed onto layered titanate and niobate supports using molecular precursors. Both systems maintain their nuclearity when Statistical analysis of HAADF-STEM images revealed that rhodium and iridium dimers had mean cluster-to-cluster distances very similar to what is expected from a random distribution of atoms over a large area, indicating that they are dispersed without aggregation. The stability of dinuclear rhodium clusters supported on titanate nanosheets was also investigated by X-ray absorption fine structure (EXAFS), DRIFTS, and first-principles calculations. Both X-ray absorption spectroscopy and HAADF-STEM simulations, guided by density functional theory (DFT)-optimized structure models, suggested that rhodium dimers adsorb onto the nanosheets in an end-on binding mode that is stable up to 100 degrees C under reducing conditions. Finally, this study highlights that crystalline nanosheets derived from layered metal oxides can be used as model supports to selectively stabilize dinuclear clusters, which could have implications for heterogeneous catalysis.

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

Energetics of polymeric carbon monoxide

The transformation of carbon monoxide (CO) from a molecular liquid to a polymeric solid under isothermal compression at room temperature is investigated using first principles theory. We report structural and thermodynamic properties from ambient density up to 2.45 g/cc obtained using density functional theory molecular dynamics simulations, including hybrid exchange corrections. The theoretical results are compared with newly obtained polymeric CO samples, synthesized in a large volume press. The explosive performance of polymeric CO is predicted and discussed. Under most favorable assumptions, it is found to be comparable to trinitrotoluene.

25 ENERGY STORAGE↗