Search NASA⌕ Search

SEARCH · Search NASA

Results for “Natural orbitals”

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

Modeling Oxidative Dehydrogenation of Propane with Supported Vanadia Catalysts Using Multireference Methods

The oxidative dehydrogenation of propane over supported vanadium oxide catalysts poses significant computational challenges due to complex electronic structure changes along the reaction coordinate, driven primarily by changes in the oxidation states of vanadium. To address these challenges, we systematically test quantum chemical methods, including multireference (MR) approaches, domain-based local pair natural orbital coupled cluster theory (DLPNO-CCSD(T)), and density functional theory (DFT). The initial C–H bond-breaking transition state requires MR treatment due to its multireference character, while subsequent steps permit efficient single-reference calculations. For the rate-limiting C–H activation step mediated by the vanadyl moiety, complete 1 active space second-order perturbation theory (CASPT2) yields an apparent activation barrier (E app 600K ) of 138 kJ/mol, consistent with experimental values (134 ± 4 kJ/mol; Gruene et al. Catal. Today 2010, 157, 137). In contrast, DLPNO-CCSD(T) overestimates this barrier (198 kJ/mol), whereas DFT predictions span 125–150 kJ/mol, depending on the functional. Our multireference investigation of this transition metal oxide-catalyzed process demonstrates that an active space that incorporates the C–H σ and V=O σ/π bonding orbitals, oxygen lone pairs, and their antibonding counterparts adequately captures electronic structure changes along the chemical transformation. Furthermore, these findings provide a general strategy for active space selection in transition metal oxide-catalyzed C/O–H bond activation reactions. The reference dataset from this work, which includes MR calculations with manually selected active spaces for all intermediates and transition states in the propane ODH reaction network, will serve as a benchmark for automating active space selection in similar systems.

Catalysts↗
DOE OSTI2021

Resonant Stimulated X-ray Raman Spectroscopy of Mixed-Valence Manganese Complexes

Resonant stimulated X-ray Raman spectroscopy of the bimetallic [Mn III Mn IV (μ-O) 2 (μ-OAC)(tacn) 2 ] 2+ manganese complex is investigated in a simulation study. A pump–probe sequence of resonant X-ray Raman excitations is employed to generate and monitor valence-electron wave packets at selected regions in the molecule, by exploiting element specific core-excited states. A two-color protocol is presented, with pump and probe pulses tuned to the Mn and N K-edges. A natural orbital decomposition allows the visualization of the electron dynamics underlying the signal.

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

The role of electron correlations in the electronic structure of putative Chern magnet TbMn6Sn6

Abstract A member of the RMn 6 Sn 6 rare-earth family materials, TbMn 6 Sn 6 , recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard U (DFT + U ) calculations to examine the electronic structure of TbMn 6 Sn 6 . To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration–interaction (CI) methods. We find that DFT + U and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT + U with Mn magnetic moments adjusted to the experiment predict the Dirac crossing in bulk to be close to the Fermi level, within ~120 meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dipole ladders with large Hubbard interaction in a moiré exciton lattice

Two-dimensional semiconductor moiré superlattices have emerged as a powerful platform for engineering correlated electronic phenomena. On the other hand, optical excitation creates charge neutral interlayer excitons with an out-of-plane electric dipole. Strong onsite dipole–dipole interaction promises the formation of correlated bosonic states, akin to the Mott states of electrons, but has not yet been demonstrated. Here we report a large interaction between excitons occupying the same moiré lattice site—characterized by the Hubbard U parameter—and consequent dipole ladders with spin- and electron-filling dependence in WSe 2 /WS 2 moiré superlattices. Photoluminescence measurements show successive peaks emerging with an energy separation of around 34 meV above the ground state as the exciton density is increased. This corresponds to the sequential injection of excitons into a single site with an energy cost to overcome the large exciton Hubbard U, forming a dipole ladder. Based on findings of local magnetic moments at two holes per moiré cell, we show that excitons can also fill a second moiré orbital, establishing the two-orbital nature of the moiré potential landscape. Our results show that the Bose–Hubbard model with possible exciton crystal phases can be investigated in interacting opto-moiré quantum matter.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Strain-tuned incompatible magnetic exchange-interaction in La2NiO4

Abstract Magnetic frustration is a route for novel ground states, including spin liquids and spin ices. Such frustration can be introduced through either lattice geometry or incompatible exchange interactions. Here, we find that epitaxial strain is an effective tool for tuning antiferromagnetic exchange interactions in a square-lattice system. By studying the magnon excitations in La 2 NiO 4 films using resonant inelastic x-ray scattering, we show that the magnon displays substantial dispersion along the antiferromagnetic zone boundary, at energies that depend on the lattice of the film’s substrate. Using first principles simulations and an effective spin model, we demonstrate that the antiferromagnetic next-nearest neighbour coupling is a consequence of the two-orbital nature of La 2 NiO 4 . Altogether, we illustrate that compressive epitaxial strain enhances this coupling and, as a result, increases the level of incompatibility between exchange interactions within a model square-lattice system.

Physics↗
DOE OSTI2020

Developing new and understanding old approximations in TDDFT

When a system has evolved far from a ground-state, the adiabatic approximations commonly used in time-dependent density functional theory calculations completely fail in some applications, while giving qualitatively good predictions in others, and sometimes even quantitative predictions. It is not clearly understood why this is so, and developing practical approximations going beyond the adiabatic approximation remains a challenge. This paper explores three different lines of investigation. First, an expression for the exact time-dependent exchange–correlation potential suggests that the accuracy of an adiabatic approximation is intimately related to the deviation between the natural orbital occupation numbers of the physical system and those of the Kohn–Sham system, and we explore this on some exactly-solvable model systems. The exact expression further suggests a path to go beyond the adiabatic approximations, and in the second part we discuss a newly proposed class of memory-dependent approximations developed in this way. Finally, we derive a new expression for the exact exchange–correlation potential from a coupling-constant path integration.

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

Improving the accuracy of the variational quantum eigensolver for molecular systems by the explicitly-correlated perturbative [2] R12 -correction

We provide an integration of the universal, perturbative explicitly correlated [2] R12 -correction in the context of the Variational Quantum Eigensolver (VQE). This approach is able to increase the accuracy of the underlying reference method significantly while requiring no additional quantum resources. The proposed approach only requires knowledge of the one- and two-particle reduced density matrices (RDMs) of the reference wavefunction; these can be measured after having reached convergence in the VQE. This computation comes at a cost that scales as the sixth power of the number of electrons. Here, we explore the performance of the VQE + [2] R12 approach using both conventional Gaussian basis sets and our recently proposed directly determined pair-natural orbitals obtained by multiresolution analysis (MRA-PNOs). Both Gaussian orbital and PNOs are investigated as a potential set of complementary basis functions in the computation of [2] R12 . In particular the combination of MRA-PNOs with [2] R12 has turned out to be very promising – persistently throughout our data, this allowed very accurate simulations at a quantum cost of a minimal basis set. Additionally, we found that the deployment of PNOs as complementary basis can greatly reduce the number of complementary basis functions that enter the computation of the correction at a complexity.

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

Energetics and mechanisms for the acetonyl radical + O 2 reaction: An important system for atmospheric and combustion chemistry

The acetonyl radical (•CH 2 COCH 3 ) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O 2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O 2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. Additionally, no submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol -1 higher in energy than acetonyl radical + O 2 . The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O 2 . Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures.

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

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

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

Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics

In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.

Chemistry↗
DOE OSTI2025

Taming the virtual space for incremental full configuration interaction

Incremental full configuration interaction (iFCI) closely approximates the FCI limit with polynomial cost through a many-body expansion of the correlation energy, providing highly accurate total energies within a given basis set. To extend iFCI beyond previous basis set limitations, this work introduces a novel natural orbital (NO) screening approach, incremental NO full configuration interaction (iNO-FCI). By consideration of the importance of virtual orbital selection in the convergence of iFCI, iNO-FCI maximizes the consistency between orbitals selected for each correlated body. iNO-FCI employs a principle of cancellation of errors and ensures that the same set of virtual NOs is used for interdependent terms. Here, this strategy significantly reduces computational cost without compromising precision. Computational savings of up to 95% are demonstrated, allowing access to larger basis sets that were previously computationally prohibitive. iNO-FCI is herein introduced and benchmarked for several difficult test cases involving double-bond dissociation, biradical systems, conjugated π systems, and the spin gap of a Cu-based transition metal complex.

Correlation energy↗
DOE OSTI2020

Electronic structure and photophysics of a supermolecular iron complex having a long MLCT-state lifetime and panchromatic absorption

Exploiting earth-abundant iron-based metal complexes as high-performance photosensitizers demands long-lived electronically excited metal-to-ligand charge-transfer (MLCT) states, but these species suffer typically from femtosecond timescale charge-transfer (CT)-state quenching by low-lying nonreactive metal-centered (MC) states. Here, we engineer supermolecular Fe(II) chromophores based on the bis(tridentate-ligand)metal(II)-ethyne-(porphinato)zinc(II) conjugated framework, previously shown to give rise to highly delocalized low-lying 3 MLCT states for other Group VIII metal (Ru, Os) complexes. Electronic spectral, potentiometric, and ultrafast pump–probe transient dynamical data demonstrate that a combination of a strong σ-donating tridentate ligand and a (porphinato)zinc(II) moiety with low-lying π*-energy levels, sufficiently destabilize MC states and stabilize supermolecular MLCT states to realize Fe(II) complexes that express 3 MLCT state photophysics reminiscent of their heavy-metal analogs. The resulting Fe(II) chromophore archetype, FeNHCPZn, features a highly polarized CT state having a profoundly extended 3MLCT lifetime (160 ps), 3 MLCT phosphorescence, and ambient environment stability. Density functional and domain-based local pair natural orbital coupled cluster [DLPNO-CCSD(T)] theory reveal triplet-state wavefunction spatial distributions consistent with electronic spectroscopic and excited-state dynamical data, further underscoring the dramatic Fe metal-to-extended ligand CT character of electronically excited FeNHCPZn. This design further prompts intense panchromatic absorptivity via redistributing high-energy absorptive oscillator strength throughout the visible spectral domain, while maintaining a substantial excited-state oxidation potential for wide-ranging photochemistry––highlighted by the ability of FeNHCPZn to photoinject charges into a SnO 2 /FTO electrode in a dye-sensitized solar cell (DSSC) architecture. Concepts enumerated herein afford opportunities for replacing traditional rare-metal–based emitters for solar-energy conversion and photoluminescence applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Optimized low-depth quantum circuits for molecular electronic structure using a separable-pair approximation

We present a classically tractable model that leads to optimized low-depth quantum circuits leveraging separable-pair approximations. The obtained circuits are well suited as a baseline circuit for emerging quantum hardware and can, in the long term, provide significantly improved initial states for quantum algorithms. The associated wave functions can be represented with linear memory requirement, which allows classical optimization of the circuits and naturally defines a minimum benchmark for quantum algorithms. In this work we employ directly determined pair-natural orbitals within a basis-set-free approach. This leads to accurate representation of the one- and many-body parts for weakly correlated systems and we explicitly illustrate how the model can be integrated into other quantum algorithms for stronger correlated systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Quasiparticle self-consistent G W band structure of CrN

The band structures of CrN are calculated using the quasiparticle self-consistent GW approach in both the low-temperature antiferromagnetic (AFM) [110] 2 Pmna orthorhombic phase and in a hypothetical ferromagnetic phase representing the paramagnetic cubic state in a saturating magnetic field. A gap of about 1 eV is found in the AFM state, while the ferromagnetic band structure is found to be half metallic. Another hypothetical AFM-1 structure [001] 1 is further considered and gives a smaller band gap. The orbital nature of the bands is studied and reveals a gap between Cr-d states, indicating a strongly correlated behavior. Here, optical dielectric functions are calculated from the interband transitions and compared to experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Constructing Hubbard models for the hydrogen chain using sliced-basis density matrix renormalization group

Sliced-basis DMRG(sb-DMRG) is used to simulate a chain of hydrogen atoms and to construct low-energy effective Hubbard-like models. The downfolding procedure first involves a change of basis to a set of atom-centered Wannier functions constructed from the natural orbitals of the exact DMRG one-particle density matrix. The Wannier function model is then reduced to a fewer-parameter Hubbard-like model, whose parameters are determined by minimizing the expectation value of the Wannier Hamiltonian in the ground state of the Hubbard Hamiltonian. This indirect variational procedure not only yields compact and simple models for the hydrogen chain, but also allows us to explore the importance of constraints in the effective Hamiltonian, such as the restricting the range of the single-particle hopping and two-particle interactions, and to assess the reliability of more conventional downfolding. The entanglement entropy for a model's ground state, cut in the middle, is an important property determining the ability of DMRG and tensor networks to simulate the model, and here we study its variation with the range of the interactions. Counterintuitively, we find that shorter ranged interactions often have larger entanglement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anomalies in the electromagnetic response and peculiarities of the domain structure in the magnetic Weyl semimetal Co 3⁢ Sn 2 ⁢S 2

In this paper, we present a comprehensive study of magnetic properties of Co 3 ⁢Sn 2 ⁢S 2 (CoSnS) crystals using macroscopic characterization and observations of the domain structure (DS) in a wide range of easy- and hard-axis fields (𝐻) and temperatures (𝑇). From macroscopic measurements, we derive temperature dependences of the basic magnetic parameters of the material [saturation magnetization 𝑀 𝑠 ⁡(𝑇), uniaxial anisotropy constant 𝐾⁡(𝑇), and exchange stiffness 𝐴⁡(𝑇)] and, based on them, calculate the thermodynamic magnetization 𝑀⁡(𝐻,𝑇) curves and DS parameters for pinning-free samples. We find important qualitative and quantitative differences between the theory and experiment defined by the efficient pinning of domain walls (DWs) in our samples. It turns out that, already at temperatures close to the ferromagnetic transition 𝑇 𝑐 , the DW dynamics falls into the heavy creep regime, resulting in a strong deferral of the DS adjustment into the equilibrium state. Here, we suggest that specific magnetic defects, responsible for pinning and causing field-tuned reproducible domain nucleation patterns observed in our high-quality samples, are pointlike. Recently found magnetic polarons formed around vacancies and impurities, which carry strong diamagnetic moments of spin-orbit nature and are intrinsic for CoSnS and expected in other Weyl semimetals, are a possible source of pinning yielding peculiar 𝑀⁡(𝐻,𝑇) response of CoSnS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Measures of complexity and entanglement in many-fermion systems

There is no unique and widely accepted definition of the complexity measure (CM) of a many-fermion wave function in the presence of interactions. The simplest many-fermion wave function is a Slater determinant. In shell-model or configuration interaction (CI) and other related methods, the state is represented as a superposition of a large number of Slater determinants, which in the case of CI calculations reaches about 20 billion terms [Johnson, arXiv:1809.07869]. Although in practice this number has been used as a CM for decades, it is ill defined: it is not unique, and it depends on the particular type and the number of single-particle wave functions used to construct the Slater determinants. Further, the canonical wave functions and/or natural orbitals [Löwdin, Adv. Phys. 5, 1 (1956); Löwdin and Shull, Phys. Rev. 101, 1730 (1956); Bardeen et al., Phys. Rev. 108, 1175 (1957); N. N. Bogoljubov, Il Nuovo Cimento 7, 794 (1958); Valatin, Il Nuovo Cimento 7, 843 (1958); de Gennes, Superconductivity of Metals and Alloys (CRC Press, Boca Raton, FL, 1999); Ring and Schuck, The Nuclear Many-Body Problem, 1st ed. (Springer-Verlag, Berlin, 2004)] and their corresponding occupation probabilities are intrinsic properties of any many-body wave function, irrespective of the representation, and they provide a unique solution to characterize the CM. The non-negative orbital entanglement entropy, which vanishes for a Slater determinant, provides the simplest CM, while a more complete measure of complexity is the entanglement spectrum. We illustrate these aspects in the case of a complex nonequilibrium time-dependent process, induced nuclear fission described within a real-time density functional theory framework extended to superfluid systems, which can describe simultaneously the long-range and the short-range correlations between fermions. The orbital entanglement entropy of the fissioning nucleus illustrates the localization mechanism of the many-body wave function in Fock and/or Hilbert space. The (minimal) number of Slater determinants required to represent such a complex many-body wave function with a well-defined number of particles in the case presented here is about 10 500 . The realistic case of the highly nonequilibrium nuclear fission process illustrated here is equivalent to a system of 23.328×10 9 interacting quantum spin-1/2 particles, a very large system for the study of quantum entanglement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab Initio Structure Factors for Spin-Dependent Dark Matter Direct Detection

We present converged ab initio calculations of structure factors for elastic spin-dependent WIMP scattering off all nuclei used in dark matter direct-detection searches: 19 F, 23 Na, 27 Al, 29 Si, 73 Ge, 12 7I, and 129,131 Xe. From a set of established two- and three-nucleon interactions derived within chiral effective field theory, we construct consistent WIMP-nucleon currents at the one-body level, including effects from axial-vector two-body currents. We then apply the in-medium similarity renormalization group to construct effective valence-space Hamiltonians and consistently transformed operators of nuclear responses. Combining the recent advances of natural orbitals with three-nucleon forces expressed in large spaces, we obtain basis-space converged structure factors even in heavy nuclei. Generally results are consistent with previous calculations but large uncertainties in 127 I highlight the need for further study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗