On Achieving High Accuracy in Quantum Chemical Calculations of 3 d Transition Metal-Containing Systems: A Comparison of Auxiliary-Field Quantum Monte Carlo with Coupled Cluster, Density Functional Theory, and Experiment for Diatomic Molecules.

The bond dissociation energies of a set of 44 3 d transition metal-containing diatomics are computed with phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing a correlated sampling technique. We investigate molecules with H, N, O, F, Cl, and S ligands, including those in the 3dMLBE20 database first compiled by Truhlar and co-workers with calculated and experimental values that have since been revised by various groups. In order to make a direct comparison of the accuracy of our ph-AFQMC calculations with previously published results from 10 DFT functionals, CCSD(T), and icMR-CCSD(T), we establish an objective selection protocol which utilizes the most recent experimental results except for a few cases with well-specified discrepancies. With the remaining set of 41 molecules, we find that ph-AFQMC gives robust agreement with experiment superior to that of all other methods, with a mean absolute error (MAE) of 1.4(4) kcal/mol and maximum error of 3(3) kcal/mol (parentheses account for reported experimental uncertainties and the statistical errors of our ph-AFQMC calculations). In comparison, CCSD(T) and B97, the best performing DFT functional considered here, have MAEs of 2.8 and 3.7 kcal/mol, respectively, and maximum errors in excess of 17 kcal/mol (for the CoS diatomic). While a larger and more diverse data set would be required to demonstrate that ph-AFQMC is truly a benchmark method for transition metal systems, our results indicate that the method has tremendous potential, exhibiting unprecedented consistency and accuracy compared to other approximate quantum chemical approaches.

[1]  Trygve Helgaker,et al.  Basis-set convergence of correlated calculations on water , 1997 .

[2]  Andreas Köhn,et al.  How To Arrive at Accurate Benchmark Values for Transition Metal Compounds: Computation or Experiment? , 2017, Journal of chemical theory and computation.

[3]  G. D. Maria,et al.  Mass-spectrometric study of gaseous molecules. ZnS, ZnSe and ZnTe , 1965 .

[4]  J. Harvey,et al.  On the accuracy of density functional theory in transition metal chemistry , 2006 .

[5]  Mario Motta,et al.  Ab initio computations of molecular systems by the auxiliary‐field quantum Monte Carlo method , 2017, 1711.02242.

[6]  F. Malone,et al.  Auxiliary-field quantum Monte Carlo calculations of the structural properties of nickel oxide. , 2018, The Journal of chemical physics.

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

[8]  D. Khomskii Transition Metal Compounds , 2014 .

[9]  Zhang,et al.  Constrained path quantum Monte Carlo method for fermion ground states. , 1995, Physical review letters.

[10]  Shiwei Zhang,et al.  Frozen-Orbital and Downfolding Calculations with Auxiliary-Field Quantum Monte Carlo. , 2013, Journal of chemical theory and computation.

[11]  C. Evans,et al.  The Pure Rotational Spectrum of AuI. , 2001, Journal of molecular spectroscopy.

[12]  D. Dixon,et al.  Prediction of Bond Dissociation Energies/Heats of Formation for Diatomic Transition Metal Compounds: CCSD(T) Works. , 2017, Journal of chemical theory and computation.

[13]  L. Tran,et al.  Spin-Unrestricted Self-Energy Embedding Theory. , 2018, The journal of physical chemistry letters.

[14]  Angela K. Wilson,et al.  Comparative Study of Single and Double Hybrid Density Functionals for the Prediction of 3d Transition Metal Thermochemistry. , 2012, Journal of chemical theory and computation.

[15]  T. Ziegler Approximate Density Functional Theory as a Practical Tool in Molecular Energetics and Dynamics , 1991 .

[16]  Krishnan Raghavachari,et al.  Assessment of Gaussian-2 and density functional theories for the computation of ionization potentials and electron affinities , 1998 .

[17]  S. Langhoff,et al.  AB Initio Studies of Transition Metal Systems , 1988 .

[18]  G. Chan,et al.  Efficient Ab Initio Auxiliary-Field Quantum Monte Carlo Calculations in Gaussian Bases via Low-Rank Tensor Decomposition. , 2018, Journal of chemical theory and computation.

[19]  T. Dunning,et al.  Valence correlation in the s2dn, sdn+1, and dn+2 states of the first‐row transition metal atoms , 1981 .

[20]  David M. Ceperley,et al.  Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods , 2017, 1705.01608.

[21]  D. Reichman,et al.  Chemical Transformations Approaching Chemical Accuracy via Correlated Sampling in Auxiliary-Field Quantum Monte Carlo. , 2017, Journal of chemical theory and computation.

[22]  M. Motta,et al.  Calculation of interatomic forces and optimization of molecular geometry with auxiliary-field quantum Monte Carlo , 2018, 1803.05599.

[23]  Shiwei Zhang,et al.  Quantum Monte Carlo Calculations in Solids with Downfolded Hamiltonians. , 2014, Physical review letters.

[24]  Jörg Meyer,et al.  Diffusion Monte Carlo for Accurate Dissociation Energies of 3d Transition Metal Containing Molecules. , 2016, Journal of chemical theory and computation.

[25]  Shiwei Zhang,et al.  Many-body study of cobalt adatoms adsorbed on graphene , 2012, 1210.6973.

[26]  Ali Alavi,et al.  Accurate Ab initio calculation of ionization potentials of the first-row transition metals with the configuration-interaction quantum Monte Carlo technique. , 2015, Physical review letters.

[27]  Werner Kutzelnigg,et al.  Quasirelativistic theory equivalent to fully relativistic theory. , 2005, The Journal of chemical physics.

[28]  D. Truhlar,et al.  Multiconfiguration Pair-Density Functional Theory and Complete Active Space Second Order Perturbation Theory. Bond Dissociation Energies of FeC, NiC, FeS, NiS, FeSe, and NiSe. , 2017, The journal of physical chemistry. A.

[29]  Donald G Truhlar,et al.  Do Practical Standard Coupled Cluster Calculations Agree Better than Kohn-Sham Calculations with Currently Available Functionals When Compared to the Best Available Experimental Data for Dissociation Energies of Bonds to 3d Transition Metals? , 2015, Journal of chemical theory and computation.

[30]  J. Beauchamp,et al.  Gas-phase studies of alkene oxidation by transition-metal oxides. Ion-beam studies of CrO+. , 1986, Journal of the American Chemical Society.

[31]  Nathan J DeYonker,et al.  Toward accurate theoretical thermochemistry of first row transition metal complexes. , 2012, The journal of physical chemistry. A.

[32]  L. Szentpály Atom-based thermochemistry: predictions of the sublimation enthalpies of group 12 chalcogenides and the formation enthalpies of their polonides. , 2008 .

[33]  A. Papakondylis Ab initio study of the ground and excited states of the zinc sulfide diatomic system, ZnS , 2011 .

[34]  Miquel Duran,et al.  How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers? , 1996 .

[35]  L. Curtiss,et al.  Assessment of Gaussian-3 and density functional theories for a larger experimental test set , 2000 .

[36]  Shiwei Zhang,et al.  Auxiliary-field quantum Monte Carlo calculations with multiple-projector pseudopotentials , 2016, 1610.06267.

[37]  Marie L. Laury,et al.  Pseudopotential-Based Correlation Consistent Composite Approach (rp-ccCA) for First- and Second-Row Transition Metal Thermochemistry. , 2015, The journal of physical chemistry. A.

[38]  K. P. Kepp,et al.  Chemical Bond Energies of 3d Transition Metals Studied by Density Functional Theory. , 2018, Journal of chemical theory and computation.

[39]  M. Morse Predissociation Measurements of Bond Dissociation Energies. , 2018, Accounts of chemical research.

[40]  Henry Krakauer,et al.  Assessing weak hydrogen binding on Ca+ centers: an accurate many-body study with large basis sets. , 2011, The Journal of chemical physics.

[41]  Stability, energetics, and magnetic states of cobalt adatoms on graphene. , 2014, Physical review letters.

[42]  Donald G Truhlar,et al.  Tests of Exchange-Correlation Functional Approximations Against Reliable Experimental Data for Average Bond Energies of 3d Transition Metal Compounds. , 2013, Journal of chemical theory and computation.

[43]  Henry Krakauer,et al.  Pressure-induced diamond to β-tin transition in bulk silicon: A quantum Monte Carlo study , 2009, 0908.4477.

[44]  Henry Krakauer,et al.  Quantum Monte Carlo method using phase-free random walks with slater determinants. , 2003, Physical review letters.

[45]  Garnet Kin-Lic Chan,et al.  Stripe order in the underdoped region of the two-dimensional Hubbard model , 2016, Science.

[46]  Markus Reiher,et al.  Exact decoupling of the relativistic Fock operator , 2012, Theoretical Chemistry Accounts.

[47]  Shiwei Zhang Ab Initio Electronic Structure Calculations by Auxiliary-Field Quantum Monte Carlo , 2018, Handbook of Materials Modeling.

[48]  Kirk A Peterson,et al.  Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc-Zn. , 2005, The Journal of chemical physics.

[49]  Kirk A. Peterson,et al.  Explicitly correlated composite thermochemistry of transition metal species. , 2013, The Journal of chemical physics.

[50]  R. Zheng,et al.  Studies on rotational spectra of the ground states of ZnO, ZnS, SiSe and SiTe isotopic molecules using an isotopic error compensation approach , 2010 .

[51]  Henry Krakauer,et al.  Auxiliary-field quantum Monte Carlo study of TiO and MnO molecules , 2005, cond-mat/0510791.

[52]  L. von Szentpály Atom-based thermochemistry: predictions of the sublimation enthalpies of group 12 chalcogenides and the formation enthalpies of their polonides. , 2008, The journal of physical chemistry. A.

[53]  D. Truhlar,et al.  Predicting bond dissociation energy and bond length for bimetallic diatomic molecules: a challenge for electronic structure theory. , 2017, Physical chemistry chemical physics : PCCP.

[54]  Donald G Truhlar,et al.  Predicting Bond Dissociation Energies of Transition-Metal Compounds by Multiconfiguration Pair-Density Functional Theory and Second-Order Perturbation Theory Based on Correlated Participating Orbitals and Separated Pairs. , 2017, Journal of chemical theory and computation.

[55]  Nathan J. DeYonker,et al.  Quantitative computational thermochemistry of transition metal species. , 2007, The journal of physical chemistry. A.

[56]  S. Peyerimhoff,et al.  MRD-CI calculations for the low-lying electronic states of scandiumhydride and titaniumhydride , 1983 .

[57]  V. Batista,et al.  The O2-Evolving Complex of Photosystem II: Recent Insights from Quantum Mechanics/Molecular Mechanics (QM/MM), Extended X-ray Absorption Fine Structure (EXAFS), and Femtosecond X-ray Crystallography Data. , 2017, Accounts of chemical research.

[58]  L. Ziurys,et al.  The Pure Rotational Spectrum of ZnS (X^{1}Σ^{+}) , 2009 .

[59]  Shiwei Zhang,et al.  Auxiliary-field quantum Monte Carlo calculations of molecular systems with a Gaussian basis. , 2006, The Journal of chemical physics.

[60]  Richard A Friesner,et al.  Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units. , 2018, Journal of chemical theory and computation.

[61]  Jie Xu,et al.  CPMC-Lab: A Matlab package for Constrained Path Monte Carlo calculations , 2014, Comput. Phys. Commun..

[62]  L. Ziurys,et al.  The pure rotational spectrum of ZnS (X1Σ , 2009 .

[63]  Henry Krakauer,et al.  Bond breaking with auxiliary-field quantum Monte Carlo. , 2007, The Journal of chemical physics.

[64]  B. Shepler,et al.  The group 12 metal chalcogenides: an accurate multireference configuration interaction and coupled cluster study , 2007 .

[65]  Yixiang Cao,et al.  A localized orbital analysis of the thermochemical errors in hybrid density functional theory: achieving chemical accuracy via a simple empirical correction scheme. , 2006, The Journal of chemical physics.

[66]  A. Trautwein Bioinorganic chemistry : transition metals in biology and their coordination chemistry : research report , 1997 .

[67]  Donald G Truhlar,et al.  Multiconfiguration pair-density functional theory: barrier heights and main group and transition metal energetics. , 2015, Journal of chemical theory and computation.

[68]  Shiwei Zhang,et al.  Auxiliary-field quantum Monte Carlo calculations of the molybdenum dimer. , 2016, The Journal of chemical physics.

[69]  D. Dixon,et al.  Extended benchmark studies of coupled cluster theory through triple excitations , 2001 .

[70]  J. Gauss,et al.  Bond Dissociation Energies for Diatomic Molecules Containing 3d Transition Metals: Benchmark Scalar-Relativistic Coupled-Cluster Calculations for 20 Molecules. , 2017, Journal of chemical theory and computation.

[71]  Shiwei Zhang,et al.  An auxiliary-field quantum Monte Carlo study of the chromium dimer. , 2014, The Journal of chemical physics.

[72]  Sandeep Sharma,et al.  PySCF: the Python‐based simulations of chemistry framework , 2018 .

[73]  J. Berkowitz,et al.  Dissociation Energies of Some Metal Sulfides , 1963 .

[74]  B. Roos The Complete Active Space Self‐Consistent Field Method and its Applications in Electronic Structure Calculations , 2007 .

[75]  R. Bartlett,et al.  Coupled-cluster theory in quantum chemistry , 2007 .

[76]  Localized orbital corrections for density functional calculations on transition metal containing systems , 2017 .

[77]  Yann Garniron,et al.  Deterministic Construction of Nodal Surfaces within Quantum Monte Carlo: The Case of FeS. , 2017, Journal of chemical theory and computation.