Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wave Functions”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Resonant inelastic X-ray scattering determination of the electronic structure of oxyhemoglobin and its model complex

Hemoglobin and myoglobin are oxygen-binding proteins with S = 0 heme {FeO 2 } 8 active sites. The electronic structure of these sites has been the subject of much debate. This study uses Fe K-edge X-ray absorption spectroscopy (XAS) and 1s2p resonant inelastic X-ray scattering (RIXS) to study oxyhemoglobin and a related heme {FeO 2 } 8 model compound, [(pfp)Fe(1-MeIm)(O 2 )] (pfp = meso-tetra(α,α,α,α- o -pivalamido-phenyl)porphyrin, or TpivPP, 1-MeIm = 1-methylimidazole) (pfpO 2 ), which was previously analyzed using L-edge XAS. The K-edge XAS and RIXS data of pfpO 2 and oxyhemoglobin are compared with the data for low-spin Fe II and Fe III [Fe(tpp)(Im) 2 ] 0/+ (tpp = tetra-phenyl porphyrin) compounds, which serve as heme references. The X-ray data demonstrate that pfpO 2 is similar to Fe II , while oxyhemoglobin is qualitatively similar to Fe III , but with significant quantitative differences. Density-functional theory (DFT) calculations show that the cjontrast between pfpO 2 and oxyhemoglobin is due to a distal histidine H bond to O 2 and the less hydrophobic environment in the protein, which lead to more backbonding into the O 2 . A valence bond configuration interaction multiplet model is used to analyze the RIXS data and show that pfpO 2 is dominantly Fe II with 6–8% Fe III character, while oxyhemoglobin has a very mixed wave function that has 50–77% Fe III character and a partially polarized Fe–O 2 π-bond.

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

Spatiotemporal analysis of a final-state shape resonance in interferometric photoemission from Cu(111) surfaces

Photoemission from solid targets includes the excitation and motion of electrons inside the substrate, followed by their propagation in vacuum and detection. It thus depends on the electronic band structure of the solid in the two distinct spectral domains of bound initial and continuum final states. While the imprint of the static (initial-state) valence electronic structure of solids on photoemission spectra is routinely examined in standard photoemission spectroscopy in the energy domain, state-of-the-art time-resolved photoelectron spectroscopy allows, in addition, the scrutiny of photoelectron propagation in the electronic continuum. Within a quantum-mechanical model for attosecond time-resolved interferometric photoelectron emission from solids, we calculated photoemission spectra as a function of the delay between the exciting primary attosecond pulse train and assisting infrared (IR) laser pulse. Accounting for final-state interactions of the photoelectron with the IR laser electric field and the periodic substrate, our numerical results for interferometric photoemission from the 3d-valence band of Cu(111) surfaces show a striking resonantly enhanced sideband yield at photoelectron kinetic energies near 24 eV, in conjunction with a pronounced increase of the photoelectron wave-function amplitude inside the solid on a length scale of a few nanometers. In this work, this resonant shift of final-state photoelectron-probability density towards the bulk can be interpreted as an increase in the photoelectron propagation time in the solid and is commensurate with the resonantly enhanced spectral sideband-phase shifts observed in recent two-pathway two-photon interference spectra by Kasmi et al.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2006

Quantum Entanglement Molecular Absorption Spectrum Simulator

Quantum Entanglement Molecular Absorption Spectrum Simulator (QE-MASS) is a computer program for simulating two photon molecular-absorption spectroscopy using quantum-entangled photons. More specifically, QE-MASS simulates the molecular absorption of two quantum-entangled photons generated by the spontaneous parametric down-conversion (SPDC) of a fixed-frequency photon from a laser. The two-photon absorption process is modeled via a combination of rovibrational and electronic single-photon transitions, using a wave-function formalism. A two-photon absorption cross section as a function of the entanglement delay time between the two photons is computed, then subjected to a fast Fourier transform to produce an energy spectrum. The program then detects peaks in the Fourier spectrum and displays the energy levels of very short-lived intermediate quantum states (or virtual states) of the molecule. Such virtual states were only previously accessible using ultra-fast (femtosecond) laser systems. However, with the use of a single-frequency continuous wave laser to produce SPDC photons, and QEMASS program, these short-lived molecular states can now be studied using much simpler laser systems. QE-MASS can also show the dependence of the Fourier spectrum on the tuning range of the entanglement time of any externally introduced optical-path delay time. QE-MASS can be extended to any molecule for which an appropriate spectroscopic database is available. It is a means of performing an a priori parametric analysis of entangled photon spectroscopy for development and implementation of emerging quantum-spectroscopic sensing techniques. QE-MASS is currently implemented using the Mathcad software package.

Nguyen, Quang-Viet↗
DOE OSTI2025

Hydrogen and deuterium tunneling in niobium

We use density functional methods to identify the atomic configurations of H and D atoms trapped by O impurities embedded in bulk Nb. The O atoms are located at the octahedral position in the Nb body-centered cubic (BCC) lattice, and H (D) atoms tunnel between two degenerate tetrahedral sites separated by a mirror plane. Using nudged elastic band (NEB) methods, we calculate the double-well potential for O-H and O-D impurities and the wave functions and tunnel splittings for H and D atoms. Our results agree with those obtained from analysis of heat capacity and neutron scattering measurements on Nb with low concentrations of O-H and O-D.

density functional theory↗
DOE OSTI2024

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Green’s function coupled cluster formulations utilizing extended inner excitations

Here, we analyze new approximations of the Green’s function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i( n, m ) method, where the excitation level of the inner auxiliary operators ( m ) used to describe the ionization potential and electron affinity effects in the N-1 and N+1 particle spaces is higher than the excitation level ( n ) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called “n + 1” rule in this category [or the GFCC-i( n, n + 1) method], which states that in order to maintain size-extensivity of the Green’s function matrix elements, the excitation level of inner auxiliary operators X p (ω) and Y q (ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N 2 and CO molecules in the inner and outer valence regimes. Lastly, in comparison with the Green’s function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.

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

Superfast quarks in deuterium

An extension to our previous study on nuclear parton distribution functions (nPDFs) [Kim and Miller, Phys. Rev. C 106, 055202 (2022)] using light-front holographic quantum chromodynamics (LFHQCD) [Brodsky, de Teramond, Dosch, and Erlich, Phys. Rep. 584, 1 (2015)] is presented. We apply the effects of nucleon motion inside the nucleus (Fermi motion/smearing) to deuterium, extending our deuterium nPDFs to the superfast , 𝑥 > 1, region [Frankfurt and Strikman, Phys. Rep. 160, 235 (1988)] where we estimate our results to be reasonable up to 𝑥 ≈ 1.7. We utilize four different deuteron wave functions (AV18, NijmI, NijmII, Nijm93). We find that our model, with no additional new parameters, shows very good agreement with deuterium EMC ratio data obtained from the BONuS experiment [Fenker et al., Nucl. Instrum. Methods Phys. Res. A 592, 273 (2008); Baillie et al. (CLAS Collaboration), Phys. Rev. Lett. 108, 142001 (2012); Phys. Rev. Lett. 108, 199902(E) (2012); Tkachenko et al. (CLAS Collaboration), Phys. Rev. C 89, 045206 (2014); Phys. Rev. C 90, 059901(E) (2014); Griffioen et al., Phys. Rev. C 92, 015211 (2015)]. Looking beyond conventional nuclear physics, and in anticipation of 12 GeV experiments at Jefferson Lab, we use a LFHQCD ansatz to predict the contributions of an exotic six-quark state to the deuteron 𝐹 2 structure function, 𝐹$^D_2$, in the superfast region. We find that the effects of using other potentials are about the same magnitude as six-quark effects—both have small effects in 𝑥 < 1, but have significant contributions at 𝑥 > 1.

Physics↗
NASA NTRS2005

Shock Waves in a Bose-Einstein Condensate

A paper presents a theoretical study of shock waves in a trapped Bose-Einstein condensate (BEC). The mathematical model of the BEC in this study is a nonlinear Schroedinger equation (NLSE) in which (1) the role of the wave function of a single particle in the traditional Schroedinger equation is played by a space- and time-dependent complex order parameter (x,t) proportional to the square root of the density of atoms and (2) the atoms engage in a repulsive interaction characterized by a potential proportional to | (x,t)|2. Equations that describe macroscopic perturbations of the BEC at zero temperature are derived from the NLSE and simplifying assumptions are made, leading to equations for the propagation of sound waves and the transformation of sound waves into shock waves. Equations for the speeds of shock waves and the relationships between jumps of velocity and density across shock fronts are derived. Similarities and differences between this theory and the classical theory of sound waves and shocks in ordinary gases are noted. The present theory is illustrated by solving the equations for the example of a shock wave propagating in a cigar-shaped BEC.

Kulikov, Igor↗
DOE OSTI2014

Validity of virial theorem in all-electron mixed basis density functional, Hartree–Fock, and GW calculations

In this paper, we calculate kinetic and potential energy contributions to the electronic ground-state total energy of several isolated atoms (He, Be, Ne, Mg, Ar, and Ca) by using the local density approximation (LDA) in density functional theory, the Hartree–Fock approximation (HFA), and the self-consistent GW approximation (GWA). To this end, we have implemented self-consistent HFA and GWA routines in our all-electron mixed basis code, TOMBO. We confirm that virial theorem is fairly well satisfied in all of these approximations, although the resulting eigenvalue of the highest occupied molecular orbital level, i.e., the negative of the ionization potential, is in excellent agreement only in the case of the GWA. We find that the wave function of the lowest unoccupied molecular orbital level of noble gas atoms is a resonating virtual bound state, and that of the GWA spreads wider than that of the LDA and thinner than that of the HFA.

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

Ballooning theory for micro-tearing mode in tokamak

This paper aims to investigate the impact of magnetic drift on the linear micro-tearing mode by using a kinetic approach to derive a reduced two-field eigen system in real space. Here, since the magnetic drift in real space has derivatives, it is more efficient to solve the mode equations in a Fourier-ballooning representation using the two-dimensional (2D) ballooning transform. The lowest-order eigen system in the Fourier-ballooning representation consists of two integral equations, which are numerically solved using the finite difference method for both eigenvalues and wave functions. The main results will be presented through graphical eigenvalue scans for each parameter. Furthermore, we present a graphical comparison between the predictions of the ballooning theory and GENE gyrokinetic code simulation in the pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Light scattering of large rough particles application to cometary grains

While the electromagnetic field scattered by a spherical particle is classically obtained by the Helmholtz equation, the general case of an arbitrary particle may be investigated in the general framework of the interaction of a wave with a scattering potential. A wave function then satisfies the Schroedinger equation. The general solution of the Schroedinger equation is given. The main disadvantage of this approach are its restriction to large particles and its scalar nature preventing the calculation of the polarization. However, Perrin and Lamy have shown how to avoid the second limitation and retrieve a vectorial description. They proved that in the case of large spheres when the ad hoc assumptions are satisfied, the expression of the scattering amplitude may be approximated by an expansion series in partial waves, i.e., on a discrete basis. The analogy may be generalized, and the ratio of the two components for a rough particle obtained by taking the ratio of the reflectivities for the two directions of polarization. These reflectivities involve the simple and double reflections calculated following the method developed by Wolff for rough surfaces. The theory is further detailed.

Lamy, Philippe L.↗
DOE OSTI2012

Communication: Density functional theory overcomes the failure of predicting intermolecular interaction energies

Density-functional theory (DFT) revolutionized the ability of computational quantum mechanics to describe properties of matter and is by far the most often used method. However, all the standard variants of DFT fail to predict intermolecular interaction energies. In recent years, a number of ways to go around this problem has been proposed. We show that some of these approaches can reproduce interaction energies with median errors of only about 5% in the complete range of intermolecular configurations. Such errors are comparable to typical uncertainties of wave-function-based methods in practical applications. Thus, these DFT methods are expected to find broad applications in modelling of condensed phases and of biomolecules.

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

Variational quantum solver employing the PDS energy functional

In our previous work (J. Chem. Phys. 2020, 153, 201102),we reported a new class of quantum algorithms that are based on the quantum computation of the connected moment expansion to find the ground and excited state energies. In particular, the Peeters-Devreese-Soldatov (PDS) formulation is found variational and bearing the potential for further combining with the existing variational quantum infrastructure. Following this direction, here we propose a variational quantum solver employing the PDS energy gradient. In comparison with the usual variational quantum eigensolver (VQE) and the original static PDS approach, the proposed variational quantum solver offers an effective approach to achieve high ac-curacy at finding the ground state and its energy through the rotation of the trial wave function of modest quality guided by the low order PDS energy gradients, thus improves the ac-curacy and efficiency of the quantum simulation. We demonstrate the performance of the proposed variational quantum solver for toy models, H2molecule, and strongly correlated planar H4system in some challenging situations. In all the case studies, the proposed variational quantum approach out-performs the usual VQE and static PDS calculations even at the lowest order.

Peng, Bo↗
DOE OSTI2015

Role of information theoretic uncertainty relations in quantum theory

Uncertainty relations based on information theory for both discrete and continuous distribution functions are briefly reviewed. We extend these results to account for (differential) Rényi entropy and its related entropy power. This allows us to find a new class of information-theoretic uncertainty relations (ITURs). The potency of such uncertainty relations in quantum mechanics is illustrated with a simple two-energy-level model where they outperform both the usual Robertson–Schrödinger uncertainty relation and Shannon entropy based uncertainty relation. In the continuous case the ensuing entropy power uncertainty relations are discussed in the context of heavy tailed wave functions and Schrödinger cat states. Again, improvement over both the Robertson–Schrödinger uncertainty principle and Shannon ITUR is demonstrated in these cases. Further salient issues such as the proof of a generalized entropy power inequality and a geometric picture of information-theoretic uncertainty relations are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Electron impact ionization in dense plasmas

Here, the distorted-wave with exchange (DWE) method is employed to evaluate electron impact ionization cross sections in dense electron–ion plasmas. Bound and continuum electron wave functions are obtained from a non-relativistic average-atom code. Plots of DWE cross sections are presented for 1s electrons in Li and Be plasmas and 2p electrons in Na and Mg plasmas. For each of these elements, cross sections are evaluated at metallic density in a range of temperatures from 10 to 100 eV and at their respective melting points. Resonances in the cross sections appear near the incident energy threshold at high temperatures. The origin of these resonances is discussed. In general, the distorted wave (DW) method without exchange is found to be a good approximation to the DWE method for electron impact ionization calculations. In the resonance region, however, exchange effects are found to be very important and cannot be neglected.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Resonant Inelastic X-ray Scattering Calculations of Transition Metal Complexes Within a Simplified Time-Dependent Density Functional Theory Framework

We present a time-dependent density functional theory (TDDFT) approach to compute the light-matter couplings between two different manifolds of excited states relative to a common ground state in the context of 4d transition metal systems. These quantities are the necessary ingredients to solve the Kramers-Heisenberg (KH) equation for resonant inelastic X-ray scattering (RIXS) and several other types of two-photon spectroscopies. The procedure is based on the pseudo-wavefunction approach, where the solutions of a TDDFT calculation can be used to construct excited-state wave- functions, and on the restricted energy window approach, where a manifold of excited states can be rigorously defined based on the energies of the occupied molecular orbitals involved in the excitation process. Thus, the present approach bypasses the need to solve the costly TDDFT quadratic-response equations. We illustrate the applicability of the method to 4d transition metal molecular complexes by calculating the 2p4d RIXS maps of three representative ruthenium complexes and comparing them to experimental results. The method can capture all the experimental features in all three complexes to allow the assignment of the experimental peaks, and relative energies correct to within 0.6 eV at the cost of two independent TDDFT calculations.

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

Multistep Quantum Master Equation Theory for Response Functions in Four Wave Mixing Electronic Spectroscopy of Multichromophoric Macromolecules

This work provides an alternative derivation of third order response functions in four wave mixing spectroscopy of multichromophoric macromolecular systems considering only single exciton states. For the case of harmonic oscillator bath linearly and diagonally coupled to exciton states, closed form expressions showing all the explicit time dependences are derived. Here, these expressions can provide more solid physical basis for understanding 2-dimensional electronic spectroscopy signals. For more general cases of system-bath coupling, the quantum master equation (QME) approach is employed for the derivation of multistep time evolution equations for Green function-like operators. Solution of these equations is feasible at the level of 2nd order non-Markovian QME, and the new approach can account for inter-exciton coupling, dephasing, relaxation, and non-Markovian effects in a consistent manner.

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

Photoelectron spectra of early 3 d -transition metal dioxide molecular anions from GW calculations

Photoelectron spectra of early 3 d -transition metal dioxide anions, Sc O 2 - , Ti O 2 - , V O 2 - , Cr O 2 - , and Mn O 2 - , are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point ( G 0 W 0 @PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV–0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew–Burke–Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G 0 W 0 @PBE0, the G -only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3 d -transition metal dioxide anions. Adding eigenvalue self-consistency also in W , on the other hand, worsens the agreement with experiment overall. Overall, our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗