Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wave Functions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Calculation of the Full Scattering Amplitude without Partial Wave Decomposition II

As is well known, the full scattering amplitude can be expressed as an integral involving the complete scattering wave function. We have shown that the integral can be simplified and used in a practical way. Initial application to electron-hydrogen scattering without exchange was highly successful. The Schrodinger equation (SE) can be reduced to a 2d partial differential equation (pde), and was solved using the finite element method. We have now included exchange by solving the resultant SE, in the static exchange approximation. The resultant equation can be reduced to a pair of coupled pde's, to which the finite element method can still be applied. The resultant scattering amplitudes, both singlet and triplet, as a function of angle can be calculated for various energies. The results are in excellent agreement with converged partial wave results.

Shertzer, J.↗
DOE OSTI2020

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born–Oppenheimer approximation

In this work, we develop a formalism to accurately account for the renormalization of the electronic structure due to quantum and thermal nuclear motions within the Born–Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from quantum Monte Carlo calculations in either the canonical or grand-canonical ensembles. The formalism applies as well to effective single electron theories such as those based on density functional theory. We show that the electronic (Bloch) crystal momentum can be restored by marginalizing the total electron–ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born–Oppenheimer approximation. Based on the Kubo–Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and, in general, differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200 K and 250 GPa and discuss the optical absorption profile of hydrogen crystals at 200 K and carbon diamond at 297 K.

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

Eikonal spin-dependent odderon and gluon Sivers function of a proton and its small-𝑥 evolution

The matrix element in the proton of the eikonal Odderon operator, with a helicity flip, has been shown to correspond to the dipole gluon Sivers function. We employ a three quark light-front model of the proton to determine the Sivers function at moderately small 𝑥 0 ∼ 0.1 and transverse momentum 𝑘 ⊥ ≲ 1 GeV. The model light-cone (LC) wave function predicts the properties of 𝑥⁢𝑓$^{⊥⁢𝑔}_{1⁢𝑇⁡}$(𝑥,𝑘 ⊥ ) such as its overall magnitude, the position of its peak in 𝑘 ⊥ , and its behavior at small 𝑘 ⊥ . We then compute numerically the Balitsky-Fadin-Kuraev-Lipatov anomalous dimension characterizing the power-law tail at 𝑘 ⊥ ≳ 1.5 GeV of the gluon Sivers function at small (but preasymptotic) LC momentum fractions 𝑥 ∼ 𝑥 0 ⁢𝑒 −1/𝛼 𝑠 : 𝑥⁢𝑓$^{⊥⁢𝑔}_{1⁢𝑇⁡}$⁡(𝑥,𝑘 ⊥ ) ∼ 𝑘$^{−3.3}_{⊥}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

New Studies of Short-Range Correlations and Their Effects on Nuclei

Over the last decade, data gathered on short range correlations (SRC) and the EMC effect have stimulated tremendous progress in our understanding of the nuclear force. In particular, the tensor force was found to be the major factor in the formation of N-N short range correlations. More recently a connection between the N-N SRC probabilities and EMC effect amplitudes was first observed. However, the real origin and characteristics of this connection have not been convincingly determined experimentally or theoretically. The project consists of several approved experiments scheduled to run at Jefferson Lab in Newport News, VA. These inclusive scattering measurements will provide a detailed map of the nature of short range correlations and its role in the in-medium modification of the nucleons. The isospin dependence in 2 N- and 3 N-SRCs will be accessed through the use of two light mirror nuclei 3 He and 3 H. An extensive A-dependence study of 2 N and 3 N-SRCs will be performed as well as a measurement of the quark distribution in short range correlations, i.e. the EMC eect in SRCs. A measurement of the deuteron tensor structure function in the deep inelastic scattering regime will verify the possible contribution of hidden color components in the deuteron wave function. This project aims to help answer the fundamental question on the nature of the short range interactions and their implications at the nucleonic and partonic levels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1998

Phase Inversion: Inferring Solar Subphotospheric Flow and Other Asphericity from the Distortion of Acoustic Waves

A method is proposed for analyzing an almost monochromatic train of waves propagating in a single direction in an inhomogeneous medium that is not otherwise changing in time. An effective phase is defined in terms of the Hilbert transform of the wave function, which is related, via the JWKB approximation, to the spatial variation of the background state against which the wave is propagating. The contaminating effect of interference between the truly monochromatic components of the train is eliminated using its propagation properties. Measurement errors, provided they are uncorrelated, are manifest as rapidly varying noise; although that noise can dominate the raw phase-processed signal, it can largely be removed by low-pass filtering. The intended purpose of the analysis is to determine the distortion of solar oscillations induced by horizontal structural variation and material flow. It should be possible to apply the method directly to sectoral modes. The horizontal phase distortion provides a measure of longitudinally averaged properties of the Sun in the vicinity of the equator, averaged also in radius down to the depth to which the modes penetrate. By combining such averages from different modes, the two-dimensional variation can be inferred by standard inversion techniques. After taking due account of horizontal refraction, it should be possible to apply the technique also to locally sectoral modes that propagate obliquely to the equator and thereby build a network of lateral averages at each radius, from which the full three-dimensional structure of the Sun can, in principle, be determined as an inverse Radon transform.

Gough, Douglas↗
DOE OSTI2019

Absolutely Localized Projection-Based Embedding for Excited States

We present a quantum embedding method that allows for calculation of local excited states embedded in a Kohn–Sham density functional theory (DFT) environment. Projection-based quantum embedding methodologies provide a rigorous framework for performing DFT-in-DFT and wave function in DFT (WF-in-DFT) calculations. The use of absolute localization, where the density of each subsystem is expanded in only the basis functions associated with the atoms of that subsystem, provide improved computationally efficiency for WF-in-DFT calculations by reducing the number of orbitals in the WF calculation. In this work, we extend absolutely localized projection-based quantum embedding to study localized excited states using EOM-CCSD-in-DFT and TDDFT-in-DFT. The embedding results are highly accurate compared to the corresponding canonical EOM-CCSD and TDDFT results on the full system, with TDDFT-in-DFT frequently more accurate than canonical TDDFT. The absolute localization method is shown to eliminate the spurious low-lying excitation energies for charge-transfer states and prevent overdelocalization of excited states. Additionally, we attempt to recover the environment response caused by the electronic excitations in the high-level subsystem using different schemes and compare their accuracy. Finally, we apply this method to the calculation of the excited-state energy of green fluorescent protein and show that we systematically converge to the full system results. Here we demonstrate how this method can be useful in understanding excited states, specifically which chemical moieties polarize to the excitation. Furthermore, this work shows absolutely localized projection-based quantum embedding can treat local electronic excitations accurately and make computationally expensive WF methods applicable to systems beyond current computational limits.

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

Ab Initio Unimolecular Reaction Kinetics of CH 2 C(＝)OCH 3 and CH 3 C(＝O)OCH 2 Radicals

The unimolecular dissociation and isomerization kinetics of the methyl ester combustion intermediates methoxycarbonylmethyl (CH 2 C(＝O)OCH 3 ) and acetyloxylmethyl (CH 3 C(＝O)OCH 2 ) are theoretically investigated using high-level ab initio methods and the Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory. Potential energy surfaces (PESs) are obtained using coupled cluster singles and doubles with perturbative triples correction (CCSD(T)), multireference singles and doubles configuration interaction (MRSDCI) with the Davidson–Silver (DS) correction, and multireference averaged coupled pair functional (MRACPF2) theory. The transition states exhibit high T1 diagnostics in coupled cluster calculations, suggesting the need for a multireference correlated wave function treatment. MRSDCI+DS and MRACPF2 capture their multiconfigurational character well, yielding lower barrier heights than CCSD(T) for these reactions. Additionally, the rate coefficients are computed using the RRKM/ME theory over a 500–2500 K temperature range and at a pressure range of 0.01 atm to the high-pressure limit. The temperature- and pressure-dependent rate coefficients are given in modified Arrhenius expressions. The β-scission of CH 2 C(＝O)OCH 3 is predicted to have a much higher barrier than the corresponding isomerization reaction and the β-scission of CH 3 C(＝O)OCH 2 . Consequently, the rate coefficients for β-scission of CH 2 C(＝O)OCH 3 are the smallest among the three reactions and the isomerization followed by decomposition to CH 3 C(＝O) and HCHO is the dominant reaction pathway for CH 2 C(＝O)OCH 3 . Furthermore, both radicals CH 2 C(＝O)OCH 3 and CH 3 C(＝O)OCH 2 are predicted to mainly decompose to CH 3 C(＝O) + HCHO rather than to the bimolecular product CH 2 C(＝O) + CH 3 O. A newly developed MA combustion mechanism, using our theoretical rate coefficients for the MA-related reactions, predicts combustion properties in good agreement with available experimental data.

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

A well-scaling natural orbital theory

Here, we introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree–Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the oneparticle density matrix.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Collins-Soper kernel and reduced soft function in lattice QCD

We evaluate the Collins-Soper kernel and the reduced soft function in lattice QCD, incorporating 𝒪⁡(𝛼 𝑠 ) matching corrections. The calculation relies on the evaluation of the quasitransverse momentum–dependent wave function with asymmetric staple-shaped quark bilinear operators and four-point meson form factors. These quantities are computed nonperturbatively using two 𝑁 𝑓 =2 + 1 + 1 twisted-mass fermion ensembles with the same lattice spacing of 𝑎 = 0.093 fm: the first ensemble has a lattice size of 24 3 × 48 and a pion mass of 346 MeV, and the second one has a lattice size of 32 3 × 64 and a pion mass of 261 MeV. The Collins-Soper kernel and the soft function are needed for the determination of the transverse momentum–dependent parton distribution functions.

Lattice QCD↗
NASA NTRS1972

Hydrogen atom in intense magnetic field.

The structure of a hydrogen atom situated in an intense magnetic field is investigaged. Three approaches are employed. An elementary Bohr picture establishes a crucial magnetic field strength, H sub a approximately equal to 5 x 10 to the 9th G. Fields in excess of H sub a are intense in that they are able to modify the characteristic atomic scales of length and binding energy. A second approach solves the Schrodinger equation by a combination of variational methods and perturbation theory. It yields analytic expressions for the wave functions and energy eigenvalues. A third approach determines the energy eigenvalues by reducing the Schrodinger equation to a one-dimensional wave equation, which is then solved numerically. Energy eigenvalues are tabulated for field strengths of 2 x 10 to the 10th G and 2 x 10 to the 12th G. It is found that at 2 x 10 to the 12th G the lowest energy eigenvalue is changed from -13.6 to about -180 eV in agreement with previous variational computations.

Canuto, V.↗
NASA NTRS2007

The Photodetachment of Ps ion and Low-Energy e(+) -H Collisions

Two calculations in the area of positron collisions are presented. The first is the calculation of the photodetachment cross section of the positronium negative ion (Ps-) using accurate variational wave functions for both the initial bound-state and the final P continuum state. The second is the calculation of partial wave cross sections for Ps(1s)-formation in ef -H(ls) collisions using the hyperspherical hidden crossing method. Since the S-wave Stiickelberg phase is close to pi, the very small S-wave Ps(1s) formation cross section can be understood in terms of destructive interference. Other examples in positron collisions are given where it is either known or expected that destructive interference is the cause of the small S-wave Ps(1s) formation cross section. In addition, examples are presented of processes in atomic physics where the Stiickelberg phase is a multiple of pi/2.

Ward, S.J.↗
DOE OSTI2022

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗
DOE OSTI2014

REVISED BIG BANG NUCLEOSYNTHESIS WITH LONG-LIVED, NEGATIVELY CHARGED MASSIVE PARTICLES: UPDATED RECOMBINATION RATES, PRIMORDIAL {sup 9}Be NUCLEOSYNTHESIS, AND IMPACT OF NEW {sup 6}Li LIMITS

We extensively reanalyze the effects of a long-lived, negatively charged massive particle, X {sup –}, on big bang nucleosynthesis (BBN). The BBN model with an X {sup –} particle was originally motivated by the discrepancy between the {sup 6,} {sup 7}Li abundances predicted in the standard BBN model and those inferred from observations of metal-poor stars. In this model, {sup 7}Be is destroyed via the recombination with an X {sup –} particle followed by radiative proton capture. We calculate precise rates for the radiative recombinations of {sup 7}Be, {sup 7}Li, {sup 9}Be, and {sup 4}He with X {sup –}. In nonresonant rates, we take into account respective partial waves of scattering states and respective bound states. The finite sizes of nuclear charge distributions cause deviations in wave functions from those of point-charge nuclei. For a heavy X {sup –} mass, m{sub X} ≳ 100 GeV, the d-wave → 2P transition is most important for {sup 7}Li and {sup 7,} {sup 9}Be, unlike recombination with electrons. Our new nonresonant rate of the {sup 7}Be recombination for m{sub X} = 1000 GeV is more than six times larger than the existing rate. Moreover, we suggest a new important reaction for {sup 9}Be production: the recombination of {sup 7}Li and X {sup –} followed by deuteron capture. We derive binding energies of X nuclei along with reaction rates and Q values. We then calculate BBN and find that the amount of {sup 7}Be destruction depends significantly on the charge distribution of {sup 7}Be. Finally, updated constraints on the initial abundance and the lifetime of the X {sup –} are derived in the context of revised upper limits to the primordial {sup 6}Li abundance. Parameter regions for the solution to the {sup 7}Li problem and the primordial {sup 9}Be abundances are revised.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗
NASA NTRS1979

Siple station experiments on wave-particle interactions in the magnetosphere

Natural and controlled whistler-mode signals have been used to study nonlinear mechanisms of wave growth and wave-wave interactions (WWI) in the magnetosphere; three general classes of WWI (triggering, suppression, and entrainment) are identified and interpreted in terms of a model based on cyclotron resonance interaction. This model is also used to estimate the wave field intensity associated with different types of WWI. A new type of triggered emission, the band-limited impulse (BLI) is interpreted in terms of the switching of phase-bunched currents. In addition, an experiment to find a threshold for the excitation of the coherent wave instability is discussed, and observed VLF wave-induced transient bursts of X-rays, light and E-region ionization enhancements are considered with regard to a step function wave interacting with all resonant electrons in a given energy range.

Helliwell, R. A.↗
DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2013

Comparison of Cluster, Slab, and Analytic Potential Models for the Dimethyl Methylphosphonate (DMMP)/TiO2 (110) Intermolecular Interaction

In a previous study (J. Phys. Chem. C 2011, 115, 12403) cluster models for the TiO2 rutile (110) surface and MP2 calculations were used to develop an analytic potential energy function for dimethyl methylphosphonate (DMMP) interacting with this surface. In the work presented here, this analytic potential and MP2 cluster models are compared with DFT "slab" calculations for DMMP interacting with the TiO2 (110) surface and with DFT cluster models for the TiO2 (110) surface. The DFT slab calculations were performed with the PW91 and PBE functionals. The analytic potential gives DMMP/ TiO2 (110) potential energy curves in excellent agreement with those obtained from the slab calculations. The cluster models for the TiO2 (110) surface, used for the MP2 calculations, were extended to DFT calculations with the B3LYP, PW91, and PBE functional. These DFT calculations do not give DMMP/TiO2 (110) interaction energies which agree with those from the DFT slab calculations. Analyses of the wave functions for these cluster models show that they do not accurately represent the HOMO and LUMO for the surface, which should be 2p and 3d orbitals, respectively, and the models also do not give an accurate band gap. The MP2 cluster models do not accurately represent the LUMO and that they give accurate DMMP/TiO2 (110) interaction energies is apparently fortuitous, arising from their highly inaccurate band gaps. Accurate cluster models, consisting of 7, 10, and 15 Ti-atoms and which have the correct HOMO and LUMO properties, are proposed. The work presented here illustrates the care that must be taken in "constructing" cluster models which accurately model surfaces.

calculations↗