A new electronic structure method for doublet states: configuration interaction in the space of ionized 1h and 2h1p determinants.

An implementation of gradient and energy calculations for configuration interaction variant of equation-of-motion coupled cluster with single and double substitutions for ionization potentials (EOM-IP-CCSD) is reported. The method (termed IP-CISD) treats the ground and excited doublet electronic states of an N-electron system as ionizing excitations from a closed-shell N+1-electron reference state. The method is naturally spin adapted, variational, and size intensive. The computational scaling is N(5), in contrast with the N(6) scaling of EOM-IP-CCSD. The performance and capabilities of the new approach are demonstrated by application to the uracil cation and water and benzene dimer cations by benchmarking IP-CISD against more accurate IP-CCSD. The equilibrium geometries, especially relative differences between different ionized states, are well reproduced. The average absolute errors and the standard deviations averaged for all bond lengths in all electronic states (58 values in total) are 0.014 and 0.007 A, respectively. IP-CISD systematically underestimates intramolecular distances and overestimates intermolecular ones, because of the underlying uncorrelated Hartree-Fock reference wave function. The IP-CISD excitation energies of the cations are of a semiquantitative value only, showing maximum errors of 0.35 eV relative to EOM-IP-CCSD. Trends in properties such as dipole moments, transition dipoles, and charge distributions are well reproduced by IP-CISD.

[1]  K. Hirao,et al.  Cluster expansion of the wave function. Electron correlations in the ground state, valence and Rydberg excited states, ionized states, and electron attached states of formaldehyde by SAC and SAC–CI theories , 1981 .

[2]  M. Ratner.  Molecular electronic-structure theory , 2000 .

[3]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

[4]  Henry F. Schaefer,et al.  On the evaluation of analytic energy derivatives for correlated wave functions , 1984 .

[5]  Anna I Krylov,et al.  Analytic gradients for the spin-conserving and spin-flipping equation-of-motion coupled-cluster models with single and double substitutions. , 2005, The Journal of chemical physics.

[6]  John F. Stanton,et al.  A simple scheme for the direct calculation of ionization potentials with coupled-cluster theory that exploits established excitation energy methods , 1999 .

[7]  Hiroshi Nakatsuji,et al.  Description of two- and many-electron processes by the SAC-CI method , 1991 .

[8]  R. Bartlett,et al.  Transformation of the Hamiltonian in excitation energy calculations: Comparison between Fock‐space multireference coupled‐cluster and equation‐of‐motion coupled‐cluster methods , 1991 .

[9]  J. VandeVondele,et al.  Electronic structure of the water dimer cation. , 2008, The journal of physical chemistry. A.

[10]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[11]  R. Chaudhuri,et al.  The eigenvalue-independent partitioning technique in Fock space: an alternative route to open-shell coupled-cluster theory for incomplete model spaces , 1989 .

[12]  P. Pieniazek,et al.  Benchmark full configuration interaction and equation-of-motion coupled-cluster model with single and double substitutions for ionized systems results for prototypical charge transfer systems: noncovalent ionized dimers. , 2007, The Journal of chemical physics.

[13]  Gradients for the partitioned equation-of-motion coupled-cluster method , 1999 .

[14]  Anna I Krylov,et al.  Double spin-flip approach within equation-of-motion coupled cluster and configuration interaction formalisms: Theory, implementation, and examples. , 2009, The Journal of chemical physics.

[15]  Anna I Krylov,et al.  Charge localization and Jahn-Teller distortions in the benzene dimer cation. , 2008, The Journal of chemical physics.

[16]  Gregory S. Tschumper,et al.  Real versus artifactual symmetry-breaking effects in Hartree-Fock, density-functional, and coupled-cluster methods. , 2004, The Journal of chemical physics.

[17]  A. Krylov,et al.  The effect of pi-stacking and H-bonding on ionization energies of a nucleobase: uracil dimer cation. , 2009, Physical chemistry chemical physics : PCCP.

[18]  P. Pieniazek,et al.  Electronic structure of the benzene dimer cation. , 2007, The Journal of chemical physics.

[19]  D. Mukherjee,et al.  A note on the direct calculation of excitation energies by quasi-degenerate MBPT and coupled-cluster theory , 1986 .

[20]  Filipp Furche,et al.  Adiabatic time-dependent density functional methods for excited state properties , 2002 .

[21]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

[22]  E. Davidson,et al.  Symmetry breaking in polyatomic molecules: real and artifactual , 1983 .

[23]  Anna I. Krylov,et al.  A spin-complete version of the spin-flip approach to bond breaking: What is the impact of obtaining spin eigenfunctions? , 2003 .

[24]  D. Mukherjee,et al.  Application of cluster expansion techniques to open shells: Calculation of difference energies , 1984 .

[25]  Mihály Kállay,et al.  Equation-of-motion coupled-cluster methods for ionized states with an approximate treatment of triple excitations. , 2005, The Journal of chemical physics.

[26]  R. Bartlett,et al.  Multireference coupled-cluster methods using an incomplete model space: Application to ionization potentials and excitation energies of formaldehyde , 1987 .

[27]  Rodney J. Bartlett,et al.  Equation of motion coupled cluster method for electron attachment , 1995 .

[28]  P. Piecuch,et al.  Active-space symmetry-adapted-cluster configuration-interaction and equation-of-motion coupled-cluster methods for high accuracy calculations of potential energy surfaces of radicals. , 2007, The Journal of chemical physics.

[29]  S. Hirata,et al.  Higher-order equation-of-motion coupled-cluster methods for ionization processes. , 2006, The Journal of chemical physics.

[30]  David Casanova,et al.  The spin-flip extended single excitation configuration interaction method. , 2008, The Journal of chemical physics.

[31]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[32]  J. Simons,et al.  Theory of electron affinities of small molecules , 1973 .

[33]  P. Szalay,et al.  Analytic energy derivatives for coupled‐cluster methods describing excited states: General formulas and comparison of computational costs , 1995 .

[34]  So Hirata,et al.  High-order determinantal equation-of-motion coupled-cluster calculations for ionized and electron-attached states , 2000 .

[35]  Kimihiko Hirao,et al.  Cluster expansion of the wavefunction. Symmetry-adapted-cluster expansion, its variational determination, and extension of open-shell orbital theory , 1978 .

[36]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[37]  J. Gauss,et al.  Analytic energy derivatives for ionized states described by the equation‐of‐motion coupled cluster method , 1994 .

[38]  Hans-Joachim Werner,et al.  Analytical energy gradients for internally contracted second-order multireference perturbation theory , 2003 .

[39]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[40]  Anna I. Krylov,et al.  Spin-flip configuration interaction: an electronic structure model that is both variational and size-consistent , 2001 .

[41]  R. Bartlett,et al.  Equation-of-motion coupled cluster method with full inclusion of connected triple excitations for electron-attached states: EA-EOM-CCSDT , 2003 .

[42]  Rodney J. Bartlett,et al.  Equation-of-motion coupled cluster method with full inclusion of the connected triple excitations for ionized states: IP-EOM-CCSDT , 2003 .