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 703 records · Page 39

Mechanistic insight into the passive NOx adsorption in the highly dispersed Pd/HBEA zeolite

Passive NOx adsorption using zeolite supported Pd catalysts has been considered one of most effective methods for low-temperature NO storage in diesel engine emission control. In the present work, first-principles density functional theory calculations were used to study NO adsorption at the highly dispersed Pd sites, i.e., monomeric PdIIO and dimeric [PdIIOPdII]2+ in the H-BEA zeolite. It has been found that each PdIIO and [PdIIOPdII]2+ sites could bind with three and four NO molecules, respectively. Upon adsorption, NO could be oxidized to NO2, resulting in the reduction of PdII to Pd0 for PdIIO and PdII to PdI, for [PdIIOPdII]2+, respectively. With increasing NO coverage at both sites, NO oxidation becomes more facile with low activation barriers of 13 and 22 kJ/mol. At the PdIIOPdII site, the binding of formed NO2 becomes weak at high NO coverage. Both PdIIO and [PdIIOPdII]2+ sites are protonated to the more stable [PdIIOH]+ and [PdIIOH]+/[PdIIOH]+ sites in the presence of H2O or neighboring BAS sites. Finally, we studied the CO effect on the NO adsorption in Pd/H-BEA. With co-fed CO in the NOx mixture, both PdIIO and [PdIIOPdII]2+ could be reduced to Pd0 and [PdI-PdI]2+ sites. Compared to the NO adsorption, our calculations show that CO adsorption is slightly stronger at the [PdIIOPdII]2+ site while it is slightly weaker at the PdIIO site. CO oxidation, rather than NO oxidation readily occurs if NO and CO co-adsorbed. The reduced PdI site and the [PdIIOH]+/[PdIIOH]+ site enhances NO adsorption.

passive NOx adsorption, density functional theory,↗

Mononuclear iridium dinitrogen complexes bonded to zeolite HY

In this study, the adsorption of N 2 on structurally well-defined dealuminated HY zeolite-supported iridium diethylene complexes was investigated. Iridium dinitrogen complexes formed when the sample was exposed to N 2 in H 2 at 298 K, as shown by infrared spectra recorded with isotopically labeled N 2 . Four supported species formed in various flowing gases: Ir(N 2 ), Ir(N 2 )(N 2 ), Ir(C 2 H 5 )(N 2 ), and Ir(H)(N 2 ). Their interconversions are summarized in a reaction network, showing, for example, that, in the presence of N 2 , Ir(N 2 ) was the predominant dinitrogen species at temperatures of 273-373 K. Ir(CO)(N 2 ) formed transiently in flowing CO, and in the presence of H 2 , rather stable iridium hydride complexes formed. Here, four structural models of each iridium complex bonded at the acidic sites of the zeolite were employed in a computational investigation, showing that the calculated vibrational frequencies agree well with experiment when full calculations are done at the level of density functional theory, independent of the size of the model of the zeolite.

36 MATERIALS SCIENCE↗

Properties of Ti/TiC Interfaces from Molecular Dynamics Simulations

Titanium carbide is used as a primary component in coating materials, thin films for electronic devices, and composites. Here, the structure of coherent and semicoherent interfaces formed between close-packed TiC (111) and Ti (0001) is investigated in classical molecular dynamics simulations. The forces on the atoms in the simulations are determined using a newly developed TiC potential under the framework of the third-generation charge optimized many-body (COMB 3 ) suite of potentials. The work of adhesion energies for the coherent interfaces is calculated and compared with the predictions of density functional theory calculations. In the case of relaxed semicoherent interfaces, a two-dimensional (2D) misfit dislocation network is predicted to form that separates the interface into different regions in which the positions of the atoms are similar to the positions at the corresponding coherent interfaces. After the interface is annealed at an elevated temperature, the climb of edge dislocations is activated which modifies the 2D misfit dislocation network and increases the work of adhesion. These findings can be used as inputs for sequential larger simulation models to understand and predict the macroscopic properties of TiC/Ti interfaces.

36 MATERIALS SCIENCE↗

Selective-Area Remote Epitaxy of ZnO Microrods Using Multilayer–Monolayer-Patterned Graphene for Transferable and Flexible Device Fabrications

Here, selective-area remote epitaxy (SA-REpi) is demonstrated for fabricating mechanically releasable position-controlled ZnO microrod (MR) arrays from donor wafers in an arrayed form. Intaglio-patterned graphene, consisting of basal single-layer graphene (SLG) overlayered with multilayer graphene (MLG) patterned with perforated holes, is transferred onto a GaN/Al 2 O 3 wafer on which the hydrothermal synthesis is performed for growing ZnO MRs. The basal SLG area exposed through the MLG pattern yields ZnO MRs, whereas the MLG plateau inhibits the growth. The noncovalent remote epitaxial heterointerface enables the release of the MR overlayer in the arrayed form, and the original source wafer is refurbished for reproducibly repeating the SA-REpi. Density-functional theory calculations suggest that localized surface charge density is induced on the surface of SLG by the underlying GaN across ultrathin SLG, which possibly provides a driving force for precursor adatoms and the following remote epitaxy of ZnO. In contrast, the induction of the charge density redistribution does not clearly occur through MLG; so, that keeps the surface of MLG nearly charge-neutral. The diameter and spacing of ZnO MRs are controlled in a designed way by changing the pattern geometries. High-resolution scanning transmission electron microscopy reveals the remote heteroepitaxial relationship at an atomic level. The remote epitaxy is expected to provide an ideal platform to transfer the addressable spatial arrays of nano- or micro-architecture semiconductor components to arbitrary target surfaces directly after the growth without the assembly procedures.

36 MATERIALS SCIENCE↗

Pressure-induced creation and annihilation of Weyl points in T d - Mo 0.5 W 0.5 Te 2 and 1 T " - Mo 0.5 W 0.5 Te 2

By means of first-principles density-functional theory calculations, we investigate the role of hydrostatic pressure in the electronic structure of the T d (Pmn2 1 ) and 1T" (Pm) phases of the Weyl semimetal Mo 0.5 W 0.5 Te 2 , which is a promising material for phase-change memory technology and superconductivity. We particularly focus on changes occurring in the distribution of the gapless Weyl points (WPs) within the 0 to 45 GPa pressure range. Here we further investigate the structural phase transition and lattice dynamics of the T d and 1T" phases within the aforementioned pressure range. Our calculations suggest that both the T d and 1T" phases of Mo 0.5 W 0.5 Te 2 host four WPs in their full Brillouin zone at zero pressure. The total number of WPs increases to 44 (36) with increasing pressure via pair creation up to 20 (15) GPa for the T d (1T") phase, and beyond this pressure pair annihilation of WPs starts occurring, leaving only 16 WPs at 45 GPa in both phases. The enthalpy versus pressure data reveal that the 1T" phase is more favorable below the critical pressure of 7.5 GPa; however, beyond this critical pressure the T d phase becomes enthalpically favorable. We also provide the calculated x-ray diffraction spectra along with the calculated Raman- and infrared-active phonon frequencies to facilitate the experimental identification of the studied phases

36 MATERIALS SCIENCE↗

CO{sub 2} binding in the (quinoline-CO{sub 2}){sup −} anionic complex

We have studied the (quinoline-CO{sub 2}){sup −} anionic complex by a combination of mass spectrometry, anion photoelectron spectroscopy, and density functional theory calculations. The (quinoline-CO{sub 2}){sup −} anionic complex has much in common with previously studied (N-heterocycle-CO{sub 2}){sup −} anionic complexes both in terms of geometric structure and covalent bonding character. Unlike the previously studied N-heterocycles, however, quinoline has a positive electron affinity, and this provided a pathway for determining the binding energy of CO{sub 2} in the (quinoline-CO{sub 2}){sup −} anionic complex. From the theoretical calculations, we found CO{sub 2} to be bound within the (quinoline-CO{sub 2}){sup −} anionic complex by 0.6 eV. We also showed that the excess electron is delocalized over the entire molecular framework. It is likely that the CO{sub 2} binding energies and excess electron delocalization profiles of the previously studied (N-heterocycle-CO{sub 2}){sup −} anionic complexes are quite similar to that of the (quinoline-CO{sub 2}){sup −} anionic complex. This class of complexes may have a role to play in CO{sub 2} activation and/or sequestration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-interaction-free electric dipole polarizabilities for atoms and their ions using the Fermi-Löwdin self-interaction correction

The static electric dipole polarizability of a system is a measure of the binding of its electrons. In density functional theory calculations, this binding is weakened by the presence of unphysical self-interaction in the density functional approximation (DFA), leading to overestimates of polarizabilities. In this work, to investigate this systematically we compare polarizabilities for the atoms from H to Ar and their anions and cations calculated in several DFAs and the corresponding self-interaction-corrected (SIC) DFAs with experiment and with high-level quantum chemistry reference values. The SIC results are obtained using the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) method. Removing self-interaction generally leads to smaller polarizabilities that agree significantly better with reference values. In conclusion, we find that FLO-SIC improves the performance of the local spin density approximation and the generalized gradient approximation (GGA) for polarizabilities to a quality that is comparable to so-called rung 4 functionals, but slightly degrades the performance of the strongly constrained and appropriately normed meta-GGA functional.

74 ATOMIC AND MOLECULAR PHYSICS↗

Ab initio molecular dynamics with nuclear quantum effects at classical cost: Ring polymer contraction for density functional theory

Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ring polymer contraction scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, this approach is as fast as a typical classical ab initio molecular dynamics simulation and 35× faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure and energetics of nanotwins in cubic boron nitrides

Recently, nanotwinned cubic boron nitrides (NT c-BN) have demonstrated extraordinary leaps in hardness. However, an understanding of the underlying mechanisms that enable nanotwins to give orders of magnitude increases in material hardness is still lacking. Here, using transmission electron microscopy, we report that the defect density of twin boundaries depends on nanotwin thickness, becoming defect-free, and hence more stable, as it decreases below 5 nm. Using ab initio density functional theory calculations, we reveal that the Shockley partials, which may dominate plastic deformation in c-BNs, show a high energetic barrier. We also report that the c-BN twin boundary has an asymmetrically charged electronic structure that would resist migration of the twin boundary under stress. These results provide important insight into possible nanotwin hardening mechanisms in c-BN, as well as how to design these nanostructured materials to reach their full potential in hardness and strength.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sooting Tendencies of Aromatic Hydrocarbons with Oxygen-Containing Side-Chains

Yield Sooting Index (YSI) was determined for 13 oxygen-containing aromatic compounds by measuring line-of-sight spectral radiance (LSSR) in a laminar non-premixed methane/air flame doped with 1000 ppm of each compound. The compounds included anisole (methoxybenzene), other aromatic ethers, and phenols. The presence of oxygen greatly reduced soot formation in some cases, but increased it in others. Thus oxygenated aromatics are potentially valuable as fuel components with intrinsically low particulate emissions, but the specific structure is important. Compounds containing an aromatic ring with an adjacent oxygen atom tended to have a lower sooting tendency than did compounds with oxygen atoms elsewhere in the side chain. This observation was quantitatively explored using Density Functional Theory (DFT) calculations. The bond dissociation energies for nine aromatic hydrocarbons and oxygenates were calculated using GGA B3LYP density functions with 6-31G(d) basis sets. When possible, unimolecular dissociation to phenoxy radical was favored and resulted in much lower YSIs than similar compounds where the oxygen atom was not adjacent to the aromatic ring. Structures that dissociated to benzyl radical had much higher sooting tendencies, sometimes by over a factor of two in comparison to phenoxyforming compounds.

09 BIOMASS FUELS↗

Physics of grain boundaries in polycrystalline photovoltaic semiconductors

Thin-film solar cells based on polycrystalline Cu(In,Ga)Se{sub 2} (CIGS) and CdTe photovoltaic semiconductors have reached remarkable laboratory efficiencies. It is surprising that these thin-film polycrystalline solar cells can reach such high efficiencies despite containing a high density of grain boundaries (GBs), which would seem likely to be nonradiative recombination centers for photo-generated carriers. In this paper, we review our atomistic theoretical understanding of the physics of grain boundaries in CIGS and CdTe absorbers. We show that intrinsic GBs with dislocation cores exhibit deep gap states in both CIGS and CdTe. However, in each solar cell device, the GBs can be chemically modified to improve their photovoltaic properties. In CIGS cells, GBs are found to be Cu-rich and contain O impurities. Density-functional theory calculations reveal that such chemical changes within GBs can remove most of the unwanted gap states. In CdTe cells, GBs are found to contain a high concentration of Cl atoms. Cl atoms donate electrons, creating n-type GBs between p-type CdTe grains, forming local p-n-p junctions along GBs. This leads to enhanced current collections. Therefore, chemical modification of GBs allows for high efficiency polycrystalline CIGS and CdTe thin-film solar cells.

36 MATERIALS SCIENCE↗

Optical properties of Eu{sup 2+}/Eu{sup 3+} mixed valence, silicon nitride based materials

Eu{sub 2}SiN{sub 3}, a mixed valence europium nitridosilicate, has been prepared via solid-state reaction synthesis and its oxidation behavior and optical properties have been determined. Furthermore, the stability of several isostructural compounds of the type M{sup 2+}L{sup 3+}SiN{sub 3} has been predicted by using the density functional theory calculations, and verified by the actual synthesis of CaLaSiN{sub 3}, CaEuSiN{sub 3} and EuLaSiN{sub 3}. The band gap of CaLaSiN{sub 3} was found around 3.2 eV giving the material its yellow color. Eu{sub 2}SiN{sub 3} on the other hand is black due to a combination of the 4f–5d absorption band of Eu{sup 2+} and the charge transfer band of Eu{sup 3+}. Thermogravimetric analysis and Raman spectroscopic study of Eu{sub 2}SiN{sub 3} revealed that oxidation of this compound in dry air takes place via a nitrogen retention complex. - Graphical abstract: Energy level scheme of Eu{sub 2}SiN{sub 3} showing the occupied N{sup 3−} 2p band (blue rectangle), unoccupied Eu{sup 2+} 5d band (white rectangle), occupied Eu{sup 2+} 4f ground states (filled red circles) and unoccupied Eu{sup 2+} ground states (open red circles). - Highlights: • Density functional theory calculations on the stability of M{sup 2+}L{sup 3+}SiN{sub 3} compounds. • Solid-state reaction synthesis of Eu{sub 2}SiN{sub 3}, CaLaSiN{sub 3}, EuLaSiN{sub 3} and CaEuSiN{sub 3}. • Determination of the Eu{sup 2+} 4f–5d and Eu{sup 3+} CT transitions in M{sup 2+}L{sup 3+}SiN{sub 3} compounds. • Oxidation of Eu{sub 2}SiN{sub 3} in dry air takes place via a nitrogen retention complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic and electronic structures of lattice mismatched Cu 2 O/TiO 2 interfaces

Heterojunction interfaces between metal oxides are often highly lattice mismatched. The atomic and electronic structures of such interfaces, however, are not well understood. Here, we have synthesized Cu 2 O/TiO 2 heterojunction thin films with 13% lattice mismatch and studied the interface via experimental methods and large-scale density function theory calculations of supercells containing 1300 atoms. We find that an interface of epitaxial quality is formed via a coincidence site lattice of 8 Cu 2 O unit cells matching 9 TiO 2 unit cells. Calculations reveal the existence of a dislocation core of the O sublattices at the interface and a random arrangement of one layer of interfacial Cu atoms. The interfacial electronic structure is found to be mostly determined by the interfacial Cu distribution, rather than by the O dislocation core. The conduction band minimum and valence band maximum states are spatially separated, and there is no strongly localized state near the core.

74 ATOMIC AND MOLECULAR PHYSICS↗

Angular Momentum in Rotating Superfluid Droplets

The angular momentum of rotating superfluid droplets originates from quantized vortices and capillary waves, the interplay between which remains to be uncovered. Here, the rotation of isolated submicrometer superfluid 4 He droplets is studied by ultrafast x-ray diffraction using a free electron laser. The diffraction patterns provide simultaneous access to the morphology of the droplets and the vortex arrays they host. In capsule-shaped droplets, vortices form a distorted triangular lattice, whereas they arrange along elliptical contours in ellipsoidal droplets. The combined action of vortices and capillary waves results in droplet shapes close to those of classical droplets rotating with the same angular velocity. The findings are corroborated by density functional theory calculations describing the velocity fields and shape deformations of a rotating superfluid cylinder.

O'Connell, Sean M.O.↗

Density functional theory study of the point defect energetics in y- LiAlO2, Li2ZrO3 and Li2TiO3 materials

The defect energetics in y-LiAlO2, Li2 TiO3, and Li2ZrO3 materials used in tritium-producing burnable absorber rods and fusion solid breeder applications was investigated using density functional theory calculations. A comprehensive analysis of the charged defects was performed for cation and anion vacancies, interstitials, antisite defects, and 3H interstitial and substitutional defects to understand the defect structures at different charge states and their stability as a function of the electron chemical potential across the electronic band gap of crystals at operating temperature and oxygen partial pressure conditions. Through investigation of the local bonding configurations and projected electronic density of states of the 3H related charged defects, including 3H interstitial (3Hint), 3Hint bound with an oxygen interstitial, and substitutional 3H defects on the Li, O, and M (M¼ Al, Ti, and Zr) sites (3HLi, 3HO, and 3HM), the 3H states of these defects were distinguished, i.e. triton (3Hþ), neutral 3H (3H0), and tritide (3H-). The two stable 3H related defects under the Li2O-LiOH rich condition at T ¼ 1000 K, between P(O2) ¼ 10^-5 and 10^-25 atm were found to be 3H interstitial (3Hint) with charge of þ1 and neutral 3HLi (substitutional 3H on the Li site) in g-LiAlO2, Li2ZrO3, and Li2TiO3 materials. For Li2TiO3, the stability of the 3HO with charge of þ1 and 3HM with charge of 3 becomes comparable to that of 3Hint atm at P(O2) ¼ 1025 atm. Overall, our results indicate the energy cost to form the stable 3H defect species (3Hint with charge of þ1 and neutral 3HLi) in g-LiAlO2 material are higher than those of Li2TiO3 and Li2ZrO3 materials at T ¼ 1000 K, between P(O2) ¼ 10^-5 and 10^-25 atm. Interactions between the point defects and the 3H species as a function of the electron chemical potential across the band gap of the materials were also examined. All the three materials exhibit strong interactions between the 3Hint and the O vacancy, Li vacancy, and O interstitial defects with the strongest interactions found in the case of y-LiAlO2 phase, indicating a higher tendency to trap 3H at such point defects in y-LiAlO2.

Charged defect analysis, Tritium-producing burnabl↗

Possible Mott transition in layered Sr 2 Mn 3 As 2 O 2 single crystals

Single crystals of Sr 2 Mn 3 As 2 O 2 have been grown for the first time, for which we show a possible layer-selective Mott insulator behavior. This compound stands out as a hybrid structure of MnO 2 and MnAs layers, analogously to the active CuO 2 and FeAs layers, respectively, in the cuprate and iron-based high-temperature superconductors. Electrical transport, neutron diffraction measurements, together with density functional theory calculations on Sr 2 Mn 3 As 2 O 2 single crystals converge toward a picture of independent magnetic order at T 1 ~ 79 K and T 2 ~ 360 K for the two Mn sublattices, with insulating behavior at odds with the metallic behavior predicted by calculations. Furthermore, our inelastic neutron-scattering studies of spin-wave dispersions for the Mn(1) sublattice reveal an effective magnetic exchange coupling of SJ ~ 3.7 meV. This is much smaller than those for the Mn(2) sublattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Influence of Nonstoichiometry on Proton Conductivity in Thin-Film Yttrium-Doped Barium Zirconate

Proton-conducting perovskites have been widely studied because of their potential application as solid electrolytes in intermediate temperature solid oxide fuel cells. Structural and chemical heterogeneities can develop during synthesis, device fabrication, or service, which can profoundly affect proton transport. Here, we use time-resolved Kelvin probe force microscopy, scanning transmission electron microscopy, atom probe tomography, and density functional theory calculations to intentionally introduce Ba-deficient planar and spherical defects and link the resultant atomic structure with proton transport behavior in both stoichiometric and nonstoichiometric epitaxial, yttrium-doped barium zirconate thin films. The defects were intentionally induced through high-temperature annealing treatment, while maintaining the epitaxial single crystalline structure of the films, with an overall relaxation in the atomic structure. As a result, the annealed samples showed smaller magnitudes of local lattice distortions because of the formation of proton polarons, thereby leading to decreased proton-trapping effect. This resulted in a decrease in the activation energy for proton transport, leading to faster proton transport.

36 MATERIALS SCIENCE↗