Configuration interaction wave functions: a seniority number approach

This work deals with the configuration interaction method when an N-electron Hamiltonian is projected on Slater determinants which are classified according to their seniority number values. We study the spin features of the wave functions and the size of the matrices required to formulate states of any spin symmetry within this treatment. Correlation energies associated with the wave functions arising from the seniority-based configuration interaction procedure are determined for three types of molecular orbital basis: canonical molecular orbitals, natural orbitals, and the orbitals resulting from minimizing the expectation value of the N-electron seniority number operator. The performance of these bases is analyzed by means of numerical results obtained from selected N-electron systems of several spin symmetries. The comparison of the results highlights the efficiency of the molecular orbital basis which minimizes the mean value of the seniority number for a state, yielding energy values closer to those provided by the full configuration interaction procedure.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada

Seniority number in spin-adapted spaces and compactness of configuration interaction wave functions

This work extends the concept of seniority number, which has been widely used for classifying N-electron Slater determinants, to wave functions of N electrons and spin S, as well as to N-electron spin-adapted Hilbert spaces. We propose a spin-free formulation of the seniority number operator and perform a study on the behavior of the expectation values of this operator under transformations of the molecular basis sets. This study leads to propose a quantitative evaluation for the convergence of the expansions of the wave functions in terms of Slater determinants. The non-invariant character of the seniority number operator expectation value of a wave function with respect to a unitary transformation of the molecular orbital basis set, allows us to search for a change of basis which minimizes that expectation value. The results found in the description of wave functions of selected atoms and molecules show that the expansions expressed in these bases exhibit a more rapid convergence than those formulated in the canonical molecular orbital bases and even in the natural orbital ones.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada

Direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions: The influence of three-index N -representability conditions

This work proposes the variational determination of two-electron reduced density matrices corresponding to the ground state of N-electron systems within the doubly occupied-configuration-interaction methodology. The P, Q, and G two-index N-representability conditions have been extended to the T1 and T2 (T2′) three-index ones and the resulting optimization problem has been addressed using a standard semidefinite program. We report results obtained from the doubly occupied-configuration-interaction method, from the two-index constraint variational procedure and from the two- and three-index constraint variational treatment. The discussion of these results along with a study of the computational cost demanded shows the usefulness of our proposal.Fil: Alcoba, Diego Ricardo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Torre, Alicia. Universidad del País Vasco; EspañaFil: Lain, Luis. Universidad del País Vasco; EspañaFil: Massaccesi, Gustavo Ernesto. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Universidad de Buenos Aires; ArgentinaFil: Oña, Ofelia Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas; ArgentinaFil: Honore, Eduardo Misael. Universidad de Buenos Aires; ArgentinaFil: Poelmans, Ward. University of Ghent; BélgicaFil: Van Neck, Dimitri. University of Ghent; BélgicaFil: Bultinck, Patrick. University of Ghent; BélgicaFil: De Baerdemacker, Stijn. University of Ghent; Bélgic

Spin contamination-free N -electron wave functions in the excitation-based configuration interaction treatment

This work deals with the spin contamination in N-electron wave functions provided by the excitation-based configuration interaction methods. We propose a procedure to ensure a suitable selection of excited N-electron Slater determinants with respect to a given reference determinant, required in these schemes. The procedure guarantees the construction of N-electron wave functions which are eigenfunctions of the spin-squared operator S 2, avoiding any spin contamination. Our treatment is based on the evaluation of the excitation level of the determinants by means of the expectation value of an excitation operator formulated in terms of spin-free replacement operators. We report numerical determinations of energies and S 2 expectation values, arising from our proposal as well as from traditional configuration interaction methods, in selected open-shell systems, in order to compare the behavior of these procedures and their computational costs.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada

Spin contamination-free N -electron wave functions in the excitation-based configuration interaction treatment

This work deals with the spin contamination in N-electron wave functions provided by the excitation-based configuration interaction methods. We propose a procedure to ensure a suitable selection of excited N-electron Slater determinants with respect to a given reference determinant, required in these schemes. The procedure guarantees the construction of N-electron wave functions which are eigenfunctions of the spin-squared operator S 2, avoiding any spin contamination. Our treatment is based on the evaluation of the excitation level of the determinants by means of the expectation value of an excitation operator formulated in terms of spin-free replacement operators. We report numerical determinations of energies and S 2 expectation values, arising from our proposal as well as from traditional configuration interaction methods, in selected open-shell systems, in order to compare the behavior of these procedures and their computational costs.Fil: Alcoba, Diego Ricardo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Torre, Alicia. Universidad del País Vasco; EspañaFil: Lain, Luis. Universidad del País Vasco; EspañaFil: Massaccesi, Gustavo Ernesto. Universidad de Buenos Aires. Ciclo Basico Comun. Departamento de Física, Química y Matemática; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Oña, Ofelia Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas; ArgentinaFil: Capuzzi, Pablo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentin

Polynomial scaling approximations and dynamic correlation corrections to doubly occupied configuration interaction wave functions

A class of polynomial scaling methods that approximate Doubly Occupied Configuration Interaction (DOCI) wave functions and improve the description of dynamic correlation is introduced. The accuracy of the resulting wave functions is analysed by comparing energies and studying the overlap between the newly developed methods and full configuration interaction wave functions, showing that a low energy does not necessarily entail a good approximation of the exact wave function. Due to the dependence of DOCI wave functions on the single-particle basis chosen, several orbital optimisation algorithms are introduced. An energy-based algorithm using the simulated annealing method is used as a benchmark. As a computationally more affordable alternative, a seniority number minimising algorithm is developed and compared to the energy based one revealing that the seniority minimising orbital set performs well. Given a well-chosen orbital basis, it is shown that the newly developed DOCI based wave functions are especially suitable for the computationally efficient description of static correlation and to lesser extent dynamic correlation.Fil: Van Raemdonck, Mario. Ghent University; BélgicaFil: Alcoba, Diego Ricardo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Poelmans, Ward. Ghent University; BélgicaFil: De Baerdemacker, Stijn. Ghent University; BélgicaFil: Torre, Alicia. Universidad del País Vasco; EspañaFil: Lain, Luis. Universidad del País Vasco; EspañaFil: Massaccesi, Gustavo Ernesto. Universidad de Barcelona. Facultad de Física. Departamento de Física Fomental; EspañaFil: Van Neck, D.. Ghent University; BélgicaFil: Bultinck, P.. Ghent University; Bélgic

Spin contamination-free N -electron wave functions in the excitation-based configuration interaction treatment

This work deals with the spin contamination in N-electron wave functions provided by the excitation-based configuration interaction methods. We propose a procedure to ensure a suitable selection of excited N-electron Slater determinants with respect to a given reference determinant, required in these schemes. The procedure guarantees the construction of N-electron wave functions which are eigenfunctions of the spin-squared operator S 2, avoiding any spin contamination. Our treatment is based on the evaluation of the excitation level of the determinants by means of the expectation value of an excitation operator formulated in terms of spin-free replacement operators. We report numerical determinations of energies and S 2 expectation values, arising from our proposal as well as from traditional configuration interaction methods, in selected open-shell systems, in order to compare the behavior of these procedures and their computational costs.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada

Picture of Bonding in Protonated Phosphorus and Arsenic Clusters. Detection of 3-Center Bonds from the Generalized Population Analysis

In this report we provide a theoretical support for the existence of nonclassical 3-center bonding recently invoked to explain the structure of protonated phosphorus and arsenic X4H+ clusters.For this purpose the electron structure of these clusters was analysed using the formalism of the generalized population analysis, in terms of which the eventual presence of multicenter bonds in a mo-lecule can be directly detected and localized. Our results confirm the presence of 3-center bonding involving X−H−X fragments in these clusters

Unrestricted treatment for the direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions

This work extends to the unrestricted orbital approach the procedure described in our previous report [Alcoba et al., J. Chem. Phys. 148, 024105 (2018)] for determining variationally the two-electron reduced density matrix arising from doubly occupied-configuration-interaction wave functions by imposing two- and three-index N-representability conditions. An analysis of the numerical results obtained in selected systems, from both restricted and unrestricted treatments, allows one to assess the performance of these methodologies as well as to show the influence of the P, Q, G, T1, and T2 positivity conditions. We highlight the satisfactory results obtained within the unrestricted scheme.Fil: Alcoba, Diego Ricardo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; ArgentinaFil: Torre, Alicia. Universidad del País Vasco; EspañaFil: Lain, Luis. Universidad del País Vasco; EspañaFil: Massaccesi, Gustavo Ernesto. Universidad de Buenos Aires. Ciclo Basico Comun. Departamento de Física, Química y Matemática; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Oña, Ofelia Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas; ArgentinaFil: Ríos, Elías Daniel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas; Argentin