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 91 records · Page 5

Non-orthogonal determinants in multi-Slater-Jastrow trial wave functions for fixed-node diffusion Monte Carlo

The accuracy and efficiency of ab initio Quantum Monte Carlo (QMC) algorithms benefit greatly from compact variational trial wave functions that accurately reproduce ground state properties of a system. We investigate the possibility of using multi-Slater-Jastrow trial wave functions with non-orthogonal determinants by optimizing identical single particle orbitals independently in separate determinants. As a test case, we compute variational and fixed-node diffusion Monte Carlo (FN-DMC) energies of a C 2 molecule. For a given multi-determinant expansion, we find that this non-orthogonal orbital optimization results in a consistent improvement in the variational energy and the FN-DMC energy on the order of a few tenths of an eV. In some cases, fewer non-orthogonal determinants are required compared to orthogonal ones in order to achieve similar accuracy in FN-DMC. Our calculations indicate that trial wave functions with non-orthogonal determinants can improve computed energies in a QMC calculation when compared to their orthogonal counterparts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Emergence of complex and spinor wave functions in scale relativity. I. Nature of scale variables

One of the main results of scale relativity as regards the foundation of quantum mechanics is its explanation of the origin of the complex nature of the wave function. The scale relativity theory introduces an explicit dependence of physical quantities on scale variables, founding itself on the theorem according to which a continuous and non-differentiable space-time is fractal (i.e., scale-divergent). In the present paper, the nature of the scale variables and their relations to resolutions and differential elements are specified in the non-relativistic case (fractal space). We show that, owing to the scale-dependence which it induces, non-differentiability involves a fundamental two-valuedness of the mean derivatives. Since, in the scale relativity framework, the wave function is a manifestation of the velocity field of fractal space-time geodesics, the two-valuedness of velocities leads to write them in terms of complex numbers, and yields therefore the complex nature of the wave function, from which the usual expression of the Schrödinger equation can be derived.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A light weight regularization for wave function parameter gradients in quantum Monte Carlo

The parameter derivative of the expectation value of the energy, ∂E/∂p, is a key ingredient in variational Monte Carlo (VMC) wave function optimization methods. In some cases, a naïve estimate of this derivative suffers from an infinite variance, which inhibits the efficiency of optimization methods that rely on a stable estimate of the derivative. In this work, we derive a simple regularization of the naïve estimator, which is trivial to implement in existing VMC codes, has finite variance, and a negligible bias, which can be extrapolated to zero bias with no extra cost. We use this estimator to construct an unbiased, finite variance estimation of ∂E/∂p for a multi-Slater–Jastrow trial wave function on the LiH molecule and in the optimization of a multi-Slater–Jastrow trial wave function on the CuO molecule. This regularized estimator is a simple and efficient estimator of ∂E/∂p for VMC optimization techniques.

36 MATERIALS SCIENCE↗

Interpolated Wave Functions for Nonadiabatic Simulations with the Fixed-Node Quantum Monte Carlo Method. In: Recent Progress in Quantum Monte Carlo, Chapter 3

Simulating nonadiabatic effects with many-body wave function approaches is an open field with many challenges. Recent interest has been driven by new algorithmic developments and improved theoretical understanding of properties unique to electron-ion wave functions. Fixed-node diffusion Monte Caro is one technique that has shown promising results for simulating electron-ion systems. In particular, we focus on the CH molecule for which previous results suggested a relatively significant contribution to the energy from nonadiabatic effects. We propose a new wave function ansatz for diatomic systems which involves interpolating the determinant coefficients calculated from configuration interaction methods. The calculated nonadiabatic contribution to the energy in the CH molecule is reduced compared to our previous results, but still remains the largest among the molecules under consideration.

Tubman, Norm M.↗

Cluster decomposition of full configuration interaction wave functions: A tool for chemical interpretation of systems with strong correlation

Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size, thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation amplitudes contain information on the complexity of the electronic wave function, but this information is contaminated by contributions from disconnected excitations, i.e., those excitations that are just products of independent lower-level excitations. The unwanted contributions can be removed via a cluster decomposition procedure, making it possible to examine the importance of connected excitations in complicated multireference molecules which are outside the reach of conventional algorithms. We present an implementation of the cluster decomposition analysis and apply it to both true FCI wave functions, as well as wave functions generated from the adaptive sampling CI algorithm. Furthermore, the cluster decomposition is useful for interpreting calculations in chemical studies, as a diagnostic for the convergence of various excitation manifolds, as well as as a guidepost for polynomially scaling electronic structure models. Applications are presented for (i) the double dissociation of water, (ii) the carbon dimer, (iii) the π space of polyacenes, and (iv) the chromium dimer. While the cluster amplitudes exhibit rapid decay with an increasing rank for the first three systems, even connected octuple excitations still appear important in Cr 2 , suggesting that spin-restricted single-reference coupled-cluster approaches may not be tractable for some problems in transition metal chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

General Forms of Wave Functions for Dipositronium, Ps2

The consequences of particle interchange symmetry for the structure of wave functions of the states of dipositronium was recently discussed by the author [I]. In the present work, the methodology is simply explained, and the wave functions are explicitly given.

Schrader, D.M.↗

Efficient antisymmetrization algorithm for the partially correlated wave functions in the free complement-local Schrödinger equation method

We propose here fast antisymmetrization procedures for the partially correlated wave functions that appear in the free complement-local Schrödinger equation (FC-LSE) method. Pre-analysis of the correlation diagram, referred to as dot analysis, combined with the determinant update technique based on the Laplace expansion, drastically reduces the orders of the antisymmetrization computations. When the complement functions include only up to single-correlated terms, the order of computations is O(N{sup 3}), which is the same as the non-correlated case. Similar acceleration is obtained for general correlated functions as a result of dot analysis. This algorithm has been successfully used in our laboratory in actual FC-LSE calculations for accurately solving the many-electron Schrödinger equations of atoms and molecules. The proposed method is general and applicable to the sampling-type methodology of other partially correlated wave functions like those in the quantum Monte Carlo and modern Hylleraas-type methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nuclear velocity perturbation theory for vibrational circular dichroism: An approach based on the exact factorization of the electron-nuclear wave function

The nuclear velocity perturbation theory (NVPT) for vibrational circular dichroism (VCD) is derived from the exact factorization of the electron-nuclear wave function. This new formalism offers an exact starting point to include correction terms to the Born-Oppenheimer (BO) form of the molecular wave function, similar to the complete-adiabatic approximation. The corrections depend on a small parameter that, in a classical treatment of the nuclei, is identified as the nuclear velocity. Apart from proposing a rigorous basis for the NVPT, we show that the rotational strengths, related to the intensity of the VCD signal, contain a new contribution beyond-BO that can be evaluated with the NVPT and that only arises when the exact factorization approach is employed. Numerical results are presented for chiral and non-chiral systems to test the validity of the approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational aspects of the continuum quaternionic wave functions for hydrogen

Over the past few years considerable attention has been given to the role played by the Hydrogen Continuum Wave Functions (HCWFs) in quantum theory. The HCWFs arise via the method of separation of variables for the time-independent Schrödinger equation in spherical coordinates. The HCWFs are composed of products of a radial part involving associated Laguerre polynomials multiplied by exponential factors and an angular part that is the spherical harmonics. In the present paper we introduce the continuum wave functions for hydrogen within quaternionic analysis ((R)QHCWFs), a result which is not available in the existing literature. In particular, the underlying functions are of three real variables and take on either values in the reduced and full quaternions (identified, respectively, with R{sup 3} and R{sup 4}). We prove that the (R)QHCWFs are orthonormal to one another. The representation of these functions in terms of the HCWFs are explicitly given, from which several recurrence formulae for fast computer implementations can be derived. A summary of fundamental properties and further computation of the hydrogen-like atom transforms of the (R)QHCWFs are also discussed. We address all the above and explore some basic facts of the arising quaternionic function theory. As an application, we provide the reader with plot simulations that demonstrate the effectiveness of our approach. (R)QHCWFs are new in the literature and have some consequences that are now under investigation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Second-Order Moller-Plesset Perturbation Theory for Molecular Dirac-Hartree-Fock Wave Functions

Moller-Plesset perturbation theory is developed to second order for a selection of Kramers restricted Dirac-Hartree-Fock closed and open-shell reference wave functions. The open-shell wave functions considered are limited to those with no more than two electrons in open shells, but include the case of a two-configuration SCF reference. Denominator shifts are included in the style of Davidson's OPT2 method. An implementation which uses unordered integrals with labels is presented, and results are given for a few test cases.

Dyall, Kenneth G.↗

Resonating quantum three-coloring wave functions for the kagome quantum antiferromagnet

Motivated by the recent discovery of a macroscopically degenerate exactly solvable point of the spin-1/2 XXZ model for J z /J = –1/2 on the kagome lattice —a result that holds for arbitrary magnetization—we develop an exact mapping between its exact ”quantum three-coloring” wave functions and the characteristic localized and topological magnons. This map, involving ”resonating two-color loops,” is developed to represent exact many-body ground state wave functions for special high magnetizations. Using this map we show that these exact ground state solutions are valid for any J z /J ≥ –1/2. This demonstrates the equivalence of the ground-state wave function of the Ising, Heisenberg, and XY regimes all the way to the J z /J = –1/2 point for these high magnetization sectors. In the hardcore bosonic language, this means that a certain class of exact many-body solutions, previously argued to hold for purely repulsive interactions (J z ≥ 0), actually hold for attractive interactions as well, up to a critical interaction strength. For the case of zero magnetization, where the ground state is not exactly known, we perform density matrix renormalization group calculations. Based on the calculation of the ground state energy and measurement of order parameters, we provide evidence for a lack of any qualitative change in the ground state on finite clusters in the Ising (J z >> J), Heisenberg (J z = J), and XY (J z = 0) regimes, continuing adiabatically to the vicinity of the macroscopically degenerate J z /J = –1/2 point. Furthermore, these findings offer a framework for recent results in the literature and also suggest that the J z /J = –1/2 point is an unconventional quantum critical point whose vicinity may contain the key to resolving the spin-1/2 kagome problem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Photoelectron Spectroscopy of CdSe Nanocrystals in the Gas Phase: A Direct Measure of the Evanescent Electron Wave Function of Quantum Dots

We present the first photoelectron spectroscopy measurements of quantum dots (semiconductor nanocrystals) in the gas phase. By coupling a nanoparticle aerosol source to a femtosecond velocity map imaging photoelectron spectrometer, we apply robust gas-phase photoelectron spectroscopy techniques to colloidal quantum dots, which typically must be studied in a liquid solvent or while bound to a surface. Working with a flowing aerosol of quantum dots offers the additional advantages of providing fresh nanoparticles for each laser shot and removing perturbations from bonding with a surface or interactions with the solvent. In this work, we perform a two-photon photoionization experiment to show that the photoelectron yield per exciton depends on the physical size of the quantum dot, increasing for smaller dots. Next, using effective mass modeling we show that the extent to which the electron wave function of the exciton extends from the quantum dot, the so-called “evanescent electron wavefunction”, increases as the size of the quantum dot decreases. We show that the photoelectron yield is dominated by the evanescent electron density due to quantum confinement effects, the difference in the density of states inside and outside of the quantum dots, and the angle-dependent transmission probability of electrons through the surface of the quantum dot. Therefore, the photoelectron yield directly reflects the fraction of evanescent electron wave function that extends outside of the quantum dot. This work shows that gas-phase photoelectron spectroscopy is a robust and general probe of the electronic structure of quantum dots, enabling the first direct measurements of the evanescent exciton wave function.

74 ATOMIC AND MOLECULAR PHYSICS↗

Many-fermion simulation from the contracted quantum eigensolver without fermionic encoding of the wave function

Quantum computers potentially have an exponential advantage over classical computers for the quantum simulation of many-fermion quantum systems. Nonetheless, fermions are more expensive to simulate than bosons due to the fermionic encoding—a mapping by which the qubits are encoded with fermion statistics. Here we generalize the contracted quantum eigensolver (CQE) to avoid fermionic encoding of the wave function. In contrast to the variational quantum eigensolver, the CQE solves for a many-fermion stationary state by minimizing the contraction (projection) of the Schrödinger equation onto two fermions. We avoid fermionic encoding of the wave function by contracting the Schrödinger equation onto an unencoded pair of particles. Solution of the resulting contracted equation by a series of unencoded two-body exponential transformations generates an unencoded wave function from which the energy and two-fermion reduced density matrix (2-RDM) can be computed. We apply the unencoded and the encoded CQE algorithms to the hydrogen fluoride molecule, the dissociation of oxygen O 2 , and a series of hydrogen chains. Both algorithms show comparable convergence towards the exact ground-state energies and 2-RDMs, but the unencoded algorithm has computational advantages in terms of state preparation and tomography.

74 ATOMIC AND MOLECULAR PHYSICS↗

Extracting Higher Central Charge from a Single Wave Function

A (2+1)D topologically ordered phase may or may not have a gappable edge, even if its chiral central charge c_ is vanishing. Recently, it was discovered that a quantity regarded as a "higher" version of chiral central charge gives a further obstruction beyond c_ to gapping out the edge. Here, in this Letter, we show that the higher central charges can be characterized by the expectation value of the partial rotation operator acting on the wave function of the topologically ordered state. This allows us to extract the higher central charge from a single wave function, which can be evaluated on a quantum computer. Our characterization of the higher central charge is analytically derived from the modular properties of edge conformal field theory, as well as the numerical results with the ν = 1/2 bosonic Laughlin state and the non-Abelian gapped phase of the Kitaev honeycomb model, which corresponds to U(1) 2 and Ising topological order, respectively. The Letter establishes a numerical method to obtain a set of obstructions to the gappable edge of (2+1)D bosonic topological order beyond c_, which enables us to completely determine if a (2+1)D bosonic Abelian topological order has a gappable edge or not. We also point out that the expectation values of the partial rotation on a single wave function put a constraint on the low-energy spectrum of the bulk-boundary system of (2+1)D bosonic topological order, reminiscent of the Lieb-Schultz-Mattis-type theorems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

General coalescence conditions for the exact wave functions: Higher-order relations for two-particle systems

We derived the necessary conditions that the non-relativistic time-independent exact wave functions for two-particle systems must satisfy at a coalescence (or cusp) point. Some of such necessary conditions are already known to be Kato's cusp condition (CC) and Rassolov and Chipman's CC. In the present study, we extended and generalized those conditions, calling them generalized coalescence conditions (GCCs). Kato's CC and Rassolov and Chipman's CC were shown to be specific cases included in the GCCs. The GCCs can be applied not only to Coulombic systems but also to any systems where the interaction between two particles is represented in a power series of the inter-particle distance. We confirmed the correctness of our derivation of these GCCs by applying the free complement wave functions of a hydrogen atom in ground and excited states, a harmonic oscillator, and a system with an interacting potential of V=r.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗