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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Altermagnetism in the layered intercalated transition metal dichalcogenide CoNb 4 Se 8

Altermagnets (AMs) are a new class of magnetic materials that combine the beneficial spintronics properties of ferromagnets and antiferromagnets, garnering significant attention recently. Here, we have identified altermagnetism in a layered intercalated transition metal diselenide, CoNb 4 Se 8 , which crystallizes with an ordered sublattice of intercalated Co atoms between NbSe 2 layers. Single crystals are synthesized, and the structural characterizations are performed using single crystal diffraction and scanning tunneling microscopy. Magnetic measurements reveal easy-axis antiferromagnetism below 168 K. Density functional theory (DFT) calculations indicate that A-type antiferromagnetic ordering with easy-axis spin direction is the ground state, which is verified through single crystal neutron diffraction experiments. Electronic band structure calculations in this magnetic state display spin-split bands, confirming altermagnetism in this compound. The layered structure of CoNb 4 Se 8 presents a promising platform for testing various predicted properties associated with altermagnetism.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of symmetry breaking and spin-orbit coupling on the band gap of halide perovskites

Halide perovskite (HP) materials have recently emerged as a class of semiconductors with immense promise for various optoelectronic applications, ranging from solar cells to light-emitting diodes. One of the unique attributes of HPs is their tunable band gaps with different factors governing their value. Furthermore, the first factor is related to relativistic corrections [“mass-Darwin,” connected to the ns 2 lone pairs, and spin-orbit coupling (SOC)] that induce an orbital shift or degeneracy splitting, resulting in a band-gap reduction. The second factor involves the structural configuration: in HPs the local symmetry of each Wyckoff position tends to be broken, inducing an opening of the band gap. Based on high-throughput density functional theory calculations, this paper systematically studies a possible self-cancelation on the band-gap correction for HPs when the polymorphous configuration—structural effects—and the SOC—electronic effects—are included. Our results indicate that the nature of interplay between SOC and symmetry breaking (SB) is that they are independent decoupling effects to describe the band-gap magnitude in halide perovskites. As a result of that, we observe a transitivity of the band-gap description; i.e., if we know the band gap of halide perovskites without SB and SOC, we can independently add the effects of band-gap reduction due to SOC and band-gap opening due to SB, regardless of the order in which these effects are considered.

14 SOLAR ENERGY↗
DOE OSTI2024

When van der Waals Met Kagome: A 2D Antimonide with a Vanadium-Kagome Network

2D materials showcase unconventional properties emerging from quantum confinement effects. In this work, a “soft chemical” route allows for the deintercalation of K + from the layered antimonide KV 6 Sb 6 , resulting in the discovery of a new metastable 2D-Kagome antimonide K 0.1(1) V 6 Sb 6 with a van der Waals gap of 3.2 Å. The structure of K 0.1(1) V 6 Sb 6 was determined via the synergistic techniques, including X-ray pair distribution function analysis, advanced transmission electron microscopy, and density functional theory calculations. The K 0.1(1) V 6 Sb 6 compound crystallizes in the monoclinic space group C2/m (a = 9.57(2) Å, b = 5.502(8) Å, c = 10.23(2) Å, β = 97.6(2)°, Z = 2). The [V 6 Sb 6 ] layers in K 0.1(1) V 6 Sb 6 are retained upon deintercalation and closely resemble the layers in the parent compound, yet deintercalation results in a relative shift of the adjacent [V 6 Sb 6 ] layers. The magnetic properties of the K 0.1(1) V 6 Sb 6 phase in the 2–300 K range are comparable to those of KV 6 Sb 6 and another Kagome antimonide KV 3 Sb 5 , consistent with nearly temperature-independent paramagnetism. Electronic band structure calculation suggests a nontrivial band topology with flat bands and opening of band crossing afforded by deintercalation. Transport property measurements reveal a metallic nature for K 0.1(1) V 6 Sb 6 and a low thermal conductivity of 0.6 W K –1 m –1 at 300 K. Additionally, ion exchange in KV 6 Sb 6 via a solvothermal route leads to a successful partial exchange of K + with A + (A = Na, Rb, and Cs). Here, this study highlights the tunability of the layered structure of the KV 6 Sb 6 compound, providing a rich playground for the realization of new 2D materials.

2D↗
DOE OSTI2026

First principles band structure of interacting phosphorus and boron/aluminum δ -doped layers in silicon

Silicon can be heavily doped with phosphorus in a single atomic layer (a δ layer), significantly altering the electronic structure of the conduction bands within the material. Recent progress has also made it possible to further dope silicon with acceptor-based δ layers using either boron or aluminum, making it feasible to create devices with interacting δ layers with opposite polarity. Using density functional theory, we calculate the electronic structure of a phosphorus-based δ layer interacting with a boron or aluminum δ layer, varying the distances between the δ layers. At separations 1 nm and smaller, the dopant potentials overlap and largely cancel each other out, leading to an electronic structure closely mimicking intrinsic silicon. At separations greater than 1 nm, the two δ layers behave independently of one another, with an equivalent electronic structure to a p–n diode with an intrinsic layer taking the place of the depletion region. One mechanism for charge transfer between δ layers at larger distances could be tunneling, where we see a tunneling probability exceeding what would be seen for a standard silicon 1.1 eV triangular barrier, indicating that the interaction between delta layers may enhance tunneling compared to a traditional junction.

delta layers↗
DOE OSTI2026

Gaunt and Breit two-electron contributions to mean-field transformations and fine structure splitting

Materials utilized by novel energy systems are often studied using weakly correlated mean-field theories. However, if these systems incorporate heavy elements, relativistic effects must be included. Therefore, a Kramers unrestricted coupled cluster with singles and doubles excitation formalism within a molecular mean-field exact two-component framework (X2C mmf ) using a four-component Dirac–Hartree–Fock (DHF) reference state is presented. The exact X2C mmf transformed normal-order Hamiltonian incorporates all one-electron and two-electron (2e) contributions from the Coulomb, Gaunt, and Breit operators and is used with the equation of motion method to calculate the excitation energies of the alkali group of elements. Using this framework, the effects of 2e Gaunt and Breit integrals are studied. Results demonstrate growing contributions from these integrals to the generated X2C mmf mean-fields and electronic fine structure calculations with increasing atomic number. Overall, this paper outlines the method, its effect within the X2C mmf approach, and lays the foundation for future theoretical development of relativistic calculations within this framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Projector‐Based Quantum Embedding Study of Iron Complexes

Projection‐based embedding theory (PBET) is used to calculate and assess the challenging spin‐crossover energies for a selection of small Fe‐containing systems by embedding the metal center into the frozen potential of the ligands. MP2, CCSD, and CCSD(T) are embedded in potentials from the SCAN and r 2 SCAN functionals and compared with the canonical values for the constituent methods and previously reported reference values. Considering the PBET calculations as a correction for the underlying DFT, the embedding calculations are able to provided improvement for most cases. In some cases, the PBET methods are able to compensate for limitations in the wave function methods and produce results similar to more rigorous calculations from the literature. For the systems with spin‐crossover energies near zero, the current methodology fails to provide consistent improvement. In conclusion, the isolated recalculation of the electronic structure around the metal center when embedded into a DFT treatment of the ligand field shows promise as a pragmatic and lower cost treatment compared to the canonical treatment of the whole system of the difficult class of spin‐crossover complexes.

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

Scaling up the transcorrelated density matrix renormalization group

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. Here, we examine systems of size up to 12×12 lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by 2.4×–14×, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

Density matrix renormalization group↗
DOE OSTI2024

Pressure-Induced Exciton Formation and Superconductivity in Platinum-Based Mineral Sperrylite

We report a comprehensive study of Sperrylite (PtAs 2 ), the main platinum source in natural minerals, as a function of applied pressures up to 150 GPa. While no structural phase transition is detected from pressure-dependent X-ray measurements, the unit cell volume shrinks monotonically with pressure following the third-order Birch–Murnaghan equation of state. The mildly semiconducting behavior found in pure synthesized crystals at ambient pressures becomes more insulating upon increasing the applied pressure before metalizing at higher pressures, giving way to the appearance of an abrupt decrease in resistance near 3 K at pressures above 92 GPa consistent with the onset of a superconducing phase. The pressure evolution of the calculated electronic band structure reveals the same physical trend as our transport measurements, with a non-monotonic evolution explained by a hole band that is pushed below the Fermi energy and an electron band that approaches it as a function of pressure, both reaching a touching point suggestive of an excitonic state. A Lifshitz transition of the electronic structure and an increase in the density of states may naturally explain the onset of superconductivity in this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Twisted Tin‐Chloride Perovskite Single‐Crystal Heterostructures

Self-assembly affords simpler synthetic routes to heterostructures compared with manual layer-by-layer stacking, yet controlling interlayer twist angles in a bulk solid remains an outstanding challenge. We report two new single-crystal heterostructures: (Sn 2 Cl 2 )(CYS) 2 SnCl 4 (CYS = + NH 3 (CH 2 ) 2 S – ; Sn_CYS) and (Sn 2 Cl 2 )(SeCYS) 2 SnCl 4 (SeCYS = + NH 3 (CH 2 ) 2 Se – ; Sn_SeCYS) synthesized in solution, with alternating perovskite and intergrowth layers. Notably, compared to the recently reported lead analog, (Pb 2 Cl 2 )(CYS) 2 PbCl 4 (Pb_CYS), the tin heterostructures feature a twist between the perovskite and intergrowth layers. We trace this twist to local distortions at the Sn centers, which change the interfacial lattice-matching requirements compared to those of the Pb analog. Electronic band structure calculations show that the striking differences in the relative energies of perovskite- and intergrowth-derived bands in Sn_CYS and Pb_CYS arise from structural and not compositional differences. The structural anisotropy of Sn_CYS is also reflected in a large in-plane photoluminescence linear anisotropy ratio. Interfacial strain further affords differential incorporation of Pb into the perovskite and intergrowth layers of the Sn heterostructures, resulting in redshifted optical absorption onsets. Thus, we posit that local structural distortions may be exploited to manipulate the twist angle and interfacial strain in bulk heterostructures, providing a new handle for tuning the band alignments of bulk quantum-well electronic structures.

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

Singly occupied 4⁢𝑓 antiferromagnetic insulators: CePO 4 and CeVO 4

Rare-earth containing wide band gap oxides, which provide spin-photon interface and narrow linewidth optical emission, are getting significant attention as the most promising candidate materials in advancing quantum transduction and memories. Here, in this work, from ab initio calculations, we identify antiferromagnetic ground states in structurally preferred monoclinic CePO 4 and tetragonal CeVO 4 exhibiting localized occupied and unoccupied Ce 4⁢𝑓 states with 4⁢𝑓−4⁢𝑓 transition characteristics. Interestingly, in CePO 4 , O 2⁢𝑝 and P 3⁢𝑝 states hybridize negligibly with Ce 4⁢𝑓 states, while in CeVO 4 , V 3⁢𝑑 and O 2⁢𝑝 states hybridize and appear as extended states in between the occupied and unoccupied Ce 4⁢𝑓 states. Here, phonon calculations and analysis identify and differentiate Raman active phonon modes along with the spin phonon coupling of Ce in both CePO 4 and CeVO 4 that ultimately lead to different 4⁢𝑓 ground-state crystal-field multiplets, which are critical to accurately describe electronic transitions for foundational quantum transduction and memories. Further, the identified 𝐶 1 site symmetry of Ce, lacking inversion symmetry in CePO 4 , is relevant for quantum memories and 𝐷 2⁢𝑑 site symmetry of Ce exhibiting inversion symmetry in CeVO 4 is relevant for quantum transduction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Electronic structure theory with molecular point group symmetries on quantum annealers

Quantum computation has the potential to revolutionize quantum chemistry through major speedups in computation times and an exponential reduction in computational resources. Here, we combine the symmetry-adapted Jordan–Wigner encoding based on the full Boolean symmetry group $\mathbb{Z}$$^{k}_{2}$ with our new implementation of the Xia–Bian–Kais (XBK) method for improving the efficiency of electronic structure theory calculations on quantum annealers, particularly by reducing the number of qubits needed to achieve the same accuracy. By providing a more extensive symmetry-adapted encoding (SAE) than previous work, we are able to simulate molecules larger than those previously reported that have been studied using methods developed for quantum annealers and without using an active space. We calculated the potential energy surfaces of H 2 , LiH, He 2 , H 2 O, O 2 , N 2 , Li 2 , F 2 , CO, BH 3 , NH 3 , and CH 4 , with the largest molecule in the STO-6G basis set requiring 16 qubits with our SAE, and compared them with full configuration interaction results. The application of SAE to the XBK method provides an exponential reduction in the size of the Hilbert space and scales well with the size of the problem. It does not introduce significant additional errors for even or large values of a key variational parameter that determines the number of ancilla qubits used in the XBK method’s Hamiltonian embedding, or for certain molecules such as He 2 and H 2 O. Here, we provide an explanation for this behavior and a recommendation on the usage of our method. In addition, we briefly discuss the potential of extracting electronic excited states from our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗
DOE OSTI2025

Green’s Function Perspective on the Nonlinear Density Response of Quantum Many-Body Systems

We derive equations of motion for higher order density response functions using the theory of thermodynamic Green’s functions. We also derive expressions for the higher order generalized dielectric functions and polarization functions. Moreover, we relate higher order response functions and higher order collision integrals within the Martin–Schwinger hierarchy. We expect our results to be highly relevant to the study of a variety of quantum many-body systems such as matter under extreme temperatures, densities, and pressures.

Density Functional Theory↗
DOE OSTI2026

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

MBX V1.2: Accelerating Data-Driven Many-Body Molecular Dynamics Simulations

The MBX software provides an advanced platform for molecular dynamics simulations, leveraging state-of-the-art MB-pol and MB-nrg data-driven many-body potential energy functions. Developed over the past decade, these potential energy functions integrate physics-based and machine-learned many-body terms trained on electronic structure data calculated at the "gold standard" coupled-cluster level of theory. Recent advancements in MBX have focused on optimizing its performance, resulting in the release of MBX v1.2. While the inherently many-body nature of MB-pol and MB-nrg ensures high accuracy, it poses computational challenges. MBX v1.2 addresses these challenges with significant performance improvements, including enhanced parallelism that fully harnesses the power of modern multicore CPUs. In conclusion, these advancements enable simulations on nanosecond time scales for condensed-phase systems, significantly expanding the scope of high-accuracy, predictive simulations of complex molecular systems powered by data-driven many-body potential energy functions.

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

Discovering Ni/Cu Single-Atom Alloy as a Highly Active and Selective Catalyst for Direct Methane Conversion to Ethylene: A First-Principles Kinetic Study

Direct methane conversion to liquid fuels or value-added chemicals is a promising technology to utilize natural resources without resorting to further petroleum extraction. However, discovering efficient catalysts for this reaction is challenging due to either coke formation or unfavorable C-H bond activation. Herein, we design single-atom alloy (SAA) catalysts to simultaneously eliminate the above two bottlenecks based on mechanism-guided strategies: (1) the active single atom enables favorable C-H bond breaking and (2) the less reactive host metal facilitates C-C coupling and thus avoids strong binding of carbonaceous species. Employing electronic structure theory calculations, we screened the stability of multiple SAAs with 3d-5d transition metals atomically dispersed on a copper surface in terms of avoiding dopant aggregation and segregation. We then evaluated reactivities of the stable SAAs as catalysts for direct methane conversion to C2 products, including methane dehydrogenation and C-C coupling mechanisms. Combining selectivity analysis with kinetic modeling, we predicted that nickel dispersed on copper, i.e., Ni/Cu SAA, is a highly active and selective catalyst that can efficiently transform methane to ethylene. This work designs efficient SAA catalysts for direct methane activation and provides chemical insights into engineering compositions of SAAs to tune their catalytic performances.

Kothakonda, Manish↗
DOE OSTI2025

Unconventional low temperature decomposition of a saturated hydrocarbon over atomically-dispersed titanium-aluminum-boron catalyst

Sonochemically-synthesized atomically-dispersed titanium-aluminum-boron nanopowder (TiAlB NP) exhibits a remarkable low-temperature catalytic activation of aliphatic C-H bonds at 750 K followed by C-C bond activation thus emerging as a potent low-cost alternative to expensive platinum group metals. Here, the model saturated hydrocarbon, exo-tetrahydrodicyclopentadiene (C10H16), undergoes catalytic decomposition on TiAlB NPs in a chemical microreactor to produce 1,3-cyclopentadiene (c-C5H6), cyclopentene (c-C5H8), and molecular hydrogen (H2) as detected in situ via isomer-selective, single-photon ionization time-of-flight mass spectrometry. Extensive electronic structure theory calculations on model clusters of the catalyst decode a unique synergy among the atomic constituents of the catalyst and chemical bonding in this stepwise, retro Diels Alder reaction: Ti, although insensitive to C-H activation in its metallic state, initiates the catalysis via chemisorption of the hydrocarbon, adjacent B centers readily abstract hydrogen atoms and store them during the catalytic cycle, while Al stabilizes the catalyst structure yet providing space for critical docking sites for the departing hydrocarbons.

Biswas, Souvick↗
DOE OSTI2025

Ground-state-based model reduction with unitary circuits

Here, we present a method to numerically obtain low-energy effective models based on a unitary transformation of the ground state. The algorithm finds a unitary circuit that transforms the ground state of the original model to a projected wavefunction with only the low-energy degrees of freedom. The effective model can then be derived using the unitary transformation encoded in the circuit. We test our method on the one-dimensional and two-dimensional square-lattice Hubbard model at half-filling, and obtain more accurate effective spin models than the standard perturbative approach.

Hubbard model↗