Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Entangled phase of simultaneous fermion and exciton condensations realized

Fermion-exciton condensates (FECs)—computationally and theoretically predicted states that simultaneously exhibit the character of superconducting states and exciton condensates—are novel quantum states whose properties may involve a hybridization of superconductivity and the dissipationless flow of energy. Here, we exploit prior investigations of superconducting states and exciton condensates on quantum devices to identify a tuneable quantum state preparation entangling the wave functions of the individual condensate states. Utilizing this state preparation, we prepare a variety of FEC states on quantum computers—realizing strongly correlated FEC states on current, noisy intermediate-scale quantum devices—and verify the presence of the dual condensate via postmeasurement analysis. This confirmation of the previously predicted condensate state on quantum devices as well as the form of its wave function motivates further theoretical and experimental exploration of the properties, applications, and stability of FECs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantum Many-Body Theory from a Solution of the N -Representability Problem

Here, in this study, we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H 8 ring.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2022

Full radiator-perturber interaction in computer simulations of hydrogenic spectral line broadening by plasmas

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating potential obtained from computer simulations. Such calculations have relied on the multipole expansion of the radiator-plasma interaction often keeping only the dipole term. With increasing density, however, higher multipoles as well as plasma perturbers overlapping the bound electron wave functions are expected to become important. For hydrogenic systems, the atomic matrix elements of the full Coulomb and screened Coulomb interactions are given by analytical formulas. Using these results, a computer simulation approach that accounts for the full radiator-plasma interaction is developed. One benefit is the removal of inherent strong collision divergences in the multipole expansion approximation. Furthermore, it yields the plasma polarization shift produced by perturbers penetrating the wave function of the radiator bound electrons. Here, the model was applied to hydrogenic argon Ly-α , Ly- β , and Ly-$\gamma$ spectral lines in a dense argon plasma at free electron densities of 10 24 or 10 25 cm ⁻3 and temperature of $800$ $eV$ relevant to plasma diagnostic techniques for inertial confinement fusion implosions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Weak Transitions in Light Nuclei

Nuclei are used for high-precision tests of the Standard Model and for studies of physics beyond the Standard Model. Without a thorough understanding of nuclei, we will not be able to meaningfully interpret the growing body of experimental data nor will we be able to disentangle new physics signals from underlying nuclear effects. This calls for accurate calculations of nuclear structure and reactions. In this work, we focus on electroweak decays in nuclei with mass number A ≤ 10 and report on ab initio Quantum Monte Carlo calculations of reduced matrix elements entering beta decays and electron captures in nuclei with mass number A ≤ 10. The many-body wave functions are calculated using selected Norfolk two- and three-nucleon potential models and associated one- and two-body axial currents at tree-level obtained from a chiral effective field theory with pions, nucleons, and Δ. The agreement with the experimental data is satisfactory except for transitions in A = 8 nuclei. In this specific case, the theory significantly underpredicts the experimental data, which indicates the need of further improvements in the corresponding nuclear wave functions. In this study, emphasis is placed on the contributions of two-body axial currents that are carefully analyzed using two-body transition densities. This allow us to study the spatial distribution and short-range behavior of two-body dynamics. In particular, the transition densities when scaled to peak at 1.0 exhibit universal short-range behavior across the considered nuclei, while they differ in the long-range tails.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Fractional Statistics

The quantum-mechanical description of assemblies of particles whose motion is confined to two (or one) spatial dimensions offers many possibilities that are distinct from bosons and fermions. We call such particles anyons. The simplest anyons are parameterized by an angular phase parameter θ. θ = 0, π correspond to bosons and fermions, respectively; at intermediate values, we say that we have fractional statistics. In two dimensions, θ describes the phase acquired by the wave function as two anyons wind around one another counterclockwise. It generates a shift in the allowed values for the relative angular momentum. Composites of localized electric charge and magnetic flux associated with an abelian U(1) gauge group realize this behavior. More complex charge-flux constructions can involve nonabelian and product groups acting on a spectrum of allowed charges and fluxes, giving rise to nonabelian and mutual statistics. Interchanges of nonabelian anyons implement unitary transformations of the wave function within an emergent space of internal states. Anyons of all kinds are described by quantum field theories that include Chern–Simons terms. The crossings of one-dimensional anyons on a ring are unidirectional, such that a fractional phase θ acquired upon interchange gives rise to fractional shifts in the relative momenta between the anyons. The quasiparticle excitations of fractional quantum Hall states have long been predicted to include anyons. Recently, the anyon behavior predicted for quasiparticles in the ν = 1/3 fractional quantum Hall state has been observed in both scattering and interferometric experiments. Excitations within designed systems, notably including superconducting circuits, can exhibit anyon behavior. Such systems are being developed for possible use in quantum information processing.

Physics↗
DOE OSTI2021

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Single-ion properties of the transverse-field Ising model material CoNb 2 O 6

We report CoNb 2 O 6 is one of the few materials that is known to approximate the one-dimensional transverse-field Ising model (1D-TFIM) near its quantum critical point. It has been inferred that Co 2+ acts as a pseudospin 1/2 with anisotropic exchange interactions that are largely Ising-like, enabling the realization of the TFIM. However, the behavior of CoNb 2 O 6 is known to diverge from the ideal TFIM under transverse magnetic fields that are far from the quantum critical point, requiring the consideration of additional anisotropic, bond-dependent (Kitaev-like) terms in the microscopic pseudospin-1/2 Hamiltonian. These terms are expected to be controlled in part by single-ion physics, namely the wave function for the pseudospin-1/2 angular momentum doublet. Here, we present the results of both inelastic neutron scattering measurements and electron paramagnetic resonance spectroscopy on CoNb 2 O 6 , which elucidate the single-ion physics of Co 2+ in CoNb 2 O 6 for the first time. We find that the system is well-described by an intermediate spin-orbit coupled Hamiltonian, and the ground state is a well-isolated Kramers doublet with an anisotropic g tensor. We provide the approximate wave functions for this doublet, which we expect will be useful in theoretical investigations of the anisotropic exchange interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Systematic comparison and cross-validation of fixed-node diffusion Monte Carlo and phaseless auxiliary-field quantum Monte Carlo in solids

Quantum Monte Carlo (QMC) methods are some of the most accurate methods for simulating correlated electronic systems. Here, we investigate the compatibility, strengths, and weaknesses of two such methods, namely, diffusion Monte Carlo (DMC) and auxiliary-field quantum Monte Carlo (AFQMC). The multideterminant trial wave functions employed in both approaches are generated using the configuration interaction using a perturbative selection made iteratively (CIPSI) technique. Complete basis-set full configuration interaction energies estimated with CIPSI are used as a reference in this comparative study between DMC and AFQMC. By focusing on a set of canonical finite–size solid-state systems, we show that both QMC methods can be made to systematically converge towards the same energy once basis-set effects and systematic biases have been removed. AFQMC shows a much smaller dependence on the trial wave function than DMC while simultaneously exhibiting a much larger basis-set dependence. We outline some of the remaining challenges and opportunities for improving these approaches.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2023

Leakage Reduces Device Coherence Demands for Pulse-Level Molecular Simulations

Quantum simulation on noisy intermediate-scale quantum devices is severely limited by short qubit coherence times. A variational pulse-shaping algorithm known as ctrl-VQE was recently proposed to address this issue by eliminating the need for parameterized quantum circuits, which lead to long state preparation times. Here, in this study, we find the fastest possible pulses that prepare target molecular wave functions for a given device Hamiltonian describing coupled transmon qubits using simulations. We find that the resulting time-optimal pulses develop a bang-bang form consistent with Pontryagin’s maximum principle. We further investigate how the minimal state preparation time is impacted by the number of energy levels active in the transmon simulations. We find that leakage outside the computational subspace (something that is usually considered problematic) speeds up the state preparation, further reducing device coherence-time demands. Our analysis reveals that this speedup is due to both an enlarged solution space of target wave functions and the appearance of additional channels connecting initial and target states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Improved quasiparticle self-consistent electronic band structure and excitons in β – LiGaO 2

The band structure of β–LiGaO 2 is calculated using the quasiparticle self-consistent QSGWˆ method where the screened Coulomb interaction Wˆ is evaluated including electron-hole interaction ladder diagrams and G is the one-electron Green's function. Improved convergence compared to previous calculations leads to a significantly larger band gap of about 7.0 eV. However, exciton binding energies are found to be large and lead to an exciton gap of about 6.0 eV if also a zero-point-motion correction of about –0.4 eV is included. Furthermore, these results are in excellent agreement with recent experimental results on the onset of absorption. Besides the excitons observed thus far, the calculations indicate the existence of a Rydberg-like series of exciton excited states, which is however modified from the classical Wannier exciton model by the anisotropies of the material and the more complex mixing of Bloch states in the excitons resulting from the Bethe-Salpeter equation. The exciton fine structure and the exciton wave functions are visualized and analyzed in various ways.

36 MATERIALS SCIENCE↗
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 OSTI2021

Twisted bilayer graphene. V. Exact analytic many-body excitations in Coulomb Hamiltonians: Charge gap, Goldstone modes, and absence of Cooper pairing

We find exact analytic expressions for the energies and wave functions of the charged and neutral excitations above the exact ground states (at rational filling per unit cell) of projected Coulomb Hamiltonians in twisted bilayer graphene. Our exact expressions are valid for any form of the Coulomb interaction and any form of A A and A B / B A tunneling. The single charge excitation energy is a convolution of the Coulomb potential with a quantum geometric tensor of the TBG bands. The neutral excitations are (high-symmetry group) magnons, and their dispersion is analytically calculated in terms of the form factors of the active bands in TBG. The two-charge excitation energy and wave functions are also obtained, and a sufficient condition on the graphene eigenstates for obtaining a Cooper pair from Coulomb interactions is obtained. For the actual TBG bands at the first magic angle, we can analytically show that the Cooper pair binding energy is zero in all such projected Coulomb models, implying that either phonons and/or nonzero kinetic energy are needed for superconductivity. Since Vafek and Kang [Phys. Rev. Lett. 125, 257602 (2020)] showed that the kinetic energy bounds on the superexchange energy are less 10 -3 in Coulomb units, the phonon mechanism becomes then very likely. If nonetheless the superconductivity is due to kinetic terms which render the bands nonflat, one prediction of our theory is that the highest T c would not occur at the highest DOS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Extracting the Pion Distribution Amplitude from Lattice QCD through Pseudo-Distributions

The Light-Cone Distribution Amplitude (LCDA) encodes the non-perturbative information of the leading Fock component of the hadron wave function, therefore required for processes including exclusive hadron production. As the Pseudo-Nambu-Goldstone boson of QCD, the nonperturbative structure of the pion is of particular interest. Progress on the Lattice QCD calculation of the pion LCDA on O(a) -improved Wilson fermion ensembles at several lattice spacings is presented. Excited-state systematics are taken into account within a Bayesian Model Averaging framework. A Renormalization-Group-Invariant (RGI) ratio of matrix elements is formed for further extraction of the pion LCDA.

Kovner, Daniel↗