Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CONFIGURATION INTERACTION”

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 127 records · Page 7

Relativistic Configuration-Interaction and Perturbation Theory Calculations for Heavy Atoms

Heavy atoms present challenges to atomic theory calculations due to the large number of electrons and their complicated interactions. Conventional approaches such as calculations based on Cowan’s code are limited and require a large number of parameters for energy agreement. One promising approach is relativistic configuration-interaction and many-body perturbation theory (CI-MBPT) methods. We present CI-MBPT results for various atomic systems where this approach can lead to reasonable agreement: La I, La II, Th I, Th II, U I, Pu II. Among atomic properties, energies, g-factors, electric dipole moments, lifetimes, hyperfine structure constants, and isotopic shifts are discussed. While in La I and La II accuracy for transitions is better than that obtained with other methods, more work is needed for actinides.

74 ATOMIC AND MOLECULAR PHYSICS↗

Complex multireference configuration interaction calculations for the K-vacancy Auger states of N{sup q+} (q = 2-5) ions

K-vacancy Auger states of N{sup q+} (q = 2-5) ions are studied by using the complex multireference single- and double-excitation configuration interaction (CMRD-CI) method. The calculated resonance parameters are in good agreement with the available experimental and theoretical data. It shows that the resonance positions and widths converge quickly with the increase of the atomic basis sets in the CMRD-CI calculations; the standard atomic basis set can be employed to describe the atomic K-vacancy Auger states well. The strong correlations between the valence and core electrons play important roles in accurately determining those resonance parameters, Rydberg electrons contribute negligibly in the calculations. Note that it is the first time that the complex scaling method has been successfully applied for the B-like nitrogen. CMRD-CI is readily extended to treat the resonance states of molecules in the near future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tutorial for Non-Orthogonal Configuration Interaction Fragment Studies Using GronOR

A brief tutorial is presented for the GronOR non-orthogonal configuration interaction application illustrating the download and build process, setting up calculations with OpenMolcas and GronOR, and a description of the file structure and input options. An overview of the computational workflow includes the generation of fragment wave functions and integrals with an adapted version of OpenMolcas, running GronOR calculations, and accounting for dynamic correlation corrections. Three end-to-end use cases are provided, representing a molecular dimer, an ion pair, and overlapping molecular fragments.

Straatsma, T.P.↗

A High Performance Block Eigensolver for Nuclear Configuration Interaction Calculations

As on-node parallelism increases and the performance gap between the processor and the memory system widens, achieving high performance in large-scale scientific applications requires an architecture-aware design of algorithms and solvers. We focus on the eigenvalue problem arising in nuclear Configuration Interaction (CI) calculations, where a few extreme eigenpairs of a sparse symmetric matrix are needed. Here, we consider a block iterative eigensolver whose main computational kernels are the multiplication of a sparse matrix with multiple vectors (SpMM), and tall-skinny matrix operations. We then present techniques to significantly improve the SpMM and the transpose operation SpMM T by using the compressed sparse blocks (CSB) format. We achieve 3-4× speedup on the requisite operations over good implementations with the commonly used compressed sparse row (CSR) format. We develop a performance model that allows us to correctly estimate the performance of our SpMM kernel implementations, and we identify cache bandwidth as a potential performance bottleneck beyond DRAM. We also analyze and optimize the performance of LOBPCG kernels (inner product and linear combinations on multiple vectors) and show up to 15× speedup over using high performance BLAS libraries for these operations. The resulting high performance LOBPCG solver achieves 1.4× to 1.8× speedup over the existing Lanczos solver on a series of CI computations on high-end multicore architectures (Intel Xeons). We also analyze the performance of our techniques on an Intel Xeon Phi Knights Corner (KNC) processor.

97 MATHEMATICS AND COMPUTING↗

Communication: Configuration interaction combined with spin-projection for strongly correlated molecular electronic structures

We present single and double particle-hole excitations in the recently revived spin-projected Hartree-Fock. Our motivation is to treat static correlation with spin-projection and recover the residual correlation, mostly dynamic in nature, with simple configuration interaction (CI). To this end, we introduce the Wick theorem for nonorthogonal determinants, which enables an efficient implementation in conjunction with the direct CI scheme. The proposed approach, termed spin-extended CI with singles and doubles, achieves a balanced treatment between dynamic and static correlations. To approximately account for the quadruple excitations, we also modify the well-known Davidson correction. We report that our approaches yield surprisingly accurate potential curves for HF, H{sub 2}O, N{sub 2}, and a hydrogen lattice, compared to traditional single reference wave function methods at the same computational scaling as regular CI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-level coupled-cluster energetics by merging moment expansions with selected configuration interaction

Inspired by our earlier semi-stochastic work aimed at converging high-level coupled-cluster (CC) energetics, we propose a novel form of the CC(P; Q) theory in which the stochastic Quantum Monte Carlo propagations, used to identify dominant higher-than-doubly excited determinants, are replaced by the selected configuration interaction (CI) approach using the perturbative selection made iteratively (CIPSI) algorithm. The advantages of the resulting CIPSI-driven CC(P; Q) methodology are illustrated by a few molecular examples, including the dissociation of F 2 and the automerization of cyclobutadiene, where we recover the electronic energies corresponding to the CC calculations with a full treatment of singles, doubles, and triples based on the information extracted from compact CI wave functions originating from relatively inexpensive Hamiltonian diagonalizations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic excitation of molecules in solution calculated using the symmetry-adapted cluster–configuration interaction method in the polarizable continuum model

The effects from solvent environment are specific to the electronic states; therefore, a computational scheme for solvent effects consistent with the electronic states is necessary to discuss electronic excitation of molecules in solution. The PCM (polarizable continuum model) SAC (symmetry-adapted cluster) and SAC-CI (configuration interaction) methods are developed for such purposes. The PCM SAC-CI adopts the state-specific (SS) solvation scheme where solvent effects are self-consistently considered for every ground and excited states. For efficient computations of many excited states, we develop a perturbative approximation for the PCM SAC-CI method, which is called corrected linear response (cLR) scheme. Our test calculations show that the cLR PCM SAC-CI is a very good approximation of the SS PCM SAC-CI method for polar and nonpolar solvents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Unbiased reduced density matrices and electronic properties from full configuration interaction quantum Monte Carlo

Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the full configuration interaction quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself) and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Isotope shifts in beryllium-, boron-, carbon-, and nitrogen-like ions from relativistic configuration interaction calculations

Energy levels, normal and specific mass shift parameters as well as electronic densities at the nucleus are reported for numerous states along the beryllium, boron, carbon, and nitrogen isoelectronic sequences. Combined with nuclear data, these electronic parameters can be used to determine values of level and transition isotope shifts. The calculation of the electronic parameters is done using first-order perturbation theory with relativistic configuration interaction wavefunctions that account for valence, core–valence, and core–core correlation effects as zero-order functions. Results are compared with experimental and other theoretical values, when available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Annual Reports in Computational Chemistry, Chapter 3: GronOR Nonorthogonal Configuration Interaction Calculations at Exascale

This chapter describes the implementation and benchmarking of a new software code for nonorthogonal configuration interaction calculations. GronOR is named after the University of Groningen and Oak Ridge National Laboratory, the two institutions that are collaborating on the development of this code for the largest massively parallel supercomputers for open science. The basic implementation approach is described and preliminary benchmarking results are given as obtained from early access to Summit, the next Oak Ridge Leadership Computing Facility's supercomputer with 200 PFlop.

Chemistry↗

Relativistic K-LL Auger spectra in the intermediate-coupling scheme with configuration interaction

Theoretical K-LL Auger spectra from relativistic Dirac-Hartree-Slater calculations in intermediate coupling with configuration interaction (ICWCI) are considered. Calculated transition rates for 25 elements with Z between 18 and 96, inclusive, are listed and compare well with experimental data. Relativistic effects are found to be important above Z equal to about 35, and ICWCI is necessary to describe the spectra for Z less than about 60.

Chen, M. H.↗

Comparison of the quadratic configuration interaction and coupled cluster approaches to electron correlation including the effect of triple excitations

The recently proposed quadratic configuration interaction (QCI) method is compared with the more rigorous coupled cluster (CC) approach for a variety of chemical systems. Some of these systems are well represented by a single-determinant reference function and others are not. The finite order singles and doubles correlation energy, the perturbational triples correlation energy, and a recently devised diagnostic for estimating the importance of multireference effects are considered. The spectroscopic constants of CuH, the equilibrium structure of cis-(NO)2 and the binding energies of Be3, Be4, Mg3, and Mg4 were calculated using both approaches. The diagnostic for estimating multireference character clearly demonstrates that the QCI method becomes less satisfactory than the CC approach as non-dynamical correlation becomes more important, in agreement with a perturbational analysis of the two methods and the numerical estimates of the triple excitation energies they yield. The results for CuH show that the differences between the two methods become more apparent as the chemical systems under investigation becomes more multireference in nature and the QCI results consequently become less reliable. Nonetheless, when the system of interest is dominated by a single reference determinant both QCI and CC give very similar results.

Taylor, Peter R.↗

Chapter Three - GronOR Nonorthogonal Configuration Interaction Calculations at Exascale

This chapter describes the implementation and benchmarking of a new software code for nonorthogonal configuration interaction calculations. GronOR is named after the University of Groningen and Oak Ridge National Laboratory, the two institutions that are collaborating on the development of this code for the largest massively parallel supercomputers for open science. The basic implementation approach is described and preliminary benchmarking results are given as obtained from early access to Summit, the next Oak Ridge Leadership Computing Facility's supercomputer with 200 PFlop.

Straatsma, T.P.↗

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↗

Relativistic Configuration-Interaction and Perturbation Theory Calculations of the Sn XV Emission Spectrum

Recently, there has been increased interest in developing advanced bright sources for lithography. Sn ions are particularly promising due to their bright emission spectrum in the required wavelength range. Cowan’s code has been used to model the emission; however, it has adjustable parameters, which limit its predictive power, and it has limited relativistic treatment. Here, we present calculations based on ab initio relativistic configuration-interaction many-body perturbation theory (CI-MBPT), with relativistic corrections included at the Dirac-Fock level and core-polarization effects with the second-order MBPT. As a proof of principle that the theory is generally applicable to other Sn ions with proper development, we focused on one ion where direct comparison with experimental observations is possible. The theory can also be used for ions of other elements to predict emissions for optimization of plasma-based bright sources.

74 ATOMIC AND MOLECULAR PHYSICS↗

Analytic derivative couplings for spin-flip configuration interaction singles and spin-flip time-dependent density functional theory

We revisit the calculation of analytic derivative couplings for configuration interaction singles (CIS), and derive and implement these couplings for its spin-flip variant for the first time. Our algorithm is closely related to the CIS analytic energy gradient algorithm and should be straightforward to implement in any quantum chemistry code that has CIS analytic energy gradients. The additional cost of evaluating the derivative couplings is small in comparison to the cost of evaluating the gradients for the two electronic states in question. Incorporation of an exchange-correlation term provides an ad hoc extension of this formalism to time-dependent density functional theory within the Tamm-Dancoff approximation, without the need to invoke quadratic response theory or evaluate third derivatives of the exchange-correlation functional. Application to several different conical intersections in ethylene demonstrates that minimum-energy crossing points along conical seams can be located at substantially reduced cost when analytic derivative couplings are employed, as compared to use of a branching-plane updating algorithm that does not require these couplings. Application to H{sub 3} near its D{sub 3h} geometry demonstrates that correct topology is obtained in the vicinity of a conical intersection involving a degenerate ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the role of dynamic electron correlation in non-orthogonal configuration interaction with fragments

Two different approaches have been implemented to include the effect of dynamic electron correlation in the Non-Orthogonal Configuration Interaction for Fragments (NOCI-F) method. The first is based on shifting the diagonal matrix elements of the NOCI matrix, while the second incorporates the dynamic correlation explicitly in the fragment wave functions used to construct the many-electron basis functions of the NOCI. The two approaches are illustrated for the calculation of the electronic coupling relevant in singlet fission and the coupling of spin moments in organic radicals. Comparison of the calculated diabatic couplings, the NOCI energies and wave functions shows that dynamic electron correlation is not only efficiently but also effectively incorporated by the shifting approach and can largely affect the coupling between electronic states. Additionally, it brings the NOCI coupling of the spin moments in close agreement with benchmark calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗