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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Atomic-Scale Understanding of Structure and Properties of Complex Pyrophosphate Crystals by First-Principles Calculations

The electronic structure and mechanical and optical properties of five pyrophosphate crystals with very complex structures are studied by first principles density functional theory calculations. The results show the complex interplay of the minor differences in specific local structures and compositions can result in large differences in reactivity and interaction that are rare in other classes of inorganic crystals. These are discussed by dividing the pyrophosphate crystals into three structural units. H 2 P 2 O 7 is the most important and dominating unit in pyrophosphates. The other two are the influential cationic group with metals and water molecules. The strongest P-O bond in P 2 O 5 is the strongest bond for crystal cohesion, but O-H and N-H bonds also play an important part. Different type of bonding between O and H atoms such as O-H, hydrogen bonding, and bridging bonds are present. Metallic cations such as Mg, Zn, and Cu form octahedral bonding with O. The water molecule provides the unique H∙∙∙O bonds, and metallic elements can influence the structure and bonding to a certain extent. The two Cu-containing phosphates show the presence of narrow metallic bands near the valence band edge. All this complex bonding affects their physical properties, indicating that fundamental understanding remains an open question.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Path integral Monte Carlo simulations of H{sub 2} adsorbed to lithium-doped benzene: A model for hydrogen storage materials

Finite temperature quantum and anharmonic effects are studied in H{sub 2}–Li{sup +}-benzene, a model hydrogen storage material, using path integral Monte Carlo (PIMC) simulations on an interpolated potential energy surface refined over the eight intermolecular degrees of freedom based upon M05-2X/6-311+G(2df,p) density functional theory calculations. Rigid-body PIMC simulations are performed at temperatures ranging from 77 K to 150 K, producing both quantum and classical probability density histograms describing the adsorbed H{sub 2}. Quantum effects broaden the histograms with respect to their classical analogues and increase the expectation values of the radial and angular polar coordinates describing the location of the center-of-mass of the H{sub 2} molecule. The rigid-body PIMC simulations also provide estimates of the change in internal energy, ΔU{sub ads}, and enthalpy, ΔH{sub ads}, for H{sub 2} adsorption onto Li{sup +}-benzene, as a function of temperature. These estimates indicate that quantum effects are important even at room temperature and classical results should be interpreted with caution. Our results also show that anharmonicity is more important in the calculation of U and H than coupling—coupling between the intermolecular degrees of freedom becomes less important as temperature increases whereas anharmonicity becomes more important. The most anharmonic motions in H{sub 2}–Li{sup +}-benzene are the “helicopter” and “ferris wheel” H{sub 2} rotations. Treating these motions as one-dimensional free and hindered rotors, respectively, provides simple corrections to standard harmonic oscillator, rigid rotor thermochemical expressions for internal energy and enthalpy that encapsulate the majority of the anharmonicity. At 150 K, our best rigid-body PIMC estimates for ΔU{sub ads} and ΔH{sub ads} are −13.3 ± 0.1 and −14.5 ± 0.1 kJ mol{sup −1}, respectively.

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

Investigation of Uranyl Perchlorate Anion Complexes in the Gas Phase via Infrared Multiphoton Dissociation and Collision‐Induced Dissociation

The infrared multiphoton dissociation spectrum of the gas-phase uranyl perchlorate anion ([UO 2 (ClO 4 ) 3 ] − ) was measured. Here, the asymmetric uranyl stretch lies somewhere between 930 and 1030 cm −1 , though it could not be deconvoluted from a closely located perchlorate stretching mode. The experimentally measured spectrum was in good agreement with spectra calculated using density functional theory with the B3LYP and TPSSh functionals and Grimme's D3 dispersion corrections with Becke–Johnson damping. The calculations suggest that the asymmetric uranyl stretch lies in the range of 970–980 cm −1 , about 20–30 cm −1 higher than for the analogous uranyl trinitrato complex and consistent with perchlorate as a weaker ligand than nitrate. Like the uranyl trinitrato anion, fragmentation of [UO 2 (ClO 4 ) 3 ] − by collision-induced dissociation resulted in loss of a ClO 3 • radical to form [UO 3 (ClO 4 ) 2 ] − . Infrared multiphoton dissociation measurements of [UO 3 (ClO 4 ) 2 ] − indicate that the structure of the product ion has an O − ligand bound to the uranyl in a T-arrangement and two perchlorate ligands, matching the findings from uranyl nitrato complexes. The uranyl stretch in this complex was slightly separated from the perchlorate modes and was measured at 906 cm −1 , slightly higher than for the analogous nitrate complex, again consistent with perchlorate as a weaker ligand than nitrate in the gas phase.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2023

Precisely doping the surface of tin-based electrocatalysts for improved CO 2 conversion to liquid chemicals

Doping tin catalysts with sulfur can improve the electrochemical CO 2 conversion into formate/formic acid, but the lack of composition-dependent activity trends hinders further catalyst development. Here, we precisely controlled the composition of sulfur-doped Sn catalysts to show that sulfur doping only improves CO 2 conversion over a very narrow composition range, achieving maximum activity at 1.4 at% S. In situ Raman spectroscopy indicted working catalysts were in a primarily metallic state (e.g. S-Sn), and we achieved some of the highest reported partial current densities in both H-cell and full-cell electrolyzer configurations. Density Functional Theory calculations predicted S atoms preferentially occupied the catalyst surface and improved CO 2 reduction by localizing charge density at the catalyst/intermediate interface, which stabilized the *OCOH intermediate and lowered the CO 2 conversion thermodynamic barrier. Our work quantifies the composition-dependent influence of S dopants on Sn-based CO 2 reduction catalysts and provides a pathway for maximizing their CO 2 conversion activity.

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

Is Cs 2 TiBr 6 a promising Pb-free perovskite for solar energy applications?

In a quest for Pb-free perovskites suitable for solar energy applications, Cs 2 TiBr 6 has recently been reported as a promising compound, with appropriate optical and electrical properties as well as high stability under environmental stresses. In this study, we pursue investigation on this compound, demonstrating phase pure Cs 2 TiBr 6 powder formation using solution synthesis and providing complementary experimental characterization and theoretical calculations. An experimental absorption onset of around 2.0 eV is extracted and a weak broad photoluminescence is measured. Density functional theory calculations predict an indirect bandgap, parity-forbidden for both the direct and indirect transitions, which explains the weak and Stokes shifted luminescence. Additionally, we highlight the strong instability of Cs 2 TiBr 6 powder in ambient atmosphere. Therefore, our experimental results supported by theoretical calculations differ from previous results and raise doubts on the suitability of Cs 2 TiBr 6 in its pristine form for solar energy applications.

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

Elastic constants and Fermi surface topology change in Calaverite AuTe{sub 2}: A density functional study

Structural, elastic, electronic and Fermi surface studies of AuTe{sub 2} have been carried out by means of first principles calculations based on density functional theory. The calculated ground state properties agree well with the experiment. Fermi surface and elastic constants are predicted for the first time and from the calculated elastic constants we find the compound to be mechanically stable satisfying the stability criteria of monoclinic structure. In addition, we also find the c-axis to be more compressible than the other two which is also speculated from the present work. The metallic behaviour of this compound is confirmed from the electronic band structure calculation as we find the bands to cross the Fermi level (E{sub F}). In addition, we also observe a FS topology change under pressure which is also explained in the present work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Cu[Ni(2,3-pyrazinedithiolate) 2 ] Metal–Organic Framework for Electrocatalytic Hydrogen Evolution

The application of metal–organic frameworks (MOFs) as electrocatalysts for small molecule activation has been an emerging topic of research. Previous studies have suggested that two-dimensional (2D) dithiolene-based MOFs are among the most active for the hydrogen evolution reaction (HER). Here, a three-dimensional (3D) dithiolene-based MOF, Cu[Ni(2,3-pyrazinedithiolate) 2 ] (1), is evaluated as an electrocatalyst for the HER. In pH 1.3 aqueous electrolyte solution, 1 exhibits a catalytic onset at –0.43 V vs the reversible hydrogen electrode (RHE), an overpotential (η 10 mA/cm 2 ) of 0.53 V to reach a current density of 10 mA/cm 2 , and a Tafel slope of 69.0 mV/dec. Interestingly, under controlled potential electrolysis, 1 undergoes an activation process that results in a more active catalyst with a 200 mV reduction in the catalytic onset and η 10 mA/cm 2 . It is proposed that the activation process is a result of the cleavage of Cu–N bonds in the presence of protons and electrons. Furthermore, this hypothesis is supported by various experimental studies and density functional theory calculations.

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

Size and structure effects of Pt{sub N} (N = 12 − 13) clusters for the oxygen reduction reaction: First-principles calculations

Size and structure effects on the oxygen reduction reaction on Pt{sub N} clusters with N = 12–13 atoms have been investigated using periodic density functional theory calculations with the generalized gradient approximation. To describe the catalytic activity, we calculated the O and OH adsorption energies on the cluster surface. The oxygen binding on the 3-fold hollow sites on stable Pt{sub 12−13} cluster models resulted more favorable for the reaction with O, compared with the Pt{sub 13}(I{sub h}) and Pt{sub 55}(I{sub h}) icosahedral particles, in which O binds strongly. However, the rate-limiting step resulted in the removal of the OH species due to strong adsorptions on the vertex sites, reducing the utility of the catalyst surface. On the other hand, the active sites of Pt{sub 12−13} clusters have been localized on the edge sites. In particular, the OH adsorption on a bilayer Pt{sub 12} cluster is the closest to the optimal target; with 0.0-0.2 eV weaker than the Pt(111) surface. However, more progress is necessary to activate the vertex sites of the clusters. The d-band center of Pt{sub N} clusters shows that the structural dependence plays a decisive factor in the cluster reactivity.

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

Temperature effects in first-principles solid state calculations of the chemical shielding tensor made simple

We study the effects of atomic vibrations on the solid-state chemical shielding tensor using first principles density functional theory calculations. At the harmonic level, we use a Monte Carlo method and a perturbative expansion. The Monte Carlo method is accurate but computationally expensive, while the perturbative method is computationally more efficient, but approximate. We find excellent agreement between the two methods for both the isotropic shift and the shielding anisotropy. The effects of zero-point quantum mechanical nuclear motion are important up to relatively high temperatures: at 500 K they still represent about half of the overall vibrational contribution. We also investigate the effects of anharmonic vibrations, finding that their contribution to the zero-point correction to the chemical shielding tensor is small. We exemplify these ideas using magnesium oxide and the molecular crystals L-alanine and β-aspartyl-L-alanine. We therefore propose as the method of choice to incorporate the effects of temperature in solid state chemical shielding tensor calculations using the perturbative expansion within the harmonic approximation. This approach is accurate and requires a computational effort that is about an order of magnitude smaller than that of dynamical or Monte Carlo approaches, so these effects might be routinely accounted for.

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

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. Here, we present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

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

Initial and secondary oxidation products on the Si(111)-(7 × 7) surface identified by atomic force microscopy and first principles calculations

We investigate the initial and secondary oxidation products on the Si(111)-(7 × 7) surface at room-temperature using atomic force microscopy (AFM) and density functional theory calculations. At the initial oxidation stages, we find that there are two types of bright spots in AFM images. One of them is identified as a Si adatom with one O atom inserted into one of the backbonds, while the other is ascribed to a Si adatom with two inserted O atoms. We observe that the latter one turns into the secondary oxidation product by a further coming O{sub 2} molecule, which appears as a more protruded bright spot. The atomic configuration of this product is identified as Si adatom whose top and all three backbonds make bonds with O atoms. The appearances of initial and secondary oxidation products are imaged as bright and dark sites by scanning tunneling microscopy, respectively. It is revealed that AFM gives us the topographic information close to the real atomic corrugation of adsorbed structures on the semiconductor surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Electronic structure of GaP/Si (001) heterojunctions and the role of hydrogen passivation

Abstract Epitaxially grown single crystal GaP on Si is of considerable interest due to being nearly lattice matched to Si, making it attractive for III‐V/Si solar cells. GaP has been used as a buffer layer for III‐V/Si solar cells and also in selective contact Si solar cells. The performance and functionality of such devices are strongly influenced by the presence of localized states at the GaP/Si interface. Here, we examine the electronic structure of GaP/Si(001) heterojunctions and the effect of hydrogen (H) passivation at the interface, in contrast to interface mixing, through density functional theory calculations. Our calculations show that due to the heterovalent mismatch nature of the GaP/Si interface, there is a high density of localized states at the abrupt GaP/Si interface due to the excess charge associated with heterovalent bonding, as reported elsewhere. We find that the addition of H leads to additional bonding at the interface, which mitigates the charge imbalance, and greatly reduces the density of localized states, leading to a nearly ideal heterojunction. A similar result is found with a completely intermixed interface (alternating cation and anion bonding) in terms of low interface state density. However, when the intermixing occurs through single‐layer or double‐layer terraces, the benefits of intermixing are lost and the interface state density reverts more to the abrupt interface case.

Meidanshahi, Reza Vatan↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Potential energy curves for the ground and low-lying excited states of CuAg

The ground and low-lying excited states of heteronuclear diatomic CuAg are examined by multi-reference configuration interaction (MRCI) method. Relativistic effects were treated and probed in two steps. Scalar terms were considered using the spin-free DKH Hamiltonian as a priori and spin-orbit coupling was calculated perturbatively via the spin-orbit terms of the Breit-Pauli Hamiltonian based on MRCI wavefunctions. Potential energy curves of the spin-free states and their corresponding Ω components correlating with the separated atom limits {sup 2}S(Cu) + {sup 2}S(Ag) and {sup 2}D(Cu) + {sup 2}S(Ag) are obtained. The results are in fine agreement with the experimental measurements and tentative conclusions for the ion-pair B0{sup +} state are confirmed by our theoretical calculations. Illustrative results are presented to reveal the relative importance and magnitude of the scalar and spin-orbit effects on the spectroscopic properties of this molecule. Time dependent density functional theory calculations, using the LDA, BLYP, B3LYP, and SAOP functionals have been carried out for CuAg and the accuracy of TD-DFT has been compared with ab initio results.

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

Transuranium Sulfide via the Boron Chalcogen Mixture Method and Reversible Water Uptake in the NaCu T S 3 Family

The behavior of 5f electrons in soft ligand environments makes actinides, and especially transuranium chalcogenides, an intriguing class of materials for fundamental studies. Due to the affinity of actinides for oxygen, however, it is a challenge to synthesize actinide chalcogenides using non-metallic reagents. Using the Boron Chalcogen Mixture (BCM) method, we achieved the synthesis of the transuranium sulfide NaCuNpS 3 starting from the oxide reagent, NpO 2 . Via the same synthetic route, the isostructural composition of NaCuUS 3 was synthesized and the material contrasted with NaCuNpS 3 . Single crystals of the U-analog, NaCuUS 3 , were found to undergo an unexpected reversible hydration process to form NaCuUS 3 ·xH 2 O (x ≈ 1.5). Here, a large combination of techniques was used to fully characterize the structure, hydration process, and electronic structures, specifically a combination of single crystal, powder, high temperature powder X-ray diffraction, extended X-ray ab-sorption fine structure, infrared, and inductively coupled plasma spectroscopies, thermogravimetric analysis, and density functional theory calculations. The outcome of these analyses enabled us to determine the composition of NaCuUS 3 ·xH 2 O and obtain a structural model that demonstrated the retention of the local structure within the [CuUS 3 ] – layers throughout the hydration-dehydration process. Band structure, density of states, and Bader charge calculations for NaCuUS 3 , NaCu-US 3 ·xH 2 O, and NaCuNpS 3 along with X-ray absorption near edge structure, UV-vis-NIR and work function measurements on ACuUS 3 (A = Na, K, and Rb) and NaCuUS 3 ·xH 2 O samples were carried out to demonstrate that electronic properties arise from the [CuTS 3 ] – layers and show surprisingly little dependence on the interlayer distance.

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

High pressure and time resolved studies of optical properties of n-type doped GaN/AlN multi-quantum wells: Experimental and theoretical analysis

High-pressure and time-resolved studies of the optical emission from n-type doped GaN/AlN multi-quantum-wells (MQWs) with various well thicknesses are analysed in comparison with ab initio calculations of the electronic (band structure, density of states) and optical (emission energies and their pressure derivatives, oscillator strength) properties. The optical properties of GaN/AlN MQWs are strongly affected by quantum confinement and polarization-induced electric fields. Thus, the photoluminescence (PL) peak energy decreases by over 1 eV with quantum well (QW) thicknesses increasing from 1 to 6 nm. Furthermore, the respective PL decay times increased from about 1 ns up to 10 μs, due to the strong built-in electric field. It was also shown that the band gap pressure coefficients are significantly reduced in MQWs as compared to bulk AlN and GaN crystals. Such coefficients are strongly dependent on the geometric factors such as the thickness of the wells and barriers. The transition energies, their oscillator strength, and pressure dependence are modeled for tetragonally strained structures of the same geometry using a full tensorial representation of the strain in the MQWs under external pressure. These MQWs were simulated directly using density functional theory calculations, taking into account two different systems: the semi-insulating QWs and the n-doped QWs with the same charge density as in the experimental samples. Such an approach allowed an assessment of the impact of n-type doping on optical properties of GaN/AlN MQWs. We find a good agreement between these two approaches and between theory and experimental results. We can therefore confirm that the nonlinear effects induced by the tetragonal strain related to the lattice mismatch between the substrates and the polar MQWs are responsible for the drastic decrease of the pressure coefficients observed experimentally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A computational framework for automation of point defect calculations

We have developed a complete and rigorously validated open-source Python framework to automate point defect calculations using density functional theory. Furthermore, the framework provides an effective and efficient method for defect structure generation, and creation of simple yet customizable workflows to analyze defect calculations. This package provides the capability to compute widely-accepted correction schemes to overcome finite-size effects, including (1) potential alignment, (2) image-charge correction, and (3) band filling correction to shallow defects. Using Si, ZnO and In2O3 as test examples, we demonstrate the package capabilities and validate the methodology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Two-Dimensional Molecular Charge Density Waves in Single-Layer-Thick Islands of a Dirac Fermion System

Charge density waves have been intensely studied in inorganic materials such as transition metal dichalcogenides however its counterpart in organic materials has yet to be explored in detail. Here we report the finding of a robust two dimensional charge density waves in molecular layers formed by α-(BEDT-TTF) 2 -I 3 on a Ag(111) surface. Low temperature scanning tunnelling microscopy images of a multi-layer thick α-(BEDT-TTF) 2 -I 3 on Ag(111) substrate reveal coexistence of 5a0 x 5a0 and √31 a0 x √31 a0 R9° charge density wave patterns commensurate with the underlying molecular lattice at 80 K. Both charge density wave patterns remain in nano-size molecular islands with just a single constituent molecular-layer thickness at 80 K and 5 K. Local tunneling spectroscopy measurements reveal the variation of the gap from 244 meV to 288 meV between the maximum and minimum charge density wave locations. Density functional theory calculations further confirm a vertical positioning of BEDT-TTF molecules in the molecular layer. While the observed charge density wave patterns are stable for the defect sites, they can be reversely switched for one molecular lattice site by means of inelastic tunnelling electron energy transfer with the electron energies exceeding 400 meV using a scanning tunneling microscope manipulation scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗