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 91 records · Page 5

Investigation of the ozone formation reaction pathway: Comparisons of full configuration interaction quantum Monte Carlo and fixed-node diffusion Monte Carlo with contracted and uncontracted MRCI

The association/dissociation reaction path for ozone (O 2 + O ↔ O 3 ) is notoriously difficult to describe accurately using ab initio electronic structure theory, due to the importance of both strong and dynamic electron correlations. Experimentally, spectroscopic studies of the highest lying recorded vibrational states combined with the observed negative temperature dependence of the kinetics of oxygen isotope exchange reactions confirm that the reaction is barrierless, consistent with the latest potential energy surfaces. Previously reported potentials based on Davidson-corrected internally contracted multireference configuration interaction (MRCI) suffer from a spurious reef feature in the entrance channel even when extrapolated towards the complete basis set limit. In this work, we report an analysis of comparisons between a variety of electronic structure methods including internally contracted and uncontracted MRCI (with and without Davidson corrections), as well as full configuration interaction quantum Monte Carlo, fixed-node diffusion Monte Carlo, and density matrix renormalization group.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-Orthogonal Configuration Interaction for Fragments

An overview is given of the non-orthogonal configuration interaction (NOCI) implementation in the open-source program GronOR. After pointing out the specifics of NOCI and its extension to ensembles of molecules NOCI-Fragments (NOCI-F), the theoretical foundations of the approaches are described. The introduction of a common molecular orbital basis, the selection of the most contributing determinants pairs, and a massively parallel and GPU-accelerated implementation have made possible the application of NOCI to medium-sized systems as illustrated in a short overview of some of the recent applications.

de Graaf, C.↗

Development of Large-Dimension Configuration-Interaction Shell-Model Code

The project period was devoted to several developments in the technical capabilities of the BIGSTIC large-dimension configuration-interaction shell-model code, written in Fortran 90. The specific computational goals for the project period were: (1) store Lanczos vectors on core in RAM to minimize I/O; (2) rewrite reorthogonalization with Lanczos vectors stored in core, consult with personnel at LLNL, LBL, ORNL, Iowa State University to maximize performance; (3) restrict creation of N-body jumps to those needed by an individual node; and (4) distribute 3-body interaction over many cores. Significant progress was made towards these goals, especially (1) and (2), although in the process they discovered intermediate tasks that had to be accomplished first. The achievements were as follows - I put into place structures and algorithms to facility fragmenting very large-dimension Lanczos intermediate vectors. Only by fragmenting the vectors can we carry out (1) and (2). In addition, I reorganized the action of the Hamiltonian matrix and created a new division of operations for MPI. Based upon earlier work, I made plans of a revised algorithm for distribution of work with MPI, with a particular eye towards breaking up the Lanczos vectors. I introduce a new derived type (opbundles) which collects the parameters for the Hamiltonian, and rewrote the application routines to use it. It has been validated and verified. I made progress towards revised MPI parallelization. Using the opbundles, I was able to compute a distribution of work over compute nodes, which should be very efficient. This new distribution is easier to derive and more efficient, in principle, than the old distribution. Furthermore, it should make applications with fragmented Lanczos vectors easier. Implementation is still in progress.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electronic spectra of DyF studied by four-component relativistic configuration interaction methods

The electronic states of the DyF molecule below 3.0 eV are studied using 4-component relativistic CI methods. Spinors generated by the average-of-configuration Hartree-Fock method with the Dirac-Coulomb Hamiltonian were used in CI calculations by the KRCI (Kramers-restricted configuration interaction) program. The CI reference space was generated by distributing 11 electrons among the 11 Kramers pairs composed mainly of Dy [4f], [6s], [6p] atomic spinors, and double excitations are allowed from this space to the virtual molecular spinors. The CI calculations indicate that the ground state has the dominant configuration (4f{sup 9})(6s{sup 2})(Ω = 7.5). Above this ground state, 4 low-lying excited states (Ω = 8.5, 7.5, 7.5, 7.5) are found with dominant configurations (4f{sup 10})(6s). These results are consistent with the experimental studies of McCarthy et al. Above these 5 states, 2 states were observed at T{sub 0} = 2.39 eV, 2.52 eV by McCarthy et al. and were named as [19.3]8.5 and [20.3]8.5. McCarthy et al. proposed that both states have dominant configurations (4f{sup 9})(6s)(6p), but these configurations are not consistent with the large R{sub e}’s (∼3.9 a.u.) estimated from the observed rotational constants. The present CI calculations provide near-degenerate states of (4f{sup 10})(6p{sub 3/2,1/2}), (4f{sup 10})(6p{sub 3/2,3/2}), and (4f{sup 9})(6s)(6p{sub 3/2,1/2}) at around 3 eV. The former two states have larger R{sub e} (3.88 a.u.) than the third, so that it is reasonable to assign (4f{sup 10})(6p{sub 3/2,1/2}) to [19.3]8.5 and (4f{sup 10})(6p{sub 3/2,3/2}) to [20.3]8.5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic excitation spectra of molecules in solution calculated using the symmetry-adapted cluster-configuration interaction method in the polarizable continuum model with perturbative approach

A perturbative approximation of the state specific polarizable continuum model (PCM) symmetry-adapted cluster-configuration interaction (SAC-CI) method is proposed for efficient calculations of the electronic excitations and absorption spectra of molecules in solutions. This first-order PCM SAC-CI method considers the solvent effects on the energies of excited states up to the first-order with using the zeroth-order wavefunctions. This method can avoid the costly iterative procedure of the self-consistent reaction field calculations. The first-order PCM SAC-CI calculations well reproduce the results obtained by the iterative method for various types of excitations of molecules in polar and nonpolar solvents. The first-order contribution is significant for the excitation energies. The results obtained by the zeroth-order PCM SAC-CI, which considers the fixed ground-state reaction field for the excited-state calculations, are deviated from the results by the iterative method about 0.1 eV, and the zeroth-order PCM SAC-CI cannot predict even the direction of solvent shifts in n-hexane for many cases. The first-order PCM SAC-CI is applied to studying the solvatochromisms of (2,2{sup ′}-bipyridine)tetracarbonyltungsten [W(CO){sub 4}(bpy), bpy = 2,2{sup ′}-bipyridine] and bis(pentacarbonyltungsten)pyrazine [(OC){sub 5}W(pyz)W(CO){sub 5}, pyz = pyrazine]. The SAC-CI calculations reveal the detailed character of the excited states and the mechanisms of solvent shifts. The energies of metal to ligand charge transfer states are significantly sensitive to solvents. The first-order PCM SAC-CI well reproduces the observed absorption spectra of the tungsten carbonyl complexes in several solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GronOR: Scalable and Accelerated Nonorthogonal Configuration Interaction for Molecular Fragment Wave Functions

GronOR is a program package for nonorthogonal configuration interaction calculations. Electronic wave functions are constructed in terms of antisymmetrized products of multiconfiguration molecular fragment wave functions. The computational complexity of the nonorthogonal methodologies implemented in GronOR applied to large molecular assemblies requires a design that takes full advantage of massively parallel supercomputer architectures and accelerator technologies. This work describes the implementation strategy and resulting performance characteristics. In addition to parallelization and acceleration, the software development strategy includes aspects of fault resiliency and heterogeneous computing. The program was designed for large-scale supercomputers but also runs effectively on small clusters and workstations for small molecular systems. GronOR is available as open source to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Communication: Approaching exact quantum chemistry by cluster analysis of full configuration interaction quantum Monte Carlo wave functions

We propose to accelerate convergence toward full configuration interaction (FCI) energetics by using the coupled-cluster approach, in which singly and doubly excited clusters, needed to determine the energy, are iterated in the presence of their three- and four-body counterparts extracted from FCI quantum Monte Carlo (FCIQMC) propagations. Finally, preliminary calculations for the water molecule at the equilibrium and stretched geometries show that we can accurately extrapolate the FCI energetics based on the early stages of FCIQMC propagations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flexible configuration-interaction shell-model many-body solver

BIGSTICK Is a flexible configuration-Interaction open-source shell-model code for the many-fermion problem In a shell model (occupation representation) framework. BIGSTICK can generate energy spectra, static and transition one-body densities, and expectation values of scalar operators. Using the built-in Lanczos algorithm one can compute transition probabflity distributions and decompose wave functions into components defined by group theory.

Johnson, CalvinW.↗

Full configuration-interaction calculations with the simplest iterative complement method: Merit of the inverse Hamiltonian

The simplest iterative complement (SIC) calculations starting from Hartree-Fock and giving full configuration interaction (CI) at convergence were performed using regular and inverse Hamiltonians. Each iteration step is variational and involves only one variable. The convergence was slow when we used the regular Hamiltonian, but became very fast when we used the inverse Hamiltonian. This difference is due to the Coulomb singularity problem inherent in the regular Hamiltonian; the inverse Hamiltonian does not have such a problem. For this reason, the merit of the inverse Hamiltonian over the regular one becomes even more dramatic when we use a better-quality basis set. This was seen by comparing the calculations due to the minimal and double-{zeta} basis sets. Similar problematic situations exist in the Krylov sequence and in the Lanczos and Arnoldi methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite-temperature vibrational full configuration interaction

Thermodynamic functions of an ideal molecular gas due to its anharmonic vibrations are evaluated in a wide range of temperature (T) by the vibrational full-configuration-interaction (FCI) method using a quartic force field and a finite number (N) of harmonic-oscillator basis functions along each normal mode. The thermodynamic functions considered are the grand potential (Ω), internal energy (U), and entropy (S). They are compared with those obtained from the Bose–Einstein theory with or without truncation of the harmonic-oscillator basis functions after quantum number N–1. The comparison reveals that the finite-basis-set errors in Ω and U are, respectively, k B Tln(k B T/Nℏω) and k B T per mode in the high-T limit, obscuring anharmonic effects when k B T > ℏω, where ω is the lowest mode frequency. Here, the benchmark data for several low-order perturbation corrections to Ω, U, and S are also obtained as the numerical derivatives of their FCI values with respect to dimensionless perturbation strength, and the domain of T and N in which these data are reliable (for the N → ∞ limits) is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Communication: An adaptive configuration interaction approach for strongly correlated electrons with tunable accuracy

We introduce a new procedure for iterative selection of determinant spaces capable of describing highly correlated systems. This adaptive configuration interaction (ACI) determines an optimal basis by an iterative procedure in which the determinant space is expanded and coarse grained until self-consistency. Two importance criteria control the selection process and tune the ACI to a user-defined level of accuracy. The ACI is shown to yield potential energy curves of N{sub 2} with nearly constant errors, and it predicts singlet-triplet splittings of acenes up to decacene that are in good agreement with the density matrix renormalization group.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin-flip configuration interaction singles with exact spin-projection: Theory and applications to strongly correlated systems

Spin-flip approaches capture static correlation with the same computational scaling as the ordinary single reference methods. Here, we extend spin-flip configuration interaction singles (SFCIS) by projecting out intrinsic spin-contamination to make it spin-complete, rather than by explicitly complementing it with spin-coupled configurations. We give a general formalism of spin-projection for SFCIS, applicable to any spin states. The proposed method is viewed as a natural unification of SFCIS and spin-projected CIS to achieve a better qualitative accuracy at a low computational cost. While our wave function ansatz is more compact than previously proposed spin-complete SF approaches, it successfully offers more general static correlation beyond biradicals without sacrificing good quantum numbers. It is also shown that our method is invariant with respect to open-shell orbital rotations, due to the uniqueness of spin-projection. We will report benchmark calculations to demonstrate its qualitative performance on strongly correlated systems, including conical intersections that appear both in ground-excited and excited-excited degeneracies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Positron-attachment to small molecules: Vibrational enhancement of positron affinities with configuration interaction level of multi-component molecular orbital approach

To theoretically demonstrate the binding of a positron to small polarized molecules, we have calculated the vibrational averaged positron affinity (PA) values along the local vibrational contribution with the configuration interaction level of multi-component molecular orbital method. This method can take the electron-positron correlation contribution into account through single electronic - single positronic excitation configurations. The PA values are enhanced by including the local vibrational contribution from vertical PA values due to the anharmonicity of the potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cluster decomposition of full configuration interaction wave functions: A tool for chemical interpretation of systems with strong correlation

Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size, thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation amplitudes contain information on the complexity of the electronic wave function, but this information is contaminated by contributions from disconnected excitations, i.e., those excitations that are just products of independent lower-level excitations. The unwanted contributions can be removed via a cluster decomposition procedure, making it possible to examine the importance of connected excitations in complicated multireference molecules which are outside the reach of conventional algorithms. We present an implementation of the cluster decomposition analysis and apply it to both true FCI wave functions, as well as wave functions generated from the adaptive sampling CI algorithm. Furthermore, the cluster decomposition is useful for interpreting calculations in chemical studies, as a diagnostic for the convergence of various excitation manifolds, as well as as a guidepost for polynomially scaling electronic structure models. Applications are presented for (i) the double dissociation of water, (ii) the carbon dimer, (iii) the π space of polyacenes, and (iv) the chromium dimer. While the cluster amplitudes exhibit rapid decay with an increasing rank for the first three systems, even connected octuple excitations still appear important in Cr 2 , suggesting that spin-restricted single-reference coupled-cluster approaches may not be tractable for some problems in transition metal chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗