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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Solid Electrolyte Interphase Formation at the Ionic Liquid Electrolyte—Lithium-Metal Interface Using an Ab Initio Molecular Dynamics Approach

An interfacial study is performed using ab initio molecular dynamics (AIMD) simulations to elucidate the electrochemical interfacial stability of an ionic liquid (IL) as component of an electrolyte in contact with a Li-metal anode when an additive, vinylene carbonate (VC), is added to the electrolyte. The IL electrolyte is composed of 1-butyl-1-methylpyrrolidinium bis(fluoro-sulfonyl)-imide (PYR14 + FSI − ) with lithium bis(trifluoromethane-sulfonyl)imide salt (Li + TFSI − ) in a 8:2 concentration ratio. Reactions observed at the interface electrolyte/Li-metal during the AIMD simulations are verified by electronic structure calculations using density functional theory calculations considering a continuum solvent model. We found that the addition of VC to the IL prevents further reduction of the counter-ion TFSI − and retards reduction of FSI − in comparison to the case when the additive VC is absent. The PYR14 + reduction is never observed with or without the additive. The 1st, 2nd, and 3rd electron affinities of the PYR14 + are −0.08, −0.36 and −3.65 eV, indicating little to large repulsive behavior against electrons. The opposite is observed for the two negative ions, FSI − and TFSI − , which yield 1st electron affinities of 1.73 and 1.83 eV, respectively. Reactions obtained under open circuit conditions. Dissociation mechanisms are presented for all IL components reduced at the IL/Li-metal interface.

Electrochemistry↗
DOE OSTI2025

Molecular Origins of Temperature-Dependent X-ray Absorption in YAG:Dy

The absorption and emission of X-rays in dysprosium-doped yttrium aluminum garnet (YAG:Dy) has produced unexpected thermographic behavior, which is investigated using a combination of finite temperature ab initio molecular dynamic simulations, structural characterization, and electronic structure calculations of X-ray characteristics. Calculated average peak X-ray absorption spectra (XAS) from simulations between 300 and 600 K result in peak intensity loss due to thermalization effects, matching experimentally measured behavior of YAG:Dy. Here, investigation of atomic snapshots indicates structural factors that correlated with the X-ray behavior, with the first Y–O coordination sphere identified as the primary structural feature unique to high XAS intensity as calculated by radial and pair distribution functions.

Ab initio molecular dynamics↗
DOE OSTI2024

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Low Thermal Conductivity in Heteroanionic Materials with Layers of Homoleptic Polyhedra

Although BiAgOSe, an analogue of a well-studied thermoelectric material BiCuOSe, is thermodynamically stable, its synthesis is complicated by the low driving force of formation from the stable binary and ternary intermediates. Here we have developed a “subtraction strategy” to suppress byproducts and produce pure phase BiAgOSe using hydrothermal methods. Electronic structure calculations and optical characterization show that BiAgOSe is an indirect bandgap semiconductor with a bandgap of 0.95 eV. The prepared sample exhibits lower lattice thermal conductivities (0.61 W·m –1 ·K –1 at room temperature and 0.35 W·m –1 ·K –1 at 650 K) than BiCuOSe. Further, lattice dynamical simulations and variable temperature diffraction measurements demonstrate that the low lattice thermal conductivity arises from both the low sound velocity and high phonon–phonon scattering rates in BiAgOSe. These in turn result primarily from the soft Ag–Se bonds in the edge-sharing AgSe 4 tetrahedra and large sublattice mismatch between the quasi-two-dimensional [Bi 2 O 2 ] 2+ and [Ag 2 Se 2 ] 2– layers. These results highlight the advantages of manipulating the chemistry of homoleptic polyhedra in heteroanionic compounds for electronic structure and phonon transport control.

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

Electronic and magnetic properties of iridium-based novel Heusler alloys

We report half-metallicity and magnetism including exchange splitting are the most significant physical parameters to predict and design a candidate material for spintronic applications. We report here an ab-initio investigation on chemical formation and dynamical stability along with electronic structure and magnetic properties of Ir 2 Cr (Si, Ge) and IrRhCr (Si, Ge) Heusler alloys. The negative formation and cohesive energies with positive phonon dispersions confirm the stabilities of these alloys. Electronic structure calculations reveal that Ir 2 Cr (Si, Ge) and IrRhCrSi alloys are half-metallic ferromagnets with unprecedented exchange splitting. In addition, IrRhCrGe also shows semi-metallic nature. All of these materials follow Slater Pauling rule with large magnetic moments and 100% spin-polarization. With Cr bearing the majority of the local magnetic moment and exchange splitting, a ferromagnetic state is more stable than a nonmagnetic state. The electronic charge distribution and population analysis confirm mixed ionic and covalent bonding. The magnetocrystalline anisotropy energy, with the easy magnetization along the [1 1 1] direction, is significantly high in Ir 2 CrGe. Elastic constants such as shear (G), bulk (B), Young’s moduli, and Poisson’s ratio indicate that the IrRhCrSi and IrRhCrGe alloys are mechanically stable, and Ir 2 CrSi and Ir 2 CrGe are mechanically unstable. The Pugh’s (B/G) and Poisson’s ratios confirm that the stable alloys are ductile.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular dynamics studies of knotted polymers

Molecular dynamics calculations have been used to explore the influence of knots on the strength of a polymer strand. In particular, the mechanism of breaking 31, 41, 51, and 52 prime knots has been studied using two very different models to represent the polymer: (1) the generic coarse-grained (CG) bead model of polymer physics and (2) a state-of-the-art machine learned atomistic neural network (NN) potential for polyethylene derived from electronic structure calculations. While there is a broad overall agreement between the results on the influence of the pulling rate on chain rupture based on the CG and atomistic NN models, for the simple 31 and 41 knots, significant differences are found for the more complex 51 and 52 knots. Notably, in the latter case, the NN model more frequently predicts that these knots can break not only at the crossings at the entrance/exit but also at one of the central crossing points. The relative smoothness of the CG potential energy surface also leads to stabilization of tighter knots compared to the more realistic NN model.

DelloStritto, Mark (ORCID:0000000206785860)↗
DOE OSTI2020

Resonant and bound states of charged defects in two-dimensional semiconductors

A detailed understanding of charged defects in two-dimensional semiconductors is needed for the development of ultrathin electronic devices. Here, we study negatively charged acceptor impurities in monolayer WS 2 using a combination of scanning tunnelling spectroscopy and large-scale atomistic electronic structure calculations. We observe several localized defect states of hydrogenic wave function character in the vicinity of the valence band edge. Some of these defect states are bound, while others are resonant. The resonant states result from the multi-valley valence band structure of WS$_2$, whereby localized states originating from the secondary valence band maximum at Γ hybridize with continuum states from the primary valence band maximum at K/K$^{\prime}$. Resonant states have important consequences for electron transport as they can trap mobile carriers for several tens of picoseconds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Lattice vacancy migration barriers in Fe-Ni alloys, and an indication as to why Ni atoms diffuse slowly: A first-principles study

Lattice vacancy migration barriers in ferromagnetic Fe 𝑥 ⁢Ni 1−𝑥 alloys (0.4 ≤ 𝑥 ≤ 0.6) are accurately quantified within the framework of ab initio electronic structure calculations using the nudged elastic band (NEB) method. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal L⁢1 0 phase—which is under consideration as a material for a rare-earth-free gap magnet for advanced engineering applications—are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni-vacancy interchanges encounter significantly higher energetic barriers than do Fe-vacancy interchanges. We contend that this aspect is a key factor in determining the differences in mobility between Fe and Ni atoms in this ferromagnetic alloy. Moreover, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarized electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarization around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. These results give atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.

density functional theory↗
DOE OSTI2020

Tunable Perovskite-Derived Bismuth Halides: Cs 3 Bi 2 (Cl 1– x I x ) 9

Bismuth-based perovskites are of interest as safer alternatives to lead-based optoelectronic materials. Prior studies have reported on the compounds Cs 3 Bi 2 Cl 9 , Cs 3 Bi 2 I 9 , and Cs 3 Bi 2 Cl 3 I 6 . Here we examine a range of compounds of the formula Cs 3 Bi 2 (Cl 1–x I x ) 9 , where x takes values from 0.09 to 0.52. Powder and single-crystal X-ray diffraction were used to determine that all of these compounds adopt the layered vacancy-ordered perovskite structure observed for Cs 3 Bi 2 Cl 3 I 6 , which is also the high-temperature phase of Cs 3 Bi 2 Cl 9 . We find that, even with very small iodine incorporation, the structure is switched to that of Cs 3 Bi 2 Cl 3 I 6 , with I atoms displaying a distinct preference for the capping sites on the BiX 6 octahedra. Optical absorption spectroscopy was employed to study the evolution of optical properties of these materials, and this is complemented by density functional theory electronic structure calculations. Finally, three main absorption features were observed for these compounds, and with increasing x, the lowest-energy features are red-shifted.

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

Noise robust detection of quantum phase transitions

Quantum computing allows for the manipulation of highly correlated states whose properties quickly go beyond the capacity of any classical method to calculate. Thus one natural problem which could lend itself to quantum advantage is the study of ground-states of condensed matter models, and the transitions between them. However, current levels of hardware noise can require extensive application of error-mitigation techniques to achieve reliable computations. In this work, we use several IBM devices to explore a finite-size spin model with multiple “phaselike” regions characterized by distinct ground-state configurations. Using preoptimized Variational Quantum Eigensolver (VQE) solutions, we demonstrate that in contrast to calculating the energy, where zero-noise extrapolation is required in order to obtain qualitatively accurate yet still unreliable results, calculations of the energy derivative, two-site spin correlation functions, and the fidelity susceptibility yield accurate behavior across multiple regions, even with minimal or no application of error-mitigation approaches. Taken together, these sets of observables could be used to identify level crossings in a simple, noise-robust manner which is agnostic to the method of ground state preparation. This work shows promising potential for near-term application to identifying quantum phase transitions, including avoided crossings and nonadiabatic conical intersections in electronic structure calculations. Published by the American Physical Society 2024

Lively, Kevin (ORCID:0000000320981494)↗
DOE OSTI2022

Combining Physics and Machine Learning for the Next Generation of Molecular Simulation

Simulating molecules and atomic systems at quantum accuracy is a grand challenge for science in the 21 st century. Quantum-accurate simulations would enable the design of new medicines and the discovery of new materials. The defining problem in this challenge is that quantum calculations on large molecules, like proteins or DNA, are fundamentally impossible with current algorithms. In this work, we explore a range of different methods that aim to make large, quantum-accurate simulations possible. We show that using advanced classical models, we can accurately simulate ion channels, an important biomolecular system. We show how advanced classical models can be implemented in an exascale-ready software package. Lastly, we show how machine learning can learn the laws of quantum mechanics from data and enable quantum electronic structure calculations on thousands of atoms, a feat that is impossible for current algorithms. Altogether, this work shows that combining advances in physics models, computing, and machine learning, we are moving closer to the reality of accurately simulating our molecular world.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Elucidating Biomass-Derived Pyrolytic Lignin Structures from Demethylation Reactions through Density Functional Theory Calculations

Pyrolytic lignin is a fraction of pyrolysis oil that contains a wide range of phenolic compounds that can be used as intermediates to produce fuels and chemicals. However, the characteristics of the raw lignin structure make it difficult to establish a pyrolysis mechanism and determine pyrolytic lignin structures. Herein this study proposes dimer, trimer, and tetramer structures based on their relative thermodynamic stability for a hardwood lignin model in pyrolysis. Different configurations of oligomers were evaluated by varying the positions of the guaiacyl (G) and syringyl (S) units and the bonds βO4 and β5 in the hardwood model lignin through electronic structure calculations. The homolytic cleavage of βO4 bonds is assumed to occur and generate two free radical fragments. These can stabilize by taking hydrogen radicals that may be in solution during the intermediate liquid (pathway 1) formation before the thermal ejection. An alternative pathway (pathway 2) could occur when the radicals use intramolecular hydrogen, turning themselves into stable products. Subsequently, a demethylation reaction can take place, thus generating a methane molecule and new oligomeric lignin-derived molecules. The most probable resulting structures were studied. We used FTIR and NMR spectra of selected model compounds to evaluate our calculation approach. Thermophysical properties were calculated using group contribution methods. The results give insights into the lignin oligomer structures and how these molecules are formed. They also provide helpful information for the design of pyrolysis oil separation and upgrading equipment.

09 BIOMASS FUELS↗
DOE OSTI2021

Angstrom-Resolved Interfacial Structure in Buried Organic-Inorganic Junctions

Charge transport processes at interfaces which are governed by complex interfacial electronic structure play a crucial role in catalytic reactions, energy storage, photovoltaics, and many biological processes. Here, the first soft X-ray second harmonic generation (SXR-SHG) interfacial spectrum of a buried interface (boron/Parylene-N) is reported. SXR-SHG shows distinct spectral features that are not observed in X-ray absorption spectra, demonstrating its extraordinary interfacial sensitivity. Comparison to electronic structure calculations indicates a boron-organic separation distance of 1.9 {\AA}, wherein changes as small as 0.1 {\AA} result in easily detectable SXR-SHG spectral shifts (ca. 100s of meV). As SXR-SHG is inherently ultrafast and sensitive to individual atomic layers, it creates the possibility to study a variety of interfacial processes, e.g. catalysis, with ultrafast time resolution and bond specificity.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2026

Chemical bond and phase stability of Ga-doped Sm2Fe17Cx magnet

Sm2Fe17C3 phase (2:17) is metastable and exhibits excellent intrinsic hard magnetic properties. Doping elements such as Ga facilitate the formation of a single-phase 2:17 structure in arc-melted Sm2Fe17Cx alloys, which opens a promising route for fabricating fully dense bulk Sm2Fe17Cx magnets via high-temperature techniques such as melting and sintering. First-principles electronic structure calculation indicates that Ga prefers to partially replace Fe at the 9d and 18h crystallographic sites in Sm2Fe17C3 and Sm2Fe17, respectively. This difference in site preference is attributed to the distinct chemical environments surrounding the Fe atoms in the two compounds. Ga substitution favors the Sm–Ga bonding formation while avoiding Ga–C interactions. Doped Ga atoms result in more negative formation energy in Sm2(Fe, Ga)17C3, indicating improved structural stability. Crystal Orbital Hamilton Population analysis reveals that carbon insertion weakens the bonding of Sm-Fe (18h) and Sm-Fe (18f) in Sm2Fe17C3. Ga doping facilitates electron redistribution across chemical bonds, thereby reinforcing Fe(18h)–Sm and Fe(18f)–Sm interactions and stabilizing the carbon-centered octahedral local structure. This synergistic effect contributes significantly to the observed enhancement in phase stability of Sm2(Fe, Ga)17Cx. These findings suggest that chemical bond engineering through the selective doping of Ga can enhance phase stability and facilitate the synthesis of Sm2Fe17C3, providing a viable strategy for developing advanced magnets.

Liu, Xubo [Critical Materials Innovation Hub, Divi↗
DOE OSTI2022

Post-Marcus Theory and Simulation of Interfacial Charge Transfer Dynamics in Organic Semiconducting Materials (Final Report)

The research program developed, validated and applies predictive computational protocols for calculating charge transfer (CT) rates in complex molecular systems, including molecular dyads and triads in liquid solution and solid-state organic semiconducting (OSC) materials. We have established a transformative computational scheme that goes beyond widely used simplifications, to achieve realistic descriptions of CT processes. The approach properly addresses the contribution of molecular environment at ambient conditions to CT processes. Our approach achieves unique insight on investigated CT processes in relevant experimental efforts. The collaborative team included three principal investigators (PIs), with complimentary expertise in classical molecular dynamics simulations and data science (Cheung), state-of-the-art electronic structure calculations (Dunietz), and cutting-edge theory and simulation techniques for modeling energy, charge and coherence transfer dynamics in molecular systems (Geva).

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

Revisiting the π → π* transition of the nitrite ion at the air/water interface: A combined experimental and theoretical study

Broadband deep ultraviolet electronic sum frequency generation spectroscopy allows measurement of electronic |χ (2) | 2 -spectra at aqueous interfaces, providing considerable improvement over deep ultraviolet electronic second harmonic generation spectroscopy that yields error-prone, pointwise spectra. Reexamination of the π → π* transition of nitrite at the air/water interface reveals that the interfacial |χ (2) | 2 -spectrum is strikingly similar to the bulk absorption spectrum. Molecular dynamics simulations and SOS-CIS(D 0 ) electronic structure calculations provide no evidence for the previously reported contact ion pair-induced red-shift of the π → π* transition at the air/water interface. Finally, our results thus revise the previous description for nitrite interfacial adsorption as a contact ion pair.

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

An efficient computational framework for charge density estimation in twisted bilayer graphene

Electronic properties such as band structure and Fermi velocity in low-angle twisted bilayer graphene (TBG) are intrinsically dependent on the atomic structure. Rigid rotation between individual graphene layers provides an approximate description of the bilayer symmetry. Upon relaxation, in-plane displacement of the atoms in low angle TBG causes a change in the symmetry through the enlargement of the AB stacking regions and the reduction in size of AA and SP stacking regions. However, the effect of this in-plane relaxation on the charge density remains unexplored, because the necessary electronic structure calculations of such large supercells of low twist angle TBG are computationally infeasible. Therefore, we develop a computationally efficient framework that enables the exploration of the charge density symmetry of the low twist angle TBG. This framework is based on the Fourier representation of the charge density which presents high intensity Bragg peaks. Here we find that with the decrease of twist angle, low intensity satellite peaks also become apparent. Our framework incorporates these satellite peaks which reveals transformation of symmetry in the charge density distribution from high to low twist angle TBG. One striking outcome is the demonstration of the electron localization in the AA region of low twist angle TBG. Our framework helps to explain the effect of the atomistic relaxation on the charge density distribution and thus, it provides information about exotic electronic properties of low twist angle TBG at a low computational expense.

36 MATERIALS SCIENCE↗