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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Homogeneous quantum electrodynamic turbulence

The electromagnetic field equations and Dirac equations for oppositely charged wave functions are numerically time-integrated using a spatial Fourier method. The numerical approach used, a spectral transform technique, is based on a continuum representation of physical space. The coupled classical field equations contain a dimensionless parameter which sets the strength of the nonlinear interaction (as the parameter increases, interaction volume decreases). For a parameter value of unity, highly nonlinear behavior in the time-evolution of an individual wave function, analogous to ideal fluid turbulence, is observed. In the truncated Fourier representation which is numerically implemented here, the quantum turbulence is homogeneous but anisotropic and manifests itself in the nonlinear evolution of equilibrium modal spatial spectra for the probability density of each particle and also for the electromagnetic energy density. The results show that nonlinearly interacting fermionic wave functions quickly approach a multi-mode, dynamic equilibrium state, and that this state can be determined by numerical means.

Shebalin, John V.↗
DOE OSTI2023

Constraint of pionless EFT using two-nucleon spectra from lattice QCD

Finite-volume pionless effective field theory (FVEFT π/ ) at next-to-leading order (NLO) is used to analyze the two-nucleon lattice QCD spectrum of Ref. [1], performed at quark masses corresponding to a pion mass of approximately 800 MeV. Specifically, the effective theory is formulated in finite volume, and variational sets of wave functions are optimized using differential programming. Using these wave functions projected to the appropriate finite-volume symmetry group, variational bounds from FVEFT π/ are obtained for the ground state, as well as excited states. By comparison with the lattice QCD GEVP spectrum, different low energy constants (LECs) are constrained. Furthermore, relativistic corrections are incorporated, allowing for the s-d-wave mixing term in the deuteron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

Development of New Open-Shell Perturbation and Coupled-Cluster Theories Based on Symmetric Spin Orbitals

A new spin orbital basis is employed in the development of efficient open-shell coupled-cluster and perturbation theories that are based on a restricted Hartree-Fock (RHF) reference function. The spin orbital basis differs from the standard one in the spin functions that are associated with the singly occupied spatial orbital. The occupied orbital (in the spin orbital basis) is assigned the delta(+) = 1/square root of 2(alpha+Beta) spin function while the unoccupied orbital is assigned the delta(-) = 1/square root of 2(alpha-Beta) spin function. The doubly occupied and unoccupied orbitals (in the reference function) are assigned the standard alpha and Beta spin functions. The coupled-cluster and perturbation theory wave functions based on this set of "symmetric spin orbitals" exhibit much more symmetry than those based on the standard spin orbital basis. This, together with interacting space arguments, leads to a dramatic reduction in the computational cost for both coupled-cluster and perturbation theory. Additionally, perturbation theory based on "symmetric spin orbitals" obeys Brillouin's theorem provided that spin and spatial excitations are both considered. Other properties of the coupled-cluster and perturbation theory wave functions and models will be discussed.

Lee, Timothy J.↗
NASA NTRS1988

Electron-impact excitation of atomic oxygen

A close-coupling approximation is employed to theoretically investigate the excitation of the 3s3S0 and 3p3P states of atomic oxygen for electron impact energies ranging from 13.87-100 eV. The target states included in the expansion of the total wave function are modeled by accurate configuration-interaction wave functions. The results are consistent with recent experimental findings for the dipole-allowed 3P-3s3S0 transition.

Tayal, S. S.↗
NASA NTRS1990

Spin bipolaron in the framework of emery model for high-T(sub c) copper oxide superconductors

The high-T(sub c) oxide compounds discovered recently exhibit a number of interesting physical properties. Two-dimensional antiferromagnetic spin order has been observed in these materials at the oxygen deficiency. This fact can be explained by strong correlation of the spins, situated on Cu sites in the conducting planes of the oxide superconductors. The doping or the oxygen deficiency lead to the occurrence of holes, occupying the oxygen p-orbitals according to the Emery model. At the small hole concentration they can move along the antiferromagnetic lattice of spins, localized on Cu sites. Researchers consider the two holes situation and describe in what way their behavior depends on the antiferromagnetic exchange interation J. It is known that in the framework of Hubbard model with strong on-site Coulomb repulsion, a single hole can form a spin polaron of the large radius. It is reasonable to admit that two holes with parallel spins (triplet) form the spin bipolaron complex owing to the hole excitations' capability to polarize Cu spin surroundings. Such an excitation was considered in the phenomenological way. Here the problem is discussed on the basis of the microscopic approach in the framework of the variational principle. A special kind of wave function is used for such a purpose. The wave function is constructed by generalizing the trial functions proposed in over two holes excitation situation (triplet) and then the region of spin bipolaron existance in the framework of Emery model is studied. In this model the Hamiltonian can be easily rewritten by forming the oxygen states transforming as the irreducible representations of the group D(sub 4).

Golub, A. A.↗
DOE OSTI2022

Quasi-deuteron model at low renormalization group resolution

The quasi-deuteron model introduced by Levinger is used to explain cross sections for knocking out high-momentum protons in photoabsorption on nuclei. This is within a framework we characterize as exhibiting high renormalization group (RG) resolution. Assuming a one-body reaction operator, the nuclear wave function must include two-body short-range correlations (SRCs) with deuteronlike quantum numbers. In Phys. Rev. C 104, 034311 (2021), we showed that SRC physics can be naturally accounted for at low RG resolution. We describe the quasi-deuteron model at low RG resolution and determine the Levinger constant, which is proportional to the ratio of nuclear photoabsorption to that for photodisintegration of a deuteron. We extract the Levinger constant based on the ratio of momentum distributions at high relative momentum. We compute momentum distributions evolved under similarity RG (SRG) transformations where the SRC physics is shifted into the operator as a universal two-body term. The short-range nature of this operator motivates using local-density approximations with uncorrelated wave functions in evaluating nuclear matrix elements, which greatly simplifies the analysis. The operator must be consistently matched to the RG scale and scheme of the interaction for a reliable extraction. We apply SRG transformations to different nucleon-nucleon (NN) interactions and use the deuteron wave functions and Weinberg eigenvalues to determine approximate matching scales. We predict the Levinger constant for several NN interactions and a wide range of nuclei comparing to experimental extractions. The predictions at low RG resolution are in good agreement with experiment when starting with a hard NN interaction and the initial operator. Similar agreement is found using soft NN interactions when the additional two-body operator induced by evolution from hard to soft is included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

GronOR: Massively Parallel and GPU-Accelerated Non-Orthogonal Configuration Interaction for Large Molecular Systems

GronOR is a program package for non-orthogonal configuration interaction calculations for an electronic wave function built in terms of anti-symmetrized products of multi-configuration molecular fragment wave functions. The two-electron integrals that have to be processed may be expressed in terms of atomic orbitals or in terms of an orbital basis determined from the molecular orbitals of the fragments. The code has been specifically designed for execution on distributed memory massively parallel and Graphics Processing Unit (GPU)-accelerated computer architectures, using an MPI+OpenACC/OpenMP programming approach. The task-based execution model used in the implementation allows for linear scaling with the number of nodes on the largest pre-exascale architectures available, provides hardware fault resiliency, and enables effective execution on systems with distinct central processing unit-only and GPU-accelerated partitions. The code interfaces with existing multi-configuration electronic structure codes that provide optimized molecular fragment orbitals, configuration interaction coefficients, and the required integrals. Algorithm and implementation details, parallel and accelerated performance benchmarks, and an analysis of the sensitivity of the accuracy of results and computational performance to thresholds used in the calculations are presented.

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

The quantum Hall effect at a weak magnetic field

Using a weak limit for the hopping integral in one direction in the Hofstadter model, we show that the fermion states in the gaps of the spectrum are determined within the Kitaev chain. The proposed approach allows us to study the behavior of Chern insulators (CI) in different classes of symmetry. We consider the Hofstadter model on the square and honeycomb lattices in the case of rational and irrational magnetic fluxes ϕ, and discuss the behavior of the Hall conductance at a weak magnetic field in a sample of finite size. We show that in the semiclassical limit at the center of the fermion spectrum, the Bloch states of fermions turn into chiral Majorana fermion liquid when the magnetic scale (1/ϕ) is equal to the sample size N. We are talking about the dielectric–metal phase transition, which is determined by the behavior of the Landau levels in 2D fermion systems in a transverse magnetic field. When a magnetic scale, which determines the wave function of fermions, exceeds the size of the sample, a jump in the longitudinal conductance occurs. The wave function describes non-localized states of fermions, the sample becomes a conductor, the system changes from the dielectric state to the metallic one. It is shown, that at 1∕ϕ>N the quantum Hall effect and the Landau levels are not realized, which makes possibility to study the behavior of CI in irrational magnetic fluxes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

The overlapping fragment approach for non-orthogonal configuration interaction with fragments

The non-orthogonal configuration interaction with fragments (NOCI-F) approach is extended opening the possibility to study intramolecular processes and materials with covalent or ionic lattices. So far, NOCI-F has been applied to study intermolecular energy and electron transfer employing ensembles of fragments that do not have atoms or bonds in common. The here presented approach divides the target system into two overlapping fragments that share one or more atoms and/or one or more bonds. After the construction of a collection of (multiconfigurational) fragment wave functions in a state specific optimization procedure, the fragment wave functions are combined to form many-electron basis functions for the non-orthogonal configuration interaction of the whole system. The orbitals in the overlapping fragment are defined by a corresponding orbital transformation of the fragment orbitals through a singular value decomposition. The overlapping fragments approach is first illustrated for a model system and then used to highlight some possible applications of NOCI with overlapping fragments. In conclusion, the results of excited state diffusion in transition metal oxide, intramolecular singlet fission and magnetic interactions in organic biradicals and ionic compounds are discussed.

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

Superluminal matter waves

The Dirac equation has resided among the greatest successes of modern physics since its emergence as the first quantum mechanical theory fully compatible with special relativity. This compatibility ensures that the expectation value of the velocity is less than the vacuum speed of light. Here, we show that the Dirac equation admits free-particle solutions where the peak amplitude of the wave function can travel at any velocity, including those exceeding the vacuum speed of light, despite having a subluminal velocity expectation value. The solutions are constructed by superposing basis functions with correlations in momentum space. These arbitrary velocity wave functions feature a near-constant profile and may impact quantum mechanical processes that are sensitive to the local value of the probability density as opposed to expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Photoionization cross-sections for atoms and ions of aluminum, silicon, and argon.

Photoionization cross sections for all levels belonging to the configurations of atoms and ions of aluminum, silicon, and argon have been calculated using Hartree-Fock bound-electron wave functions and close-coupling approximation free-electron wave functions. The results are presented in the form of a computationally convenient interpolation formula and should find wide astrophysical application.

Chapman, R. D.↗
NASA NTRS1977

Line strengths, collision strengths and excitation rates for multiply-charged silicon ions

In the present paper, the line strengths, collision strengths, and rate coefficients are calculated for a variety of transitions in multiply charged silicon ions from Si(VI) to Si(XIV). The line strengths are obtained by using Clementi wave functions for the ground-state configuration, and excited-state wave functions generated by a semiempirical method. The collision strengths are calculated in an LS coupling scheme in the distorted-wave approximation, neglecting exchange except for the helium-like transitions. These results are then integrated over a Maxwellian velocity distribution function to yield rate coefficients. The rates are presented graphically and also in terms of a two-parameter fit.

Davis, J.↗
DOE OSTI2025

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
NASA NTRS1973

Photoionization from excited states of helium

The cross sections for photoionization from the 2 1S, 2 3S, 2 1P and 2 3P excited states of helium are calculated for photoelectron energies below the n = 2 threshold of He(+) using Hylleraas bound state wave functions and 1s-2s-2p close coupling final state wave functions. The resonant structures associated with the lowest-lying 1S, 1P, 3P, and 1D autoionizing states of helium are found to be characterized by large values of the line profile parameter q. The cross sections and the photoelectron angular distribution asymmetry parameters for the P-states are calculated for various polarization states of the target atom and the incident photon. Experiments which would lead to the separate determinations of the S- and D- wave partial photoionization cross sections are discussed.

Jacobs, V. L.↗
NASA NTRS1995

Studies of radar backscatter as a function of wave properties and the winds in the turbulent marine atmosphere

The research on model functions for ADEOS and ERS-1 are summarized and an analysis of the differences between the three kinds of models is provided in this final report. The success of the AMI on ERS-1 obtained at GSFC and NMC is highlighted. The problem of wind stress description is reviewed within and the scatterometer model being developed for high winds monitoring for the AMI on ERS-1 and ERS-2 is described.

Pierson, Willard J., Jr.↗
DOE OSTI2021

The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions

Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density ap- proximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vor- tices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.

Jin, Shi↗
DOE OSTI2025

SchrödingerNet: A Universal Neural Network Solver for the Schrödinger Equation

Recent advances in machine learning have facilitated numerically accurate solution of the electronic Schrödinger equation (SE) by integrating various neural network (NN)-based wave function ansatzes with variational Monte Carlo methods. Nevertheless, such NN-based methods are all based on the Born–Oppenheimer approximation (BOA) and require computationally expensive training for each nuclear configuration. In this work, we propose a novel NN architecture, SchrödingerNet, to solve the full electronic-nuclear SE by defining a loss function designed to equalize local energies across the system. This approach is based on a translationally, rotationally and permutationally symmetry-adapted total wave function ansatz that includes both nuclear and electronic coordinates. Furthermore, this strategy not only allows for an efficient and accurate generation of a continuous potential energy surface at any geometry within the well-sampled nuclear configuration space, but also incorporates non-BOA corrections, through a single training process. Comparison with benchmarks of atomic and small molecular systems demonstrates its accuracy and efficiency.

Chemical calculations↗
DOE OSTI2023

A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo

Abstract We present a discussion of recent progress in excited‐state‐specific quantum chemistry and quantum Monte Carlo alongside a demonstration of how a combination of methods from these two fields can offer reliably accurate excited state predictions across singly excited, doubly excited, and charge transfer states. Both of these fields have seen important advances supporting excited state simulation in recent years, including the introduction of more effective excited‐state‐specific optimization methods, improved handling of complicated wave function forms, and ways of explicitly balancing the quality of wave functions for ground and excited states. To emphasize the promise that exists at this intersection, we provide demonstrations using a combination of excited‐state‐specific complete active space self‐consistent field theory, selected configuration interaction, and state‐specific variance minimization. These demonstrations show that combining excited‐state‐specific quantum chemistry and variational Monte Carlo can be more reliably accurate than either equation of motion coupled cluster theory or multi‐reference perturbation theory, and that it can offer new clarity in cases where existing high‐level methods do not agree. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry

Otis, Leon↗