Search NASA⌕ Search

SEARCH · Search NASA

Results for “SLATER METHOD”

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 19 records

Slater transition methods for core-level electron binding energies

Methods for computing core-level ionization energies using self-consistent field (SCF) calculations are evaluated and benchmarked. These include a “full core hole” (or “ΔSCF”) approach that fully accounts for orbital relaxation upon ionization, but also methods based on Slater’s transition concept in which the binding energy is estimated from an orbital energy level that is obtained from a fractional-occupancy SCF calculation. A generalization that uses two different fractional-occupancy SCF calculations is also considered. The best of the Slater-type methods afford mean errors of 0.3–0.4 eV with respect to experiment for a dataset of K-shell ionization energies, a level of accuracy that is competitive with more expensive many-body techniques. An empirical shifting procedure with one adjustable parameter reduces the average error below 0.2 eV. Here, this shifted Slater transition method is a simple and practical way to compute core-level binding energies using only initial-state Kohn–Sham eigenvalues. It requires no more computational effort than ΔSCF and may be especially useful for simulating transient x-ray experiments where core-level spectroscopy is used to probe an excited electronic state, for which the ΔSCF approach requires a tedious state-by-state calculation of the spectrum. As an example, we use Slater-type methods to model x-ray emission spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Excited State Specific Multi-Slater Jastrow Wave Functions

We combine recent advances in excited state variational principles, fast multi-Slater Jastrow methods and selective configuration interaction, to create multi-Slater Jastrow wave function approximations that are optimized for individual excited states. In addition to the Jastrow variables and linear expansion coefficients, this optimization includes state-specific orbital relaxations in order to avoid the compromises necessary in state-averaged approaches. Here, we demonstrate that, when combined with variance matching to help balance the quality of the approximation across different states, this approach delivers accurate excitation energies even when very modest multi-Slater expansions are used. Intriguingly, this accuracy is maintained even when we study a difficult chlorine-anion-to-π* charge transfer in which traditional state-averaged multireference methods must contend with different states that require drastically different orbital relaxations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio calculations in atoms, molecules, and solids, treating spin–orbit coupling and electron interaction on an equal footing

We incorporate explicit, non-perturbative treatment of spin–orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin–orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br 2 and I 2 , and solid Bi.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Monte Carlo calculation of the binding energy of the beryllium dimer

The accurate calculation of the binding energy of the beryllium dimer is a challenging theoretical problem. In this study, the binding energy of Be{sub 2} is calculated using the diffusion Monte Carlo (DMC) method, using single Slater determinant and multiconfigurational trial functions. DMC calculations using single-determinant trial wave functions of orbitals obtained from density functional theory calculations overestimate the binding energy, while DMC calculations using Hartree-Fock or CAS(4,8), complete active space trial functions significantly underestimate the binding energy. In order to obtain an accurate value of the binding energy of Be{sub 2} from DMC calculations, it is necessary to employ trial functions that include excitations outside the valence space. Our best estimate DMC result for the binding energy of Be{sub 2}, obtained by using configuration interaction trial functions and extrapolating in the threshold for the configurations retained in the trial function, is 908 cm{sup −1}, only slightly below the 935 cm{sup −1} value derived from experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A quantum Monte Carlo study of the molybdenum dimer (Mo 2 )

In this report we have studied the molybdenum dimer (Mo 2 ) system. The binding energy was calculated by means of the fixed-node DMC (FN-DMC) method. The Slater part of the trial wave function was constructed by the Selected-CI method by using the orbitals generated by the KS-DFT method with a hybrid meta-GGA exchange and correlation functional, TPSSh. We also carried out CCSD(T) calculations which were subsequently extrapolated to the complete basis set (CBS) limit. The results are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optical characterization of ZnO nanomaterial with praseodymium ions

ZnO nanomaterial with praseodymium ions was prepared by chemical synthesis method. The ZnO nanomaterial was characterized by XRD, SEM and TEM. Their absorption in UV-VIS/NIR regions was measured at room temperature. The experimental oscillator strengths were calculated from the areas under the absorption bands. Eight absorption bands have been observed. From these spectral data various energy interaction parameters like Slater–Condon, Lande, Racah, Nephelauxetic ratio and bonding parameters have been computed. Judd-Ofelt analysis has been carried out using the absorption spectra to evaluate the radiative properties for luminescent levels of the praseodymium ion and discussed. The observed nano particle size is 2nm.

36 MATERIALS SCIENCE↗

Non-orthogonal multi-Slater determinant expansions in auxiliary field quantum Monte Carlo

We investigate here the use of non-orthogonal multi-Slater determinant (NOMSD) expansions as trial wavefunctions in auxiliary field quantum Monte Carlo simulations of molecular systems. We show that NOMSD trial wavefunctions with as few as twenty determinants are sufficient in order to achieve chemical accuracy across most of the G1 molecular test set. We also show that NOMSD trial wavefunctions are useful for more challenging strongly correlated systems by computing relative energies along the isomerization path of the [Cu 2 O 2 ] 2+ molecule. Our results for [Cu 2 O 2 ] 2+ compare favorably with other accurate quantum many-body methods, including density matrix renormalization group and completely renormalized coupled cluster methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Approximate singly excited states from a two-component Hartree-Fock reference

For many molecules, relaxing the spin symmetry constraint on the wave function results in the lowest energy mean-field solution. The two-component Hartree-Fock (2cHF) method relaxes all spin symmetry constraints, and the wave function is no longer an eigenfunction of the total spin, spin projection, or time-reversal symmetry operators. For ground state energies, 2cHF is a superior mean-field method for describing spin-frustrated molecules. For excited states, the utility of 2cHF is uncertain. Here, we implement the 2cHF extensions of two single-reference excited state methods, the two-component configuration interaction singles and time-dependent Hartree-Fock. We compare the results to the analogous methods based off of the unrestricted Hartree-Fock approximation, as well as the full configuration interaction for three small molecules with distinct 2cHF solutions, and discuss the nature of the 2cHF excited state solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Properties of coupled-cluster equations originating in excitation sub-algebras

In this paper we discuss properties of single-reference coupled cluster (CC) equations associated with the existence of nilpotent sub-algebras of excitations that allow one to represent CC equations in a hybrid fashion where the cluster amplitudes associated with these sub-algebras can be obtained by solving corresponding eigenvalue problem. For closed-shell formulations analyzed in this paper, the hybrid representation of CC equations provides a natural way for extending concepts behind actives-spaces and seniority numbers to provide more accurate description of electron correlation effects. Moreover, new representation can be utilized to re-define iterative algorithms used to solve CC equations, especially for though cases defined by the presence of strong static and dynamic correlation effects. We will also explore invariance properties associated with the nilpotent sub-algebras to define a new class of CC approximations referred in this paper to as the nilpotent-flow-based CC methods. We illustrate the performance of these methods on the example of ground- and excited-state calculations for commonly used small benchmark systems.

Kowalski, Karol↗

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations

In this research, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cooper minima in high-Z atoms: effects of correlation and relativity on np photoionization

Abstract Photoionization dipole transition matrix elements pass through a zero or attain a minimum that leaves imprints on photoionization parameters like the cross-section, angular distribution asymmetry parameter, phase shift, and photoionization time delay. This minimum is commonly known as the ‘Cooper minimum’ (CM). The CM, in general, is strongly affected by relativistic and correlation effects. Previous works investigated CM in the 6pand 5psubshell photoionization up toZ= 100 using the single-particle Dirac-Slater (DS) method. The present work extends the earlier work toZup to 120 using more accurate methods; Dirac–Hartree–Fock (DHF) which includes the relativistic effects and exchange correlations, and the relativistic random phase approximation (RRPA) which includes both initial and final state electron-electron correlations along with relativistic effects. In addition to the study of photoionization from the 6pand 5psubshells, the 4psubshell has also been investigated in the present work. To demonstrate the prominent effects in the high-Z atoms, Rn (Z= 86), Ra (Z= 88), No (Z= 102), Cn (Z = 112), Og (Z= 118), and Ubn (Z= 120) are investigated.

Optics↗

Assessment of multireference approaches to explicitly correlated full configuration interaction quantum Monte Carlo

The Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method has proved able to provide near-exact solutions to the electronic Schrödinger equation within a finite orbital basis set, without relying on an expansion about a reference state. However, a drawback to the approach is that being based on an expansion of Slater determinants, the FCIQMC method suffers from a basis set incompleteness error that decays very slowly with the size of the employed single particle basis. The FCIQMC results obtained in a small basis set can be improved significantly with explicitly correlated techniques. Here, we present a study that assesses and compares two contrasting “universal” explicitly correlated approaches that fit into the FCIQMC framework: the [2]{sub R12} method of Kong and Valeev [J. Chem. Phys. 135, 214105 (2011)] and the explicitly correlated canonical transcorrelation approach of Yanai and Shiozaki [J. Chem. Phys. 136, 084107 (2012)]. The former is an a posteriori internally contracted perturbative approach, while the latter transforms the Hamiltonian prior to the FCIQMC simulation. These comparisons are made across the 55 molecules of the G1 standard set. We found that both methods consistently reduce the basis set incompleteness, for accurate atomization energies in small basis sets, reducing the error from 28 mE{sub h} to 3-4 mE{sub h}. While many of the conclusions hold in general for any combination of multireference approaches with these methodologies, we also consider FCIQMC-specific advantages of each approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pressure induced structural phase transition in IB transition metal nitrides compounds

Transition metal mononitrides are known as refractory compounds, and they have, relatively, high hardness, brittleness, melting point, and superconducting transition temperature, and they also have interesting optical, electronic, catalytic, and magnetic properties. Evolution of structural properties would be an important step towards realizing the potential technological scenario of this material of class. In the present study, an effective interionic interaction potential (EIOP) is developed to investigate the pressure induced phase transitions in IB transition metal nitrides TMN [TM = Cu, Ag, and Au] compounds. The long range Coulomb, van der Waals (vdW) interaction and the short-range repulsive interaction upto second-neighbor ions within the Hafemeister and Flygare approach with modified ionic charge are properly incorporated in the EIOP. The vdW coefficients are computed following the Slater-Kirkwood variational method, as both the ions are polarizable. The estimated value of the phase transition pressure (Pt) and the magnitude of the discontinuity in volume at the transition pressure are consistent as compared to the reported data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exchange correlation potentials from full configuration interaction in a Slater orbital basis

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. As a result, examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

74 ATOMIC AND MOLECULAR PHYSICS↗

State Preparation of Antisymmetrized Geminal Power on a Quantum Computer without Number Projection

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen–Cooper–Schrieffer (PBCS) wave function. It is also an elementary symmetric polynomial (ESP) state. We generalize previous research on deterministically implementing the Dicke state to a state preparation algorithm for an ESP state, or equivalently AGP, on a quantum computer. Our method is deterministic and has polynomial cost, and it does not rely on number symmetry breaking and restoration. We also show that our circuit is equivalent to a disentangled unitary paired coupled cluster operator and a layer of unitary Jastrow operator acting on a single Slater determinant. Here, the method presented herein highlights the ability of disentangled unitary coupled cluster to capture nontrivial entanglement properties that are hardly accessible with traditional Hartree–Fock based electronic structure methods.

97 MATHEMATICS AND COMPUTING↗

Geminal embedding scheme for optimal atomic basis set construction in correlated calculations

We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-dependent projected Hartree-Fock

Projected Hartree-Fock (PHF) has recently emerged as an alternative approach to describing degenerate systems where static correlation is abundant, when the spin-symmetry is projected. Here, we derive a set of linearized time-dependent equations for PHF in order to be able to access excited states. The close connection of such linear-response time-dependent PHF (TDPHF) to the stability condition of a PHF wave function is discussed. Expanding this analysis also makes it possible to give analytical expressions for the projected coupling terms of Hamiltonian and overlaps between excited Slater determinants. TDPHF with spin-projection (TDSUHF) and its Tamm-Dancoff approximation are benchmarked for several electronically degenerate molecules including the dissociating H{sub 2}, F{sub 2} and O{sub 3} at equilibrium, and the distorted ethylene. It is shown that they give consistently better descriptions of excited states than does time-dependent HF (TDHF). Furthermore, we demonstrate that they offer not only singly but also doubly excited states, which naturally arise upon spin-projection. We also address the thermodynamic limit of TDSUHF, using non-interacting He gas. While TDPHF singly excited states tend to converge to those of HF with the size of the system due to the lack of size-extensivity of PHF, doubly excited states remain reasonable even at the thermodynamic limit. We find that the overall performance of our method is systematically better than the regular TDHF in many cases at the same computational scaling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A hybrid configuration interaction treatment based on seniority number and excitation schemes

We present a configuration interaction method in which the Hamiltonian of an N-electron system is projected on Slater determinants selected according to the seniority-number criterion along with the traditional excitation-based procedure. This proposed method is especially useful to describe systems which exhibit dynamic (weak) correlation at determined geometric arrangements (where the excitation-based procedure is more suitable) but show static (strong) correlation at other arrangements (where the seniority-number technique is preferred). The hybrid method amends the shortcomings of both individual determinant selection procedures, yielding correct shapes of potential energy curves with results closer to those provided by the full configuration interaction method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗