Search NASA⌕ Search

SEARCH · Search NASA

Results for “electronic wave function”

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

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↗

Ab Initio Calculations on the Ground and Excited Electronic States of Thorium–Ammonia, Thorium–Aza-Crown, and Thorium–Crown Ether Complexes

Positively charged metal–ammonia complexes are known to host peripheral, diffuse electrons around their molecular skeleton. The resulting neutral species form materials known as expanded or liquid metals. Alkali, alkaline earth, and transition metals have been investigated previously in experimental and theoretical studies of both the gas and condensed phase. This work is the first ab initio exploration of an f-block metal–ammonia complex. The ground and excited states are calculated for Th0–3+ complexes with ammonia, crown ethers, and aza-crown ethers. For Th3+ complexes, the one valence electron Th populates the metal’s 6d or 7f orbitals. For Th0–2+, the additional electrons prefer occupation of the outer s- and p-type orbitals of the complex, except Th(NH3)10, which uniquely places all four electrons in outer orbitals of the complex. Although thorium coordinates up to ten ammonia ligands, octa-coordinated complexes are more stable. Crown ether complexes have a similar electronic spectrum to ammonia complexes, but excitations of electrons in the outer orbitals of the complex are higher in energy. Aza-crown ethers disfavor the orbitals perpendicular to the crowns, attributed to the N-H bonds pointing along the plane of the crowns.

74 ATOMIC AND MOLECULAR PHYSICS↗

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-principles calculation of excited states of diatomic molecules: a benchmark for the Gutzwiller conjugate gradient minimisation method

We recently proposed the Gutzwiller conjugate gradient minimisation (GCGM) method for efficient and accurate calculation of the ground state total energy of molecular and bulk systems. The GCGM method is developed under the framework of Gutzwiller wave function but goes beyond the commonly adopted Gutzwiller approximation to improve the accuracy and flexibility in treating the correlation effects. In this conference proceeding, we benchmark the GCGM method with the calculation of excited state potential energy curves of three diatomic molecules, namely H 2 , N 2 , and O 2 . Our calculations demonstrate the flexibility and reasonable accuracy of the method.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dynamical electron correlation and the chemical bond. I. Covalent bonds in AH and AF (A = B–F)

Dynamical electron correlation has a major impact on the computed values of molecular properties and the energetics of molecular processes. This study focused on the effect of dynamical electron correlation on the spectroscopic constants ($R$ e , $ω$ e , $D$ e ), and potential energy curves, Δ$E$($R$), of the covalently bound AH and AF molecules, A = B–F. The changes in the spectroscopic constants (Δ$R$ e , Δ$ω$ e , Δ$D$ e ) caused by dynamical correlation are erratic and, at times, even surprising. These changes can be understood based on the dependence of the dynamical electron correlation energies of the AH and AF molecules as a function of the bond distance, i.e., Δ$E$ DEC ($R$). At large $R$, the magnitude of Δ$E$ DEC ($R$) increases nearly exponentially with decreasing $R$, but this increase slows as $R$ continues to decrease and, in many cases, even reverses at very short $R$. The changes in Δ$E$ DEC ($R$) in the region around $R$ e were as unexpected as they were surprising, e.g., distinct minima and maxima were found in the curves of Δ$E$ DEC ($R$) for the most polar molecules. The variations in Δ$E$ DEC ($R$) for $R$ ≲ $R$ e are directly correlated with major changes in the electronic structure of the molecules as revealed by a detailed analysis of the spin-coupled generalized valence bond wave function. In conclusion, the results reported here indicate that we have much to learn about the nature of dynamical electron correlation and its effect on chemical bonds and molecular properties and processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Classical Preoptimization Approach for ADAPT-VQE: Maximizing the Potential of High-Performance Computing Resources to Improve Quantum Simulation of Chemical Applications

The ADAPT-VQE algorithm is a promising method for generating a compact ansatz based on derivatives of the underlying cost function, and it yields accurate predictions of electronic energies for molecules. In this work, we report the implementation and performance of ADAPT-VQE with our recently developed sparse wave function circuit solver (SWCS) in terms of accuracy and efficiency for molecular systems with up to 52 spin orbitals. The SWCS can be tuned to balance computational cost and accuracy, which extends the application of ADAPT-VQE for molecular electronic structure calculations to larger basis sets and a larger number of qubits. Using this tunable feature of the SWCS, we propose an alternative optimization procedure for ADAPT-VQE to reduce the computational cost of the optimization. Furthermore, by preoptimizing a quantum simulation with a parametrized ansatz generated with ADAPT-VQE/SWCS, we aim to utilize the power of classical high-performance computing in order to minimize the work required on noisy intermediate-scale quantum hardware, which offers a promising path toward demonstrating quantum advantage for chemical applications.

ADAPT-VQE↗

Generation and Control of Self-Organized Nonlinear Kinetic Structures in High Energy Density Plasmas in the Presence of Intense Magnetic Fields and Ultrashort Laser Pulses

Goals were to study the interplay between electron plasma waves (EPW), KEEN waves and externally generated magnetic fields. In particular, the Weibel instability B field generation process and its interrelationship with the existence of nonlinear EPWs in high energy density plasmas. We focused on a number of models of how Kinetic, Nonlinear Electron Plasma Waves, KNL-EPW and KEEN waves create anisotropic electron velocity distribution functions, e- VDF, and how these anisotropic e- VDFs in turn drive the Weibel instability and generate B fields. Our goal is to control the SRS and SKEENS processes that generate the KNL-EPW, control the anisotropy, and thus also control the dynamics of the resulting B fields, their influence on the transport coefficients and heat transport that results, their modification of SRS itself and the reinforced anisotropy driven loop gain.

(Kinetic electrostatic electron nonlinear) KEEN wa↗

Measurement of Thomson-scattering spectra with continuous angular resolution (invited)

Here, a novel Thomson-scattering diagnostic with continuous angular resolution over a span of 120° was developed for the characterization of plasmas produced at the Omega Laser Facility. Spectrally resolving light scattered from electron plasma wave features as a function of emission angle provides a means to efficiently probe a large range of plasma frequencies and k vectors. Together, these spectra contain critical constraints on the plasma-physics models used to interpret the data and enable experimental measurements of the electron-velocity distribution function over several orders of magnitude without assumptions about its mathematical form. Major components of the instrument include (1) a reflective collection objective that gathers light over a range of 120° × 12°; (2) a spatial-filter image relay for measurement localization; (3) cylindrical optics for producing a line image of the collection aperture; (4) a transmission grating spectrometer; and (5) a time-gated, image-intensified camera. Thomson-scattered light collected from an ~(50 – μM) 3 volume of plasma is recorded with 0.8-nm spectral and 1° angular resolution. Initial experiments examined the properties of the electron-velocity distribution in gas-jet-produced plasmas in the presence of heating via inverse bremsstrahlung absorption.

47 OTHER INSTRUMENTATION↗

Average-atom calculations of bound-free and free-free cross sections in dense plasmas

Cross sections for photo-induced bound-free and free-free transitions in plasmas are evaluated in the average-atom approximation and applied to determine opacities of dense plasmas of light elements. Parameters characterizing the plasmas (chemical potential, average ionic charge, free electron density, bound and continuum wave functions, and occupation numbers) are obtained from the average-atom model. Lowest-order calculations of the free-free cross sections, which diverge in the low-frequency limit, are regularized by accounting for the finite electron-ion relaxation time. The resulting analysis provides the basis for average-atom studies of plasma opacities. Such studies are presented for dense lithium, beryllium, boron, and carbon. Applications are given to Rosseland mean opacities of dense hydrogen and deuterium plasmas and to comparisons of free-free to bound-free opacities in shock-compressed plasmas. Average-atom cross section and opacity calculations are extended to plasmas consisting of more than one ionic species, boron nitride, polystyrene, and a composite H, He, C plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Data for "Plasma shape effects on the Alfvén eigenmode spectrum through Alfvén slow-magnetosonic wave coupling"

electron density profile as a function of normalized poloidal flux: psinorm ne(10^20/m^3) dne/dpsiN Pressure profile as a function of normalized poloidal flux: psinorm ptot(KPa) dptot/dpsiN scaled to oftain beta=5% at magnetic axis Magnetic safety factor profile as a function of normalized poloidal flux: psinorm q_profile All needed to perform the NOVA calculations.

Alfven waves↗

Nonclassical Spin Transfer Effects in an Antiferromagnet

We simulate scattering of electrons by a chain of antiferromagnetically coupled quantum Heisenberg spins, to analyze spin-transfer effects not described by the classical models of magnetism. Our simulations demonstrate efficient excitation of dynamical states that would be forbidden by the semiclassical symmetries, such as generation of multiple magnetic excitation quanta by a single electron. Furthermore, quantum interference of spin wave functions enables generation of magnetization dynamics with amplitudes exceeding the transferred magnetic moment. Here, the efficiency of excitation is almost independent of the electron spin polarization, and is governed mainly by the transfer of energy. Nonclassical spin transfer may thus enable efficient electronic control of antiferromagnets not limited by the classical constraints.

36 MATERIALS SCIENCE↗

Nanoscale Imaging of Magnetotransport around a Circular 𝑝−𝑛 Junction in Graphene

Magnetoresistance studies of 2D systems are often shaped by the motion of electrons that occupy spatially confined wave functions, such as topological edge modes and disorder-induced bound states. Directly probing how such states form and behave in situ , under applied currents, provides a clear way of connecting microscopic physics to the macroscopic transport response. Here, in this Letter, scanning tunneling potentiometry is used to probe the local, current-induced electrochemical potential of carriers in graphene near circular 𝑝−𝑛 junctions in an out-of-plane magnetic field ranging from 0 to 1.4 T. These measurements provide detailed information about the motion of carriers at the nanometer scale, revealing how it evolves with increasing field. The electrochemical potential displays distinct patterns, such as dipoles, spirals, and concentric disks in weak, moderate, and high fields, respectively. The size and orientation of these patterns can be used to understand how local carrier dynamics change in different transport regimes as well as to directly extract physical parameters such as the electron mean free path and cyclotron diameter.

36 MATERIALS SCIENCE↗

On the nonadiabatic collisional quenching of OH(A) by H 2 : a four coupled quasi-diabatic state description

A four-state diabatic potential energy matrix (DPEM), H d , for the description of the nonadiabatic quenching of OH(A 2 Σ + ) by collisions with H 2 is reported. The DPEM is constructed as a fit to adiabatic energies, energy gradients, and derivative couplings obtained exclusively from multireference configuration interaction wave functions. A four-adiabatic-electronic-state representation is used in order to describe all energetically accessible regions of the nuclear coordinate space. Partial permutation-inversion symmetry is incorporated into the representation. Here, the fit is based on electronic structure data at 42 882 points, described by over 1.6 million least squares equations with a root mean square (mean unsigned) error of 178(83) cm –1 . Comparison of ab initio and H d determined minima, saddle points, and energy minimized points on C 2v , C s , C ∞v , and C 1 (noncoplanar) portions of two conical intersection seams are used to establish the accuracy of the H d .

74 ATOMIC AND MOLECULAR PHYSICS↗

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code was efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Ground state wave functions for single-band Hubbard models from the Gutzwiller conjugate gradient minimisation theory

The Gutzwiller conjugate gradient minimisation (GCGM) theory is an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. Previous applications of GCGM provides satisfying accuracy of ground-state energy of Hubbard models. In the current work, we address the problem of whether the obtained wave function is a good approximation for the true ground state by comparing the correlation functions with the benchmark data. Additionally, our results confirms the accuracy of the reproduced ground state of the regular Hubbard model, but with some exception for the frustrated Hubbard model.

74 ATOMIC AND MOLECULAR PHYSICS↗

Application of the Screened Hydrogenic Model to Light Atoms

The purpose of this work was to develop a Screened Hydrogenic Model (SHM) to accurately calculate the electron energies for light atoms and ions with up to ten electrons for atomic numbers up to 18. The total energy of an atom or ion was calculated with effective nuclear charges and screening parameters for each electron type (1s, 2s and 2p) within a specific electron configuration. Multiple energy states, centered at the total energy, were calculated for elec-tron configurations that have Russell Saunders coupling. The energy of each electron included its relativistic energy, EREL, but close overall agreement between the calculated and experimental energies of multi-electron configura-tions required that the one-electron expression for EREL be modified in a simple manner. In the present work, 98% of the 587 calculated energies for light atoms/ions have a relative error within ±0.1% of the corresponding experimental energies. The effective nuclear charges described in this work allow hydrogen-like wave functions to be defined for the electrons within a multi-electron configuration. The SHM, described in this work, is available for future calculations involving light atoms and ions.

36 MATERIALS SCIENCE↗