Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional approximations”

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 919 records · Page 51

Orientation-dependent structural properties during growth and growth mechanism of CoO films

The orientation-dependent local structural properties of CoO films on sapphire substrates during growth were investigated through linearly-polarized extended X-ray absorption fine structure (EXAFS) measurements. Specifically, CoO(1 1 1) and (100) crystals with a rock-salt structure (Fm3m) were epitaxially grown on alpha-Al2O3(0001) and (1012) substrates, respectively, at 700 degrees C using a radio-frequency sputtering system. The local structural properties of CoO films in the in-plane and out-of-plane orientations were quantitatively determined using linearly-polarized EXAFS at the Co K-edge during growth. The EXAFS analysis revealed that during the initial stages of growth, the local structural properties exhibit significant differences compared to thick films, with short atomic distances and large (small) Debye-Waller factors observed in the out-of-plane (in-plane) orientations. The local structural strain mostly diminished when approximately 20 CoO layers accumulated on the substrate. Density functional theory (DFT) calculations further supported these findings by confirming that cobalt atoms initially form stable bonds with the sapphire surface, leading to the simultaneous growth of oxygen and cobalt layers in a coordinated manner through layer-by-layer growth.

36 MATERIALS SCIENCE↗

Periodic domain boundary ordering in a dense molecular adlayer: Sub-saturation carbon monoxide on Pd(111)

Here, we describe a previously unreported ordered phase of carbon monoxide adsorbed on the (111) facet of single crystal palladium at near-saturation coverage. The adlayer superstructure is identified from low energy electron diffraction to be c(16×2) with respect to the underlying Pd(111) surface net. The ideal coverage is determined to be 0.6875 ML, approximately 92% of the 0.75–ML saturation coverage. Density functional theory calculations support a model for the molecular packing characterized by strips of locally-saturated (2×2) regions, with the CO bound near high-symmetry surface sites, separated by antiphase domain boundaries. The structure exists in a narrow coverage range and is prepared by heating the saturated adlayer to desorb a small fraction of the CO. Comparison of the c(16×2) domain-boundary structure with structural motifs at lower coverages suggests that between 0.6 and 0.6875 ML the adlayer order may be more strongly influenced by interadsorbate repulsion than by adsorption-site-specific interactions. The system is an example of the structural complexity that results from the compromise between adsorbate–substrate and adsorbate–adsorbate interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-Time Time-Dependent Electronic Structure Theory

Real-time electronic structure methods provide an unprecedented view of electron dynamics and ultrafast spectroscopy on the atto- and femtosecond time scale with vast potential to yield new insights into the electronic behavior of molecules and materials. Here, we discuss the fundamental theory underlying various real-time electronic structure methods as well as advantages and disadvantages of each. We give an overview of the numerical techniques that are widely used for real-time propagation of the quantum electron dynamics with an emphasis on Gaussian basis set methods. We also showcase many of the chemical applications and scientific advances made by using real-time electronic structure calculations and provide an outlook of possible new directions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into Nonelectroactive C–C Bond Formation on Cu(100) during Electrochemical CO 2 Reduction from Multiconfigurational Wavefunction Theory

Carbon–carbon (C–C) bond formation is necessary for hydrocarbon (and oxygenate) synthesis beyond methane (and formate/formic acid) during electrochemical CO and CO 2 reduction (ECOR and ECO 2 R). Cu has notable ability to form hydrocarbons compared to other pure metals. In particular, the (100) facet of face-centered cubic Cu forms ethylene competitively with H 2 and methane during both ECOR and ECO 2 R. Past simulations based on density functional theory (DFT) with standard exchange-correlation functional approximations predict fast nonelectroactive C–C bond formation channels involving adsorbed (*) CO together with another *CO, formyl (*CHO), or hydroxymethylidyne (*COH), forming OC*–*CO, OC*–CHO*, and OC*–*COH, respectively. Such simulations support the prevailing hypothesis that emergence of C 2 products is kinetically determined at the early stages of the reduction chemistry. Here we show, via simulations with more accurate many-body, i.e., “correlated”, wavefunction theory (enabled by an embedding scheme), that the coupling of *CO with a *CO or a *COH (previously predicted at the same level of theory to kinetically dominate over *CHO as the one-electron reduction product of *CO) is highly activated (kinetically impeded), with free energy barriers >1 eV, in contradiction to previous DFT-based simulations. Intriguingly, we find that the coupling of two adjacent *COHs incurs only a small barrier (<0.3 eV) and is exoergic (< –1 eV); however, given the predicted low surface mobility of *COH, the emergence of HOC*–*COH is also improbable, at least at low *COH coverages. We therefore conclude that it is highly unlikely for *CO to participate in nonelectroactive C–C bond formation on pristine Cu(100), contrary to conventional wisdom, and that the energetically favorable *COH dimerization may occur only after substantial buildup of *COH on the surface.

adsorption↗

Metal Halide Perovskites Demonstrate Radiation Hardness and Defect Healing in Vacuum

Herein, we subject formamidinium lead iodide films to oxygen-containing gases (flowing O 2 or free diffusion of lab atmosphere), inert gases (flowing He, Ar, or N 2 ), and vacuum. Our films are irradiated by Cu Ka X-rays and held at 75 °C while X-ray diffraction is recorded. Under all gas conditions, we observe a reproducible 1.1 ± 0.5 A 3 perovskite lattice contraction from an initial unit cell volume of 256.5 ± 0.8 A 3 concurrent with continuous perovskite loss and lead iodide growth. Oxygen-containing gases increase the reaction rates without materially altering perovskite structural changes. Under the same temperature and irradiation conditions in vacuo, a self-healing reaction is observed, exhibited by a reproducible (0.9 ± 0.3 A 3 ) lattice expansion and stabilization of the perovskite. Furthermore, interactions between the perovskite, defects, and minority phases are simulated by generalized gradient approximation Perdew-Burke-Ernzerhof (GGA-PBE) density functional theory. Lattice contraction indicates an increase in the concentration of Schottky defects-pairs of formamidinium and iodine vacancies. Under irradiation in every atmospheric condition, a solid solution of Schottky defects with a concentration of several percent diffuses and precipitates forming lead iodide and consuming the defects. In the presence of ionized gases, this framework is modified to include the continual loss of formamidinium and iodine ions from the perovskite forming Schottky defects.

14 SOLAR ENERGY↗

Two-Component Relativistic Calculations of Electric-Field Gradients Using Exact Decoupling Methods: Spin–orbit and Picture-Change Effects

Electric field gradients (EFGs) for heavy nuclei in diatomic molecules (hydrogen halides and group-13 iodides), uranyl (UO 2 2+ ), and a uranyl carbonate complex, were computed with different formally exactly (or nearly exactly) decoupled two-component relativistic Hamiltonians: the “exact two-component” (X2C) method, the Barysz–Sadlej–Snijders (BSS) approach, and up to 35th order in the Douglas–Kroll–Hess (DKH) expansion, utilizing a new implementation in the open-source NWChem quantum chemistry package. Results from two-component Hartree–Fock and density functional theory (DFT) calculations at the scalar relativistic approximation as well as including spin–orbit coupling are reported. Picture-change corrected EFGs obtained with X2C, BSS, and high-order DKH are shown to be numerically equivalent. Effects from spin–orbit coupling on the EFGs tend to be moderate but should not be excluded for reliable predictions. Picture-change effects on the EFG of a heavy atom can be as large as the correct picture-change corrected value, as is shown for the uranium EFG in the uranyl systems. The results demonstrate the necessity of simultaneously including spin–orbit and picture-change effects in calculations of heavy atom EFGs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charge Density Wave and Narrow Energy Gap at Room Temperature in 2D Pb 3– x Sb 1+ x S 4 Te 2-δ with Square Te Sheets

We report a new two-dimensional compound Pb 3– x Sb 1+ x S 4 Te 2-δ has a charge density wave (CDW) at room temperature. The CDW is incommensurate with q-vector of 0.248(6)a* + 0.246(8)b* + 0.387(9)c* for x = 0.29(2) and d = 0.37(3) due to positional and occupational long range ordering of Te atoms in the sheets. The modulated structure was refined from the single crystal X-ray diffraction data with a superspace group P$\bar{1}$(αβγ)0 using (3 + 1)-dimensional crystallography. The resistivity increases with decreasing temperature, suggesting semiconducting behavior. The transition temperature (T CDW ) of the CDW is ~ 345 K above which the Te square sheets become disordered with no q-vector. Lastly, first-principles density functional theory calculations on the undistorted structure and an approximate commensurate supercell reveal that the gap is due to the structure modulation.

36 MATERIALS SCIENCE↗

Anisotropic in-plane thermal conductivity of black phosphorus nanoribbons at temperatures higher than 100 K

Black phosphorus attracts enormous attention as a promising layered material for electronic, optoelectronic and thermoelectric applications. Here we report large anisotropy in in-plane thermal conductivity of single-crystal black phosphorus nanoribbons along the zigzag and armchair lattice directions at variable temperatures. Thermal conductivity measurements were carried out under the condition of steady-state longitudinal heat flow using suspended-pad micro-devices. We discovered increasing thermal conductivity anisotropy, up to a factor of two, with temperatures above 100 K. A size effect in thermal conductivity was also observed in which thinner nanoribbons show lower thermal conductivity. Analysed with the relaxation time approximation model using phonon dispersions obtained based on density function perturbation theory, the high anisotropy is attributed mainly to direction-dependent phonon dispersion and partially to phonon–phonon scattering. Lastly, our results revealing the intrinsic, orientation-dependent thermal conductivity of black phosphorus are useful for designing devices, as well as understanding fundamental physical properties of layered materials.

36 MATERIALS SCIENCE↗

Magic in twisted transition metal dichalcogenide bilayers

The long-wavelength moiré superlattices in twisted 2D structures have emerged as a highly tunable platform for strongly correlated electron physics. We study the moiré bands in twisted transition metal dichalcogenide homobilayers, focusing on WSe 2 , at small twist angles using a combination of first principles density functional theory, continuum modeling, and Hartree-Fock approximation. We reveal the rich physics at small twist angles θ < 4°, and identify a particular magic angle at which the top valence moiré band achieves almost perfect flatness. In the vicinity of this magic angle, we predict the realization of a generalized Kane-Mele model with a topological flat band, interaction-driven Haldane insulator, and Mott insulators at the filling of one hole per moiré unit cell. The combination of flat dispersion and uniformity of Berry curvature near the magic angle holds promise for realizing fractional quantum anomalous Hall effect at fractional filling. We also identify twist angles favorable for quantum spin Hall insulators and interaction-induced quantum anomalous Hall insulators at other integer fillings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Full-scale computation for all the thermoelectric property parameters of half-Heusler compounds

The thermoelectric performance of materials relies substantially on the band structures that determine the electronic and phononic transports, while the transport behaviors compete and counter-act for the power factor PF and figure-of-merit ZT. These issues make a full-scale computation of the whole set of thermoelectric parameters particularly attractive, while a calculation scheme of the electronic and phononic contributions to thermal conductivity remains yet challenging. In this work, we present a full-scale computation scheme based on the first-principles calculations by choosing a set of doped half- Heusler compounds as examples for illustration. The electronic structure is computed using the WIEN2k code and the carrier relaxation times for electrons and holes are calculated using the Bardeen and Shockley’s deformation potential (DP) theory. The finite-temperature electronic transport is evaluated within the framework of Boltzmann transport theory. In sequence, the density functional perturbation combined with the quasi-harmonic approximation and the Klemens’ equation is implemented for calculating the lattice thermal conductivity of carrier-doped thermoelectric materials such as Tidoped NbFeSb compounds without losing a generality. The calculated results show good agreement with experimental data. Lastly, the present methodology represents an effective and powerful approach to calculate the whole set of thermoelectric properties for thermoelectric materials.

36 MATERIALS SCIENCE↗

Near-edge band structures and band gaps of Cu-based semiconductors predicted by the modified Becke-Johnson potential plus an on-site Coulomb U

Diamond-like Cu-based multinary semiconductors are a rich family of materials that hold promise in a wide range of applications. Unfortunately, accurate theoretical understanding of the electronic properties of these materials is hindered by the involvement of Cu d electrons. Density functional theory (DFT) based calculations using the local density approximation or generalized gradient approximation often give qualitative wrong electronic properties of these materials, especially for narrow-gap systems. The modified Becke-Johnson (mBJ) method has been shown to be a promising alternative to more elaborate theory such as the GW approximation for fast materials screening and predictions. However, straightforward applications of the mBJ method to these materials still encounter significant difficulties because of the insufficient treatment of the localized d electrons. We show that combining the promise of mBJ potential and the spirit of the well-established DFT + U method leads to a much improved description of the electronic structures, including the most challenging narrow-gap systems. A survey of the band gaps of about 20 Cu-based semiconductors calculated using the mBJ + U method shows that the results agree with reliable values to within ±0.2 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Infrared and Raman spectra, DFT-calculations and spectral assignments of germacyclohexane

Raman spectra of germacyclohexane in liquid and solid states were recorded and depolarization data obtained. Infrared absorption spectra of the vapor and liquid have been studied. The wavenumbers of the vibrational modes were derived in the harmonic and anharmonic approximation in B3LYP/ccpVTZ calculations. According to the calculations, germacyclohexane exists in the stable chair conformation, whereas a possible twist form should have more than 15 kJ·mol{sup -1} higher enthalpy of formation what makes this conformer experimentally not observable. The 27 A' and 21 A'' fundamentals were assigned on the basis of the calculations and infrared and Raman band intensities, contours of gas phase infrared spectral bands and Raman depolarization measurements. An average discrepancy of ca. 0.77 % was found between the observed and the calculated anharmonic wavenumbers for the 48 modes. Substitution of carbon atom with Ge atom in the cyclohexane ring is reasoning flattening of the ring.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab initio study of structural, electronic, magnetic alloys: XTiSb (X = Co, Ni and Fe)

Structural, electronic and magnetic properties of three semi-Heusler compounds of CoTiSb, NiTiSb and FeTiSb were calculated by the method (FP-LAPW) which is based on the DFT code WIEN2k. We used the generalized gradient approximation (GGA (06)) for the term of the potential exchange and correlation (XC) to calculate structural properties, electronic properties and magnetic properties. Structural properties obtained as the lattice parameter are in good agreement with the experimental results available for the electronic and magnetic properties was that: CoTiSb is a semiconductor NiTiSb is a metal and FeTiSb is a half-metal ferromagnetic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hydrogen molecule on lithium adsorbed graphene: A DFT study

Electronic structure calculations for the adsorption of molecular hydrogen on lithium (Li) decorated and pristine graphene have been studied systematically using SIESTA code [1] within the framework of the first-principle DFT under the Perdew-Burke-Ernzerhof (PBE) form of the generalized gradient approximation (GGA)[2], including spin polarization. The energy of adsorption of hydrogen molecule on graphene is always enhanced by the presence of co-adsorbed lithium. The most efficient adsorption configuration is when H{sub 2} is lying parallel to lithium adsorbed graphene which is in contrast to its adsorption on pristine graphene (PG) where it prefers perpendicular orientation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Effect of doping on electronic properties of HgSe

First principle study of electronic properties of pure and doped HgSe have been performed using all electron Full Potential Linearized Augmented Plane Wave (FP-LAPW) method using ELK code. The electronic exchange and co-relations are considered using Generalized Gradient Approach (GGA). Lattice parameter, Density of States (DOS) and Band structure calculations have been performed. The total energy curve (Energy vs Lattice parameter), DOS and band structure calculations are in good agreement with the experimental values and those obtained using other DFT codes. The doped material is studied within the Virtual Crystal Approximation (VCA) with doping levels of 10% to 25% of electrons (hole) per unit cell. Results predict zero band gap in undopedHgSe and bands meet at Fermi level near the symmetry point Γ. For doped HgSe, we found that by electron (hole) doping, the point where conduction and valence bands meet can be shifted below (above) the fermi level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Blind test of density-functional-based methods on intermolecular interaction energies

In the past decade, a number of approaches have been developed to fix the failure of (semi)local density-functional theory (DFT) in describing intermolecular interactions. In this work, the performance of several such approaches with respect to highly accurate benchmarks is compared on a set of separation-dependent interaction energies for ten dimers. Since the benchmarks were unknown before the DFT-based results were collected, this comparison constitutes a blind test of these methods.

36 MATERIALS SCIENCE↗

Coupled electron-ion Monte Carlo simulation of hydrogen molecular crystals

We performed simulations for solid molecular hydrogen at high pressures (250 GPa ≤ P ≤ 500 GPa) along two isotherms at T = 200 K (phase III) and at T = 414 K (phase IV). At T = 200 K, we considered likely candidates for phase III, the C2c and Cmca12 structures, while at T = 414 K in phase IV, we studied the Pc48 structure. We employed both Coupled Electron-Ion Monte Carlo (CEIMC) and Path Integral Molecular Dynamics (PIMD). The latter is based on Density Functional Theory (DFT) with the van der Waals approximation (vdW-DF). The comparison between the two methods allows us to address the question of the accuracy of the exchange-correlation approximation of DFT for thermal and quantum protons without recurring to perturbation theories. In general, we find that atomic and molecular fluctuations in PIMD are larger than in CEIMC which suggests that the potential energy surface from vdW-DF is less structured than the one from quantum Monte Carlo. We find qualitatively different behaviors for systems prepared in the C2c structure for increasing pressure. Within PIMD, the C2c structure is dynamically partially stable for P ≤ 250 GPa only: it retains the symmetry of the molecular centers but not the molecular orientation; at intermediate pressures, it develops layered structures like Pbcn or Ibam and transforms to the metallic Cmca-4 structure at P ≥ 450 GPa. Instead, within CEIMC, the C2c structure is found to be dynamically stable at least up to 450 GPa; at increasing pressure, the molecular bond length increases and the nuclear correlation decreases. For the other two structures, the two methods are in qualitative agreement although quantitative differences remain. Finally, we discuss various structural properties and the electrical conductivity. Here, we find that these structures become conducting around 350 GPa but the metallic Drude-like behavior is reached only at around 500 GPa, consistent with recent experimental claims.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ferroelectric barium titanate derivatives containing Mo and Mg for transparent photovoltaic applications

Using first-principles methods, we investigate the electronic properties of the [Ba(Mo 1/2 ,Mg 1/2 )O 3 ] x -[BaTiO 3 ] 1-x solid solution derived from barium titanate as a potential candidate to be used in photovoltaic devices. Focusing on the bandgap and its origin, we study the effect of different possible Mo and Mg contents, arrangements, and phases of [Ba(Mo 1/2 ,Mg 1/2 )O 3 ]x-[BaTiO 3 ] 1-x . We find that [Ba(Mo 1/2 ,Mg 1/2 )O 3 ] 0.25 -[BaTiO 3 ] 0.75 is a viable candidate for use in transparent photovoltaics due to its energy bandgap of 2.6 eV in the rhombohedral phase. In all cases, [Ba(Mo 1/2 ,Mg 1/2 )O 3 ]x-[BaTiO 3 ] 1-x materials exhibit spontaneous polarization that allows the exploitation of the bulk photovoltaic effect and in principle may allow high power conversion efficiency exceeding the Shockley-Queisser limit for these materials.

14 SOLAR ENERGY↗