Numerically accurate linear response-properties in the configuration-interaction singles (CIS) approximation.

In the present work, we report an efficient implementation of configuration interaction singles (CIS) excitation energies and oscillator strengths using the multi-resolution analysis (MRA) framework to address the basis-set convergence of excited state computations. In MRA (ground-state) orbitals, excited states are constructed adaptively guaranteeing an overall precision. Thus not only valence but also, in particular, low-lying Rydberg states can be computed with consistent quality at the basis set limit a priori, or without special treatments, which is demonstrated using a small test set of organic molecules, basis sets, and states. We find that the new implementation of MRA-CIS excitation energy calculations is competitive with conventional LCAO calculations when the basis-set limit of medium-sized molecules is sought, which requires large, diffuse basis sets. This becomes particularly important if accurate calculations of molecular electronic absorption spectra with respect to basis-set incompleteness are required, in which both valence as well as Rydberg excitations can contribute to the molecule's UV/VIS fingerprint.

[1]  Hans-Joachim Werner,et al.  Systematically convergent basis sets for explicitly correlated wavefunctions: the atoms H, He, B-Ne, and Al-Ar. , 2008, The Journal of chemical physics.

[2]  R. Ahlrichs,et al.  Treatment of electronic excitations within the adiabatic approximation of time dependent density functional theory , 1996 .

[3]  P. Dirac.  The Quantum Theory of the Emission and Absorption of Radiation , 1927 .

[4]  Dmitrij Rappoport,et al.  Density functional methods for excited states: equilibrium structure and electronic spectra , 2005 .

[5]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[6]  Nicholas C. Handy,et al.  Multiresolution quantum chemistry in multiwavelet bases: time-dependent density functional theory with asymptotically corrected potentials in local density and generalized gradient approximations , 2005 .

[7]  Dmitrij Rappoport,et al.  Property-optimized gaussian basis sets for molecular response calculations. , 2010, The Journal of chemical physics.

[8]  Robert J. Harrison,et al.  Krylov subspace accelerated inexact Newton method for linear and nonlinear equations , 2004, J. Comput. Chem..

[9]  Kirk A Peterson,et al.  Molecular core-valence correlation effects involving the post-d elements Ga-Rn: benchmarks and new pseudopotential-based correlation consistent basis sets. , 2010, The Journal of chemical physics.

[10]  Hideo Sekino,et al.  Basis set limit Hartree-Fock and density functional theory response property evaluation by multiresolution multiwavelet basis. , 2008, The Journal of chemical physics.

[11]  David E. Woon,et al.  Gaussian basis sets for use in correlated molecular calculations. IV. Calculation of static electrical response properties , 1994 .

[12]  Andreas Dreuw,et al.  Single-reference ab initio methods for the calculation of excited states of large molecules. , 2005, Chemical reviews.

[13]  H. Koch,et al.  Gauge invariance of the coupled cluster oscillator strength , 1998 .

[14]  M. Jungen,et al.  Universal Gaussian basis sets for an optimum representation of Rydberg and continuum wavefunctions , 1989 .

[15]  Erwin Schrödinger,et al.  Quantisierung als Eigenwertproblem , 1925 .

[16]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[17]  J. Olsen,et al.  Triplet excitation energies in full configuration interaction and coupled-cluster theory , 2001 .

[18]  H. Sekino,et al.  Time‐dependent solution of molecular quantum systems using multiwavelet , 2011 .

[19]  Erwin Schrödinger,et al.  Über das Verhältnis der Heisenberg‐Born‐Jordanschen Quantenmechanik zu der meinem , 1926 .

[20]  Poul Jørgensen,et al.  Response functions from Fourier component variational perturbation theory applied to a time-averaged quasienergy , 1998 .

[21]  Benjamin T. Miller,et al.  A parallel implementation of the analytic nuclear gradient for time-dependent density functional theory within the Tamm–Dancoff approximation , 1999 .

[22]  Edward F. Valeev,et al.  Computing molecular correlation energies with guaranteed precision. , 2013, The Journal of chemical physics.

[23]  Gregory Beylkin,et al.  Multiresolution quantum chemistry: basic theory and initial applications. , 2004, The Journal of chemical physics.

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

[25]  Luca Frediani,et al.  Fully adaptive algorithms for multivariate integral equations using the non-standard form and multiwavelets with applications to the Poisson and bound-state Helmholtz kernels in three dimensions , 2013 .

[26]  Robert G. Harrison,et al.  Attosecond electron dynamics: A multiresolution approach , 2012 .

[27]  Robert J. Harrison,et al.  MADNESS: A Multiresolution, Adaptive Numerical Environment for Scientific Simulation , 2015, SIAM J. Sci. Comput..

[28]  Martin J. Mohlenkamp,et al.  Algorithms for Numerical Analysis in High Dimensions , 2005, SIAM J. Sci. Comput..

[29]  P. Jørgensen,et al.  Large-scale calculations of excitation energies in coupled cluster theory: The singlet excited states of benzene , 1996 .

[30]  H. Sekino,et al.  Basis set limit computation of dynamic polarizability at near‐resonance region , 2013 .