Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo

Abstract

Transition metal complexes are ubiquitous in biology and chemical catalysis, yet they remain difficult to accurately describe with ab initio methods due to the presence of a large degree of dynamic electron correlation, and, in some cases, strong static correlation which results from a manifold of low-lying states. Progress has been hindered by a scarcity of high quality gas-phase experimental data, while exact ab initio predictions are usually computationally unaffordable due to the large size of the systems. In this work, we present a data set of 34 3d metal-containing complexes with gas-phase ligand-dissociation energies that have reported uncertainties of $\le$ 2 kcal/mol. We perform all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multi-determinant trial wavefunctions selected by a blackbox procedure. We compare the results with those from DFT with various functionals, and DLPNO-CCSD(T). We find MAE of 1.09 $\pm$ 0.28 kcal/mol for our best ph-AFQMC method, vs 2.89 kcal/mol for DLPNO-CCSD(T) and 1.57 - 3.87 kcal/mol for DFT. We find maximum errors of 2.96 $\pm$ 1.71 kcal/mol for our best ph-AFQMC method, vs 9.15 kcal/mol for DLPNO-CCSD(T) and 5.98 - 13.69 kcal/mol for DFT. The reasonable performance of several functionals is in stark contrast to the much poorer accuracy previously demonstrated for diatomics, suggesting a moderation in electron correlation due to ligand coordination. However, the unpredictably large errors for a small subset of cases with both DFT and DLPNO-CCSD(T) leave cause for concern, especially due to the unreliability of common multi-reference indicators. In contrast, the robust and, in principle, systematically improvable results of ph-AFQMC for these realistic complexes establish it as a useful tool for elucidating the electronic structure of transition metal-containing complexes and predicting their gas-phase properties.

arXiv comment37 pages, 8 figures

References in corpus (14)

• Correlated Electrons in High Temperature Superconductors
• Intrinsic atomic orbitals: An unbiased bridge between quantum theory and chemical concepts
• Quantum Monte Carlo method using phase-free random walks with Slater determinants
• Automated construction of molecular active spaces from atomic valence orbitals
• Ab initio computations of molecular systems by the auxiliary-field quantum Monte Carlo method
• Auxiliary-field quantum Monte Carlo calculations of molecular systems with a Gaussian basis
• On Achieving High Accuracy in Quantum Chemical Calculations of 3d Transition Metal Systems: A Comparison of Auxiliary-Field Quantum Monte Carlo with Coupled Cluster, Density Functional Theory, and Experiment for Diatomic Molecules
• Direct comparison of many-body methods for realistic electronic Hamiltonians
• Calculation of Ligand Dissociation Energies in Large Transition-Metal Complexes
• Distinguishing Artificial and Essential Symmetry Breaking in a Single Determinant: Approach and Application to the C${}_{60}$, C${}_{36}$, and C${}_{20}$ Fullerenes
• Singlet-Triplet Energy Gaps of Organic Biradicals and Polyacenes with Auxiliary-Field Quantum Monte Carlo
• Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units
• Chemical Transformations Approaching Chemical Accuracy via Correlated Sampling in Auxiliary-Field Quantum Monte Carlo
• Spin-unrestricted self-energy embedding theory

Cited by in corpus (6)

• Some Recent Developments in Auxiliary-Field Quantum Monte Carlo for Real Materials
• Correlation-consistent Gaussian basis sets for solids made simple
• Optimizing the Regularization in Size-Consistent Second-Order Brillouin-Wigner Perturbation Theory
• Scalable Ab Initio Electronic Structure Methods with Near Chemical Accuracy for Main Group Chemistry
• The Role of High-Order Electron Correlation Effects in a Model System for Non-valence Correlation-bound Anions
• Non-perturbative Many-Body Treatment of Molecular Magnets