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

Seniority eigenstate configuration interaction

Zero-seniority methods have shown great promise for the description of strongly correlated electronic systems. Other seniority sectors have been much less explored, and in particular, the maximal seniority sector and zero seniority have the same underlying algebraic structure. We introduce a seniority eigenstate configuration interaction in which the wave function is constrained to have good fixed local seniority for each paired orbital, by which we mean we partition orbitals into a pairing set with seniority zero, and a spin set with seniority one. Here, we show how to build the effective Hamiltonian for this ansatz, and demonstrate that high-seniority wave functions have unexpectedly excellent accuracy for strongly correlated fermionic systems, with accuracy competitive with or better than seniority zero for the Hubbard model and for the dissociation of the nitrogen molecule.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Single-Determinant Ground State in Ce 4+ Imidophosphorane Complexes

X-ray spectroscopy techniques are critical in the electronic structure analysis of high-valent lanthanides. The interpretation of multipeaked features at the lanthanide L 3 -edge has remained a challenging question, as it is observed across a range of material classes. A series of structurally related Ce 4+ complexes were prepared to probe the potential ligand field perturbation of the ground state within the series. The tuning of relative 4f and ligand orbital energies in homoleptic and heteroleptic tetravalent Ce imidophosphorane complexes is achieved through ligand derivatization and is clearly demonstrated by UV−vis spectroscopy and electrochemically measured redox potentials. However, Ce L 3 - edge high-energy resolution fluorescence-detected (HERFD) X-ray absorption near-edge structure (XANES) spectra present features at consistent numbers and energies across the range of complexes. Resonant inelastic X-ray scattering (RIXS) is employed to visualize the observed features in the HERFD-XANES spectra. Large complete active space configuration interaction singles and doubles (CASCISD) calculations demonstrate that the ground-state wave function of all complexes can be described employing a single determinant. As a result, the multielectron feature at the L 3 -edge observed in this study for the Ce 4+ imidophosphorane complexes is described as excited-state multiconfigurational behavior that is independent of ligand variation, in systems where the ground state is described using the simple single determinant wave function.

Lanthanides↗
DOE OSTI2025

Applications of the modified Hulthén-Kohn method for bound and scattering states

We adapt the Hulthén–Kohn method suggested by Efros [Phys. Rev. C 99, 034620 (2019)] for calculating various observables in the continuum and discrete spectrum using two-body interactions in single- and coupled channel systems. We explore the convergence of phase shifts and wave functions as well as the location of S-matrix poles which enables obtaining both resonance and bound state parameters. We find that employing a harmonic oscillator basis, together with an interaction smoothing scheme introduced by Gyarmati et al. [Nucl. Phys. A 326, 119 (1979)], and adopting approximate bound-state solutions for the short-range components of basis wave functions lead to good convergence even with restricted oscillator quanta accessible for modern no-core shell model codes. The adapted Efros method will facilitate ab initio many-body nuclear structure applications.

Nuclear reactions↗
DOE OSTI2024

Spatial imaging of polarized deuterons at the Electron-Ion Collider

We study diffractive vector meson production at small-x in the collision of electrons and polarized deuterons e + d $\uparrow$ . We consider the polarization dependence of the nuclear wave function of the deuteron, which results in an azimuthal angular dependence of the produced vector meson when the deuteron is transversely polarized. The Fourier coefficients extracted from the azimuthal angular dependence of the vector meson differential cross-section exhibit notable differences between longitudinally and transversely polarized deuterons. The angular dependence of the extracted effective deuteron radius provides direct insight into the structure of the polarized deuteron wave function. Furthermore, we observe slightly increased gluon saturation effects when the deuteron is longitudinally polarized compared to the transversely polarized case. The small-x observables studied in this work will be accessible at the future Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Observable measurement-induced transitions

One of the main postulates of quantum mechanics is that measurements destroy quantum coherence (wave function collapse). Recently it was discovered that in a many-body system dilute local measurements still preserve some coherence across the entire system. As the measurement density is increased, a phase transition occurs that is characterized by the disentanglement of different parts of the system. Unfortunately, this transition is impossible to observe experimentally for macroscopic systems because it requires an exponentially costly full tomography of the many-body wave function or a comparison with the simulation on an oracle classical computer. In this work we report the discovery of another measurement-induced phase transition that can be observed experimentally if quantum dynamics can be reversed. On one side of this phase transition the quantum information encoded in some part of the Hilbert space is fully recovered after the time inversion. On the other side, all quantum information is corrupted. This transition also manifests itself as the change in the behavior of the probability to observe the same measurement outcome in the process that consists of identical blocks repeated many times. In each block the unitary evolution is followed by the measurement. On one side of the transition the probability decreases exponentially with the number of repetitions, on the other it tends to a constant as the number of repetitions is increased. We confirm the existence of the proposed phase transition through numerical simulations of realistic quantum circuits and analytical calculations using an effective random-matrix theory model.

Measurement-induced phase transitions↗
DOE OSTI2024

A new-old approach to composite scalars with chiral fermion constituents

We develop a dynamical, Lorentz invariant theory of composite scalars in configuration space consisting of chiral fermions, interacting by the perturbative exchange of a massive “gluon” of coupling g 0 and mass M 0 2 (the coloron model). The formalism is inspired by, but goes beyond, old ideas of Yukawa and the Nambu-Jona-Lasinio (NJL) model. It yields a non-pointlike internal wave-function of the bound state, ϕ ( r ) , which satisfies a Schrödinger-Klein-Gordon (SKG) equation with eigenvalue μ 2 . For super-critical coupling, g 0 2 > g 0 c 2 , we have μ 2 < 0 leading to spontaneous symmetry breaking. The binding of chiral fermions is semiclassical, not loop-level as in NJL. The mass scale is determined by the interaction as in NJL. We mainly focus on the short-distance, large M 0 2 limit, yielding an NJL pointlike interaction, but the bound state internal wave-function, ϕ ( r → ) , remains spatially extended and dilutes ϕ ( 0 ) . This leads to power-law suppression of the induced Yukawa and quartic couplings and requires radically less fine-tuning of a hierarchy than does the NJL model. We include a discussion of loop corrections of the theory. A realistic top condensation model appears possible.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Reproducibility of fixed-node diffusion Monte Carlo across diverse community codes: The case of water–methane dimer

Fixed-node diffusion quantum Monte Carlo (FN-DMC) is a widely trusted many-body method for solving the Schrödinger equation, known for its reliable predictions of material and molecular properties. Furthermore, its excellent scalability with system complexity and near-perfect utilization of computational power make FN-DMC ideally positioned to leverage new advances in computing to address increasingly complex scientific problems. Even though the method is widely used as a computational gold standard, reproducibility across the numerous FN-DMC code implementations has yet to be demonstrated. This difficulty stems from the diverse array of DMC algorithms and trial wave functions, compounded by the method’s inherent stochastic nature. Here, this study represents a community-wide effort to assess the reproducibility of the method, affirming that yes, FN-DMC is reproducible (when handled with care). Using the water–methane dimer as the canonical test case, we compare results from eleven different FN-DMC codes and show that the approximations to treat the non-locality of pseudopotentials are the primary source of the discrepancies between them. In particular, we demonstrate that, for the same choice of determinantal component in the trial wave function, reliable and reproducible predictions can be achieved by employing the T-move, the determinant locality approximation, or the determinant T-move schemes, while the older locality approximation leads to considerable variability in results. These findings demonstrate that, with appropriate choices of algorithmic details, fixed-node DMC is reproducible across diverse community codes—highlighting the maturity and robustness of the method as a tool for open and reliable computational science.

Della Pia, Flaviano [Univ. of Cambridge (United Ki↗
DOE OSTI2025

Quantum magic and multipartite entanglement in the structure of nuclei

Motivated by the Gottesman-Knill theorem, we present a detailed study of the quantum complexity of p -shell and s d -shell nuclei. Valence-space nuclear shell-model wave functions generated by the code are mapped to qubit registers using the Jordan-Wigner mapping (12 qubits for the p shell and 24 qubits for the s d shell), from which measures of the many-body entanglement ( n -tangles) and magic (nonstabilizerness) are determined. While exact evaluations of these measures are possible for nuclei with a modest number of active nucleons, Monte Carlo simulations are required for the more complex nuclei. The broadly applicable Pauli-string I ̂ Z ̂ exact (PSIZe) Markov chain Monte Carlo (MCMC) technique is introduced to accelerate the evaluation of measures of magic in deformed nuclei (with hierarchical wave functions), by factors of ≈ 8 for some nuclei. Significant multinucleon entanglement is found in the s d shell, dominated by proton-neutron configurations, along with significant measures of magic. This is evident not only for the deformed states, but also for nuclei on the path to instability via regions of shape coexistence and level inversion. These results indicate that quantum-computing resources will accelerate precision simulations of such nuclei and beyond. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantum entanglement in nuclear fission

Nuclear fission presents a unique example of quantum entanglement in strongly interacting many-body systems. A heavy nucleus can split into hundreds of combinations of two complementary fragments in the fission process. The entanglement of fragment wave functions is persistent even after separation and impacts the partition of particles and energies between fragments. Based on microscopic dynamical calculations of the fission of 240 Pu, this work finds that dynamical quantum entanglement is indispensable in the appearance of sawtooth distributions of average excitation energies of fragments and thus neutron multiplicities, but not in average neutron excess of fragments. Both sawtooth slopes from particle-number projections are found to be steep – a feature which can be alleviated by random fluctuations. The persistent entanglement is mainly due to non-adiabatic dynamics since the final splitting is so fast that the non-localization of wave functions is kept during the separation. These findings may impact the understanding of quantum entanglement more broadly in mesoscopic systems

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Fragment-based initialization for quantum subspace methods

Here, we present a novel quantum-classical algorithm called LAS-QKSD for multireference systems, by combining a classical localized active space (LAS) fragment-based multireference algorithm with the quantum Krylov subspace diagonalization (QKSD) method for quantum computers. The algorithm uses wave function information from a LAS self-consistent field (LASSCF) calculation to prepare an initial state with better overlap with the target ground state than the Hartree-Fock state. This is coupled with the use of QKSD to ultimately converge to the exact energy, providing faster convergence than starting from the Hartree-Fock state. Fragmentation has the two-fold benefit of fewer configurations on the classical side of the algorithm as well as fewer state preparation gates on the quantum side. First, we compare the LAS-QKSD method to the classical LASSCF method and to QKSD with a Hartree-Fock initial state. We then examine ways to load the LASSCF wave function using direct initialization and a QKSD-motivated spectral filtering approach. Finally, using a bimetallic complex, we show that the LAS-QKSD method is an efficient alternative to highly expensive complete active space SCF (CASSCF) calculations on strongly correlated systems.

D'Cunha, Ruhee↗
DOE OSTI2025

Unexpected origin of the quantum efficiency reduction in long-wavelength (In,Ga)⁢N light-emitting diodes

Differential carrier lifetime (DCL) measurements were performed on c-plane In 𝑥 ⁢Ga 1−𝑥⁢ N/Ga⁢N single-quantum-well (QW) light-emitting diodes (LEDs) with varying indium-content QWs (x = 13.5%, 16%, 22%), emitting with violet, blue, and green wavelengths. The recombination lifetimes of LEDs were found to increase with increasing indium composition, resulting in increased carrier densities n measured by DCL. Extraction of the A coefficients, which are assumed to not vary with n, and of the effective B(n) and C(n) coefficients of the ABC model of the internal quantum efficiency (IQE) of QWs showed no significant changes in the A coefficient with increasing indium content [In] in the In 𝑥 ⁢Ga 1−𝑥 ⁢N QW, and a reduction in the B(n) and C(n) coefficients with increasing [In]. When looking at the Shockley-Read-Hall (SRH) recombination rate An, the radiative recombination rate B(n)𝑛 2 , and the Auger-Meitner (AM) recombination rate C(n)𝑛 3 , we observed that at any given n, the SRH rate is the same for the three [In] measured, while the radiative rate decreases (by up to approximately 9 times) and the AM rate decreases (by up to approximately 7 times) with increasing indium content. This shows that the larger reduction in radiative recombination rates relative to nonradiative recombination rates with increasing [In] is the largest contributor to the decreased quantum efficiencies of LEDs (the “green gap”) at any given n. While some of the reduction in the effective recombination coefficients B(n) and C(n) can be explained by the reduced wave-function overlaps of the QW with increasing indium content, the larger reduction of B(n) relative to C(n) motivates further study of the intrinsic recombination coefficients 𝐵 0 and 𝐶 0 of bulk In 𝑥 ⁢Ga 1−𝑥 ⁢N alloys and how they are related to the effective recombination coefficients B(n) and C(n) in QWs. The relative contributions of nonradiative recombination are enhanced at any given current density due to the sublinear relationship between the carrier density n and the current density. Thus, the solution to the green gap, up to any intrinsic limitations set by the indium content of the QW, from an IQE perspective, requires the design of long-wavelength (In,Ga)⁢N LEDs with improved wave-function overlaps that can operate with a lower carrier density at any given current density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2024

Dynamical model of J/ψ photoproduction on the nucleon

Here, a dynamical model based on a phenomenological charm quark-nucleon (c – N) potential v cN and the Pomeronexchange mechanism is constructed to investigate the J/ψ photoproduction on the nucleon from threshold to invariant mass W = 300 GeV. The J/ψ – N potential, V J/ψN (r), is constructed by folding v cN into the wave function φ J/ψ ($c\bar{c}$) of J/ψ within a constituent quark model (CQM) of Segovia et al. [Int. J. Mod. Phys. E 22, 1330026 (2013)]. A photoproduction amplitude is also generated by v cN by a $c\bar{c}$–loop integration over the γ → $c\bar{c}$ vertex function and φ J/ψ ($c\bar{c}$). No commonly used vector meson dominance assumption is used to define this photoproduction amplitude which is needed to describe the data near the threshold. The c – N potential v cN (r) is parameterized in a form such that the predicted V J/ψN (r) at large distances has the same Yukawa potential form extracted from a lattice QCD (LQCD) calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)]. The parameters of vcN are determined by fitting the total cross-section data of Jefferson Laboratory (JLab) by performing calculations that include J/ψ – N final-state interactions (FSI). The resulting differential cross sections dσ /dt are found in good agreements with the data. It is shown that the FSI effects dominate the cross section in the very near-threshold region, allowing for sensitive testing of the predicted J/ψ – N scattering amplitudes. By imposing the constraints of J/ψ – N potential extracted from the LQCD calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)], we have obtained three J/ψ – N potentials which fit the JLab data equally well. The resulting J/ψ – N scattering lengths are in the range of a = [-0.05, -0.25] fm. With the determined v cN (r) and the wave functions generated from the same CQM, the constructed model is used to predict the cross sections of photoproduction of η c (1S) and ψ(2S) mesons for future experimental tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Muon capture on Li 6 , C 12 , and O 16 from ab initio nuclear theory

Muon capture on nuclei is one of the most promising probes of the nuclear electroweak current driving the yet-hypothetical neutrinoless double-beta (0νββ) decay. Both processes involve vector and axial-vector currents at finite momentum transfer, q ~ 100 MeV, as well as the induced pseu doscalar and weak-magnetism currents. Comparing measured muon-capture rates with reliable ab initio nuclear-theory predictions could help us validate these currents. To this end, we compute partial muon-capture rates for 6 Li, 12 C and 16 O, feeding the ground and excited states in 6 He, 12 B and 16 N, using ab initio no-core shell model with two- and three-nucleon chiral interactions. Here, we remove the spurious center-of-mass motion by introducing translationally invariant operators and approximate the effect of hadronic two-body currents by Fermi-gas model. We solve the bound-muon wave function from the Dirac wave equations in the Coulomb field created by a finite nucleus. We find that the computed rates to the low-lying states in the final nuclei are in good agreement with the measured counterparts. We highlight sensitivity of some of the transitions to the sub-leading three-nucleon interaction terms. We also compare summed rates to several tens of final states with the measured total capture rates and note that we slightly underestimate the total rate with this simple approach due to limited range of excitation energies

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Attosecond X-Ray Core-Level Chronoscopy of Aromatic Molecules

Attosecond photoemission or photoionization delays are a unique probe of the structure and the electronic dynamics of matter. However, the spectral congestion of valence photoelectron spectra sets fundamental limits to the complexity of systems that can be studied, and the delocalization of valence electron wave functions blurs the spatial origin of the photoelectron wave packet. Using attosecond x-ray pulses from LCLS, we demonstrate the key advantages of measuring core-level delays: The photoelectron spectra remain atomlike, the measurements become element specific, and the observed scattering dynamics originate from a pointlike source when multicenter interference effects are negligible. We exploit these unique features to reveal the effects of changing functional groups (C-H vs N) and symmetry on attosecond scattering dynamics by measuring and calculating the photoionization delays between N−1⁢𝑠 and C−1⁢𝑠 core shells of a series of aromatic azabenzene molecules. Remarkably, the delays increase with the number of nitrogen atoms in the molecule and reveal multiple resonances. We identify two previously unknown mechanisms regulating the associated attosecond dynamics, namely the enhanced confinement of the trapped wave function with the replacement of C-H groups by N atoms and the decrease of the coupling strength among the photoemitted partial waves with increasing symmetry. This study demonstrates the unique opportunities opened by measurements of core-level photoionization delays for unraveling attosecond electron dynamics in complex matter.

Ji, Jia-Bao [Eidgenoessische Technische Hochschule↗