Search NASA⌕ Search

SEARCH · Search NASA

Results for “trial 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cohesion and excitations of diamond-structure silicon by quantum Monte Carlo: Benchmarks and control of systematic biases

We have carried out quantum Monte Carlo (QMC) calculations of silicon crystal focusing on the accuracy and systematic biases that affect the electronic structure characteristics. The results show that 64 and 216 atom supercells provide an excellent consistency for extrapolated energies per atom in the thermodynamic limit for ground, excited, and ionized states. We have calculated the ground state cohesion energy with both systematic and statistical errors below ≈ 0.05 eV. The ground state exhibits a fixed-node error of only 1.3(2)% of the correlation energy, suggesting an unusually high accuracy of the corresponding single-reference trial wave function. We obtain a very good agreement between optical and quasiparticle gaps that affirms the marginal impact of excitonic effects. Our most accurate results for band gaps differ from the experiments by about 0.2 eV. This difference is assigned to a combination of residual finite-size and fixed-node errors. In conclusion, we have estimated the crystal Fermi level referenced to vacuum that enabled us to calculate the edges of valence and conduction bands in agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Nuclear states projected from a pair condensate

Atomic nuclei exhibit deformation, pairing correlations, and rotational symmetries. To meet these competing demands in a computationally tractable formalism, we revisit the use of general pair condensates with good particle number as a trial wave function for even-even nuclei. After minimizing the energy of the condensate, allowing for general triaxial deformations, we project out states with good angular momentum with a fast projection technique. To show applicability, we present example calculations from pair condensates in several model spaces and compare against angular-momentum projected Hartree-Fock and full configuration-interaction shell-model calculations. This approach successfully generates spherical, vibrational, and rotational spectra, demonstrating potential for modeling medium- to heavy-mass nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Interlayer coupling in two-dimensional MoS 2 and phosphorene bilayers: Benchmark quantum Monte Carlo study of interaction energies and quasiparticle band gaps

Using high-accuracy many-body quantum Monte Carlo (QMC) methods, we study the effect of interlayer coupling on the properties of two-dimensional freestanding bilayers (BLs) of MoS 2 and phosphorene. The properties of the two BL-materials are very different and largely determined by the interlayer interaction, which is purely van der Waals in MoS 2 and partially electronic/chemical in phosphorene, resulting in a modest layer-dependent property modulation in MoS 2 and strong modulation in phosphorene. Multireference and symmetry considerations are used to construct state-of-the-art accuracy QMC trial wave functions. We determine the quasiparticle band gaps for both materials, $Δ^{\textrm{qp}}_{Γ→\textrm{K}}$ = 2.45 ± 0.05 eV in BL-MoS 2 and $Δ^{\textrm{qp}}_{Γ→Γ}$=1.59 ±0.1eV in BL-phosphorene. In conclusion, these benchmark band-gap values make it possible to consolidate the interpretation of the widely scattered experimental and theory data.

Huang, Yongda [Slovak Academy of Sciences (SAS), B↗
DOE OSTI2022

VQE method: a short survey and recent developments

Abstract The variational quantum eigensolver (VQE) is a method that uses a hybrid quantum-classical computational approach to find eigenvalues of a Hamiltonian. VQE has been proposed as an alternative to fully quantum algorithms such as quantum phase estimation (QPE) because fully quantum algorithms require quantum hardware that will not be accessible in the near future. VQE has been successfully applied to solve the electronic Schrödinger equation for a variety of small molecules. However, the scalability of this method is limited by two factors: the complexity of the quantum circuits and the complexity of the classical optimization problem. Both of these factors are affected by the choice of the variational ansatz used to represent the trial wave function. Hence, the construction of an efficient ansatz is an active area of research. Put another way, modern quantum computers are not capable of executing deep quantum circuits produced by using currently available ansatzes for problems that map onto more than several qubits. In this review, we present recent developments in the field of designing efficient ansatzes that fall into two categories—chemistry–inspired and hardware–efficient—that produce quantum circuits that are easier to run on modern hardware. We discuss the shortfalls of ansatzes originally formulated for VQE simulations, how they are addressed in more sophisticated methods, and the potential ways for further improvements.

Fedorov, Dmitry A. (ORCID:0000000316598580)↗
NASA NTRS1972

D autoionization states of He and ionic H

Positions of the lowest 1,3De autoionization states of He and H(-) below the n = 2 level of the He(+) and H were calculated variationally, using Feshbach's Q-operator formalism. The trial wave function is of the Hylleraas-type with appropriate angular momentum factors. The widths and the shifts of the states have also been calculated. The shifts are found to be positive for all the states calculated here. The results with 112 terms for most states are lower than any previously calculated. The calculated lowest autoionization states of the He and H(-) (relative to the ground states of He and H respectively) are 59.902 eV and 10.1185 eV, in good agreement with the observed values of 59.9 eV and 10.13 + or 0.015 eV.

Bhatia, A. K.↗
NASA NTRS1986

Variational calculations of subbands in a quantum well with uniform electric field - Gram-Schmidt orthogonalization approach

Variational calculations of subband eigenstates in an infinite quantum well with an applied electric field using Gram-Schmidt orthogonalized trial wave functions are presented. The results agree very well with the exact numerical solutions even up to 1200 kV/cm. It is also shown that, for increasing electric fields, the energy of the ground state decreases, while that of higher subband states increases slightly up to 1000 kV/cm and then decreases for a well size of 100 A.

Ahn, Doyeol↗
DOE OSTI2024

Ab initio description of bcc iron with correlation matrix renormalization theory

We applied the ab initio spin-polarized correlation matrix renormalization theory to the ferromagnetic state of the bulk bcc iron. We showed that it was capable of reproducing the equilibrium physical properties and the pressure-volume curve in good comparison with experiments. We then focused on the analysis of its local electronic correlations. By exploiting different local fluctuation-related physical quantities as measures of electronic correlation within target orbits, we elucidated the different roles of t 2g and e g states in both spin channels and presented compelling evidence to showcase this distinction in their electronic correlation.

3-dimensional systems↗
NASA NTRS1986

Variational divergence in wave scattering theory with Kirchhoffean trial functions

In a recent study of variational improvement of the Kirchhoff approximation for electromagnetic scattering by rough surfaces, a key ingredient in the variational principle was found to diverge for important configurations (e.g., backscatter) if the polarization had any vertical component. The cause and a cure of this divergence are discussed here. The divergence is demonstrated to occur for arbitrary perfectly conducting scatterers and its universal characterstics are determined, by means of a general divergence criterion that is derived. A variational cure for the divergence is prescribed, and it is tested successfully on a standard scattering model.

Bird, J. F.↗
DOE OSTI2025

Systematic improvement of trial states in phaseless auxiliary-field quantum Monte Carlo

We extend the use of coupled cluster (CC) trial states in the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method beyond single and double excitations to include both triple and quadruple excitations. With this AFQMC/CC hierarchy, we are able to systematically benchmark the method's performance on molecular systems as the quality of the trial is improved. Our results show that the phaseless AFQMC energy improves systematically and is typically significantly more accurate than the energy of the underlying trial state. However, the relative improvement compared to the trial CC energy decreases as we ascend the CC hierarchy. As the CC wavefunction is usually further approximated when used as an AFQMC trial, we also explore the relationship between the components of the CC wavefunction and the resulting AFQMC/CC error. Our results suggest that improving the representation of the CC wave function in the AFQMC trial does not always lower the bias even when it increases the fidelity of the trial with the exact ground state.

Chemical Physics (physics.chem-ph)↗
NASA NTRS1975

Positronium-hydrogen elastic scattering

Elastic scattering of positronium atoms by hydrogen atoms below the excitation threshold at 5.1 eV is investigated with the investigation limited to the spin-zero (singlet) state of the two electrons. Nonresonant S-wave phase shifts have been obtained from normalizable trial functions by extracting the center-of-mass wave function and examining its spatial behavior. An effective-range expansion is fitted to the phase shifts and the previously obtained positronium-hydride bound-state energy to yield a scattering length approximately equal to 5.3 Bohr radii as compared with Fraser's (1961) corrected exchange result of 7.28 Bohr radii. The stabilization method suggests that a resonance occurs at a scattering energy of 4.5 eV, about 0.6 eV below the excitation threshold. Using the complex-rotation method, evidence is obtained which confirms the existence of the resonance, and its width is estimated to be about 0.06 eV.

Drachman, R. J.↗
NASA NTRS1977

Correlated wave functions for three-particle systems with Coulomb interaction - The muonic helium atom

A computational procedure for calculating correlated wave functions is proposed for three-particle systems interacting through Coulomb forces. Calculations are carried out for the muonic helium atom. Variational wave functions which explicitly contain interparticle coordinates are presented for the ground and excited states. General Hylleraas-type trial functions are used as the basis for the correlated wave functions. Excited-state energies of the muonic helium atom computed from 1- and 35-term wave functions are listed for four states.

Huang, K.-N.↗
NASA NTRS1986

Electromagnetic scattering theory

Electromagnetic scattering theory is discussed with emphasis on the general stochastic variational principle (SVP) and its applications. The stochastic version of the Schwinger-type variational principle is presented, and explicit expressions for its integrals are considered. Results are summarized for scalar wave scattering from a classic rough-surface model and for vector wave scattering from a random dielectric-body model. Also considered are the selection of trial functions and the variational improvement of the Kirchhoff short-wave approximation appropriate to large size-parameters. Other applications of vector field theory discussed include a general vision theory and the analysis of hydromagnetism induced by ocean motion across the geomagnetic field. Levitational force-torque in the magnetic suspension of the disturbance compensation system (DISCOS), now deployed in NOVA satellites, is also analyzed using the developed theory.

Bird, J. F.↗
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 OSTI2021

Fast & accurate emulation of two-body scattering observables without wave functions

We combine Newton's variational method with ideas from eigenvector continuation to construct a fast & accurate emulator for two-body scattering observables. The emulator will facilitate the application of rigorous statistical methods for interactions that depend smoothly on a set of free parameters. Our approach begins with a trial K or T matrix constructed from a small number of exact solutions to the Lippmann–Schwinger equation. Subsequent emulation only requires operations on small matrices. We provide several applications to short-range potentials with and without the Coulomb interaction and partial-wave coupling. It is shown that the emulator can accurately extrapolate far from the support of the training data. When used to emulate the neutron-proton cross section with a modern chiral interaction as a function of 26 free parameters, it reproduces the exact calculation with negligible error and provides an over 300x improvement in CPU time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

Excited electronic states of O2/-/.

Configuration interaction (CI) variational trial functions (which assure formally correct dissociation) and the single-configuration self-consistent-field approximation are applied to the calculation of the excited electronic states of the negatively charged oxygen molecule. The energy curves and wave functions of these states are analyzed. All the excited states are found to have equilibrium-internuclear separations at least 1 a.u. larger than the ground state. The two lowest energy states are characterized, respectively, as shape and valence Feshbach resonances.

Krauss, M.↗