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 91 records · Page 5

One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography

Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reproducibility of calculations on Li species with correlation-consistent basis sets

Here we have found that two different basis set series, both labeled cc-pVxZ, are widely available and widely used for calculation on Li compounds such that calculations that should be the same are different, causing a reproducibility problem. Here we study the magnitude of the disagreement of these nominally-the-same basis sets for wave function and density functional calculations of Li clusters. We find deviations among the basis sets of as large as 22 millihartrees. In most cases one obtains lower energies with the noncanonical choices. Therefore, it is important to specify which version of the basis set is used for calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measurement of the $\textit{B}(\textit{E}2 ↑)$ strengths of 36 Ca and 38 Ca

In this work, the $B(E2, 0^+_1 → 2^+_1)$ strengths of 36 Ca and 38 Ca are measured to be 131(20) $e^2\text{fm}^4$ and 101(11) $e^2\text{fm}^4$, respectively. The $\textit{B}(\textit{E}2)$ value for 36 Ca required a measurement of the p/γ branching ratio because the 2 + state is proton unbound. This branching ratio is $B_p$=0.087(8). These $\textit{B}(\textit{E}2)$ and branching-ratio values can be reproduced in the shell-model with the ZMB2 interaction, an interaction that predicts the $\textit{Z} = 20 \textit{sd}$-shell closure is incomplete with large proton pf-shell occupancies in the ground state. These occupancies are at odds with other shell-model and energy-density-functional calculations of 36 Ca. New data are used to provide an update on constraints of the density dependence of the symmetry energy through mirror charge-radii differences as well as to help reduce uncertainties of the astrophysical important 35 K(p,γ) reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hybrid Functional Calculations for Antimony Doping in CdTe

CdTe-based solar cells are leading thin-film photovoltaic technology, with efficiencies over 22%, but still much lower than the theoretical maximum of 29%. Further improvements will rely on increasing the open-circuit voltage Voc, which, in turn, depends on carrier density and lifetime. Using hybrid density functional calculations, we investigate Sb doping of CdTe, focusing on its limitation as shallow acceptor and the formation of compensating AX center. Paying special attention to supercell size and effects of spin-orbit coupling, we predict an ionization energy of 116 meV in the dilute limit, much closer to recent experimental value of 103 meV from temperature-dependent Hall measurements in bulk single crystals, and in contrast to much larger values from previous calculations. We also find that the Sb-related AX centers are not major compensation centers in Sb-doped CdTe.

14 SOLAR ENERGY↗

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations

In this research, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast optical melting of trimer superstructure in layered 1T′-TaTe2

Abstract Quasi-two-dimensional transition-metal dichalcogenides are a key platform for exploring emergent nanoscale phenomena arising from complex interactions. Access to the underlying degrees-of-freedom on their natural time scales motivates the use of advanced ultrafast probes sensitive to self-organised atomic-scale patterns. Here, we report the ultrafast investigation of TaTe 2 , which exhibits unique charge and lattice trimer order characterised by a transition upon cooling from stripe-like chains into a (3 × 3) superstructure of trimer clusters. Utilising MeV-scale ultrafast electron diffraction, we capture the photo-induced TaTe 2 structural dynamics – exposing a rapid ≈ 1.4 ps melting of its low-temperature ordered state followed by recovery via thermalisation into a hot cluster superstructure. Density-functional calculations indicate that the initial quench is triggered by intra-trimer Ta charge transfer which destabilises the clusters, unlike melting of charge density waves in other TaX 2 compounds. Our work paves the way for further exploration and ultimately rapid optical and electronic manipulation of trimer superstructures.

36 MATERIALS SCIENCE↗

The deep-acceptor nature of the chalcogen vacancies in 2D transition-metal dichalcogenides

Abstract Chalcogen vacancies in the semiconducting monolayer transition-metal dichalcogenides (TMDs) have frequently been invoked to explain a wide range of phenomena, including both unintentional p-type and n-type conductivity, as well as sub-band gap defect levels measured via tunneling or optical spectroscopy. These conflicting interpretations of the deep versus shallow nature of the chalcogen vacancies are due in part to shortcomings in prior first-principles calculations of defects in the semiconducting two-dimensional TMDs that have been used to explain experimental observations. Here we report results of hybrid density functional calculations for the chalcogen vacancy in a series of monolayer TMDs, correctly referencing the thermodynamic charge transition levels to the fundamental band gap (as opposed to the optical band gap). We find that the chalcogen vacancies are deep acceptors and cannot lead to n-type or p-type conductivity. Both the (0/−1) and (−1/−2) transition levels occur in the gap, leading to paramagnetic charge states S = 1 / 2 and S = 1, respectively, in a collinear-spin representation. We discuss trends in terms of the band alignments between the TMDs, which can serve as a guide to future experimental studies of vacancy behavior.

36 MATERIALS SCIENCE↗

Microscopic Origin of Polarization Charges at Ga N / ( Al , Ga ) N Interfaces

Ga⁢N/(Al,Ga)⁢N heterojunctions are at the heart of high-electron-mobility transistors that are being adopted for high-power and high-frequency applications. The strong polarization fields present at this interface significantly enhance the density of the two-dimensional electron gas (2DEG) that is confined on the Ga⁢ N side of the junction. The microscopic origin of these electrons has been debated over the years: after excluding that they would be contributed by bulk donors, a model that identifies surface states on the (Al, Ga)⁢N surface as the source of electrons has become widely adopted. Recently, it has become clear, however, that the measured density of surface states is insufficient to account for the high electron density in the 2DEG. Here we demonstrate, based on state-of-the-art first-principles calculations, that the source of electrons is intrinsic to the overall structure and that the negative charge in the 2DEG is balanced by fixed charge on the surface. We perform a rigorous study of polarization, using our recently developed methodology for quantifying polarization fields within the finite-sized systems that can be addressed with density-functional calculations. The results show that the electrons that appear in the 2DEG originate locally at the interface, and that the net charge at the interface is predominantly compensated by fixed charge on the surface, rather than surface states. We elucidate the source of this fixed charge and associate it with surface reconstructions or the presence of heterovalent impurities (such as oxygen). Our results force a reassessment of the impact of surface states on the density of the 2DEG: rather than serving as a supply of electrons, the surface states mainly act to pin the Fermi level. Our conclusions allow a fresh interpretation of experimental observations and allow devising guidelines for optimizing carrier densities in the 2DEG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Theoretical Evaluation of Metal–Ligand Bonding in Neptunium Compounds in Relation to 237 Np Mössbauer Spectroscopy

The 237 Np Mössbauer isomer shift and quadrupole splitting (QS) are powerful probes for the metal–ligand bonding of neptunium, a 5f-element of vital importance in the nuclear fuel cycle. A large set of Np compounds with different oxidation states (III) to (VII) is studied to investigate, by first-principles calculations, isomer shifts and the QS trends in relation to the Np oxidation state. Natural Bond Orbital analysis reveals that in addition to donation bonding to the 5f shell, participation of the 6d and 7s neptunium shells in covalent (donation) bonding substantially impacts the isomer shifts. The isomer shift cannot be interpreted solely by the 5f shell electron count. The isomer shift for Np(II) compounds is estimated to be in the range of 31–34 mm/s, less positive than for Np(III) compounds. For the QS, density functional calculations fail to reproduce the quadrupole splitting for some Np(VI) ionic solids. A multiconfigurational wave function approach reproduces the observed QS trends. The calculations give a semiquantitative interpretation of the trends for Np oxidation states (V) to (VII). Furthermore, the contrasting QS for standard and “reverse” neptunyl(VI), at the opposite extremes of the observed QS scale, arises predominantly from the different crystal environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Monte Carlo Calculations of Chemical Binding and Reactions

The auxiliary field quantum Monte Carlo method developed by the PIs has been shown to provide the most accurate description of strongly correlated electronic systems, from molecules to solids. Unlike other explicitly many‐body approaches, the quantum Monte Carlo method scales as a low order polynomial of systems size, similar to mean‐field methods such as density functional theory. However, the auxiliary field quantum Monte Carlo algorithm is significantly more expensive than traditional density functional calculations. This creates a bottleneck for applications to extended systems, such as large molecules and solids. One principal objective of this proposal was to develop new auxiliary field quantum Monte Carlo computational strategies to achieve improved scaling with system size, using downfolding and localization schemes, without sacrificing the predictive power of the calculations. A second goal is to extend the reach of auxiliary field quantum Monte Carlo to calculate excited states. This final report summarizes what has been achieved during the course the project toward these goals.

97 MATHEMATICS AND COMPUTING↗

Electronic structure Fermi liquid theory of high T(sub c) superconductors: Comparison with experiments

For years, there has been controversy on whether the normal state of the Cu-oxide superconductors is a Fermi liquid or some other exotic ground state. However, some experimentalists are clarifying the nature of the normal state of the high T(sub c) superconductors by surmounting the experimental difficulties in producing clean, well characterized surfaces so as to obtain meaningful high resolved photoemission data, which agrees with earlier positron-annihilation experiments. The experimental work on high resolution angle resolved photoemission by Campuzano et al. and positron-annihilation studies by Smedskjaer et al. has verified the calculated Fermi surfaces in YBa2Cu3O7 superconductors and has provided evidence for the validity of the energy band approach. Similar good agreement was found for Bi2Sr2CaCu2O8 by Olson et al. As a Fermi liquid (metallic) nature of the normal state of the high T(sub c) superconductors becomes evident, these experimental observations have served to confirm the predictions of the local density functional calculations and hence the energy band approach as a valid natural starting point for further studies of their superconductivity.

Freeman, A. J.↗

Energy and mechanical properties predictions in Fe-Ni binary system by ab initio calculations

Here, the fcc random solid solutions in Fe-Ni binary system were investigated using ab initio calculations. Density functional theory (DFT) was applied along with supercells generated by the Special Quasi-Random Structure (SQS) method to predict the enthalpy of mixing and magnetic and elastic properties. The results were compared with relevant experimental and computational data produced by other groups. The effect of SQS supercell size and relaxation procedure on the predicted properties of the alloys were also studied. Our calculations show that 16-atoms SQS model provides adequate approach to the mixing energy of the alloy. However, calculation of mechanical properties requires larger supercell. Excellent agreement between experimental and calculated elastic properties was found in fcc region of Fe-Ni system.

36 MATERIALS SCIENCE↗

Operando Elucidation on the Working State of Immobilized Fluorinated Iron Porphyrin for Selective Aqueous Electroreduction of CO 2 to CO

Iron porphyrin-based molecular catalysts can electrocatalyze CO 2 reduction to CO at nearly 100% selectivity in water. Nevertheless, the associated active sites and reaction mechanisms remain debatable, impeding the establishment of design guidelines for effective catalysts. This study reports coupling in operando experiments and theoretical calculations for immobilized 5,10,15,20-tetrakis(pentafluorophenyl) porphyrin Fe(III) chloride (FeF 20 TPP) for electrocatalytic CO 2 reduction in an aqueous phase. In operando UV–vis and X-ray absorption near-edge structure spectra indicated the persisting presence of Fe(II) species during the cathodic reaction, acting as catalytic sites that accommodate CO as Fe(II)–CO adducts. Consistently, the density functional calculations pointed out that the ligand-reduced state with oxidized Fe, namely, [Fe(II)F 20 (TPP • )] - , prevails in the catalytic cycle prior to the rate-controlling step. This work provides the conclusive representation related to the working states of Fe-based molecular catalysts under reaction conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unoccupied electronic resonances of Sc adsorbed on W(001) by k-resolved inverse photoemission

Scandium adsorbed on the (001) face of tungsten has been studied using Auger-electron spectroscopy, low-energy electron diffraction, k-resolved inverse photoelectron spectroscopy, work-function measurements, and relativistic-electronic-structure calculations. We find that the work function of W(001) does not decrease monotonically as a function of scandium adsorption but reaches a minimum value of Phi = 3.25 eV at 50% coverage. For a complete monolayer, an order 1 X 1 Sc overlayer is formed. Inverse photoemission in the isochromat mode was used to map the unoccupied energy levels of this overlayer along the bar-Gamma(bar-Mu) symmetry direction of the surface Brillouin zone. Local-density-functional calculations using the muffin-tin orbital method were performed for a 1 X 1 Sc overlayer on W and are compared with the experimental two-dimensional band structure.

Lamouri, A.↗

Isotope effect and critical magnetic fields of superconducting YH 6 : A Migdal-Eliashberg theory approach

The emergence of near-ambient temperature superconductivity under pressure in metal hydride systems has motivated a desire to further understand such remarkable properties, specifically critical magnetic fields. Here, YH 6 is suggested to be a departure from conventional superconductivity, due to apparent anomalous behavior. Using density functional calculations in conjunction with Migdal-Eliashberg theory we show that in YH 6 the critical temperature and the isotope effect under pressure, as well as the high critical fields, are consistent with strong-coupling conventional superconductivity, a property anticipated to extend to other related systems. Furthermore, strong-coupling corrections occur to the expected BCS values for the isotope effect coefficient (α), Ginzburg-Landau parameter [κ 1 ⁡(T)], London penetration depth [λ L (T)], electromagnetic coherence length [ξ⁡(T)], and the energy gap (Δ 0 ).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Magnetic interactions and excitations in SrMnSb 2

The magnetic interactions in the antiferromagnetic (AFM) Dirac semimetal candidate SrMnSb 2 are investigated using ab initio linear response theory and inelastic neutron scattering (INS). Our calculations reveal that the first two nearest in-plane couplings (J 1 and J 2 ) are both AFM in nature, indicating a significant degree of spin frustration, which aligns with experimental observations. The orbital resolution of exchange interactions shows that J 1 and J 2 are dominated by direct and superexchange, respectively. In a broader context, a rigid-band model suggests that electron doping fills the minority spin channel and results in a decrease in the AFM coupling strength for both J 1 and J 2 . To better compare with INS measurements, we calculate the spin-wave spectra within a linear spin-wave theory, utilizing the computed exchange parameters. Although the calculated spin-wave spectra somewhat overestimate the magnon bandwidth, they exhibit overall good agreement with measurements from INS experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coexisting paramagnetic spins and long-range magnetic order in Ba4⁢(Ru0.92⁢Ir0.08)3⁢O10

We investigate the effect of dilute Ir substitution on the magnetism of the trimer-based ruthenate Ba4⁢Ru3⁢O10 using neutron diffraction, magnetic susceptibility measurements, first-principles calculations, and atomistic simulations. Neutron diffraction shows that Ir doping preserves the zigzag antiferromagnetic structure and the ordered-moment magnitude of the parent compound, in which the moments reside exclusively on the two outer Ru(2) sites of each Ru3⁢O12 trimer, while the central Ru(1) site remains nonmagnetic. The Néel temperature is reduced from ∼105 to 84.0(1) K upon 8% Ir substitution, while magnetic susceptibility reveals a pronounced low-temperature Curie-like upturn, indicating the coexistence of paramagnetic spins with long-range antiferromagnetic order. Density-functional calculations show that Ir preferentially occupies the central Ru(1) site, which disrupts the exchange pathways of the affected Ru3⁢O12 trimer state. Atomistic simulations incorporating this paramagnetic dilution reproduce the suppressed ordering temperature and the coexistence of ordered and paramagnetic components.

Islam, Farhan [Ames Laboratory]↗

Fine-tuning catalytic selectivity by modulating catalyst-environment interactions: CO 2 hydrogenation over Pd-based catalysts

Capturing catalytic behaviors under operational conditions is pivotal to gaining a mechanistic understanding and promoting the design of robust catalysts. The challenge lies in the difficulty of monitoring real-time surface dynamics driven by catalyst-environment interactions. Here, in this work, we introduce a framework based on density functional calculations and kinetic modeling. This framework significantly improves the accuracy of theoretical models’ descriptions of experimental observations by quantifying environmental impacts on surface phases and active sites. CO 2 hydrogenation over Pd-based catalysts is taken as a showcase. The observed selectivity variations of Pd and Pd-M bimetallic catalysts strongly correlate with hydrogen coverage maintained under typical CO 2 hydrogenation conditions. By reducing the amount of surface hydrogen, the selectivity tuned effectively from formic acid toward CO and methanol. This study not only deepens the comprehension of dynamics of active sites under active chemical conditions but also introduces an alternative opportunity for catalytic tuning by modulating catalyst-environment interactions.

03 NATURAL GAS↗