Development of Monte Carlo configuration interaction: natural orbitals and second-order perturbation theory.

Approximate natural orbitals are investigated as a way to improve a Monte Carlo configuration interaction (MCCI) calculation. We introduce a way to approximate the natural orbitals in MCCI and test these and approximate natural orbitals from Møller-Plesset perturbation theory and quadratic configuration interaction with single and double substitutions in MCCI calculations of single-point energies. The efficiency and accuracy of approximate natural orbitals in MCCI potential curve calculations for the double hydrogen dissociation of water, the dissociation of carbon monoxide, and the dissociation of the nitrogen molecule are then considered in comparison with standard MCCI when using full configuration interaction as a benchmark. We also use the method to produce a potential curve for water in an aug-cc-pVTZ basis. A new way to quantify the accuracy of a potential curve is put forward that takes into account all of the points and that the curve can be shifted by a constant. We adapt a second-order perturbation scheme to work with MCCI (MCCIPT2) and improve the efficiency of the removal of duplicate states in the method. MCCIPT2 is tested in the calculation of a potential curve for the dissociation of nitrogen using both Slater determinants and configuration state functions.

[1]  J. Greer.  Monte Carlo Configuration Interaction , 1998 .

[2]  C. David Sherrill,et al.  Natural orbitals as substitutes for optimized orbitals in complete active space wavefunctions , 2004 .

[3]  Daniel J. Taylor,et al.  Monte Carlo configuration interaction applied to multipole moments, ionization energies, and electron affinities , 2013, J. Comput. Chem..

[4]  J. Greer.  Consistent treatment of correlation effects in molecular dissociation studies using randomly chosen configurations , 1995 .

[5]  P. Löwdin.  Quantum Theory of Many-Particle Systems. I. Physical Interpretations by Means of Density Matrices, Natural Spin-Orbitals, and Convergence Problems in the Method of Configurational Interaction , 1955 .

[6]  S. Matsika,et al.  High-Multiplicity Natural Orbitals in Multireference Configuration Interaction for Excited States. , 2012, Journal of chemical theory and computation.

[7]  Peter J. Knowles,et al.  A determinant based full configuration interaction program , 1989 .

[8]  Jürgen Gauss,et al.  State-of-the-art density matrix renormalization group and coupled cluster theory studies of the nitrogen binding curve. , 2004, The Journal of chemical physics.

[9]  Martin Head-Gordon,et al.  A perturbative correction to the quadratic coupled-cluster doubles method for higher excitations , 2002 .

[10]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[11]  R. Bartlett,et al.  Monte Carlo configuration interaction predictions for the electronic spectra of Ne, CH2, C2, N2, and H2O compared to full configuration interaction calculations. , 2008, The Journal of chemical physics.

[12]  Harrison Shull,et al.  NATURAL ORBITALS IN THE QUANTUM THEORY OF TWO-ELECTRON SYSTEMS , 1956 .

[13]  C. Q. Lee,et al.  The Computer Journal , 1958, Nature.

[14]  M. Nolan,et al.  A basis set study for the calculation of electronic excitations using Monte Carlo configuration interaction , 2001 .

[15]  Michael W. Schmidt,et al.  A natural orbital diagnostic for multiconfigurational character in correlated wave functions , 1999 .

[16]  Jeppe Olsen,et al.  Full configuration interaction benchmarking of coupled-cluster models for the lowest singlet energy surfaces of N2 , 2000 .

[17]  M. Nolan,et al.  A Monte Carlo configuration generation computer program for the calculation of electronic states of atoms, molecules, and quantum dots , 2000 .

[18]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

[19]  Peter J. Knowles,et al.  A new determinant-based full configuration interaction method , 1984 .

[20]  J. P. Coe,et al.  Calculations of potential energy surfaces using Monte Carlo configuration interaction. , 2012, The Journal of chemical physics.

[21]  Klaus Ruedenberg,et al.  Split-localized orbitals can yield stronger configuration interaction convergence than natural orbitals , 2003 .

[22]  Martin Head-Gordon,et al.  Quadratic configuration interaction. A general technique for determining electron correlation energies , 1987 .

[23]  J. Muenter,et al.  Electric dipole moment of carbon monoxide , 1975 .

[24]  J. Paldus,et al.  Comparison of the open‐shell state‐universal and state‐selective coupled‐cluster theories: H4 and H8 models , 1995 .

[25]  Robert J. Harrison,et al.  Approximating full configuration interaction with selected configuration interaction and perturbation theory , 1991 .

[26]  James C. Greer,et al.  Estimating full configuration interaction limits from a Monte Carlo selection of the expansion space , 1995 .