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 811 records · Page 45

Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework: Theory, implementation, and validation against multi-reference methods

In this work, we report an implementation of non-adiabatic coupling (NAC) forces within the equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) framework via the summed-state approach. Using illustrative examples, we compare NAC forces computed with EOM-CCSD and multi-reference (MR) wave functions (for selected cases, we also consider configuration interaction singles). In addition to the magnitude of the NAC vectors, we analyze their direction, which is important for the calculations of the rate of non-adiabatic transitions. Our benchmark set comprises three doublet radical-cations (hexatriene, cyclohexadiene, and uracil), neutral uracil, and sodium-doped ammonia clusters. When the characters of the states agree among different methods, we observe good agreement between the respective NAC vectors, both in the Franck-Condon region and away. In the cases of large discrepancies between the methods, the disagreement can be attributed to the difference in the states’ character, which, in some cases, is very sensitive to electron correlation, both within single-reference and multi-reference frameworks. Finally, the numeric results confirm that the accuracy of NAC vectors depends critically on the quality of the underlying wave functions. Within their domain of applicability, EOM-CC methods provide a viable alternative to MR approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effective matrix model for nuclear physics on a quantum computer

In this project we studied a simplified potential from [1] describing nuclear physics from quark and gluon fields, where spatial degrees of freedom have been integrated over. This potential deals with the phase structure of the strong interaction. We created a matrix representation of the Hamiltonian for this theory using discrete quantum mechanics [2]. We then feed this Hamiltonian to IBM's Qiskit package [3]. We run on IBM's 32-qubit quantum simulator, which approximates the ground state energy using the variational quantum eigensolver (VQE) algorithm [4]. This result is then compared with the exact ground state energy, which we can calculate classically since the matrices we use are still small enough. The bosonic Hamiltonian for a 6-qubit representation took over 1000 Pauli factors to represent fully. To improve our performance, we can look for a more efficient representation of this model's Hamiltonian. We will develop a better variational trial wave function and initial state, as well as include the effect of fermions and heavy quarks on the potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗

Entanglement rearrangement in self-consistent nuclear structure calculations

Background: Entanglement plays a central role in a diverse array of increasingly important research areas, including quantum computation, simulation, measurement, sensing, and communication. Extensive suites of investigations have been performed to better understand entanglement in atomic and molecular quantum many-body systems, while the exploration of entanglement in the structure of nuclei and their reactions is presently in its infancy. Purpose: The goal of this work is to begin investigating the entanglement properties of nuclei from first-principles nuclear many-body calculations. We attempt to identify common features and emergent structures of entanglement that could ultimately lead to new and natural many-body schemes. With an eye toward quantum accelerators in future hybrid-supercomputers, criteria for partitioning nuclear many-body calculations into quantum and classical components may provide advantages in future large-scale computations. Along the way we look for explanations of the relative success of phenomenological models such as the nuclear shell model, and for better ways to match to low-energy nuclear effective field theories and lattice QCD calculations to nuclear many-body techniques that are based upon entanglement. Method: We explore the entanglement between single-particle states in 4 He and 6 He. The patterns of entanglement emerging from different single-particle bases are compared, and possible links with the convergence of observables are explored, in particular, ground-state energies. The nuclear wave functions are obtained by performing active-space no-core configuration-interaction calculations using a two-body nucleon-nucleon interaction derived from chiral effective field theory. Entanglement measures within single-particle bases exhibiting different degrees of complexity are determined, in particular, harmonic oscillator (HO), Hartree-Fock (HF), natural (NAT) and variational natural (VNAT) bases. Specifically, single-orbital entanglement entropy, two-orbital mutual information, and negativity are studied. Results: The entanglement structures in 4 He and 6 He are found to be more localized within NAT and VNAT bases than within a HO basis for the optimal HO parameters we have worked with. In particular a core-valence structure clearly emerges from the full no-core calculation of 6 He. Here, the two-nucleon mutual information shows that the VNAT basis, which typically exhibits good convergence properties, effectively decouples the active and inactive spaces. Conclusions: Measures of one- and two-nucleon entanglement are found to be useful in analyzing the structure of nuclear wave functions, in particular the efficacy of basis states, and may provide useful metrics toward developing more efficient schemes for ab initio computations of the structure and reactions of nuclei, and quantum many-body systems more generally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamics of a qubit in a high-impedance transmission line from a bath perspective

We investigate the quantum dynamics of a generic model of light-matter interaction in the context of high-impedance waveguides, focusing on the behavior of the photonic states generated in the waveguide. The model treated consists simply of a two-level system coupled to a bosonic bath (the Ohmic spin-boson model). Quantum quenches as well as scattering of an incident coherent pulse are studied using two complementary methods. First, we develop an approximate ansatz for the electromagnetic waves based on a single multimode coherent state wave function; formally, this approach combines in a single framework ideas from adiabatic renormalization, the Born-Markov approximation, and input-output theory. Second, we present numerically exact results for scattering of a weak intensity pulse by using numerical renormalization group (NRG) calculations. NRG provides a benchmark for any linear response property throughout the ultrastrong-coupling regime. We find that in a sudden quantum quench, the coherent state approach produces physical artifacts, such as improper relaxation to the steady state. These previously unnoticed problems are related to the simplified form of the ansatz that generates spurious correlations within the bath. In the scattering problem, NRG is used to find the transmission and reflection of a single photon, as well as the inelastic scattering of that single photon. Simple analytical formulas are established and tested against the NRG data that predict quantitatively the transport coefficients for up to moderate environmental impedance. These formulas resolve pending issues regarding the presence of inelastic losses in the spin-boson model near absorption resonances, and could be used for comparison to experiments in Josephson waveguide quantum electrodynamics. Finally, the scattering results using the coherent state wave-function approach are compared favorably to the NRG results for very weak incident intensity. Finally, we end our study by presenting results at higher power where the response of the system is nonlinear.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Proximity-enhanced valley Zeeman splitting at the WS 2 /graphene interface

The valley Zeeman physics of excitons in monolayer transition metal dichalcogenides provides valuable insight into the spin and orbital degrees of freedom inherent to these materials. Being atomically-thin materials, these degrees of freedom can be influenced by the presence of adjacent layers, due to proximity interactions that arise from wave function overlap across the 2D interface. Here, we report 60 T magnetoreflection spectroscopy of the A- and B- excitons in monolayer WS 2 , systematically encapsulated in monolayer graphene. While the observed variations of the valley Zeeman effect for the A- exciton are qualitatively in accord with expectations from the bandgap reduction and modification of the exciton binding energy due to the graphene-induced dielectric screening, the valley Zeeman effect for the B- exciton behaves markedly different. We investigate prototypical WS 2 /graphene stacks employing first-principles calculations and find that the lower conduction band of WS 2 at the <?CDATA $K/K^{^{\prime}}$?> K / K ′ valleys (the <:tex-math><?CDATA $\mathrm{CB}^-$?> C B − band) is strongly influenced by the graphene layer on the orbital level. Specifically, our detailed microscopic analysis reveals that the conduction band at the Q point of WS 2 mediates the coupling between <?CDATA $\mathrm{CB}^-$?> C B − and graphene due to resonant energy conditions and strong coupling to the Dirac cone. This leads to variations in the valley Zeeman physics of the B- exciton, consistent with the experimental observations. Our results therefore expand the consequences of proximity effects in multilayer semiconductor stacks, showing that wave function hybridization can be a multi-step energetically resonant process, with different bands mediating the interlayer interactions. Such effects can be further exploited to resonantly engineer the spin-valley degrees of freedom in van der Waals and moiré heterostructures.

36 MATERIALS SCIENCE↗

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↗

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↗

Nucleon and nucleon-pair momentum distributions in A≤12 nuclei

We report variational Monte Carlo calculations of single-nucleon momentum distributions for A≤12 nuclei and nucleon-pair and nucleon-cluster momentum distributions for A≤8. The wave functions have been generated for a Hamiltonian containing the Argonne ν 18 two-nucleon and Urbana X three-nucleon potentials. The single-nucleon and nucleon-pair momentum distributions exhibit universal features attributable to the one-pion-exchange tensor interaction The single-nucleon distributions are broken down into proton and neutron components and spin-up and spin-down components where appropriate. The nucleon-pair momentum distributions are given separately for pp and pn pairs. The nucleon-cluster momentum distributions include dp in 3 He, tp and dd in S 4 He, αd in 6 Li,αt in 7 Li, and αα in 8 Be. Detailed tables are provided on-line for download.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Intrinsic transverse momentum and dynamical chiral symmetry breaking

We study the effect of QCD vacuum structure on the intrinsic transverse momentum distribution of partons in the nucleon at a low scale. The dynamical breaking of chiral symmetry is caused by non-perturbative interactions at distances of the order rho ~ 0.2 - 0.3 fm, much smaller than the typical nucleon size R ~ 1 fm, resulting in a two-scale picture of nucleon structure. Using an effective dynamical model based on chiral constituent quark degrees of freedom and the 1/N_c expansion (chiral quark-soliton model), we calculate the transverse momentum distribution of quarks and antiquarks at a low scale. The distribution of valence quarks is localized at p_T ~ 1/R. The distribution of flavor-singlet unpolarized sea quarks exhibits a power-like tail extending up to the chiral-symmetry-breaking scale 1/{rho}. A similar tail is present in the flavor-nonsinglet polarized sea. These features are model-independent and represent the imprint of the QCD vacuum on the nucleon&#x27;s partonic structure. At the level of the nucleon&#x27;s light-cone wave function, we show that sea quarks partly exist in correlated pairs of transverse size {rho} << R, analogous to short-range NN correlations in nuclei. We discuss the implications of our findings for the transverse momentum distributions in hard scattering processes (semi-inclusive DIS, Drell-Yan pair production) and possible experimental tests of the non-perturbative parton correlations induced by QCD vacuum structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Invariant currents and scattering off locally symmetric potential landscapes

We study the effect of discrete symmetry breaking in inhomogeneous scattering media within the framework of generic wave propagation. Our focus is on one-dimensional scattering potentials exhibiting local symmetries. We find a class of spatially invariant nonlocal currents, emerging when the corresponding generalized potential exhibits symmetries in arbitrary spatial domains. These invariants characterize the wave propagation and provide a spatial mapping of the wave function between any symmetry related domains. This generalizes the Bloch and parity theorems for broken reflection and translational symmetries, respectively. Their nonvanishing values indicate the symmetry breaking, whereas a zero value denotes the restoration of the global symmetry where the well-known forms of the two theorems are recovered. These invariants allow for a systematic treatment of systems with any local symmetry combination, providing a tool for the investigation of the scattering properties of aperiodic but locally symmetric systems. To this aim we express the transfer matrix of a locally symmetric potential unit via the corresponding invariants and derive quantities characterizing the complete scattering device which serve as key elements for the investigation of transmission spectra and particularly of perfect transmission resonances. -- Highlights: •We show that local discrete symmetries yield invariant currents. •Bloch and parity theorems are generalized when the associated symmetries are broken. •Formulation of scattering via the symmetry-induced invariant currents. •We provide sum rules for the invariant currents characterizing perfect transmission.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diversity of Chemical Bonding and Oxidation States in MS 4 Molecules of Group 8 Elements

The geometric and electronic ground-state structures of six MS 4 molecules (M = group-8 metals Fe, Ru, Os, Hs, Sm, and Pu) have been studied by using quantum-chemical density-functional and correlated wave-function approaches. The MS 4 species are compared to analogous MO 4 species recently investi-gated (Inorg. Chem. 2016, 55: 4616). Metal oxidation state (MOS) of high value VIII appears in low- spin singlet Td geometric species (Os,Hs)S 4 and (Ru,Os,Hs)O 4 , whereas low MOS=II appears in high- spin septet D 2d species Fe(S 2 ) 2 and (slightly excited) metastable Fe(O 2 ) 2 . The ground states of all other molecules have intermediate MOS values, containing S 2- , S 2 2- , S2 1- (and resp. O 2- -, O 1- , O 2 2- , O 2 1- ) ligands, bonded by ionic, covalent and correlative contributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytical solutions of a generalized non-central potential in N-dimensions

We present that N-dimensional non-relativistic wave equation for the generalized non-central potential with arbitrary angular momentum is analytically solvable in the hyperspherical coordinates. Asymptotic iteration method as a different approach is applied to obtain N-dimensional energy eigenvalues and the corresponding eigenfunctions. In hyperspherical coordinates, the wave function solutions are obtained in terms of hypergeometric functions and Jacobi polynomials. The bound states of quantum systems under consideration for some special cases, such as Hartmann and Makarov potentials, have been discussed in N-dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

On the bond distance in methane

The equilibrium bond distance in methane was optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valence correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.

Bowen-Jenkins, Philippa↗

On the bond distance in methane

The equilibrium bond distance in methane has been optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valance correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.

Bowen-Jenkins, Philippa↗

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↗