Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electronic Structure 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Atomic-Layer-Deposited Aluminum Oxide Thin Films Probed with X-ray Scattering and Compared to Molecular Dynamics and Density Functional Theory Models

A better understanding of amorphous aluminum oxide’s structure and electronic properties is obtained through combined experimental and computational approaches. Grazing incidence X-ray scattering measurements were carried out on aluminum oxide thin films grown using thermal atomic layer deposition. The corresponding pair distribution functions (PDFs) showed structures similar to previously reported PDFs of solid-state amorphous alumina and molten alumina. Structural models based on crystalline alumina polymorphs (PDFgui) and amorphous alumina (molecular dynamics, MD) were examined for structural comparisons to the experimental PDF data. Smaller MD models were optimized and verified against larger models to allow for quantum chemical electronic structure calculations. The electronic structure of the amorphous alumina models yields additional insight into the band structure and electronic defects present in amorphous alumina that are not present in crystalline samples.

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

Rare Earth and Transition Metal Containing Glasses

Transition metal (TM) and rare earth (RE) ions have been incorporated into many glass systems such as silicate, phosphate, and borosilicate-based oxide glasses, as well as in halide and chalcogenide glasses, that find applications ranging from optical, photonic, and magnetic devices, solid-state battery, to nuclear waste disposal. Understanding the structural role of RE and TM in these glasses can help to develop glass compositions for targeted applications with either high-optical emission efficiency, electrical conductivity, or chemical durability. In this chapter, we first provide a general introduction of the applications and structural features of RE and TM in glasses, then the critical aspects of molecular dynamics (MD) simulations of these glasses such as interatomic potentials, structural analysis tools to study RE and TM ions in glasses and their clustering behaviors, Quantitative Structure–Property Analysis (QSPR), diffusion and dynamic property calculations, and electronic structure calculations to understand electronic defects such as charge trapping and radiation effects are introduced. Three representative case studies are presented: the first one is on MD simulations of erbium- and europium-doped silica and silicate glasses, as well as cerium doped aluminophosphate glasses, that revealed the effect of glass composition on RE ion local structure and clustering behavior. Electronic structure calculations of cerium-doped glass show how the existence of multioxidation states help to mediate radiation-induced damages caused by excited electron–hole pairs was also discussed. The second one focuses on alkali vanadophosphate glasses where the existence of two vanadium oxidation states help to provide electronic conduction in the glasses while alkali ions provide ionic conduction. MD simulations were used to understand vanadium environments and other structural aspects in the phosphate glasses, as well as the ionic transport behaviors of alkali ions. The third case study is on zirconium-containing borosilicate and aluminosilicate glasses that find wide applications in nuclear waste disposal. MD simulations help to provide structural details of zirconium ions that are validated by diffraction and EXAFS spectra. The structural information was used to interpret changes of mechanical properties and chemical durability by using QSPR and other analyses-based MD-generated structure models.

Du, Jincheng↗
DOE OSTI2022

A nearsighted force-training approach to systematically generate training data for the machine learning of large atomic structures

A challenge of atomistic machine-learning (ML) methods is ensuring that the training data are suitable for the system being simulated, which is particularly challenging for systems with large numbers of atoms. Most atomistic ML approaches rely on the nearsightedness principle (“all chemistry is local”), using information about the position of an atom’s neighbors to predict a per-atom energy. Here, in this work, we develop a framework that exploits the nearsighted nature of ML models to systematically produce an appropriate training set for large structures. We use a per-atom uncertainty estimate to identify the most uncertain atoms and extract chunks centered around these atoms. It is crucial that these small chunks are both large enough to satisfy the ML’s nearsighted principle (that is, filling the cutoff radius) and are large enough to be converged with respect to the electronic structure calculation. We present data indicating when the electronic structure calculations are converged with respect to the structure size, which fundamentally limits the accuracy of any nearsighted ML calculator. These new atomic chunks are calculated in electronic structures, and crucially, only a single force—that of the central atom—is added to the growing training set, preventing the noisy and irrelevant information from the piece’s boundary from interfering with ML training. The resulting ML potentials are robust, despite requiring single-point calculations on only small reference structures and never seeing large training structures. We demonstrated our approach via structure optimization of a 260-atom structure and extended the approach to clusters with up to 1415 atoms.

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

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2023

Roadmap on electronic structure codes in the exascale era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ThC 2 @C 82 versus Th@C 84 : unexpected formation of triangular thorium carbide cluster inside fullerenes

Synthesis of the first thorium-containing clusterfullerenes, ThC 2 @C s (6)–C 82 and ThC 2 @C 2 (5)–C 82 , is reported. These two novel actinide fullerene compounds were characterized by mass spectrometry, single-crystal X-ray diffraction crystallography, UV–vis–NIR spectroscopy, and theoretical calculations. Crystallographic studies reveal that the encapsulated ThC 2 clusters in both C s (6)–C 82 and C 2 (5)–C 82 feature a novel bonding structure with one thorium metal center connected by a C≡C unit, forming an isosceles triangular configuration, which has not been hitherto observed for endohedral fullerenes or for solid phase thorium carbides. Electronic structure calculations assign a formal electronic structure of [Th 4+ (C 2 ) 2- ] 2+ @[C 82 ] 2- , with pronounced donation bonding from (C 2 ) 2- to Th 4+ , secondary backbonding from the fullerene to thorium and Th–C double bond character in both compounds. This work presents a new family of endohedral fullerenes, MC 2 @C 2n-2 , being unexpected isomers of MC 2n , and provides broader understanding of thorium bonding.

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

Heteroanionic Ruddlesden-Popper ferroelectrics from anion order and octahedral tilts

We describe a strategy to design ferroelectric heteroanionic materials based on the coupling of anion order and octahedral tilts in n = 1 Ruddlesden-Popper structures, which allows noncentrosymmetric and polar compounds to arise from centrosymmetric anion-ordered structures. Here we investigate the relative phase stabilities of the polymorphs of Sr 2 ScO 3 F and Ca 2 ScO 3 F derived from this coupling using electronic structure calculations. We find that large degrees of octahedral tilting can stabilize different anion orders derived from assembly of [ScO 5 F] 8– octahedra. To further understand the link between octahedral tilting and anion order, we quantitatively separate the contributions of electrostatics and covalent interactions to the stability of tilts. We find that the tilts are driven primarily by covalent interactions and that local out-of-plane polar displacements, induced from the anion order, further stabilize the octahedral tilts through the pseudo-Jahn-Teller effect.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Flexibility of the factorized form of the unitary coupled cluster Ansatz

The factorized form of the unitary coupled cluster Ansatz is a popular state preparation Ansatz for electronic structure calculations of molecules on quantum computers. It is often viewed as an approximation (based on the Trotter product formula) for the conventional unitary coupled cluster operator. In this work, we show that the factorized form is quite flexible, allowing one to range from a conventional configuration interaction, to conventional unitary coupled cluster, to efficient approximations that lie in between these two. The variational minimization of the energy often allows simpler factorized unitary coupled cluster approximations to achieve high accuracy, even if they do not accurately approximate the Trotter product formula. This is similar to how quantum approximate optimization algorithms can achieve high accuracy with a small number of levels.

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

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Orbital hybridisation effects in B2 phase Cr doped Co 2 MnAl

Here we investigate here the magnetic, transport, local structural and electronic properties of Co 2 Mn 1−x Cr x Al (x = 0, 0.05, 0.1 and 0.2). Our results show that all the compounds stabilise in B2 phase and are ferromagnets. The results reveal disorder at the structural and magnetic levels. X-ray absorption near edge structure (XANES) analysis reveal signature of antisite disorder between Mn and Al atoms with equal ratio. The electronic structure calculations suggest enhancement in the half metallicity, localisation of electrons at the Fermi level and an increment in density of states with doping. The combined results of electronic structure calculations and XANES studies suggest transfer of electrons to the Co site. The results of high temperature resistivity measurements suggest the conduction electrons are undergoing transition from delocalisation to weak localisation to activated behaviour with Cr doping. The extended x-ray absorption spectroscopic analysis shows that the local structure around Mn atom is different from the global structure as obtained from the x-ray diffraction results. The behaviour of the edge region is in line with the trend as obtained from the compositional analysis. We observe link between the hybridisation of 3d like states at the Mn, Cr sites with that at the Co site and the transport properties. This could help in understanding the unusual decrement in the lattice parameter with doping. These results reveal the role of local structure in understanding the physical properties of such systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗
DOE OSTI2023

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2026

Physical properties investigation of tetragonal Ba T 2 P 2 (T = Ru, Pd)

The structural, thermal, electrical transport, and magnetic properties of ternary compounds BaT 2 P 2 (T = Ru, Pd) are reported together with electronic structure calculations. Our results show that BaRu 2 P 2 crystallises in the layered tetragonal ThCr 2 Si 2 -type structure (space group I4/mmm) whereas BaPd 2 P 2 adopts the primitive tetragonal CeMg 2 Si 2 -type structure (space group P4/mmm). The combined experimental results and electronic structure calculations suggest metallic and diamagnetic ground states for both materials. Electrical transport measurements further indicate non-Fermi liquid behaviour at low temperatures. Notably, BaRu 2 P 2 exhibits an anomalous feature in the magnetic susceptibility at T ~ 260 K, similar to that observed in BaRu 2 As 2 . A comparison between the bare density of states obtained from electronic structure calculations and the experimentally estimated values suggests the presence of strong correlation effects in the material. Furthermore, the heat capacity of BaPd 2 P 2 , calculated using the full phonon dispersion relation, is in good agreement with the experimental data, particularly at low temperatures, highlighting the significant contribution of optical phonon branches to C P in this temperature range.

flat bands↗
DOE OSTI2021

Electronic Structure Theory Calculations Using Modern Architectures: KNL vs Haswell

The time to solution and parallel efficiency of several commonly used electronic structure methods (Hartree–Fock, density functional theory, second order perturbation theory, resolution of the identity second order perturbation theory, coupled cluster) are evaluated on both the Intel Xeon Haswell and the Intel Xeon Phi Knights Landing (KNL) architectures. The Haswell completes the benchmark calculations with a faster time to solution than the KNL for all molecules and methods tested. While the Haswell exhibits an average speedup of at least 3.5 relative to the KNL for all nonthreaded computations, the KNL has a better parallel efficiency than the Haswell with increasing core counts. The architectures are further tested using a more computationally costly coupled cluster method on a transition state reaction. The Haswell appears to be the best choice to minimize the time to solution, though for very large systems and high levels of theory that require memory intensive processes the superior memory hierarchy and larger on node memory of the KNL can make it a better choice. These results are used to showcase aspects of novel architectures that will increase efficiency for quantum chemistry applications.

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

Hybrid functional calculations of electronic structure and carrier densities in rare-earth monopnictides

The structural parameters and electronic structure of rare-earth pnictides are calculated using density functional theory (DFT) with the Heyd, Scuseria, and Ernzerhof (HSE06) screened hybrid functional. Specifically, we focus on RE-V compounds, with RE=La, Gd, Er, and Lu, and V=As, Sb, and Bi, and analyze the effects of spin-orbit coupling and treating the RE4$\textit{f}$ electrons as valence electrons in the projector augmented wave approach. The results of HSE06 calculations are compared with DFT within the generalized gradient approximation (DFT-GGA) and other previous calculations. We find that all these RE-V compounds are semimetals with electron pockets at the $\textit{X}$ point and hole pockets at Γ. Whereas in DFT-GGA the carrier density is significantly overestimated, the computed carrier densities using HSE06 is in good agreement with the available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Doping topological Dirac semimetal with magnetic impurities: Electronic structure of Mn-doped Cd 3 ⁢As 2

The prospect of transforming a Dirac topological semimetal (TSM) into a Weyl TSM phase, following doping by magnetic impurities, is central to TSM applications. The magnetic field from polarized 𝑑 levels of magnetic impurities produces a field with a sharp local structure. To what extent magnetic impurities act in the same manner as an applied field and what are the effects of such a field on the electronic structure of a Dirac TSM are the subject of this paper. We present electronic structure calculations of bulk Cd 3 ⁢As 2 with substitutional doping of Mn impurities in the dilute alloy range. Quasiparticle 𝐺⁡𝑊 (QS⁢𝐺⁡𝑊) ab initio electronic structure calculations are used in conjunction with 𝑘 · 𝑝 model Hamiltonian calculations. As expected, we observe the splitting of the Dirac points into pairs of Weyl points following the doping with Mn. We also show that the electronic structure of Mn-doped Cd 3⁢ As 2 can be emulated by the electronic structure of pristine Cd 3⁢ As 2 with an appropriate external magnetic field. Some properties of the conductivity of bulk Cd 3 ⁢As 2 for different magnetic field orientations are also investigated. Our results inform future opportunities for unique device functionality based on band structure tuning not found in conventional magnetic Weyl TSM.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Charge Ordering and Low-Temperature Lattice Distortion in the {beta}'-(BDT-TTF){sub 2}CF{sub 3}CF{sub 2}SO{sub 3} Dime Mott Insulator.

We present single-crystal x-ray diffraction measurements, optical investigations, and electronic structure calculations for the organic charge-transfer salt beta'-(BEDT-TTF)(2)CF3CF2SO3 synthesized by electrocrystallization. Electronic structure calculations confirm the quasi-one-dimensional behavior of the compound and optical conductivity measurements reveal the dimer-Mott insulating nature of the system. The splitting of the charge-sensitive nu(2) mode in Raman spectra demonstrates the onset of an interlayer charge-ordered phase below 25 K, also suggested by the crystal structure considerations. This transition is accompanied by clear signatures of a lattice distortion in the BEDT-TTF donor layer, as shown by a splitting of the vibrational nu(3) mode in infrared spectra. At the same time, the sharp redshift of the nu(1) mode involving the BEDT-TTF ethylene groups strongly suggests a significant modification of the hydrogen-type bonding present between the BEDT-TTF donor layer and the CF3CF2SO3- anion layer. These observations point to a subtle interplay of charge and lattice degrees of freedom at the phase transition.

Olejniczak, I.↗