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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A density functional theory model of mechanically activated silyl ester hydrolysis

To elucidate the mechanism of the mechanically activated dissociation of chemical bonds between carboxymethylated amylose (CMA) and silane functionalized silicon dioxide, we have investigated the dissociation kinetics of the bonds connecting CMA to silicon oxide surfaces with density functional calculations including the effects of force, solvent polarizability, and pH. We have determined the activation energies, the pre-exponential factors, and the reaction rate constants of candidate reactions. The weakest bond was found to be the silyl ester bond between the silicon and the alkoxy oxygen atom. Under acidic conditions, spontaneous proton addition occurs close to the silyl ester such that neutral reactions become insignificant. Upon proton addition at the most favored position, the activation energy for bond hydrolysis becomes 31 kJ mol{sup −1}, which agrees very well with experimental observation. Heterolytic bond scission in the protonated molecule has a much higher activation energy. The experimentally observed bi-exponential rupture kinetics can be explained by different side groups attached to the silicon atom of the silyl ester. The fact that different side groups lead to different dissociation kinetics provides an opportunity to deliberately modify and tune the kinetic parameters of mechanically activated bond dissociation of silyl esters.

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

D{sup 0} magnetism in Ca doped narrow carbon nanotubes: First principle chirality effect study

Curvature has always had crucial effects on the physical properties of narrow carbon nanotubes (CNTs) and here spin-polarized density functional calculations were employed to study electronic and magnetic properties of calcium-decorated narrow (5,5) and (9,0)CNTs with close diameters (∼7 Å) and different chiralities. Our results showed that chirality had great impact on the electronic structure and magnetization of the doped CNTs. In addition, internally or externally doping of the calcium atoms was studied comparatively and although for the (9,0)CNT the internal doping was the most stable configuration, which involves a novel kind of spin-polarization originated from Ca-4s electrons, but for the (5,5)tube the external doping was the most stable one without any spin-polarization. On the other hand, calcium doping in the center of the (5,5)CNT was an endothermic process and led to the spin-polarization of unoccupied Ca-3d orbitals via direct exchange interaction between adjacent Ca atoms. In the considered systems, the existence of magnetization in the absence of any transition-metal elements was an example of valuable d{sup 0} magnetism title.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

Thermodynamic ground state of MgB{sub 6} predicted from first principles structure search methods

Crystalline structures of magnesium hexaboride, MgB{sub 6}, were investigated using unbiased structure searching methods combined with first principles density functional calculations. An orthorhombic Cmcm structure was predicted as the thermodynamic ground state of MgB{sub 6}. The energy of the Cmcm structure is significantly lower than the theoretical MgB{sub 6} models previously considered based on a primitive cubic arrangement of boron octahedra. The Cmcm structure is stable against the decomposition to elemental magnesium and boron solids at atmospheric pressure and high pressures up to 18.3 GPa. A unique feature of the predicted Cmcm structure is that the boron atoms are clustered into two forms: localized B{sub 6} octahedra and extended B{sub ∞} ribbons. Within the boron ribbons, the electrons are delocalized and this leads to a metallic ground state with vanished electric dipoles. The present prediction is in contrast to the previous proposal that the crystalline MgB{sub 6} maintains a semiconducting state with permanent dipole moments. MgB{sub 6} is estimated to have much weaker electron-phonon coupling compared with that of MgB{sub 2}, and therefore it is not expected to be able to sustain superconductivity at high temperatures.

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

Effect of chlorination on the TlBr band edges for improved room temperature radiation detectors: Effect of chlorination on the TlBr band edges for radiation detectors

Thallium bromide (TlBr) crystals subjected to hydrochloric acid (HCl) chemical treatments have been shown to advantageously affect device performance and longevity in TlBr-based room temperature radiation detectors, yet the exact mechanisms of the improvements remain poorly understood. Here in this paper, we investigate the influence of several HCl chemical treatments on device-grade TlBr and describe the changes in the composition and electronic structure of the surface. Composition analysis and depth profiles obtained from secondary ion mass spectrometry (SIMS) identify the extent to which each HCl etch condition affects the detector surface region and forms of a graded TlBr/TlBr 1-x CL x surface heterojunction. Using a combination of X-ray photoemission spectroscopy (XPS) and hybrid density functional calculations, we are able to determine the valence band offsets, band gaps, and conduction band offsets as a function of Cl content over the entire composition range of TIBr 1-x C1 X . This study establishes a strong correlation between device process conditions, surface chemistry, and electronic structure with the goal of further optimizing the long-term stability and radiation response of TlBr-based detectors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies

We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke- Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Surface facet dependence of competing alloying mechanisms

Metal alloys are ubiquitous in many branches of heterogeneous catalysis, and it is now fairly well established that the local atomic structure of an alloy can have a profound influence on its chemical reactivity. While these effects can be difficult to probe in nanoparticle catalysts, model studies using well defined single crystal surfaces alloyed with dopants enable these structure–function correlations to be drawn. The first step in this approach involves understanding the alloying mechanism and the type of ensembles formed. In this study, we examined the atomic structure of RhCu single-atom alloys formed on Cu(111), Cu(100), and Cu(110) surfaces. Our results show a striking difference between Rh atoms alloying in Cu(111) vs the more open Cu(100) and Cu(110) surface facets. Unlike Cu(111) on which Rh atoms preferentially place-exchange with Cu atoms in the local regions above step edges leaving the majority of the Cu surface free of Rh, highly dispersed, homogeneous alloys are formed on the Cu(100) and (110) surfaces. These dramatically different alloying mechanisms are understood by quantifying the energetic barriers for atomic hopping, exchange, swapping, and vacancy filling events for Rh atoms on different Cu surfaces through theoretical calculations. Density functional theory results indicate that the observed differences in the alloying mechanism can be attributed to a faster hopping rate, relatively high atomic exchange barriers, and stronger binding of Rh atoms in the vicinity of step edges on Cu(111) compared to Cu(110) and Cu(100). Furthermore, these model systems will serve as useful platforms for examining structure sensitive chemistry on single-atom alloys.

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

Defects responsible for abnormal n-type conductivity in Ag-excess doped PbTe thermoelectrics

Density functional calculations have been performed to investigate the role of Ag defects in PbTe thermoelectric materials. Ag-defects can be either donor, acceptor, or isovalent neutral defect. When Ag is heavily doped in PbTe, the neutral (Ag-Ag) dimer defect at Pb-site is formed and the environment changes to the Pb-rich/Te-poor condition. Under Pb-rich condition, the ionized Ag-interstitial defect (Ag{sub I}{sup +}) becomes the major donor. The formation energy of Ag{sub I}{sup +} is smaller than other native and Ag-related defects. Also it is found that Ag{sub I}{sup +} is an effective dopant. There is no additional impurity state near the band gap and the conduction band minimum. The charge state of Ag{sub I}{sup +} defect is maintained even when the Fermi level is located above the conduction band minimum. The diffusion constant of Ag{sub I}{sup +} is calculated based on the temperature dependent Fermi level, formation energy, and migration energy. When T > 550 K, the diffusion length of Ag within a few minutes is comparable to the grain size of the polycrystalline PbTe, implying that Ag is dissolved into PbTe and this donor defect is distributed over the whole lattice in Ag-excess doped polycrystalline PbTe. The predicted solubility of Ag{sub I}{sup +} well explains the increased electron carrier concentration and electrical conductivity reported in Ag-excess doped polycrystalline PbTe at T = 450–750 K [Pei et al., Adv. Energy Mater. 1, 291 (2011)]. In addition, we suggest that this abnormal doping behavior is also found for Au-doped PbTe.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Transferability and accuracy by combining dispersionless density functional and incremental post-Hartree-Fock theories: Noble gases adsorption on coronene/graphene/graphite surfaces

The accuracy and transferability of the electronic structure approach combining dispersionless density functional theory (DFT) [K. Pernal et al., Phys. Rev. Lett. 103, 263201 (2009)] with the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)], are validated for the interaction between the noble-gas Ne, Ar, Kr, and Xe atoms and coronene/graphene/graphite surfaces. This approach uses the method of increments for surface cluster models to extract intermonomer dispersion-like (2- and 3-body) correlation terms at coupled cluster singles and doubles and perturbative triples level, while periodic dispersionless density functionals calculations are performed to estimate the sum of Hartree-Fock and intramonomer correlation contributions. Dispersion energy contributions are also obtained using DFT-based symmetry-adapted perturbation theory [SAPT(DFT)]. An analysis of the structure of the X/surface (X = Ne, Ar, Kr, and Xe) interaction energies shows the excellent transferability properties of the leading intermonomer correlation contributions across the sequence of noble-gas atoms, which are also discussed using the Drude oscillator model. We further compare these results with van der Waals-(vdW)-corrected DFT-based approaches. As a test of accuracy, the energies of the low-lying nuclear bound states supported by the laterally averaged X/graphite potentials (X = {sup 3}He, {sup 4}He, Ne, Ar, Kr, and Xe) are calculated and compared with the best estimations from experimental measurements and an atom-bond potential model using the ab initio-assisted fine-tuning of semiempirical parameters. The bound-state energies determined differ by less than 6–7 meV (6%) from the atom-bond potential model. The crucial importance of including incremental 3-body dispersion-type terms is clearly demonstrated, showing that the SAPT(DFT) approach effectively account for these terms. With the deviations from the best experimental-based estimations smaller than 2.3 meV (1.9%), the accuracy of the combined DFT and post-HF incremental scheme is established for all the noble-gas atoms. With relative deviations smaller than 4% and 11%, good agreement is also achieved by applying the vdW-corrected DFT treatments PBE-D3 and vdW-DF2 for noble-gas atoms heavier than neon.

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

Experimental and theoretical studies of the structural and electronic properties of vanadium–benzene sandwich clusters and their anions: V{sub n}Bz{sub n}{sup 0/−} (n = 1–5) and V{sub n}Bz{sub n−1}{sup 0/−} (n = 2–5)

One end open V{sub n}Bz{sub n}{sup −} (n = 1–5; Bz = benzene) and both ends open V{sub n}Bz{sub n−1}{sup −} (n = 2–5) vanadium–benzene cluster anions were studied using anion photoelectron spectroscopy and density functional calculations. The smaller (n ≤ 3) V{sub n}Bz{sub n} and V{sub n}Bz{sub n−1} clusters and corresponding anions were found to have structural isomers, whereas full-sandwiched V{sub n}Bz{sub n+1} clusters preferred to form multiple-decker sandwich structures. Several isomeric V{sub 2}Bz{sub 2} structures were identified theoretically and the anion photoelectron spectra of V{sub 2}Bz{sub 2}{sup 0/−} were explained well by the coexistence of two isomeric structures: (1) a V{sub 2}-core structure sandwiched between benzene molecules and (2) an alternating sandwich structure with the spin state strongly dependent on the structure. The adiabatic electron affinity of both V{sub n}Bz{sub n} and V{sub n}Bz{sub n−1} was found to increase with the cluster size at larger sizes (n = 4 or 5) and approaches to that of V{sub n}Bz{sub n+1}. The evolution of the structural and electronic properties of V{sub n}Bz{sub m} and V{sub n}Bz{sub m}{sup −} (m = n and n − 1) with size is discussed in comparison with V{sub n}Bz{sub n+1} and V{sub n}Bz{sub n+1}{sup −}.

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

Equation of state of boron carbide B 4 ⁢C

We present the results of recent experiments conducted on the Sandia Z machine and a new tabular equation of state for B 4 ⁢C. The equation of state was calibrated to a combination of density functional calculations reported here and fits to preexisting data. It was constructed partly to recover the effects of a shock-driven, polymorphic phase transition of unknown character beginning at particle velocities of just under 3 km/s (shock pressures of 95 GPa). Some of the Z experiments included sound speeds determined by the overtaking rarefaction method, from which we calculate the Grüneisen parameter and compare with previous experiments conducted at the OMEGA laser [Fratanduono et al ., Phys. Rev. B 94 , 184107 (2016)], our own first principles calculations, and another recent tabular equation of state [Zhang et al ., Phys. Rev. E 102 , 053203 (2020)]. We also compare our results with previous static compression, thermophysical, and melt studies, finding mixed consistency. We predict the onset and completion of shock melting at 225 and 265 GPa, respectively, and predict a melt curve that is largely flat to pressures of several hundred GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Tuning the carrier concentration to improve the thermoelectric performance of CuInTe{sub 2} compound

The electronic and transport properties of CuInTe{sub 2} chalcopyrite are investigated using density functional calculations combined with Boltzmann theory. The band gap predicted from hybrid functional is 0.92 eV, which agrees well with experimental data and leads to relatively larger Seebeck coefficient compared with those of narrow-gap thermoelectric materials. By fine tuning the carrier concentration, the electrical conductivity and power factor of the system can be significantly optimized. Together with the inherent low thermal conductivity, the ZT values of CuInTe{sub 2} compound can be enhanced to as high as 1.72 at 850 K, which is obviously larger than those measured experimentally and suggests there is still room to improve the thermoelectric performance of this chalcopyrite compound.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Li 8 MnO 6 : A Novel Cathode Material with Only Anionic Redox

In Li-excess transition-metal-oxide cathode materials, anionic oxygen redox can offer high capacity and high voltages, although peroxo and superoxo species may cause oxygen loss, poor cycling performance, and capacity fading. Previous work showed that undesirable formation of peroxide and superoxide bonds is controlled to some extent by Mn substitution, and the present work uses density functional calculations to examine the reasons for this by studying the anionic redox mechanism Li 8 MnO 6 . This material is obtained by substituting Mn for Sn in Li 8 SnO 6 or for Zr in Li 8 ZrO 6 , and we also compare to previous work on those materials. The calculations predict that Li 8 MnO 6 is stable at room temperature (with a band gap of 3.19 eV as calculated HSE06 and 1.82 eV as calculated with the less reliable with PBE+U), and they elucidate the chemical and structural effects involved in the inhibition of oxygen release in this cathode. Throughout the whole delithiation process, only O 2- ions are oxidized. The directional Mn-O bonds formed from unfilled 3d orbitals effectively inhibit the formation of O-O bonds, and the layered structure is maintained even after removing 3 Li per Li 8 MnO 6 formula unit. The calculated average voltage for removal of 3 Li is 3.69 V by HSE06, and the corresponding capacity is 389 mAh/g. The high voltage of oxygen anionic redox and the high capacity result in a high energy density of 1436 Wh/kg. The Li-ion diffusion barrier for the dominant interlayer diffusion path along the c-axis is 0.57 eV by PBE+U. Finally, these results help us to understand the oxygen redox mechanism in a new lithium-rich Li 8 MnO 6 cathode material and contribute to the design of high-energy-density lithium-ion-battery cathode materials with favorable electrochemical properties based on anionic oxygen redox.

Li8MnO6↗
DOE OSTI2016

First-principle optimal local pseudopotentials construction via optimized effective potential method

The local pseudopotential (LPP) is an important component of orbital-free density functional theory, a promising large-scale simulation method that can maintain information on a material’s electron state. The LPP is usually extracted from solid-state density functional theory calculations, thereby it is difficult to assess its transferability to cases involving very different chemical environments. Here, we reveal a fundamental relation between the first-principles norm-conserving pseudopotential (NCPP) and the LPP. On the basis of this relationship, we demonstrate that the LPP can be constructed optimally from the NCPP for a large number of elements using the optimized effective potential method. Specially, our method provides a unified scheme for constructing and assessing the LPP within the framework of first-principles pseudopotentials. Our practice reveals that the existence of a valid LPP with high transferability may strongly depend on the element.

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

Structural and electronic properties of Si{sub 1–x}Ge{sub x} alloy nanowires

We present first-principles density-functional calculations of Si{sub 1–x}Ge{sub x} alloy nanowires. We show that given the composition of the alloy, the structural properties of the nanowires can be predicted with great accuracy by means of Vegard's law, linearly interpolating the values of a pure Si and a pure Ge nanowire of the same diameter. The same holds, to some extent, also for electronic properties such as the band-gap. We also assess to what extend the band-gap varies as a function of disorder, i.e., how it changes for different random realization of a given concentration. These results make possible to tailor the desired properties of SiGe alloy nanowires starting directly from the data relative to the pristine wires.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Structural, electronic and magnetic properties of carbon doped boron nitride nanowire: Ab initio study

Using spin-polarized density functional theory calculations, we demonstrated that carbon doped boron nitride nanowire (C-doped BNNW) has diverse electronic and magnetic properties depending on position of carbon atoms and their percentages. Our results show that only when one carbon atom is situated on the edge of the nanowire, C-doped BNNW is transformed into half-metal. The calculated electronic structure of the C-doped BNNW suggests that doping carbon can induce localized edge states around the Fermi level, and the interaction among localized edge states leads to semiconductor to half-metal transition. Overall, the bond reconstruction causes of appearance of different electronic behavior such as semiconducting, half-metallicity, nonmagnetic metallic, and ferromagnetic metallic characters. The formation energy of the system shows that when a C atom is doped on surface boron site, system is more stable than the other positions of carbon impurity. Our calculations show that C-doped BNNW may offer unique opportunities for developing nanoscale spintronic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Origin of the magnetic and orbital ordering in α-Sr 2 CrO 4

Motivated by recent experimental progress in transition metal oxides with the K 2 NiF 4 structure, in this study we investigate the magnetic and orbital ordering in α-Sr 2 CrO 4 . Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the ab initio band structure near the Fermi level. The unique reverse splitting of t 2g orbitals in α-Sr 2 CrO 4 , with the 3d 2 electronic configuration for the Cr 4+ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground state phase diagram for the two dimensional compound α-Sr 2 CrO 4 . We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of α-Sr 2 CrO 4 , we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound α-Sr 2 CrO 4 , in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with d xy driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the d xz and d yz orbitals driving the orbital order.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Oxygen disorder, a way to accommodate large epitaxial strains in oxides

Density-functional calculations (total-energy comparisons) and checks for negative-frequency phonon modes are employed as a stability indicator to show that, in rutile- and fluorite-structure oxides, e.g., zirconia strained by a strontium titanate substrate, oxygen-sublattice disorder can be the energetically preferred way to accommodate strain.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Dry Reforming of Methane on a Highly‐Active Ni‐CeO 2 Catalyst: Effects of Metal‐Support Interactions on C−H Bond Breaking

Abstract Ni‐CeO 2 is a highly efficient, stable and non‐expensive catalyst for methane dry reforming at relative low temperatures (700 K). The active phase of the catalyst consists of small nanoparticles of nickel dispersed on partially reduced ceria. Experiments of ambient pressure XPS indicate that methane dissociates on Ni/CeO 2 at temperatures as low as 300 K, generating CH x and CO x species on the surface of the catalyst. Strong metal–support interactions activate Ni for the dissociation of methane. The results of density‐functional calculations show a drop in the effective barrier for methane activation from 0.9 eV on Ni(111) to only 0.15 eV on Ni/CeO 2− x (111). At 700 K, under methane dry reforming conditions, no signals for adsorbed CH x or C species are detected in the C 1s XPS region. The reforming of methane proceeds in a clean and efficient way.

Liu, Zongyuan↗