SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

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

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 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 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 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 OSTI2024

Photoluminescence of a Uranium(IV) Alkoxide Complex

In this report, we describe the photoluminescence of a homoleptic uranium(IV) alkoxide complex. Excitation of [Li(THF)] 2 [U IV (O t Bu) 6 ] leads to the first example of photoluminescence from a well-defined actinide complex originating from an f–f excitation, supported by second order multiconfigurational electronic structure calculations including spin–orbit coupling. These calculations show strong spin–orbit coupling between the excited triplet and singlet states for the 5f-orbital manifold, which leads to a long-lived excited state lifetime of 0.85 s at low temperature. The photophysical properties of homoleptic uranium(V) and uranium(VI) tertbutoxide complexes are also presented; we find that oxidation of the uranium(IV) alkoxide results in quenching of luminescence in [Li(THF)][U V (O t Bu) 6 ] and [U VI (O t Bu) 6 ]. This is attributed to competing ligand to metal charge transfer absorption processes shifted to lower energy upon oxidation of the actinide center, which mask the relevant f–f transitions in the visible region of the electronic absorption spectrum.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

New layered quaternary Zintl pnictide oxides Ba 2 Zn 2 Pn 2 O ( Pn = Sb, Bi): Discovery, crystal structures, band engineering, and transport properties

Three new heteroanionic oxypnictides, Ba 2 Zn 2 Sb 2 O, Ba 2 Zn 2 Bi 2 O, and the solid solution Ba 2 Zn 2 Sb 2−x Bi x O (x ≈ 1.1–1.6), have been synthesized and structurally characterized. They are isostructural with their Mn-bearing analog, adopting the Ba 2 Mn 2 Sb 2 O-type structure (space group P6 3 /mmc, No. 194), and feature a double-layered 2D $^{2}_{∞}$ [Zn 2 Pn 2 O] 2- substructure (Pn = Sb, Bi, Sb/Bi) composed of corner-sharing, distorted tetrahedral ZnPn 3 O units. Electronic structure calculations reveal a systematic progression from semiconducting Ba 2 Zn 2 Sb 2 O to metallic Ba 2 Zn 2 Bi 2 O as Bi content increases. These trends are corroborated by transport property measurements, with Ba 2 Zn 2 Sb 0.9(1) Bi 1.1 O exhibiting relatively low electrical resistivity, high Hall mobilities of ∼160 cm 2 /V·s, and large Seebeck coefficients from 69 to 132 μV K −1 over the 300–600 K temperature range. Comparison with structurally related Zintl pnictides, such as SrIn 2 As 2 and PrZn 3 As 3 phases, situates Ba 2 Zn 2 Pn 2 O (Pn = Sb, Bi) within a broader family of heteroanionic oxypnictide Zintl compounds, highlighting their structural flexibility and amenability to band engineering. Finally, electronic structure and bonding considerations point to tunable semiconducting behavior and underscore the relevance of these materials for thermoelectric and topological applications.

Band engineering↗

DOE OSTI2026

Anisotropic Ferromagnetism in CrAu 3 Sb 6

The crystal structure and properties of CrAu 3 Sb 6 are presented, determined by measurements on single crystal and polycrystalline samples and first-principles calculations. The trigonal structure (space group P$\overline{3}$1m) comprises a CdI 2 -like sublattice of AuSb 2 with Cr occupying octahedral holes in a fully ordered triangular array. It can be viewed as a variation of the interesting and well-known families of partially intercalated transition metal dichalcogenides, but with stronger interactions along the stacking direction evidenced by short Cr–Au distances. The compound is metallic and ferromagnetic with a Curie temperature of 164 K. A strong anomalous contribution to the Hall effect is seen in the ferromagnetic state, and quantum oscillations are observed in magnetization at 2 K. Magnetization measurements show that the ordered moments of 1.5 μ B per Cr are oriented along the c-axis with relatively strong magnetocrystalline anisotropy. Electronic structure calculations confirm this uniaxial anisotropy and the important role of spin–orbit coupling in CrAu 3 Sb 6 and reveal strongly favored ferromagnetic ground state consistent with the measured Curie temperature. Through combined experiment and theory, this work provides a detailed picture of the basic properties and behaviors of this uniquely structured, Cr-based, anisotropic ferromagnet.

Crystal structure↗

DOE OSTI2026

Complex and Unusual Excited-State Relaxation Dynamics of 9,9′-Bifluorenylidene Revealed by Comprehensive Time-Resolved Spectroscopy and MRSF-TDDFT Calculations

Design of novel photochemical molecular motors often requires molecular building blocks that exhibit rather unusual photoactivity, for which conventional analyses of spectroscopic data can lead to conflicting interpretations. We here systematically investigated the excited-state relaxation dynamics of one such molecule, 9,9′-bifluorenylidene (BF), through comprehensive and complementary integration of ultrafast transient absorption (TA) and femtosecond stimulated Raman (FSRS) spectroscopies and first-principles mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT). TA and FSRS identified two sequentially formed transients following photoexcitation. The decay kinetics of the two intermediates differ in response to excitation wavelengths and the viscosity/polarity of solvents. MRSF-TDDFT calculations reveal a direct, barrierless internal-conversion pathway from the bright Franck−Condon state to a dark S1 minimum, where the excited-state population is transiently trapped, accounting for the first transient species observed in spectroscopic experiments. Further tracking down along the PES with MRSF-TDDFT mapped out two nonradiative relaxation pathways via conical intersections that connect the dark S 1 state to three configurations in the ground-state manifolds, within which a ring structure with a C8−C8′ bond and the vibrationally excited ground-state BF were identified from spectroscopic and kinetic data. The complexity of relaxation kinetics was attributed to the flexible torsional and twisting motions about the C9−C9′ bridge bond enabled by the diradical character of the S 1 state. These findings clarify unusual photoactive relaxation dynamics stemming from a novel correlation between structural flexibility and shifting electronic characteristics, and they demonstrate the importance of integrating spectroscopic and advanced electronic structure calculation studies for judicious clarification of complex, competing relaxation pathways of excited states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗

DOE OSTI2024

Metadynamics investigation of lanthanide solvation free energy landscapes and insights into separations energetics

Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N′-dimethyl, N,N′-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.

Wang, Xiaoyu↗

DOE OSTI2025

Fragme∩t: An Open‐Source Framework for Multiscale Quantum Chemistry Based on Fragmentation

Fragment-based quantum chemistry offers a means to circumvent the nonlinear computational scaling of conventional electronic structure calculations, by partitioning a large calculation into smaller subsystems then considering the many-body interactions between them. Variants of this approach have been used to parameterize classical force fields and machine learning potentials, applications that benefit from interoperability between quantum chemistry codes. However, there is a dearth of software that provides interoperability yet is purpose-built to handle the combinatorial complexity of fragment-based calculations. To fill this void we introduce “Fragme∩t”, an open-source software application that provides a tool for community validation of fragment-based methods, a platform for developing new approximations, and a framework for analyzing many-body interactions. Fragme∩t includes algorithms for automatic fragment generation and structure modification, and for distance- and energy-based screening of the requisite subsystems. Checkpointing, database management, and parallelization are handled internally and results are archived in a portable database. Interfaces to various quantum chemistry engines are easy to write and exist already for Q-Chem, PySCF, xTB, Orca, CP2K, MRCC, Psi4, NWChem, GAMESS, and MOPAC. Applications reported here demonstrate parallel efficiencies around 96% on more than 1000 processors but also showcase that the code can handle large-scale protein fragmentation using only workstation hardware, all with a codebase that is designed to be usable by non-experts. Fragme∩t conforms to modern software engineering best practices and is built upon well established technologies including Python, SQLite, and Ray. The source code is available under the Apache 2.0 license.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Strong, Yet Split Hydrogen Bonding with Ice Rules in Delafossite (H/D)RhO 2

Despite remaining enigmatic, strong hydrogen bonding provides an advanced design handle for tailoring the properties of functional materials. Here, in this study, 3 R –(H/D)RhO 2 delafossites (prepared by ion exchange of Na + from NaRhO 2 ) contain H/D in linear coordination with O, linking Rh III O 2 layers. Bragg and real-space X-ray and neutron scattering analysis, vibrational and solid-state NMR spectroscopy, and density functional theory (DFT)–based electronic structure calculations have been employed to understand the nature of the hydrogen bonding. Despite short distances between H/D and the two O to which they are bonded, a clear double-minimum corresponding to a shorter and longer (H/D)–O distance is established. The triangular lattices formed by H/D appear to display ice-like disorder, corroborated by low-temperature heat capacity measurements.

NMR spectroscopy↗

DOE OSTI2024

Iodide Double Perovskites and the Limits of their Structural Stability

AbstractWhile halide double perovskites A2M(I)M(III)X6 have attracted significant attention, examples involving iodides are rare. We examine the limits of the structural stability of iodide double perovskites, presenting the synthesis and single‐crystal structures of Cs2NaScI6 and Cs2NaYI6. Bypassing the common expectation that iodides have small band gaps, these compounds display optical gaps of 3.10 eV [M(III) = Sc] and 3.65 eV [M(III) = Y]. Cs2NaScI6 is the only iodide double perovskite to exhibit a cubic crystal structure at room temperature. Density functional theory‐based electronic structure calculations help understand the role of competing Cs3M(III)2I9 (3 : 2 : 9) phases and provide possible reasons for why iodide double perovskites based around In(III), Sb(III) and Bi(III) cations have proved elusive. We confirm design rules for halide double perovskites based around concepts of the tolerance factor and the radius ratio of the smaller, trivalent ion, but also point to situations such as what is observed for Cs2NaScI6 where a double perovskite can be trapped in a metastable structure.

Mulligan, Anya S.↗

DOE OSTI2024

Runtime performance of a GAMESS quantum chemistry application offloaded to GPUs

Summary Computational chemistry is at the forefront of solving urgent societal problems, such as polymer upcycling and carbon capture. The complexity of modeling these processes at appropriate length and time scales is mainly manifested in the number and types of chemical species involved in the reactions and may require models of several thousand atoms and large basis sets to accurately capture the chemical complexity and heterogeneity in the physical and chemical processes. The quantum chemistry package General Atomic and Molecular Electronic Structure System (GAMESS) has a wide array of methods that can efficiently and accurately treat complex chemical systems. In this work, we have used the GAMESS Effective Fragment Molecule Orbital (EFMO) method for electronic structure calculation of a challenging mesoporous silica nanoparticle (MSN) model surrounded by about 4700 water molecules to investigate the strong scaling and GPU offloading on hybrid CPU‐GPU nodes. Experiments were performed on the Perlmutter platform at the National Energy Research Scientific Computing Center. Good strong scaling and load balancing have been observed on up to 88 hybrid nodes for different settings of the execution parameters for the calculation considered here. When GPUs are oversubscribed by offloading work from multiple CPU processes, using the NVIDIA multi‐process service (MPS) has consistently reduced time to solution and energy consumed. Additionally, for some configuration parameter settings, oversubscription with MPS improved performance by up to 5.8% over the case without oversubscription.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Influence of Peripheral Substituents on Fe(II) Spin State in Complexes with Tridentate Schiff‐Base Ligands

Abstract Six Fe(II) complexes with substituted 6‐((quinoline‐8‐ylimino)methyl)phenolates (R 2 qsal) were synthesized. X‐ray crystal structure determination for [Fe(R 2 qsal) 2 ], where R=F ( 1 ), Cl ( 2 ), Br ( 3 ), I ( 4 ), CF 3 /H ( 5 ), or CN/H ( 6 ), revealed mononuclear complexes that crystallize either solvent‐free in case of 1 , 5 , and 6 or as CH 2 Cl 2 solvates in the case of 2 , 3 , and 4 . Complexes 1 – 4 , which contain π‐donating substituents, exhibit only high‐spin (HS) state, as shown by crystal structure analysis and magnetic measurements. In contrast, complexes 5 and 6 , which contain electron‐withdrawing substituents, exist in the low‐spin (LS) state at lower temperatures but exhibit crossover to the HS state above 300 K. Electronic structure calculations show a good correlation between the relative energies of the LS and HS states and the experimentally observed behavior.

Gakiya‐Teruya, Miguel↗

DOE OSTI2025

Driving rapid atomic order in MnAl via low-magnitude magnetic field annealing

Application of a mild (60 mT), uniform magnetic field during short-term thermal treatment of kinetically retained, atomically disordered (paramagnetic) ε-MnAl was found to deliver a significant ~50 % increase in the formation of L1 0 atomically ordered (ferromagnetic) τ-MnAl product phase, compared to that produced by conventional (i.e., zero-field) annealing under identical thermal conditions. The magnetic field, applied in a passive closed-circuit configuration during annealing, induced significant changes in the structural, magnetic, and phase evolution of the material. Computational results based on electronic structure calculations demonstrate that the effective magnetic susceptibility of τ-MnAl is sensitive to the orientation, rather than the magnitude, of an applied magnetic field in the vicinity of the Curie temperature. The uniaxial magnetocrystalline anisotropy of the L1 0 structure is proposed to act as a filter for selective propagation of the population of τ-MnAl variants that are favorably aligned with the applied field. In this manner, crystallographic “gridlock” is alleviated that would otherwise arise from the coexistence of multiple, energetically equivalent τ-phase variants within the parent ε-phase matrix. These results confirm that static, low-magnitude magnetic field annealing is able to accelerate L1 0 atomic ordering in the MnAl system and likely can exert similar influences in relevant magnetic systems, facilitating efficient tailoring of structure-sensitive magnetic properties for the manufacture of magnetic materials.

36 MATERIALS SCIENCE↗