Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Near-Infrared Emission from Sulfur Heteroatom Defects in Single-Walled Carbon Nanotubes

Quantum defects introduced into single-walled carbon nanotubes (SWCNTs) that support near-infrared emission are promising candidates for many applications. Here, we present a type of quantum defect in which sulfur atoms are incorporated into the crystal lattice of SWCNTs through their reactions with a mild reducing agent, sodium dithionite. When radical anions generated by sodium dithionites react with (6,5)-SWCNTs, two new emission features emerge at 1120 and 1250 nm. Comparing these spectral features with density functional theory calculation results, we attribute them to defects with episulfide or thioether structures. Here, by investigating the influences of reaction conditions including reactant concentration, light irradiation, and solvent type, we identify the reaction mechanism between SWCNTs and sodium dithionite and achieve good control over defect densities. Due to the intrinsic reactivity of sulfur species, the sulfur-containing quantum defects discovered here may facilitate the introduction of more complex and functionally diverse moieties into SWCNTs.

Anions↗
DOE OSTI2022

In-situ synthesized N-doped ZnO for enhanced CO 2 sensing: Experiments and DFT calculations

Chemiresistive CO 2 sensing is attractive due to low cost and ease of chip-level integration. Our previous studies (Yong Xia, 2021) showed the well-developed ZnO material fabricated by in-situ annealing exhibited good CO 2 sensing performance. Here, we have expanded on those studies, including CO 2 cyclic tests under both dry air and N 2 background whereby a much higher response to CO 2 in N 2 background was observed. In this study, detailed density functional theory calculations were conducted to understand the behavior. The results indicated nitrogen doping is mainly responsible for the observed response. In the presence of pre-adsorbed O 2 , N-doped ZnO can no longer interact with CO 2 , which agrees well with the observation of higher response in N 2 background. Furthermore, density of states analysis showed N sp 2 hybridized orbital and N 2p orbital of the N dopant mixed with sp 2 hybridized orbital of C atom and 2p orbitals of C/O atoms in CO 2 to form σ and π bonds, respectively. However, they mixed with O 2s/2p orbitals of O atom in O 2 when pre-adsorbed O 2 was present, hindering CO 2 interaction with N-doped ZnO, and resulting in limited response in air. The illustrated mechanism does not only further the understanding of metal oxide-based CO 2 sensing, but also guide the design of new functional materials for CO 2 sensing or capture.

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

Calculations of CO Oxidation over a Au/TiO 2 Catalyst: A Study of Active Sites, Catalyst Deactivation, and Moisture Effects

The reaction mechanism of CO oxidation on Au/TiO 2 catalysts remains elusive. Here, we employ density functional theory calculations to gain an understanding of several important aspects of the system, including CO adsorption, the active oxygen species, catalyst deactivation, and the promoting effect of moisture on catalytic activity. Distinct from previous theoretical studies, which tend to address these questions individually, here we construct a model of the catalytic system which can address all of the issues mentioned. For this, we have considered the complex interactions among reactants, products, and catalysts under reaction conditions. The main findings of our present study are (1) the Au/TiO 2 interface boundary can be easily oxidized, (2) CO adsorption on oxidized Au results in the formation of O–Au–CO species, (3) surface lattice oxygen on a TiO 2 support is the active oxygen species, (4) CO 2 binds strongly on the O top /Ti 5c site, forming carbonate that blocks the active site, and finally (5) water can accelerate O 2 dissociation and carbonate decomposition. The results of our theoretical model are compared with existing experimental observations and found to be largely consistent with them.

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

Ab Initio Calculations of New Full Heusler Alloys Rh{sub 2}ZrX (X = Al, Ga, In, Si, Ge, Sn)

In recent years, Heusler alloys have attracted extensive attention because of their usefulness in spintronic devices. In this paper, the structural, electronic, and magnetic properties of full Heusler alloys Rh{sub 2}ZrX (X = Al, Ga, In, Si, Ge, Sn) have been studied using Quantum Espresso software package based on density functional theory. Calculations have been done by GGA and GGA+U methods. The investigations on the structural properties of the alloys showed that the Rh{sub 2}ZrSi alloy had the most stable structure (E{sub coh} = − 1.807 eV) among the six compounds reviewed. The results of GGA calculation on electronic properties of Rh{sub 2}ZrX alloys show that while the Rh{sub 2}ZrX (X = Al, Ga, In) alloys had metallic nature due to overlapping between conduction and valance band at Fermi level for both spin directions, the Rh{sub 2}ZrX (X = Si, Ge, Sn) alloys had half-metal characteristic. Also, Rh{sub 2}ZrSi and Rh{sub 2}ZrGe alloys had 100% spin polarization. By studying the magnetic properties, it was found that the Rh{sub 2}ZrAl and Rh{sub 2}ZrSi alloys were antiferromagnetic, while the Rh{sub 2}ZrGa, Rh{sub 2}ZrIn, Rh{sub 2}ZrGe, and Rh{sub 2}ZrSn alloys were ferromagnetic materials. Also, the total spin magnetic moment of alloys was found in accordance with the Slater-Pauling rule. The results obtained using GGA+U calculation are slightly different from the results of GGA calculation.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Octahedral tilt transitions in the relaxor ferroelectric Na{sub 1/2}Bi{sub 1/2}TiO{sub 3}

The kinetics of octahedral tilt transitions in the lead-free relaxor material sodium bismuth titanate Na{sub 1/2}Bi{sub 1/2}TiO{sub 3} (NBT) is investigated by electronic structure calculations within density functional theory. Energy barriers for transitions between tetragonal, rhombohedral and orthorhombic tilts in cation configurations with [001]- and [111]-order on the A-sites are determined by nudged elastic band calculations. By tilting entire layers of octahedra simultaneously we find that the activation energy is lower for structures with 001-order compared to such with 111-order. The energetic coupling between differently tilted layers is, however, negligibly small. By introducing a single octahedral defect we create local tilt disorder and find that the deformation energy of the neighboring octahedra is less in a rhombohedral than in a tetragonal structure. By successively increasing the size of clusters of orthorhombic defects in a rhombohedral matrix with 001-order, we determine a critical cluster size of about 40 Å . Thus groups of about ten octahedra can be considered as nuclei for polar nanoregions, which are the cause of the experimentally observed relaxor behavior of NBT. - Graphical abstract: Nine orthorhombic oxygen octahedral tilt defects in a rhombohedral tilt configuration. - Highlights: • Chemical order influences energy barriers of octahedral tilt transitions. • The octahedral deformation energy is lower in rhombohedral phases. • Tilt defect clusters are more likely in rhombohedral structures. • Tilt defect clusters can act as nuclei for polar nanoregions.

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

Using first-principles calculations to screen for fragile magnetism: Case study of LaCrGe 3 and LaCrSb 3

In this study, we present a coupled experimental/theoretical investigation of pressure effect on the ferromagnetism of LaCrGe 3 and LaCrSb 3 compounds. The magnetic, electronic, elastic, and mechanical properties of LaCrGe 3 and LaCrSb 3 at ambient condition are studied by first-principles density-functional theory calculations. The pressure dependences of the magnetic properties of LaCrGe 3 and LaCrSb 3 are also investigated. The ferromagnetism in LaCrGe 3 is rather fragile, with a ferro- to paramagnetic transition at a relatively small pressure (around 7 GPa from our calculations, and 2 GPa in experiments). The key parameter controlling the magnetic properties of LaCrGe 3 is found to be the proximity of the peak of Cr density of states to the Fermi level, a proximity that is strongly correlated with the distance between Cr atoms along the c axis, suggesting that there would be a simple way to suppress magnetism in systems with one-dimensional arrangement of magnetic atoms. By contrast, the ferromagnetism in LaCrSb 3 is not fragile. In conclusion, our theoretical results are consistent with our experimental results and demonstrate the feasibility of using first-principles calculations to aid experimental explorations in screening for materials with fragile magnetism.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Massively-parallel electron dynamics calculations in real-time and real-space: Toward applications to nanostructures of more than ten-nanometers in size

A highly efficient program of massively parallel calculations for electron dynamics has been developed in an effort to apply the method to optical response of nanostructures of more than ten-nanometers in size. The approach is based on time-dependent density functional theory calculations in real-time and real-space. The computational code is implemented by using simple algorithms with a finite-difference method in space derivative and Taylor expansion in time-propagation. Since the computational program is free from the algorithms of eigenvalue problems and fast-Fourier-transformation, which are usually implemented in conventional quantum chemistry or band structure calculations, it is highly suitable for massively parallel calculations. Benchmark calculations using the K computer at RIKEN demonstrate that the parallel efficiency of the program is very high on more than 60 000 CPU cores. The method is applied to optical response of arrays of C{sub 60} orderly nanostructures of more than 10 nm in size. The computed absorption spectrum is in good agreement with the experimental observation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Insights into the Stability of Zeolitic Imidazolate Frameworks in Humid Acidic Environments from First-Principles Calculations

Understanding the degradation mechanisms of zeolitic imidazolate frameworks (ZIFs) is crucial to improve their chemical stability and realize their potential industrial applications. In this work, we conduct a comprehensive study utilizing dispersion-corrected density functional theory calculations to investigate the chemical stability of bulk ZIFs and their external surfaces under conditions of acid-gas exposure. We examine the influence of steric factors such as topology and ligand functionalization on the relative chemical stability of prototypical ZIFs (ZIF-2 and ZIF-8), including their hypothetical polymorphs. We find that defect formation is more thermodynamically and kinetically favorable at ZIF external surfaces versus the bulk, and that both topology and ligand functionalization impact defect formation. In addition, we provide a detailed mechanism for the reaction of ZIFs with sulfurous and sulfuric acids, of which the latter serves as a catalyst in potential degradation reactions of ZIFs. We also provide information about the adsorption strength of a range of acid gases to defective ZIF structures, which can inform potential strategies to regenerate ZIFs and/or achieve defect engineering in these materials.

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

Native defects in antiferromagnetic topological insulator MnBi 2 Te 4

Using scanning tunneling microscopy and spectroscopy, in this study we visualized the native defects in antiferromagnetic topological insulator MnBi 2 Te 4 . Two native defects, Mn Bi and Bi Te antisites, can be well resolved in the topographic images. Mn Bi tend to suppress the density of states at the conduction band edge. Spectroscopy imaging reveals a localized peaklike local density of state at ~ 80 meV below the Fermi energy. A careful inspection of topographic and spectroscopic images, combined with density functional theory calculation, suggests this results from Bi Mn antisites at Mn sites. The random distribution of Mn Bi and Bi Mn antisites results in spatial fluctuation of local density of states near the Fermi level in MnBi 2 Te 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

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

Single pair of Weyl nodes in the spin-canted structure of EuCd 2 As 2

We report time reversal symmetry breaking Weyl semimetals are unique among Weyl materials in allowing the minimal number of Weyl points, thus offering the clearest signatures of the associated physics. Here we present neutron diffraction, density-functional theory, and transport measurement results which indicate that EuCd 2 As 2 , under ambient field, strain, and pressure, is such a material with a single pair of Weyl points. Our work reveals a magnetic structure (magnetic space group C2'/m') with Eu moments pointing along the [210] direction in plane and canted ~ 30° out of plane. Density functional theory calculations using this structure show that the observed canting drastically alters the relevant electronic bands, relative to the in-plane order, leading to a single set of well-defined Weyl points. Furthermore, we find the canting angle can tune the distance of the Weyl points above the Fermi level, with the smallest distance at low canting angles. Finally, transport measurements of the anomalous Hall effect and longitudinal magnetoresistance exhibit properties indicative of a chiral anomaly, thus supporting the neutron scattering and DFT results suggesting EuCd 2 As 2 is close to the ideal situation of the Weyl hydrogen atom.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Electron Thermal Conductivity of Pu and Zr Substituted Gamma-Uranium

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (k) and fissile density; however, the effects of alloying elements on k remain unclear. Here, the electron thermal conductivity (k_e) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the k_e of ?-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate k_e; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the method that applies the electron relaxation time of ?-U is best, while if the elements are dissimilar, a mixed method that mixes several parameters associated with k_e from each element in the alloy is best. The introduction of all alloying elements decreases k_e; however, in binary compounds, Pu and Zr have different effects. Pu generally flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases the k_e more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel with without sacrificing k. We also note that these k_e calculation methods can be applied to non-fuel alloys that require k_e predictions, such as cladding, heat exchanger, and structural materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

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

Mirror symmetry breaking and lateral stacking shifts in twisted trilayer graphene

Here, we construct a continuum model of twisted trilayer graphene using ab initio density-functional-theory calculations, and we apply it to address twisted trilayer electronic structure. Our model accounts for moiré variation in site energies, hopping between outside layers, and hopping within layers. We focus on the role of a mirror symmetry present in ABA graphene trilayers with a middle layer twist. The mirror symmetry is lost intentionally when a displacement field is applied between layers, and unintentionally when the top layer is shifted laterally relative to the bottom layer. We calculate two band-structure characteristics that are directly relevant to transport measurements, namely the Drude weight and the weak-field Hall conductivity. We relate them via the Hall density to assess the influence of the accidental lateral stacking shifts currently present in all experimental devices on electronic properties. Finally, we comment on the role of the possible importance of accidental lateral stacking shifts for superconductivity in twisted trilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Regulating Top-Surface Multilayer/Single-Crystal Graphene Growth by “Gettering” Carbon Diffusion at Backside of the Copper Foil

A unique strategy is reported to constrain the nucleation centers for multilayer graphene (MLG) and, later, single-crystal graphene domains by gettering carbon source on backside of the flat Cu foil, during chemical vapor deposition. Hitherto, for a flat Cu foil, the top-surface-based growth mechanism is emphasized, while overlooking the graphene on the backside. However, the systematic experimental findings indicate a strong correlation between the backside graphene and the nucleation centers on the top-surface, governed by the carbon diffusion through the bulk Cu. This understanding steers to devise a strategy to mitigate the carbon diffusion to the top-surface by using a carbon “getter” substrate, such as nickel, on the backside of the Cu foil. Depth profiling of the nickel substrate, along with the density functional theory calculations, verifies the gettering role of the nickel support. The implementation of the backside carbon gettering approach on single-crystal graphene growth results in lowering the nucleation density by two orders of magnitude. This enables the single-crystal domains to grow by 6 mm laterally on the untreated Cu foil. Finally, the growth of large-area polycrystalline single layer graphene, free of unwanted MLG domains, with significantly improved field-effect mobility of ≈6800 cm 2 V -1 s -1 is demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Structure and Dynamics of Polysulfide Clusters in a Nonaqueous Solvent Mixture of 1,3-Dioxolane and 1,2-Dimethoxyethane

Molecular clustering and associated dynamic processes of lithium polysulfide species were unraveled using classical molecular dynamics and ab initio metadynamics calculations. The spectroscopic signatures of polysulfide clusters were analyzed using a multimodal analysis including experimental and computational nuclear magnetic resonance (NMR) and X-ray absorption spectroscopies. Lithium polysulfide solutes (Li 2 S 4 , Li 2 S 6 , and Li 2 S 8 ) and their mixtures in a 1,3-dioxolane and 1,2-dimethoxyethane (DOL/DME) solvent undergo aggregation driven by intramolecular lithium–sulfur (Li–S) interactions, leading to distributions of cluster sizes, which could critically influence the functioning of lithium–sulfur batteries. Representative polysulfide clusters with systematic increases in molecular size were extracted from the classical molecular dynamics trajectories for subsequent structural and spectroscopic property calculations using density functional theory analysis. Structural analysis of these clusters reveals progressively decreasing solvent involvement in Li+ coordination varying from Li 2 S 4 to Li 2 S 8 , with more pronounced variation and changes in DME compared with those in DOL. These observations are reflected in the analysis of the experimental and theoretical 7 Li and 17 O NMR chemical shifts and pulsed field gradient-NMR diffusion measurements. A comparison of experimental and theoretical S K-edge X-ray absorption near edge spectra shows that relatively large lithium sulfide chain clusters are likely to occur in the DOL/DME-solvated lithium sulfide systems. Ab initio metadynamics simulations and NMR analysis indicate that Li + solvated by only the solvent can occur through Li + dissociation from sulfide chains. However, the occurrence of “sulfide-free” Li + is a minor mechanism compared with the dynamic aggregation and shuttling processes of polysulfide solvates in DOL/DME-based electrolytes of Li–S batteries. Overall, atomistic insights gained about clustering and lithium exchange dynamics will be critical for the predictive understanding of the polysulfide shuttling and nucleation process that dictates the Li–S battery performance.

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

Nanoscale modulation of flat bands via controllable charge density wave defects in 4⁢𝐻⁡𝑏−Ta⁢S 2

Electron correlation is a main driver of exotic quantum phases and their interplay. The 4H b-TaS 2 system, possessing an intrinsic heterostructure of 1T- and 1H -TaS 2 monolayers, offers a unique opportunity to control electron correlation by distorting the atomic lattice or tuning interlayer coupling. Here, we investigated intrinsically deformed charge-density waves (CDWs) in the 1T layer of 4H b-TaS 2 to elucidate and control their effects on flat bands using scanning tunneling microscopy and spectroscopy (STM/S) combined with first-principles calculations. We identified two types of CDW defects: Type 1 has structural distortion and locally suppressed flat bands, while Type 2 features an increased flat band filling factor of intact CDW structure. Density functional theory calculations indicate that a sulfur vacancy in the 1T layer distorts the CDW structure and gives rise to a Type 1, whereas a sulfur vacancy in the 1H layer reduces the interlayer charge transfer and leads to a Type 2. Furthermore, we demonstrated creating and erasing individual CDW defects via STM manipulation. Here, our findings provide a pathway to not only tune flat bands but also selectively manipulate the interaction between CDW, the atomic lattice, and interlayer coupling in strongly correlated systems with atomic precision.

Charge density waves↗
DOE OSTI2021

Strengths of plasmon satellites in XPS: Real-time cumulant approach

Plasmon satellites observed in x-ray photoemission spectra (XPS) reflect many-body effects beyond the single-electron approximation. However, theoretical calculations are challenging and simplified theories based on electron gas models are not always satisfactory. Here, a cumulant Green’s function approach with a real-time time-dependent density functional theory calculation of the cumulant is described for calculations of these effects, without the need for adjustable parameters. Here, the approach yields good agreement with experimental XPS results, as illustrated for Ge and Al. Corrections due to extrinsic losses and interference are also briefly discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗