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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Extremely Slow Trap-Mediated Hole Relaxation in Room-Temperature InP/ZnSe/ZnS Quantum Dots

Time-resolved photoluminescence (PL) and transient absorption (TA) spectroscopy have been used to elucidate the hole relaxation dynamics in high-quality InP/ZnSe/ZnS quantum dots (QDs) at room temperature. In QDs having an In/P ratio close to 1:1, hole cooling is fast compared to the ~40 ps time-correlated photon-counting instrument response. However, in QDs that luminesce at wavelengths greater than around 600 nm and have excess indium, a significant fraction of the PL exhibits a slower rise time. In this case, the PL rises over a range of timescales varying from close to the instrument response to approximately 300–400 ps. The time constant of the slowest component increases with increasing InP core size and the fraction of the PL that has the slow rise time increases with excess indium in the particle. In all cases, a slow rise time and any spectral evolution on that timescale are absent in the TA results, indicating that the slow dynamics may be assigned to relaxation in the valence band. QDs emitting at wavelengths below 560 nm show very different dynamics, exhibiting a barely detectable <50 ps rise time. Based on these results and measured elemental compositions, the slow hole relaxation is assigned to the transient population of indium-based traps in the ZnSe shell that are above the valence band edge only for the larger, red-emitting QDs. Here, the relative trap to band edge relaxation rates can be understood in terms of effective mass approximation hole wave function calculations, which show that the relaxation rate is largely determined by the extent to which the hole wave function in the lowest exciton penetrates into the ZnSe shell, overlapping the trapped-hole wave functions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2021

Localized Active Space Pair-Density Functional Theory

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

On the spin separation of algebraic two-component relativistic Hamiltonians: Molecular properties

The idea for separating the algebraic exact two-component (X2C) relativistic Hamiltonians into spin-free (sf) and spin-dependent terms [Z. Li, Y. Xiao, and W. Liu, J. Chem. Phys. 137, 154114 (2012)] is extended to both electric and magnetic molecular properties. Taking the spin-free terms (which are correct to infinite order in α ≈ 1/137) as zeroth order, the spin-dependent terms can be treated to any desired order via analytic derivative technique. This is further facilitated by unified Sylvester equations for the response of the decoupling and renormalization matrices to single or multiple perturbations. For practical purposes, explicit expressions of order α{sup 2} in spin are also given for electric and magnetic properties, as well as two-electron spin-orbit couplings. At this order, the response of the decoupling and renormalization matrices is not required, such that the expressions are very compact and completely parallel to those based on the Breit-Pauli (BP) Hamiltonian. However, the former employ sf-X2C wave functions, whereas the latter can only use nonrelativistic wave functions. As the sf-X2C terms can readily be interfaced with any nonrelativistic program, the implementation of the O(α{sup 2}) spin-orbit corrections to sf-X2C properties requires only marginal revisions of the routines for evaluating the BP type of corrections.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1980

Two-photon excitation cross-section in light and intermediate atoms

The method of explicit summation over the intermediate states is used along with LS coupling to derive an expression for two-photon absorption cross section in light and intermediate atoms in terms of integrals over radial wave functions. Two selection rules, one exact and one approximate, are also derived. In evaluating the radial integrals, for low-lying levels, the Hartree-Fock wave functions, and for high-lying levels, hydrogenic wave functions obtained by the quantum defect method are used. A relationship between the cross section and the oscillator strengths is derived. Cross sections due to selected transitions in nitrogen, oxygen, and chlorine are given. The expression for the cross section is useful in calculating the two-photon absorption in light and intermediate atoms.

Omidvar, K.↗

DOE OSTI2016

Delayed Update Algorithms for Quantum Monte Carlo Simulation on GPU: Extended Abstract

QMCPACK is open source scientific software designed to perform Quantum Monte Carlo (QMC) simulation, a first-principles method for describing many-body systems. The evaluation of each Monte Carlo move requires finding the determinant of a dense matrix in the wave functions. This calculation forms a key computational kernel in QMCPACK. After each accepted event, the wave function matrix undergoes a rank-one update to represent a single particle move within the system. The matrix inverse is updated via the Sherman-Morrison formula; occasionally, the explicit inverse must be recomputed to maintain numerical stability.; We propose an alternate approach to this kernel that utilizes QR factorization to maintain stability without refactorization. In addition, algorithms based on a novel delayed update scheme are explored in this effort. This strategy enables probability evaluation for multiple successive Monte Carlo moves, with application of accepted moves to the wave function matrix delayed until an event is denied, or a predetermined limit on acceptances p is reached. Accepted events grouped in this manner are then applied to the matrix en bloc with enhanced arithmetic intensity and computational efficiency. Delayed update Sherman-Morrison probability evaluation algorithms with CPU and GPU-accelerated implementations are tested, profiled, and analyzed. A QR based delayed update algorithm remains in development. Results are evaluated against existing methods for numerical stability and efficiency; emphasis is placed on large systems, for which acceleration is critical.

Computer Science↗

DOE OSTI2026

Fokker-Planck Equation Governing the Distribution of Walkers in Auxiliary-Field Quantum Monte Carlo

Auxiliary-field quantum Monte Carlo (AFQMC) is typically formulated as an open-ended random walk in an overcomplete space of Slater determinants, implemented through a Langevin equation. However, the explicit form of the underlying Fokker-Planck equation governing the walker population distribution has remained unknown. Here, in this Letter, we derive the Fokker-Planck equation for AFQMC and propose a novel numerical scheme to solve it. The solution of the Fokker-Planck equation reveals the wave function actually sampled by the AFQMC algorithm. Interestingly, we find that even when the exact ground state is used as a guiding wave function in constrained path AFQMC, contrary to the common assumption, the wave function sampled by AFQMC is not exact. Beyond clarifying several fundamental aspects of AFQMC, the availability of a Fokker-Planck equation formulation opens new avenues for systematically improving its accuracy, which we outline in this Letter.

Monte Carlo methods↗

NASA NTRS1978

Electron impact excitation of highly charged sodium-like ions

Optical transition probabilities and electron collision strengths for Ca X, Fe XVI, Zn XX, Kr XXVI and Mo XXXII are calculated for transitions between n equal to 3 and n equal to 4 levels. The calculations neglect relativistic effects on the radial functions. A semi-empirical approach provides wave functions of the excited states; a distorted wave function without exchange is employed to obtain the excitation cross sections. The density dependence of the relative intensities of certain emission lines in the sodium isoelectronic sequence is also discussed.

Blaha, M.↗

DOE OSTI2017

Flux Rule, U(1) Instanton, and Superconductivity

We examine some consequences of the duality that a U(1) phase factor added on a wave function describes a whole system motion and also plays the role of a U(1) gauge potential. First, we show that the duality solves a long-standing puzzling problem that the ‘flux rule’ (the Faraday’s induction formula) and the Lorentz force calculation for an emf emerging in an electron system moving in a magnetic field give the same result (Feynman et al. 1963). Next, we examine a U(1) phase factor induced on the wave function for an electron system due to the single-valuedness requirement of the wave function with respect to the electron coordinates, and its consequential appearance of a U(1) instanton. This instanton explains the Meissner effect, supercurrent generation, flux quantization in the units of (h/(2e)), and the voltage quantization in the units of ((hf)/(2e)) across the Josephson junction in the presence of a radiation field with frequency f. In the experiment, a radiation field must be present to have a finite voltage across the Josephson junction; but a clear explanation for it has been lacking. The present work provides an explanation for it, and also explains the high precision of the quantized voltage as due to a topological effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2017

Excitation variance matching with limited configuration interaction expansions in variational Monte Carlo

In the regime where traditional approaches to electronic structure cannot afford to achieve accurate energy differences via exhaustive wave function flexibility, rigorous approaches to balancing different states’ accuracies become desirable. As a direct measure of a wave function’s accuracy, the energy variance offers one route to achieving such a balance. Here, we develop and test a variance matching approach for predicting excitation energies within the context of variational Monte Carlo and selective configuration interaction. In a series of tests on small but difficult molecules, we demonstrate that the approach it is effective at delivering accurate excitation energies when the wave function is far from the exhaustive flexibility limit. Results in C3, where we combine this approach with variational Monte Carlo orbital optimization, are especially encouraging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Effect of Configuration Interaction of Rare-Earth Ion States on the Intensity of Intermultiplet Transitions

We present the analysis of wave functions and the effect of configuration interaction on the intensity of absorption bands for Pr{sup 3+}, Tm{sup 3+}, and Er{sup 3+} ions. It is shown that due to the strong spin–orbit interaction the wave functions of rare-earth ions are multicomponent superpositions of ″pure″ {sup 2S+1}L{sub J} multiplets with different values of S and L and the same values of the angular momentum J. For this reason, the effect of configurational interaction depends not only on the energy gap up to excited configurations, but also on the presence of high-energy {sup 2S+1}L{sub J} components. The presence of such components yields a strong influence on the configuration interaction even at the deep layers. Another consequence of the multicomponent composition of the wave functions is the formation of groups of strongly interrelated states. When the configuration interaction is not fully taken into account, the error in describing the oscillator strengths is distributed among all the states of the group, and the low overall accuracy of the description is not perceived as a contradiction to the theory. To confirm the proposed assumptions, the oscillator strengths of the absorption transitions of laser materials activated by Pr{sup 3+}, Tm{sup 3+}, and Er{sup 3+} ions are described using the modified Judd–Ofelt theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Anti-symmetric barron functions and their approximation with sums of determinants

A fundamental problem in quantum physics is to encode functions that are completely anti-symmetric under permutations of identical particles. The architecture of neural network models for the electron wave function typically comprises an equivariant component followed by a summation of determinants. The recently introduced Generic Antisymmetric (GA) block is designed to enhance the expressivity of such neural wave functions, and it was found that the 2-layer GA block achieved more accurate energies than the corresponding single-determinant FermiNet architecure, suggesting its promise as a way to improve the expressivity of neural wave functions. In this paper we show how the function expressed by the 2-layer GA block can be decomposed into a sum of determinants. We formalize this result by defining the antisymmetric Barron space as a generalized version of the 2-layer GA block and providing an appromation theorem for this function class. This result can be viewed as a negative result showing that the 2-layer GA block is not more expressive than using multiple determinants.

Abrahamsen, Nilin↗

DOE OSTI2013

Intrinsic transverse momentum and parton correlations from dynamical chiral symmetry breaking

The dynamical breaking of chiral symmetry in QCD is caused by nonperturbative interactions on a distance scale rho ~ 0.3 fm, much smaller than the typical hadronic size R ~ 1 fm. These short-distance interactions influence the intrinsic transverse momentum distributions of partons and their correlations at a low normalization point. We study this phenomenon in an effective description of the low-energy dynamics in terms of chiral constituent quark degrees of freedom, which refers to the large-N_c limit of QCD. The nucleon is obtained as a system of constituent quarks and antiquarks moving in a self-consistent classical chiral field (relativistic mean-field approximation, or chiral quark-soliton model). The calculated transverse momentum distributions of constituent quarks and antiquarks are matched with QCD quarks, antiquarks and gluons at the chiral symmetry--breaking scale rho^{-2}. We find that the transverse momentum distribution of valence quarks is localized at p_T^2 ~ R^{-2} and roughly of Gaussian shape. The distribution of unpolarized sea quarks exhibits a would-be power-like tail ~1/p_T^2 extending up to the chiral symmetry-breaking scale. Similar behavior is observed in the flavor-nonsinglet polarized sea. The high-momentum tails are the result of short-range correlations between sea quarks in the nucleon's light-cone wave function, which are analogous to short-range NN correlations in nuclei. We show that the nucleon's light-cone wave function contains correlated pairs of transverse size rho << R with scalar-isoscalar (Sigma) and pseudoscalar-isovector (Pi) quantum numbers, whose internal wave functions have a distinctive spin structure and become identical at p_T^2 ~ rho^{-2} (restoration of chiral symmetry). These features are model-independent and represent an effect of dynamical chiral symmetry breaking on the nucleon's partonic structure. Our results have numerous implications for the transverse momentum distributions of particles produced in hard scattering processes. Under certain conditions the nonperturbative parton correlations predicted here could be observed in particle correlations between the current and target fragmentation regions of deep-inelastic scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2013

Phononic crystals of spherical particles: A tight binding approach

The vibrational dynamics of a fcc phononic crystal of spheres is studied and compared with that of a single free sphere, modelled either by a continuous homogeneous medium or by a finite cluster of atoms. For weak interaction among the spheres, the vibrational dynamics of the phononic crystal is described by shallow bands, with low degree of dispersion, corresponding to the acoustic spheroidal and torsional modes of the single sphere. The phonon displacements are therefore related to the vibrations of a sphere, as the electron wave functions in a crystal are related to the atomic wave functions in a tight binding model. Important dispersion is found for the two lowest phonon bands, which correspond to zero frequency free translation and rotation of a free sphere. Brillouin scattering spectra are calculated at some values of the exchanged wavevectors of the light, and compared with those of a single sphere. With weak interaction between particles, given the high acoustic impedance mismatch in dry systems, the density of phonon states consist of sharp bands separated by large gaps, which can be well accounted for by a single particle model. Based on the width of the frequency gaps, tunable with the particle size, and on the small number of dispersive acoustic phonons, such systems may provide excellent materials for application as sound or heat filters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Machine learning the dynamics of quantum kicked rotor

Highlights: • Quantum kicked rotor. • Quantum phase transition. • Machine learning. • LSTM network. • Convolutional neural network. Using the multilayer convolutional neural network (CNN), we can detect the quantum phases in random electron systems, and phase diagrams of two and higher dimensional Anderson transitions and quantum percolations as well as disordered topological systems have been obtained. Here, instead of using CNN to analyze the wave functions, we analyze the dynamics of wave packets via long short-term memory network (LSTM). We adopt the quasi-periodic quantum kicked rotors, which simulate the three and four dimensional Anderson transitions. By supervised training, we let LSTM extract the features of the time series of wave packet displacements in localized and delocalized phases. We then simulate the wave packets in unknown phases and let LSTM classify the time series to localized and delocalized phases. We compare the phase diagrams obtained by LSTM and those obtained by CNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Interference in nonlinear Compton scattering using a Schrödinger-equation approach

The interference between Compton scattering and nonlinear Compton scattering from a two-color field in the X-ray regime is theoretically examined for bound electrons. Here, the underlying phase shifts are analyzed using a perturbative approach in the incoming classical field. The perturbative approach is bench marked with a non-perturbative approach in the classical field. The interference for different combinations of linear polarization of the two fields is examined when the Compton and the nonlinear Compton scattered waves have the same wave vector and polarization. Only two cases exhibit interference. When there is interference, the calculations reveal an intrinsic phase difference between the Compton scattered wave function and the nonlinear Compton scattered wave function of either 0 or pi depending on the scattering angle.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2016

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/t from 2 to 8, for a range of densities including half-filling and n = 0.3,0.5,0.6, 0.75, and 0.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in this paper in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to 16 × 16 are also provided for benchmark purposes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1974

Radiative transitions involving the (2p2)(3 Pe) metastable autodetaching of H(-)

The absorption coefficient for the free-bound transition H (ls) + e(-)+ h omega yields H(-)(2 sq p,(3)P(e)) is calculated (together with the differential emission rate for the inverse process) using ls - 2s - 2p close coupling continuum wave functions and a Hylleraas bound state wave function. A maximum in the absorption and emission spectra is found to occur at a photon wavelength of 1219.5 A, which is 2 A closer to the Lyman alpha line than predicted by the calculations of Drake, and is in closer agreement with the stellar absorption feature identified by Heap and Stecher. The free-bound absorption process appears to be a significant source of continuous ultraviolet opacity.

Jacobs, V. L.↗

NASA NTRS1994

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗