Search NASA⌕ Search

SEARCH · Search NASA

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 109 records · Page 6

Toward Full Configuration Interaction for Transition-Metal Complexes

In this work, an efficacious approximation to full configuration interaction (FCI) is adapted to calculate singlet-triplet gaps for transition-metal complexes. This strategy, incremental FCI (iFCI), uses a many-body expansion to systematically add correlation to a simple reference wave function and therefore achieves greatly reduced computational costs compared to FCI. iFCI through the 3-body expansion is demonstrated on four model transition-metal complexes involving the metals Zn, V, and Cu. Screening techniques to increase the computational efficiency of iFCI are proposed and tested, showing reduction in the number of 3-body terms by more than 90% with controlled errors. The largest complex treated by iFCI has 142 valence electrons, all of which are correlated among the full set of 444 active orbitals. Computed spin gaps approach experimental results for the four complexes, though room for improvement remains.

74 ATOMIC AND MOLECULAR PHYSICS↗

Systematic calculations of energy levels and transition rates of BE-like ions with Z = 10–30 using a combined configuration interaction and many-body perturbation theory approach

We report calculations of energy levels and radiative rates for transitions among the lowest 116 fine-structure levels arising from the n⩽5 configurations in Be-like ions with Z = 10–30. The wavelengths, oscillator strengths, line strengths, and radiative rates for all possible electric dipole, magnetic dipole, electric quadrupole, and magnetic quadrupole transitions among the 116 levels have been calculated using the combined configuration interaction and many-body perturbation method. The accuracy of the results is determined through extensive comparisons with existing laboratory measurements and theoretical results. The present complete set of results should be of great help in line identification and the interpretation of spectra, as well as in the modeling and diagnostics of astrophysical and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hyperfine structures and Landé g{sub J}-factors for n=2 states in beryllium-, boron-, carbon-, and nitrogen-like ions from relativistic configuration interaction calculations

Energy levels, hyperfine interaction constants, and Landé g{sub J}-factors are reported for n=2 states in beryllium-, boron-, carbon-, and nitrogen-like ions from relativistic configuration interaction calculations. Valence, core–valence, and core–core correlation effects are taken into account through single and double-excitations from multireference expansions to increasing sets of active orbitals. A systematic comparison of the calculated hyperfine interaction constants is made with values from the available literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Trajectory-guided configuration interaction simulations of multidimensional quantum dynamics

We propose an approach to modelling multidimensional quantum systems which uses direct-dynamics trajectories to guide wavefunction propagation. First, trajectory simulations are used to generate a sample of dynamically relevant configurations on the potential energy surface (PES). Second, the sampled configurations are used to construct an n-mode representation of the PES using a greedy algorithm. Finally, the time-dependent Schroedinger equation is solved using a configuration interaction expansion of the wavefunction, with individual basis functions derived directly from the 1-mode contributions to the n-mode PES. This approach is successfully demonstrated by application to a 20-dimensional benchmark problem describing tunnelling in the presence of coupled degrees of freedom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin-flip, tensor equation-of-motion configuration interaction with a density-functional correction: A spin-complete method for exploring excited-state potential energy surfaces

We revisit the formalism of the spin-adapted, spin-flip (SA-SF) configuration-interaction singles (CIS) method based on a tensor equation-of-motion formalism that affords proper spin eigenstates without sacrificing single-reference simplicity. Matrix elements for SA-SF-CIS are then modified in a manner similar to collinear spin-flip time-dependent density functional theory (SF-TDDFT), to include a DFT exchange-correlation correction. The performance of this method, which we call SA-SF-DFT, is evaluated numerically and we find that it systematically improves the energies of electronic states that exhibit significant spin contamination within the conventional SF-TDDFT approach. The new method cures the state assignment problem that plagues geometry optimizations and ab initio molecular dynamics simulations using traditional SF-TDDFT, without sacrificing computational efficiency, and furthermore provides correct topology at conical intersections, including those that involve the ground state, unlike conventional TDDFT. As such, SA-SF-DFT appears to be a promising method for generating excited-state potential energy surfaces at DFT cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fully Quantum Embedding with Density Functional Theory for Full Configuration Interaction Quantum Monte Carlo

We here develop a fully quantum embedded version of initiator full configuration interaction quantum Monte Carlo (i–FCIQMC) and apply it to study an ionic bond (lithium hydride, LiH) and a covalent bond (hydrogen flouride, HF) physisorbed to a benzene molecule. The embedding is performed using a recently developed Huzinaga projection operator approach, which affords good synergy with i–FCIQMC by minimizing the number of orbitals in the calculation. When considering the dissociation energy of these bonds into closed-shell ionic fragments, we find that i–FCIQMC embedded in density functional theory (i–FCIQMC-in-DFT) delivers comparable accuracy with coupled cluster singles and doubles with perturbative triples embedded in density functional theory (CCSD(T)-in-DFT). In treating the bond dissociation energy curve of HF, i–FCIQMC-in-DFT has improved accuracy over CCSD(T)-in-DFT due to the presence of strong correlation. Lastly, we discuss the implications of the new i–FCIQMC-in-DFT method as applied to bond breaking in catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The role of Rydberg and continuum levels in computing high harmonic generation spectra of the hydrogen atom using time-dependent configuration interaction

We study the role of Rydberg bound-states and continuum levels in the field-induced electronic dynamics associated with the High-Harmonic Generation (HHG) spectroscopy of the hydrogen atom. Time-dependent configuration-interaction (TD-CI) is used with very large atomic orbital (AO) expansions (up to L= 4 with sextuple augmentation and off-center functions) to describe the bound Rydberg levels, and some continuum levels. To address the lack of ionization losses in TD-CI with finite AO basis sets, we employed a heuristic lifetime for energy levels above the ionization potential. The heuristic lifetime model is compared against the conventional atomic orbital treatment (infinite lifetimes), and a third approximation which is TD-CI using only the bound levels (continuum lifetimes go to zero). The results suggest that spectra calculated using conventional TD-CI do not converge with increasing AO basis set size, while the zero lifetime and heuristic lifetime models converge to qualitatively similar spectra, with implications for how best to apply bound state electronic structure methods to simulate HHG. The origin of HHG spectral features including the cutoff and extent of interference between peaks is uncovered by separating field-induced coupling between different types of levels (ground state, bound Rydberg levels, and continuum) in the simulated electronic dynamics. Thus the origin of deviations between the predictions of the semi-classical three step model and the full simulation can be associated with particular physical contributions, which helps to explain both the successes and the limitations of the three step model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fe{sup 15+} dielectronic recombination and the effects of configuration interaction between resonances with different captured electron principal quantum numbers

Dielectronic recombination (DR) of Na-like Fe{sup 15+} forming Mg-like Fe{sup 14+} via excitation of a 2l core electron has been investigated. We find that configuration interaction (CI) between DR resonances with different captured electron principal quantum numbers n can lead to a significant reduction in resonance strengths for n{>=}5. Previous theoretical work for this system has not considered this form of CI. Including it accounts for most of the discrepancy between previous theoretical and experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Configuration-interaction many-body-perturbation-theory energy levels of four-valent Si i

The mixed configuration-interaction (CI) many-body-perturbation-theory method is accurate in divalent atoms. In more complex atoms, with the number of valence electrons it becomes progressively more difficult to saturate CI space. Here in this study, a four-valence electron atom, Si i, is considered. It is found that by using a relatively small cavity of 30 a.u. and by choosing carefully configuration space, it is possible to obtain quite accurate agreement between the theory and experiment. After subtraction of systematic shifts of 481 and -426 cm -1 for the lowest even and odd states, respectively, the deviation between theory and experiment becomes at the level of 100 cm -1 . Finally, this agreement is comparable to that in divalent atoms where the CI saturation has been achieved. It is anticipated that the approach can also give good results for atoms with more valence electrons to be considered in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpretation of the photoelectron, ultraviolet, and vacuum ultraviolet photoabsorption spectra of bromobenzene by ab initio configuration interaction and DFT computations

New photoelectron, ultraviolet (UV), and vacuum UV (VUV) spectra have been obtained for bromobenzene by synchrotron study with higher sensitivity and resolution than previous work. This, together with use of ab initio calculations with both configuration interaction and time dependent density functional theoretical methods, has led to major advances in interpretation. The VUV spectrum has led to identification of a considerable number of Rydberg states for the first time. The Franck-Condon (FC) analyses including both hot and cold bands lead to identification of the vibrational structure of both ionic and electronically excited states including two Rydberg states. The UV onset has been interpreted in some detail, and an interpretation based on the superposition of FC and Herzberg-Teller contributions has been performed. In a similar way, the 6 eV absorption band which is poorly resolved is analysed in terms of the presence of two ππ* states of {sup 1}A{sub 1} (higher oscillator strength) and {sup 1}B{sub 2} (lower oscillator strength) symmetries, respectively. The detailed analysis of the vibrational structure of the 2{sup 2}B{sub 1} ionic state is particularly challenging, and the best interpretation is based on equation-of-motion-coupled cluster with singles and doubles computations. A number of equilibrium structures of the ionic and singlet excited states show that the molecular structure is less subject to variation than corresponding studies for iodobenzene. The equilibrium structures of the 3b{sub 1}3s and 6b{sub 2}3s (valence shell numbering) Rydberg states have been obtained and compared with the corresponding ionic limit structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗