Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quantifying Large Lattice Relaxations in Photovoltaic Devices

Temporal variations of Cu(In,Ga)Se 2 photovoltaic device properties during light exposure at various temperatures and voltage biases for times up to 100 h are analyzed using the kinetic theory of large lattice relaxations. Open-circuit voltage and p-type doping increased with charge injection and decreased with temperature at low injection conditions. Furthermore, lattice relaxation can account for both trends and activation energies extracted from the data are approximately 0.9 and 1.2 eV for devices with lower and higher sodium content, respectively. In these devices, increased sodium content resulted in higher initial p-type doping with greater stability. First-principles calculations providing revised activation energies for the (V Se -V Cu ) complex suggest that this defect does not account for the metastability observed here.

14 SOLAR ENERGY↗
DOE OSTI2020

Ti- and Fe-related charge transition levels in $β$ – Ga 2 O 3

Here, deep-level transient spectroscopy measurements on $β$ – Ga 2 O 3 crystals reveal the presence of three defect signatures labeled $E$ 2a , $E$ 2b , and $E$ 3 with activation energies at around 0.66 eV, 0.73 eV, and 0.95 eV below the conduction band edge. Using secondary ion mass spectrometry, a correlation between the defect concentration associated with $E$ 3 and the Ti concentration present in the samples was found. Particularly, it is found that $E$ 3 is the dominant Ti-related defect in $β$ – Ga 2 O 3 and is associated with a single Ti atom. This finding is further corroborated by hybrid functional calculations that predict Ti substituting on an octahedral Ga site, denoted as Ti GaII , to be a good candidate for $E$ 3 . Moreover, the deep level transient spectroscopy results show that the level previously labeled $E$ 2 and attributed to Fe substituting on a gallium site (FeGa) consists of two overlapping signatures labeled E2a and E2b. We tentatively assign $E$ 2a and $E$ 2b to Fe substituting for Ga on a tetrahedral or an octahedral site, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Coordination Behaviour of Cu I Photosensitizers Bearing Multidentate Ligands Investigated by X-ray Absorption Spectroscopy

A systematic series of four novel homo- and heteroleptic CuI photosensitizers based on tetradentate 1,10-phenanthroline ligands of the type X^N^N^X containing two additional donor moieties in the 2,9-position (X=SMe or OMe) were designed. Their solid-state structures were assessed by X-ray diffraction. Cyclic voltammetry, UV-vis absorption, emission and X-ray absorption spectroscopy were then used to determine their electrochemical, photophysical and structural features in solution. Following, time-resolved X-ray absorption spectroscopy in the picosecond time scale, coupled with time-dependent density functional theory calculations, provided in-depth information on the excited state electron configurations. For the first time, a significant shortening of the Cu–X distance and a change in the coordination mode to a pentacoordinated geometry is shown in the excited states of the two homoleptic complexes. These findings are important with respect to a precise understanding of the excited state structures and a further stabilization of this type of photosensitizers.

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

Vibrational Entropy Contribution to Mixing Free Energy of Ni-Rich LiNi 1- y Co y O 2

First principles density-functional-theory calculations are performed of the mixing free energy of the layered oxide LiNi 1-y Co y O 2 , a model Li-ion-battery cathode material, at 1000 K. Analysis of internal energy, LN exchange, configurational entropy and vibrational entropy are combined to obtain the net mixing free energy. The vibrational entropy requires a more indirect treatment than the other terms; analysis based on heat capacity is presented here, to show that the vibrational entropy is of comparable magnitude to the configurational entropy, and it tends to promote solid-solution behavior in LiNi 1-y Co y O 2 .

25 ENERGY STORAGE↗
DOE OSTI2022

A metal-supported single-atom catalytic site enables carbon dioxide hydrogenation

Nitrogen-doped graphene-supported single atoms convert CO 2 to CO, but fail to provide further hydrogenation to methane – a finding attributable to the weak adsorption of CO intermediates. To regulate the adsorption energy, here we investigate the metal-supported single atoms to enable CO 2 hydrogenation. We find a copper-supported iron-single-atom catalyst producing a high-rate methane. Density functional theory calculations and in-situ Raman spectroscopy show that the iron atoms attract surrounding intermediates and carry out hydrogenation to generate methane. The catalyst is realized by assembling iron phthalocyanine on the copper surface, followed by in-situ formation of single iron atoms during electrocatalysis, identified using operando X-ray absorption spectroscopy. The copper-supported iron-single-atom catalyst exhibits a CO 2 -to-methane Faradaic efficiency of 64% and a partial current density of 128 mA cm –2 , while the nitrogen-doped graphene-supported one produces only CO. The activity is 32 times higher than a pristine copper under the same conditions of electrolyte and bias.

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

The striking influence of oxophilicity differences in heterometallic Mo–Mn oxide cluster reactions with water

Mixed-metal oxides have proven to be effective catalysts for the hydrogen evolution reaction, often outperforming either of the binary metal oxides. The reactivity of Mn x MoO y - (x = 1, 2; y = 3, 4) clusters toward H 2 O was investigated via time-of-flight mass spectrometry with clear evidence of cluster oxidation and corresponding H 2 production, specifically for Mn x MoO 3 - (x = 1, 2) clusters. Unlike previously studied MoxOy- clusters, which assumed a broad distribution of stoichiometries (typically x ≤ y ≤ 3x), both MnMoOy- and Mn 2 MoO y - preferentially formed y = 3 and 4 compositions in significant quantities under our source conditions. The electronic and molecular structures of the Mn x MoO y (x = 1, 2; y = 3, 4) anion and neutral clusters were probed with anion photoelectron spectroscopy and analyzed with supporting density functional theory calculations. Our studies suggest that both metal centers are involved in initial cluster–water complex formation, while Mo is the center that undergoes oxidation; hence, reactivity terminates when Mo is saturated in its highest oxidation state of +6. Across these four clusters, Mn remains relatively reduced and is stable in a high-spin electronic configuration. The preferential reactivity of water molecules toward the Mo center rather than Mn is rationalized by the much lower relative oxophilicity of Mn.

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

Chemical control of the Rashba spin splitting size of α-GeTe(111) surface states by adjusting the potential at the topmost atomic layer

As control of the Rashba spin splitting size is highly desirable for spintronic devices, intensive studies have been performed to vary the splitting size by, for example, applying an electric field or designing novel heterostructures. However, direct observation of Rashba spin splitting size change via spectroscopic measurements has not been done so far. Here, we report results of angle-resolved photoemission studies on ferroelectric α-GeTe(111). We observe that the Rashba splitting size of α-GeTe(111) surface states is reduced upon dosing with potassium (K) which has a very low electron affinity. Based on density functional theory calculations, we find that the electric potential energy at the topmost atomic layer and the surface potential energy barrier is reduced upon K dosing. Furthermore, this change in the surface potential induces both delocalization of the surface states and reduction of the splitting size. We expect to increase the splitting size by dosing other elements or molecules with high electron affinity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Evidence for nitrogen binding to surface defects for topological insulator Bi 2 Se 3

Using scanning tunneling spectroscopy and theoretical simulations we have studied the effects of nitrogen gas exposure on the electronic density of states of Bi 2 Se 3 , a well-studied topological insulator. In carefully controlled measurements, Bi 2 Se 3 crystals were initially cleaved in a helium gas environment and then exposed to a 22 SCFH flow of ultra-high purity N 2 gas. We observe a resulting change in the spectral curves, with the exposure effect saturating after approximately 50 min, ultimately bringing the Dirac point about 50 meV closer to the Fermi level. Furthermore, these results are compared to density functional theory calculations, which support a picture of N 2 molecules physisorbing near Se vacancies. Furthermore, ab initio molecular dynamics simulations aided by a Blue Moon ensemble method reveal the dissociative adsorption of N 2 molecules which then bind strongly to Se vacancies at the surface. In this scenario, the binding of the N atom to a Se vacancy site removes the surface defect state created by the vacancy and changes the position of the Fermi energy with respect to the Dirac point.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impacts of hydrogen bonding interactions with Np(v/vi)O 2 Cl 4 complexes: vibrational spectroscopy, redox behavior, and computational analysis

The neptunyl (Np(V)O 2 + /Np(VI)O 2 2+ ) cation is the dominant form of 237 Np in acidic aqueous solutions and the stability of the Np(V) and Np(VI) species is driven by the specific chemical constituents present in the system. Hydrogen bonding with the oxo group may impact the stability of these species, but there is limited understanding of how these intermolecular interactions influence the behavior of both solution and solid-state species. In the current study, we systematically evaluate the interactions between the neptunyl tetrachloride species and hydrogen donors in coordination complexes and in the related aqueous solutions. Both Np(V) compounds (N 2 C 4 H 12 ) 2 [Np(V)O 2 Cl 4 ]Cl (Np(V)pipz) and (NOC 4 H 10 ) 3 [Np(V)O 2 Cl 4 ] (Np(V)morph) exhibit directional hydrogen bonding to the neptunyl oxo group while Np(VI) compounds (NC 5 H 6 ) 2 [Np(VI)O 2 Cl 4 ] (Np(VI)pyr) and (NOC 4 H 10 ) 4 [Np(VI)O 2 Cl 4 ]·2Cl (Np(VI)morph) assemble via halogen interactions. Here, the Raman spectra of the solid-state phases indicate the activation of vibrational bands when there is asymmetry of the neptunyl bond, while these spectral features are not observed within the related solution phase spectra. Density functional theory calculations of the Np(V)pipz system suggest that activation of the ν 3 asymmetric stretch and other combination modes lead to additional complexity within the solid-state spectra. Electrochemical analyses of complexes in the solution phases are consistent with the results of the crystallization experiments as the voltammetric potentials of Np(V)/Np(VI) complexes in the presence of protonated heterocycles differ from the potentials of pure Np(V) and may correlate with the hydrogen bonding interactions.

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

Superconducting and antiferromagnetic properties of dual-phase V 3 Ga

The binary compound V 3 Ga can exhibit two near-equilibrium phases, the A15 structure that is superconducting and the Heusler D0 3 structure that is semiconducting and antiferromagnetic. Density functional theory calculations show that these two phases are nearly degenerate, being separated in energy by only ±10 meV/atom. Our magnetization measurements on bulk-grown samples show antiferromagnetism and superconducting behavior below 14 K. Overall, these results indicate the possibility of using V 3 Ga for quantum technology devices exploiting the co-existence of superconductivity and antiferromagnetism in a dual-phase material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electronic structure and magnetic tendencies of trilayer La 4 Ni 3 O 10 under pressure: Structural transition, molecular orbitals, and layer differentiation

Motivated by the recent observation of superconductivity in the pressurized trilayer Ruddlesden-Popper (RP) nickelate La 4 Ni 3 O 10 , we explore its structural, electronic, and magnetic properties as a function of hydrostatic pressure from first-principles calculations. We find that an orthorhombic (monoclinic)-to-tetragonal transition under pressure takes place concomitantly with the onset of superconductivity. The electronic structure of La 4 Ni 3 O 10 can be understood using a molecular trimer basis wherein n molecular subbands arise as the d z 2 orbitals hybridize strongly along the c axis within the trilayer. The magnetic tendencies indicate that the ground state at ambient pressure is formed by nonmagnetic inner planes and stripe-ordered outer planes that are antiferromagnetically coupled along the c axis, resulting in an unusual ↑, 0, ↓ stacking that is consistent with the spin density wave model previously suggested by neutron diffraction. Such a state is destabilized at the pressure where superconductivity arises. Despite the presence of d z 2 states at the Fermi level, the d x 2 –y 2 orbitals also play a key role in the electronic structure of La4Ni3O10. Finally, this active role of the d x 2 –y 2 states in the low-energy physics of the trilayer RP nickelate, together with the distinct electronic behavior of the inner and outer planes, resembles the physics of multilayer cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. Here, we have focused our workflow on TiH 2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

Realization of Organocerium-Based Fullerene Molecular Materials Showing Mott Insulator-Type Behavior

Abstract Electron-rich organocerium complexes (C 5 Me 4 H) 3 Ce and [(C 5 Me 5 ) 2 Ce(ortho-oxa)], with redox potentials E 1/2 = –0.82 V and E 1/2 = –0.86 V versus Fc/Fc + respectively, were reacted with fullerene (C 60 ) in different stoichiometries to obtain molecular materials. Structurally characterized co-crystals: [(C 5 Me 4 H) 3 Ce] 2 •C 60 (1) and [(C 5 Me 5 ) 2 Ce(ortho-oxa)] 3 •C 60 (2) of C 60 with cerium-based rare earth molecular precursors are reported for the first time. The extent of charge transfer in 1 and 2 was evaluated using a series of physical measurements: FT-IR, Raman, solidstate UV-vis-NIR spectroscopy, X-ray absorption near edge structure (XANES) spectroscopy, and magnetic susceptibility measurements. The physical measurements indicate that 1 and 2 comprise the cerium(III) oxidation state with formally neutral C 60 as a co-crystal in both cases. Pressure-dependent periodic density functional theory calculations were performed to study the electronic structure of 1. Inclusion of a Hubbard-U parameter removes Ce f states from the Fermi level, opens up a band gap, and stabilizes FM/AFM magnetic solutions that are isoenergetic because of the large distances between the Ce(III) cations. Furthermore, the electronic structure of this strongly correlated Mott insulator-type system is reminiscent of the well-studied Ce 2 O 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of Heat Treatments in H 2 and Ar on the E 1 center in β -Gs 2 O 3

In this report the influence of heat treating n-type bulk β-Ga 2 O 3 in hydrogen (H2) and argon (Ar) gases on the presence of the defect level commonly labeled as E 1 was studied. Fourier transform-infrared spectroscopy confirms that hydrogen (H) is incorporated into β-Ga 2 O 3 during H 2 annealing at 900 °C. Deep-level transient spectroscopy measurements reveal that the concentration of the E 1 level is promoted by the introduction of H, in contrast to what is observed in samples heat-treated in an Ar flow. We further find the E 1 level to be stable against heat treatments at 650 K, both with and without an applied reverse-bias voltage. Potential candidates for the defect origin of E 1 are investigated using hybrid-functional calculations, and three types of defect complexes involving H are found to exhibit charge-state transition levels compatible with E 1 , including substitutional H at one of the threefold coordinated O sites, Ga-substitutional shallow donor impurities passivated by H, and certain configurations of singly hydrogenated Ga–O divacancies. Among these types, only the latter exhibit H binding energies that are consistent with the observed thermal stability of E 1 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational catalyst discovery: Active classification through myopic multiscale sampling

We report the recent boom in computational chemistry has enabled several projects aimed at discovering useful materials or catalysts. We acknowledge and address two recurring issues in the field of computational catalyst discovery. First, calculating macro-scale catalyst properties is not straightforward when using ensembles of atomic-scale calculations [e.g., density functional theory (DFT)]. We attempt to address this issue by creating a multi-scale model that estimates bulk catalyst activity using adsorption energy predictions from both DFT and machine learning models. The second issue is that many catalyst discovery efforts seek to optimize catalyst properties, but optimization is an inherently exploitative objective that is in tension with the explorative nature of early-stage discovery projects. In other words, why invest so much time finding a “best” catalyst when it is likely to fail for some other, unforeseen problem? We address this issue by relaxing the catalyst discovery goal into a classification problem: “What is the set of catalysts that is worth testing experimentally?” Here, we present a catalyst discovery method called myopic multiscale sampling, which combines multiscale modeling with automated selection of DFT calculations. It is an active classification strategy that seeks to classify catalysts as “worth investigating” or “not worth investigating” experimentally. Our results show an ~7–16 times speedup in catalyst classification relative to random sampling. These results were based on offline simulations of our algorithm on two different datasets: a larger, synthesized dataset and a smaller, real dataset.

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

Ultra-large-scale ab initio quantum chemical computation of bio-molecular systems: The case of spike protein of SARS-CoV-2 virus

The COVID-19 pandemic poses a severe threat to human health with an unprecedented social and economic disruption. Spike (S) glycoprotein of the SARS-CoV-2 virus is pivotal in understanding the virus anatomy, since it initiates the first contact with the ACE2 receptor in the human cell. We report results of ab initio computation of the spike protein, the largest ab initio quantum chemical computation to date on any bio-molecular system, using a divide and conquer strategy by focusing on individual structural domains. In this approach we divided the S-protein into seven structural domains: N-terminal domain (NTD), receptor binding domain (RBD), subdomain 1 (SD1), subdomain 2 (SD2), fusion peptide (FP), heptad repeat 1 with central helix (HR1-CH) and connector domain (CD). The entire Chain A has 14,488 atoms including the hydrogen atoms but excluding the amino acids with missing coordinates based on the PDB data (ID: 6VSB). The results include structural refinement, ab initio calculation of intramolecular bonding mechanism, 3- dimensional non-local inter-amino acid interaction with implications for the inter-domain interaction. Details of the electronic structure, interatomic bonding, partial charge distribution and the role played by hydrogen bond network are discussed. In the interaction among structural domains, we present new insights for crucial hinge-like movement and fusion process. Extension of such calculation to the interface between the S-protein binding domain and ACE2 receptor can provide a pathway for computational understanding of mutations and the design of therapeutic drugs to combat the COVID-19 pandemic.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Mechanistic Insights and Rational Design of Ca-Doped CeO 2 Catalyst for Acetic Acid Ketonization

Carboxylic acid ketonization has recently gained significant attention to produce biomass-derived hydrocarbon fuels as it not only removes the highly reactive carboxylic functional group but also increases the size of the carbon chain. In this work, Ca-doped CeO 2 -based catalysts were investigated for acetic acid ketonization using a combined experimental and computational approach. Acetic acid conversion was performed across a range of temperatures including higher temperatures relevant to catalytic hot gas filtration (450 °C). Ca addition slightly decreases overall acetic acid ketonization reactivity yet stabilizes the catalyst at the higher temperatures necessary for catalytic hot gas filtration. From density functional theory calculations of the ketonization reaction mechanism, the C–C coupling and water formation steps are identified as two of the most energy-consuming steps on a CeO 2 surface with a proximal oxygen vacancy and the presence of a Ca dopant stabilizes the key intermediates. Calculations predict an optimal structure comprising three Ca ensembles to minimize the reaction free energies for C–C coupling and water formation steps. These findings provide a priori information to guide future experiments for ketonization catalyst design and development.

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

Density functional theory study on the coupling and reactions of diferuloylputrescine as a lignin monomer

Diferuloylputrescine has been found in a variety of plant species, and recent work has provided evidence of its covalent bonding into lignin. Results from nuclear magnetic resonance spectroscopy revealed the presence of bonding patterns consistent with homo-coupling of diferuloylputrescine and the possibility of cross-coupling with lignin. In the present work, density functional theory calculations have been applied to assess the energetics associated with radical coupling, rearomatization, and dehydrogenation for possible homo-coupled dimers of diferuloylputrescine and cross-coupled dimers of diferuloylputrescine and coniferyl alcohol. Furthermore, the values obtained for these reaction energetics are consistent with those reported for monolignols and other novel lignin monomers. As such, this study shows that there would be no thermodynamic impediment to the incorporation of diferuloylputrescine into the lignin polymer and its addition to the growing list of non-canonical lignin monomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗