Search NASA⌕ Search

SEARCH · Search NASA

Results for “wavefunction analysis”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

libwfa: Wavefunction analysis tools for excited and open‐shell electronic states

Abstract An open‐source software library for wavefunction analysis, libwfa, provides a comprehensive and flexible toolbox for post‐processing excited‐state calculations, featuring a hierarchy of interconnected visual and quantitative analysis methods. These tools afford compact graphical representations of various excited‐state processes, provide detailed insight into electronic structure, and are suitable for automated processing of large data sets. The analysis is based on reduced quantities, such as state and transition density matrices (DMs), and allows one to distill simple molecular orbital pictures of physical phenomena from intricate correlated wavefunctions. The implemented descriptors provide a rigorous link between many‐body wavefunctions and intuitive physical and chemical models, for example, exciton binding, double excitations, orbital relaxation, and polyradical character. A broad range of quantum‐chemical methods is interfaced with libwfa via a uniform interface layer in the form of DMs. This contribution reviews the structure of libwfa and highlights its capabilities by several representative use cases. This article is categorized under: Software > Quantum Chemistry Theoretical and Physical Chemistry > Spectroscopy

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

Unveiling in-gap states and Majorana zero modes in superconductor-topological insulator bilayer model

Interfaces between topological insulators and superconductors (SCs) are promising platforms for realizing Majorana zero modes (MZMs) via the superconducting proximity effect. We introduce a bilayer model consisting of the surface states of a three-dimensional topological insulator (3DTI) coupled to an 𝑠-wave superconductor and systematically study the role of interlayer tunneling strength (𝑡⊥) motivated by the recent growth of the Fe⁡(Te,Se)/Bi 2 ⁢Te 3 heterostructure. We find that increasing 𝑡⊥ shifts the proximity-induced (PrI) gap minimum away from the Γ point, giving rise to momentum-selective interference patterns that manifest as spatial oscillations in the in-gap states. We introduce an antidot with a magnetic vortex in the model and investigate the resulting in-gap states, including MZMs and Caroli–de Gennes–Matricon (CdGM) modes. With increasing hybridization strength, the energy separation between MZMs and CdGM states increases, enhancing the isolation of MZMs. Importantly, in the strong hybridization limit, the leading CdGM separation remains large in spite of the decrease in the PrI gap. Spin- and spatially resolved wavefunction analysis reveals angular momentum asymmetries absent in conventional 𝑠-wave systems. A direct comparison with a standalone 𝑠-wave superconductor confirms the emergence of distinct 𝑝-wave-like features in the bilayer geometry. Our results provide experimentally relevant predictions for tuning the stability of MZMs and their differentiation from the CdGM modes in SC-3DTI heterostructures and offer a theoretical framework for probing unconventional superconductivity in engineered topological systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

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

Resonant double-core excitations with ultrafast, intense X-ray pulses

Intense few-to-sub-femtosecond soft X-ray pulses can produce neutral, two-site excited double-core-hole states by promoting two core electrons to the same unoccupied molecular orbital. We theoretically investigate double nitrogen K-edge excitations of nitrous oxide (N$_2$O) with multiconfigurational electronic structure calculations. We show that the second core-excitation energy is reduced with respect to its ground state value. A site-selective double core-excitation mechanism using intense few-femtosecond x-rays is investigated using time-dependent Schrodinger equation (TDSE) simulations. The subsequent two-step Auger-Meitner and two-electrons-one-electron decay spectra of the double core-excited states are analysed using a Mulliken population analysis of the multiconfirational wavefunctions. The change in the electron emission lineshape between the absorption of 1 or 2 photons in the resonant core-excitation is predicted by combining this approach with the TDSE simulations. Furthermore, we examine the possibility of resolving the double core-excited states with X-ray pump-probe techniques by calculating the chemical shifts of the core-electron binding energy of the core-excited states and decay products.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Quantum-classical electron distributions in atoms and atomic ions

A quantum-classical approach is used to obtain the velocity distributions in atoms and positive and negative ions in both ground and excited states. In the analysis, Hartree-Fock electronic wavefunctions are used to determine the radial electron distributions, and the central-field approximation is used to study the the dynamic properties of the localized electrons. The distributions for the outer and inner shells are found to agree well with exact results obtained by numerical calculations.

Kunc, Joseph A.↗
DOE OSTI2023

Quantum textures of the many-body wavefunctions in magic-angle graphene

Interactions among electrons create novel many-body quantum phases of matter with wavefunctions that often reflect electronic correlation effects, broken symmetries, and novel collective excitations. A wide range of quantum phases has been discovered in MATBG, including correlated insulating, unconventional superconducting, and magnetic topological phases. The lack of microscopic information, including precise knowledge of possible broken symmetries, has thus far hampered our understanding of MATBG's correlated phases. Here we use high-resolution scanning tunneling microscopy to directly probe the wavefunctions of the correlated phases in MATBG. The squares of the wavefunctions of gapped phases, including those of the correlated insulators, pseudogap, and superconducting phases, show distinct patterns of broken symmetry with a √3 x √3 super-periodicity on the graphene atomic lattice that has a complex spatial dependence on the moiré superlattice scale. We introduce a symmetry-based analysis to describe our measurements of the wavefunctions of MATBG's correlated phases with a set of complex-valued local order parameters. For the correlated insulators in MATBG, at fillings of ±2 electrons per moiré unit cell relative to charge neutrality, we compare the observed quantum textures to those expected for proposed theoretical ground states. In typical MATBG devices, the textures of correlated insulators' wavefunctions closely match those of the theoretically proposed IKS order, while in ultra-low-strain samples our data has local symmetries like those of a T-IVC phase. Here, we also study the wavefunction of MATBG's superconducting state, revealing strong signatures of intervalley coherence that can only be distinguished from those of the insulator with our phase-sensitive measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Illuminating Ligand Field Contributions to Molecular Qubit Spin Relaxation via T 1 Anisotropy

Electron spin relaxation in paramagnetic transition metal complexes constitutes a key limitation on the growth of molecular quantum information science. However, there exist very few experimental observables for probing spin relaxation mechanisms, leading to a proliferation of inconsistent theoretical models. Here we demonstrate that spin relaxation anisotropy in pulsed electron paramagnetic resonance is a powerful spectroscopic probe for molecular spin dynamics across a library of highly coherent Cu(II) and V(IV) complexes. Here, neither the static spin Hamiltonian anisotropy nor contemporary computational models of spin relaxation are able to account for the experimental T 1 anisotropy. Through analysis of the spin-orbit coupled wavefunctions, we derive an analytical theory for the T 1 anisotropy that accurately reproduces the average experimental anisotropy of 2.5. Furthermore, compound-by-compound deviations from the average anisotropy provide a promising approach for observing specific ligand field and vibronic excited state coupling effects on spin relaxation. Finally, we present a simple density functional theory workflow for computationally predicting T 1 anisotropy. Analysis of spin relaxation anisotropy leads to deeper fundamental understanding of spin-phonon coupling and relaxation mechanisms, promising to complement temperature-dependent relaxation rates as a key metric for understanding molecular spin qubits.

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

Chemical information from XPS: Theory and experiment for Ni(OH) 2

For this work, the features and the electronic character of the states for the Ni 2p x-ray photoelectron spectroscopy (XPS) of Ni(OH) 2 were analyzed. This detailed analysis is based on ab initio molecular orbital wavefunctions for a cluster model of Ni(OH) 2 . The theory is validated by comparison with experiment. Then, advanced methods are used to explain and contrast the properties of different groups of ionic states. An important conclusion is that in most cases, the ionic states cannot be described with a single configuration or determinant. Despite this essential many-body character of the XPS, we demonstrate that it is possible to understand the origin of the main and satellite XPS features in terms of their orbital character.

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

Exponential acceleration of macroscopic quantum tunneling in a Floquet Ising model

The exponential suppression of macroscopic quantum tunneling (MQT) in the number of elements to be reconfigured is an essential element of broken symmetry phases. This suppression is also a core bottleneck in quantum algorithms, such as traversing an energy landscape in optimization, and adiabatic state preparation more generally. In this work, we demonstrate exponential acceleration of MQT through Floquet engineering with the application of a uniform, high frequency transverse drive field. Using the ferromagnetic phase of the transverse field Ising model in one and two dimensions as a prototypical example, we identify three phenomenological regimes as a function of drive strength. For weak drives, the system exhibits exponentially decaying tunneling rates but robust magnetic order; in the crossover regime at intermediate drive strength, we find polynomial decay of tunnelling alongside vanishing magnetic order; and at very strong drive strengths both the Rabi frequency and time-averaged magnetic order are approximately constant with increasing system size. We support these claims with extensive full wavefunction and tensor network numerical simulations, and theoretical analysis. An experimental test of these results presents a technologically important and novel scientific question accessible on NISQ-era quantum computers.

Grattan, George↗
DOE OSTI2022

Simulating Effective QED on Quantum Computers

In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of T -gates needed to perform such simulations on a 3 D lattice of n s sites scales at worst as O ( n s 3 / ϵ ) 1 + o ( 1 ) in the thermodynamic limit for position basis simulations and O ( n s 4 + 2 / 3 / ϵ ) 1 + o ( 1 ) in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of O ~ ( n s 2 + 2 / 3 / ϵ ) for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of O ~ ( n s 5 + 2 / 3 / ϵ ) . We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than 10 13 non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Benchmarking third-order cluster perturbation theory for electronically excited states

In this study, we investigate the reliability of cluster perturbation (CP) theory applied to the calculation of electronically excited states through a comprehensive benchmark. In CP theory, perturbative corrections are added to the properties of a parent excitation space, which converge toward the properties of a target excitation space. For the CPS(D-n) model, perturbative corrections through order n are added to the coupled cluster singles (CCS) excitation energies to target the coupled cluster singles and doubles (CCSD) excitation energies. Through a comparative analysis of excitation energy calculations across a diverse set of molecules and wavefunction methods, we present a comprehensive evaluation of the accuracy of the third-order CPS(D) model, CPS(D-3), in calculating excitation energies. Further, our findings demonstrate that CPS(D-3) is a reliable alternative to established methods, particularly CCSD, while systematically overestimating the excitation energies compared to high-level coupled cluster methods such as CC3. These results highlight the strengths and limitations of CPS(D-3), as well as the promising directions for its future development.

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

Charting C–C coupling pathways in electrochemical CO 2 reduction on Cu(111) using embedded correlated wavefunction theory

The electrochemical CO 2 reduction reaction (CO 2 RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C 2+ ) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO 2 RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C 2+ product formation via CO 2 RR electrocatalysts.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Partial Wave Analysis of the omega pi^- Final State Photoproduced at GlueX

The meson has traditionally been understood as the bound state of a quark and an anti-quark, but there is experimental evidence for exotic meson states, which have properties that cannot be produced by this quark-antiquark model. One candidate for such exotic systems is hybrid meson, predicted by lattice QCD, which includes gluonicexcitation in its wavefunction. One of the prominent decay modes of the lightest predicted exotic hybrid is expected to be b1pi. THis dissertation presernts an analysis of the gamma rho -> omega pi- Delta ++ channel produced at the GlueX experiment, in the mass range where the b1(1235) meson is prominent. An amplitude analysis is performed on the decay of the omega pi- system, in order to extract the resonance parameters of the b-1, as well as analyze the resonant and non-resonant background.

Schertz, Amy↗
NASA NTRS1992

Negative-energy states in the Dirac-Hartree-Fock problem - The effect of omission of two-electron integrals involving the small component

The effect of omission of two-electron integrals involving basis functions for the small component of the wavefunction on the eigenvalue spectrum in the Dirac-Hartree-Fock problem is studied. From an analysis of the Fock matrix it is shown that omission of these integrals moves the negative-energy states down, not up. Their complete omission does not give rise to intruder states. The appearance of intruder states occurs when only some of the core integrals are omitted, due to the nature of particular contraction schemes used for the core basis functions. Use of radially localized functions rather than atomic functions alleviates the intruder state problem.

Dyall, Kenneth G.↗
DOE OSTI2023

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Despite decades of practice, finite-size errors in many widely used electronic structure theories for periodic systems remain poorly understood. For periodic systems using a general Monkhorst-Pack grid, there has been no comprehensive and rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order Møller-Plesset perturbation theory (MP2), which are the simplest wavefunction based method, and the simplest post-Hartree-Fock method, respectively. Such calculations can be viewed as a multi-dimensional integral discretized with certain trapezoidal rules. Due to the Coulomb singularity, the integrand has many points of discontinuity in general, and standard error analysis based on the Euler-Maclaurin formula gives overly pessimistic results. The lack of analytic understanding of finite-size errors also impedes the development of effective finite-size correction schemes. We propose a unified analysis to obtain sharp convergence rates of finite-size errors for the periodic HF and MP2 theories. Our main technical advancement is a generalization of the result of Lyness [Math. Comp. 30 (1976), pp. 1–23] for obtaining sharp convergence rates of the trapezoidal rule for a class of non-smooth integrands. Our result is applicable to three-dimensional bulk systems as well as low dimensional systems (such as nanowires and 2D materials). Our unified analysis also allows us to prove the effectiveness of the Madelung-constant correction to the Fock exchange energy, and the effectiveness of a recently proposed staggered mesh method for periodic MP2 calculations (see X. Xing, X. Li, and L. Lin [J. Chem. Theory Comput. 17 (2021), pp. 4733–4745]). In conclusion, our analysis connects the effectiveness of the staggered mesh method with integrands with removable singularities, and suggests a new staggered mesh method for reducing finite-size errors of periodic HF calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Manipulation of chiral interface states in a moiré quantum anomalous Hall insulator

Moiré systems made from stacked two-dimensional materials host correlated and topological states that can be electrically controlled with applied gate voltages. One prevalent form of topological state that can occur are Chern insulators that display a quantum anomalous Hall effect. Here we manipulate Chern domains in an interaction-driven quantum anomalous Hall insulator made from twisted monolayer–bilayer graphene and observe chiral interface states at the boundary between different domains. By tuning the carrier concentration, we stabilize neighbouring domains of opposite Chern number that then provide topological interfaces devoid of any structural boundaries. This allows the wavefunction of chiral interface states to be directly imaged using a scanning tunnelling microscope. In conclusion, our theoretical analysis confirms the chiral nature of observed interface states and allows us to determine the characteristic length scale of valley polarization reversal across neighbouring Chern domains.

36 MATERIALS SCIENCE↗