Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave function”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electronic Structure of Excitons in Hematite Fe 2 O 3

We present a computational study of the electronic structure of excitons in hematite Fe 2 O 3 using density functional theory (DFT) and time-dependent DFT theory. Excitons resulting from solar light absorption are precursors to photogenerated charge carriers, which are the active species in solar-to-fuel redox conversions. Upon vertical excitation, the lowest energy exciton is found to have a multideterminantal wave function, describing an electron excited from O 2p states to Fe 3d states of the nearest basal plane, possibly supporting a recent XUV-derived exciton structure with a radius of a single Fe–O bond. Subsequently, the exciton self-traps into stable electron–hole pair structures whose wave functions are well-approximated by single excitations conveniently labeled “HOMO–LUMO”. Self-trapping yields electron–hole pair structures whereby the electron state (electron polaron-like) is separated from the hole state (hole polaron-like) by 3, 5, 7, 9, ... basal planes in structures referred to as Exc-3, Exc-5, Exc-7, Exc-9, ... The natural transition orbitals (electron–hole states) exhibit a strong “localized” character. The hole state is best described as a (FeO 6 ) + octahedral moiety carrying ~0.65 h + charge, and overall ~70% of the hole is assigned to O 2p atomic states, mostly on the four (4) equatorial O atoms of the moiety. The electron state is best described as a (FeFe) – moiety carrying ~0.80 e – charge in Fe 3d atomic states. The lattice distortions around the hole site exhibit Fe–O bond shortening due to the removal of electron density from an orbital state with O–Fe–O antibonding character. The electron site shows Fe–Fe distance shortening due to excess electron density added to an orbital state with in-phase Fe–Fe interactions. Furthermore, these excitonic structures with increasing electron–hole separation can be viewed as the onset of exciton separation into isolated electrons and hole polarons.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Graphically Contracted Function Construction with the Recursive Pairwise Merge Algorithm

A new implementation of a recursive pairwise merge algorithm to construct a GCF from a list of CSF expansion coefficients is presented. The essential new feature is the preallocation of some work arrays used within the intermediate steps of the merge procedure. This results in roughly an order of magnitude improvement in overall efficiency and also approximately eliminates a factor of n, the molecular orbital dimension, from the original implementation. Initial application of this merge procedure to a series of H m molecules shows that the GCF wave functions can be represented well both with delocalized canonical Hartree-Fock orbitals and with localized molecular orbitals. Finally, for a given wave function complexity, as measured by the average facet count, $\bar{\text {f}}$, the delocalized Hartree-Fock orbitals show smaller errors for small $\bar{\text {f}}$ values, while the localized orbitals show smaller errors for larger $\bar{\text {f}}$ values.

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

Holographic View of Mixed-State Symmetry-Protected Topological Phases in Open Quantum Systems

We establish a holographic duality between d -dimensional mixed-state symmetry-protected topological (mSPT) phases and ( d + 1 ) -dimensional subsystem symmetry-protected topological (SSPT) states. Specifically, we show that the reduced density matrix of the boundary layer of a ( d + 1 ) -dimensional SSPT state with subsystem symmetry S and global symmetry G corresponds to a d -dimensional mSPT phase with strong S and weak G symmetries. Conversely, we demonstrate that the wave function of an SSPT state can be constructed by replicating the density matrix of the corresponding lower-dimensional mSPT phase. This mapping links the density matrix in lower dimensions to the entanglement properties of higher-dimensional wave functions, providing an approach for analyzing nonlinear quantities and quantum information metrics in mixed-state systems. Our duality offers a new perspective for studying intrinsic mSPT phases that are unique to open quantum systems, without pure-state analogues. We show that strange correlators and twisted Rényi- N correlators can diagnose these nontrivial phases and we explore their connection to strange correlators in pure-state SSPT phases. Furthermore, we discuss several implications of this holographic duality, including a method for preparing intrinsic mSPT phases through the duality. Published by the American Physical Society 2025

Sun, Shijun (ORCID:0000000168880030)↗
DOE OSTI2022

On the role of dynamic electron correlation in non-orthogonal configuration interaction with fragments

Two different approaches have been implemented to include the effect of dynamic electron correlation in the Non-Orthogonal Configuration Interaction for Fragments (NOCI-F) method. The first is based on shifting the diagonal matrix elements of the NOCI matrix, while the second incorporates the dynamic correlation explicitly in the fragment wave functions used to construct the many-electron basis functions of the NOCI. The two approaches are illustrated for the calculation of the electronic coupling relevant in singlet fission and the coupling of spin moments in organic radicals. Comparison of the calculated diabatic couplings, the NOCI energies and wave functions shows that dynamic electron correlation is not only efficiently but also effectively incorporated by the shifting approach and can largely affect the coupling between electronic states. Additionally, it brings the NOCI coupling of the spin moments in close agreement with benchmark calculations.

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

Pentaquarks made of light quarks and their admixture to baryons

This paper is a continuation of our studies of multiquark hadrons. The antisymmetrization of their wave functions required by Fermi statistics is nontrivial, as it mixes orbital, color, spin, and flavor structures. In our previous papers we developed a method to find them based on the representations of the permutation group, and derived the explicit wave functions for baryons excited to the first and second shells (L = 1, 2), tetraquarks $qq$$\overline{q}$$\overline{q}$ and hexaquarks (6q). Now we apply it to light pentaquarks ($qqq$$\overline{q}$), in the S- and P-shells (L = 0, 1). Using Jacobi coordinates, one can use the hyperdistance approximation in 12-dimensional space. We further address the issue of “unquenching” of baryons, by considering their mixing with pentaquarks, via two channels, through the addition of σ-like or π-like $\overline{q}$$q$ pairs. This mixing is central for understanding of the observed flavor asymmetry of the antiquark sea, the amount of orbital motion issue as well as other nucleon properties.

Baryons↗
DOE OSTI2025

Diffusion Monte Carlo Study of the Structure and Spectroscopy of H 3 O –

A potential energy surface for H 3 O – has been developed based on the NN+(MOB-ML) approach we developed for studies of complexes of OH – with two or three water molecules. Unlike those systems, H 3 O – has two low-energy isomers, H – ·H 2 O and OH – ·H 2 , which differ in energy by less than 2.5 kcal mol –1 , and which are separated by a barrier of roughly 4.5 kcal mol –1 . We find that by training the NN+(MOB-ML) model using structures based on diffusion Monte Carlo (DMC) simulations initiated in the two potential minima, we are able to obtain a potential surface that describes both isomers. Using these potentials, the structure and spectra of H 3 O – and its deuterated analogues are investigated using DMC. These calculations show that the ground state wave function for H 3 O – is mainly localized in the H – ·H 2 O minimum in the potential, with a small amount of the probability amplitude (<5%) in the region of the OH – ·H 2 minimum. The delocalization of the wave function into the secondary minimum is lowered by deuteration of the water molecule, while replacing H – with D – increases the isomerization due to the shortening of the average distance between the hydride ion and the hydrogen atom in water to which it is bound. Introducing one quantum of excitation in the H – ···H 2 O stretching vibration increases the amount of isomerization, while the isomerization decreases with additional excitation of this mode. Excitation of the free OH stretch in the water molecule also increases the amount of isomerization, while excitation of the out-of-plane bending vibration suppresses the isomerization. Furthermore, the effects of partial deuteration on the frequencies for these vibrations are also explored.

Chemical structure↗
NASA NTRS1978

Study of diatomic molecules. 2: Intensities

The theory of perturbations, giving the diatomic effective Hamiltonian, is used for calculating actual molecular wave functions and intensity factors involved in transitions between states arising from Hund's coupling cases a,b, intermediate a-b, and c tendency. The Herman and Wallis corrections are derived, without any knowledge of the analytical expressions of the wave functions, and generalized to transitions between electronic states with whatever symmetry and multiplicity. A general method for studying perturbed intensities is presented using primarily modern spectroscopic numerical approaches. The method is used in the study of the ScO optical emission spectrum.

Femenias, J. L.↗
DOE OSTI2025

Impact of hole polaron formation on excitonic transitions in MgO from first principles

Here, we present a first-principles investigation of the excitonic properties of magnesia (MgO), an ionic insulator known to host hole polarons. We combine a density functional theory-based approach for structural relaxation in the presence of the hole and many-body perturbation theory to describe the excitonic properties. We determine that the hole polaron introduces new in-gap occupied states 0.6–0.8 eV above the valence band maximum that lead to two low-energy peaks in the optical spectrum. The predicted redshift of the lowest-energy transition due to polaron formation of 0.8 eV agrees well with the experimental Stokes shift of 0.8–0.9 eV. Analysis of the exciton wave function indicates that the electron-hole pair consists of a localized hole and delocalized electron, but that the wave function retains its Wannier-Mott character even in the presence of the hole polaron. Our study demonstrates that combining these previously established methods allows for a relatively computationally inexpensive approach to studying the exciton polaron in materials where only one charge carrier forms a polaron.

electronic structure↗
DOE OSTI2020

Localized excitons in defective monolayer germanium selenide

Germanium selenide (GeSe) is a van der Waals–bonded layered material with promising optoelectronic properties, which has been experimentally synthesized for two-dimensional (2D) semiconductor applications. In the monolayer, due to reduced dimensionality and, thus, screening environment, perturbations such as the presence of defects have a significant impact on its properties. We apply density functional theory and many-body perturbation theory to understand the electronic and optical properties of GeSe containing a single selenium vacancy in the –2 charge state. We predict that the vacancy results in midgap “trap states” that strongly localize the electron and hole density and lead to sharp, low-energy optical absorption peaks below the predicted pristine optical gap. Analysis of the exciton wave function reveals that the 2D Wannier-Mott exciton of the pristine material monolayer is strongly modified by the presence of the defect. The lowest-energy exciton wave function is highly localized around the defect, with the Bohr radius reduced by a factor of 4 when compared to the pristine material and a dipole moment along the out-of-plane axis due to the defect-induced symmetry breaking. Altogether, these results suggest that the vacancy significantly modifies the excitonic properties of the system, demonstrating the importance of considering defects in the context of material design.

2-dimensional systems↗
DOE OSTI2020

QCD hidden-color hexadiquark in the core of nuclei

Hidden-color configurations are a key prediction of QCD with important physical consequences. In this work we examine a QCD color-singlet configuration in nuclei formed by combining six scalar [ud] diquarks in a strongly bound SU (3) C channel. The resulting hexadiquark state is a charge-2, spin-0, baryon number-4, isospin-0, color-singlet state. It contributes to alpha clustering in light nuclei and to the additional binding energy not saturated by ordinary nuclear forces in as well as the alpha-nuclei sequence of interest for nuclear astrophysics. We show that the strongly bound combination of six scalar isospin-0 [ud] diquarks within the nuclear wave function - relative to free nucleons - provides a natural explanation of the EMC effect measured by the CLAS collaboration's comparison of nuclear parton distribution function ratios for a large range of nuclei. Overall, these experiments confirmed that the EMC effect; i.e., the distortion of quark distributions within nuclei, is dominantly identified with the dynamics of neutron-proton (“isophobic”) short-range correlations within the nuclear wave function rather than proton-proton or neutron-neutron correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantum Simulation of Molecular Dynamics Processes─A Benchmark Study Using a Classical Simulator and Present-Day Quantum Hardware

Here, we explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM’s superconducting qubits and IonQ’s trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Electron-Hydrogen Elastic Scattering

Scattering by single-electron systems is always of interest because the wave function of the target is known exactly. Various approximations have been employed to take into account distortion produced in the target. Among them are the method of polarized orbitals and the close coupling approximation. Recently, e-H and e-He+ S-wave scattering in the elastic region has been studied using the Feshbach projection operator formalism. In this approach, the usual Hartree-Fock and exchange potentials are augmented by an optical potential and the resulting phase shifts have rigorous lower bounds. Now this method is being applied to the e-H P-wave scattering in the elastic region. The number of terms in the Hylleraas-type wave function for the 1,3 P phase shifts is 84 and the resulting phase shifts (preliminary) are given. The results have been given up to five digits because to that accuracy they are rigorous lower bounds. They are in general agreement with the variational (VAR) results of Armstead, and those obtained from the intermediate energy R-matrix method (RM) of Scholz et al., and the finite element method (FEM) of Botero and Shertzer. The later two methods do not provide any bounds on phase shifts.

Bhatia, A. K.↗
DOE OSTI2023

Efficient solutions of fermionic systems using artificial neural networks

In this study, we explore the similarities and differences between variational Monte Carlo techniques that employ conventional and artificial neural network representations of the ground-state wave function for fermionic systems. Our primary focus is on shallow neural network architectures, specifically the restricted Boltzmann machine, and we examine unsupervised learning algorithms that are appropriate for modeling complex many-body correlations. We assess the advantages and drawbacks of conventional and neural network wave functions by applying them to a range of circular quantum dot systems. Our findings, which include results for systems containing up to 90 electrons, emphasize the efficient implementation of these methods on both homogeneous and heterogeneous high-performance computing facilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Bardeen-Cooper-Schrieffer pairing of composite fermions

Topological pairing of composite fermions has led to remarkable ideas, such as excitations obeying non-Abelian braid statistics and topological quantum computation. Here, we construct a p-wave paired Bardeen-Cooper-Schrieffer (BCS) wave function for composite fermions in the torus geometry, which is a convenient geometry for formulating momentum space pairing as well as for revealing the underlying composite-fermion Fermi sea. Following the standard BCS approach, we minimize the Coulomb interaction energy at half filling in the lowest and the second Landau levels, which correspond to filling factors ν = 1/2 and ν = 5/2 in GaAs quantum wells, by optimizing two variational parameters that are analogous to the gap and the Debye cut-off energy of the BCS theory. Our results show no evidence for pairing at ν = 1/2 but a clear evidence for pairing at ν = 5/2. To a good approximation, the highest overlap between the exact Coulomb ground state at ν = 5/2 and the BCS state is obtained for parameters that minimize the energy of the latter, thereby providing support for the physics of composite-fermion pairing as the mechanism for the 5/2 fractional quantum Hall effect. We discuss the issue of modular covariance of the composite-fermion BCS wave function, and calculate its Hall viscosity and pair correlation function. By similar methods, we look for but do not find an instability to s-wave pairing for a spin-singlet composite-fermion Fermi sea at half-filled lowest Landau level in a system where the Zeeman splitting has been set to zero.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Charge ordering in Ir dimers in the ground state of Ba 5 AlIr 2 O 11

It has been well established experimentally that the interplay of electronic correlations and spin-orbit interactions in Ir 4+ and Ir 5+ oxides results in insulating J eff = 1/2 and J eff = 0 ground states, respectively. However, in compounds where the structural dimerization of iridium ions is favorable, the direct Ir d-d hybridization can be significant and takes a key role. Here, in this study, we investigate the effects of direct Ir d-d hybridization in comparison with electronic correlations and spin-orbit coupling in Ba 5 AlIr 2 O 11 , a compound with Ir dimers. Using a combination of ab initio many-body wave-function quantum chemistry calculations and resonant inelastic x-ray scattering experiments, we elucidate the electronic structure of Ba 5 AlIr 2 O 11 . We find excellent agreement between the calculated and the measured spin-orbit excitations. Contrary to expectations, the analysis of the many-body wave function shows that the two Ir (Ir 4+ and Ir 5+ ) ions in the Ir 2 O 9 dimer unit in this compound preserve their local J eff character close to 1/2 and 0, respectively. The local point group symmetry at each of the Ir ions plays an important role, significantly limiting the direct d-d hybridization. Our results emphasize that minute details in the local crystal field environment can lead to dramatic differences in the electronic states in iridates and 5d oxides in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

Relativistic calculation of atomic M-shell ionization by protons

Relativistic plane-wave Born-approximation calculations of cross sections for M-shell ionization of Ho-67, Au-79, U-92 by protons with incident energies from 0.05 to 1 MeV are reported. Dirac-Hartree-Slater wave functions were employed and binding-energy change and Coulomb deflection were taken into account. Associated X-ray production cross sections were also computed. Results are compared with previous theoretical predictions and with experimental data. Definite improvement in the theory has been attained by the use of realistic wave functions and consistent inclusion of the effects of relativity.

Chen, M. H.↗
NASA NTRS1987

Oscillator strengths for S I, S II, and S III

A series of calculations for atomic data of various sulfur and oxygen ions is examined. Recent observations of the Io plasma torus obtained with the Voyager UV Spectrometer, the IUE satellite short wavelength spectrograph, and the rocket-borne faint object telescope are discussed. The calculation of oscillator strengths for S II, the P I sequence, S I, and S III in terms of configuration interaction effects is described. The derivation of orbital wave functions is considered. The use of the close coupling method to estimate collision strengths is studied. The accuracy of these calculations depends on: (1) the number of states used in the close coupling expansion; (2) resonance contributions to the thermally averaged collision strength; and (3) the quality of the target state wave functions. Tables of the derived oscillator strengths are presented.

Ho, Y. K.↗