Full Wave Function Optimization with Quantum Monte Carlo - A study

14 hours ago - Full Wave Function Optimization with Quantum Monte Carlo - A study of the Dissociation Energies of ZnO, FeO, FeH, and CrS ...
0 downloads 0 Views 430KB Size
Subscriber access provided by CAL STATE UNIV SAN FRANCISCO

Quantum Electronic Structure

Full Wave Function Optimization with Quantum Monte Carlo - A study of the Dissociation Energies of ZnO, FeO, FeH, and CrS Jil Ludovicy, Kaveh Haghighi Mood, and Arne Luechow J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.9b00241 • Publication Date (Web): 21 Aug 2019 Downloaded from pubs.acs.org on August 24, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Full Wave Function Optimization with Quantum Monte Carlo - A study of the Dissociation Energies of ZnO, FeO, FeH, and CrS Jil Ludovicy, Kaveh Haghighi Mood, and Arne Lüchow∗ Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, 52062 Aachen, Germany E-mail: [email protected] Phone: +49 241 80 94748

Abstract The dissociation energies of four transition metal dimers are determined using diffusion Monte Carlo. The Jastrow, CI, and molecular orbital parameters of the wave function are both partially and fully optimized with respect to the variational energy. The pivotal role is thereby ascribable to the optimization of the molecular orbital parameters of a complete active space wave function in the presence of a Jastrow correlation function. Excellent results are obtained for ZnO, FeO, FeH, and CrS. In addition, potential energy curves are computed for the first three compounds at multi-reference diffusion Monte Carlo (MR-DMC) level, from which spectroscopic constants such as the equilibrium bond distance, the harmonic frequency, and the anharmonicity are extracted. All of those quantities agree well with the experiment. Furthermore, it is shown for CrS that a restricted active space calculation can yield improved initial orbitals by including

1

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

single and double excitations from the original active space into a set of virtual orbitals. We demonstrated in this study that the fixed-node error in DMC can be systematically reduced for multi-reference systems by orbital optimization in compact active spaces. While DMC calculations with a large number of determinants are possible and very accurate, our results demonstrate that compact wave functions may be sufficient in order to obtain accurate nodal surfaces, which determine the accuracy of DMC, even in the case of transition metal compounds.

Introduction Transition metals and compounds thereof constitute chemical systems that are of great interest in catalytic processes, in electrochemistry as well as in biochemistry. 1–4 Transition metals are highly interesting systems since they compromise open d-shells, they exhibit several oxidation states, and they often show magnetic properties. The precise understanding of how catalysts work is of major importance when it comes to elucidating and predicting catalytic processes. The bond-breaking between transition metals and main group elements plays the primary role in these processes, which justifies the study of transition metal compounds by means of high-level quantum chemical methods. Although the ultimate goal may be to analyze catalytic reactions and other bulk material properties, the investigation of transition metal dimers is the first important step toward an accurate theoretical description of bond breaking processes in bulk materials. The study of these small compounds already poses a considerable challenge due to their strong correlation and complicated electronic structures, thus leading to the calculation of these systems proving rather complex. For the late transition metals, being of great interest especially for catalysis, the static electron correlation yields an important contribution to their properties. The development of efficient methods that are able to capture this part of the correlation thus constitutes a highly active field of research. Progress has been made by designing suitable density functional theory (DFT) functionals 5,6 as well as by further developing multi-reference wave 2

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

function based methods, such as multi-reference coupled cluster (MRCC) 7 . The use of quantum Monte Carlo (QMC) methods is justified because not only are they highly parallelizable but also because they scale, in general, favorably with the number of electrons. 8,9 It remains a challenge to retain the low-order polynomial scaling in QMC when accounting for static correlation. The variational (VMC) and diffusion (DMC) quantum Monte Carlo approaches 10–12 are the most widely spread stochastic methods applied in chemistry and physics to determine the properties of systems when an accuracy beyond mean-field theory is required. Transition metal compounds have successfully been studied with VMC and DMC in the past. Wagner and Mitas 13 investigated transition metal oxides with fixed-node DMC (FN-DMC). Petz and Lüchow 14 reported accurate dissociation energies and ionization potentials for sulfide compounds with FN-DMC. Diedrich et al. 15 studied transition metal carbonyls with regard to their dissociation energies with FN-DMC. Horváthová et al. 16 presented energetics for transition metal organometallics employing QMC methods. Recently, Doblhoff-Dier et al. 17 published dissociation energies of 3d transition metal compounds calculated with DMC. However, for several transition metal compounds, a rather large discrepancy between the theoretical and the experimental dissociation energies can be observed. These compounds are believed to exhibit prominent multi-reference character, the single-determinant approach is not capable of describing these systems correctly. Recently, two of the present authors have shown that a multi-reference ansatz in combination with the optimization of the orbital parameters is necessary to predict the right ground state for and to reproduce the dissociation energy of FeS. 18 In their study, the authors found that the re-optimization of the molecular orbitals (MOs) in the presence of a Jastrow correlation function was the key to obtaining accurate results. This multi-reference DMC (MR-DMC) approach with optimization of the orbital parameters is therefore applied to systems for which the single-reference FN-DMC results show significant deviations to the experimental dissociation energies. Caffarel and coworkers introduced the use of CIPSI wave functions (configuration inter-

3

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

action using perturbative selection made iteratively) as guide functions in DMC. 19 CIPSI wave functions are large but efficient approximations of the Full-CI wave function, where the dynamic correlation is described through an expansion into determinants. The CIPSI approach was applied to the computation of the F2 potential energy curve 20 as well as to the calculation of atomic systems 19,21 . In both cases, impressive results were obtained. Scemama et al. 22 recently applied the CISPI-DMC approach to FeS and confirmed our results. In contrast to the CIPSI-DMC method, we attempt, in our approach, to describe the systems by a small, physically motivated CAS with both orbitals and CI parameters optimized in the presence of a Jastrow correlation function, which accurately describes the dynamic correlation. The accuracy of the DMC method is determined by the accuracy of the nodal surface of the trial wave function. When selected CI wave functions, such as CIPSI, are used, the trial wave function is found to converge toward the Full-CI solution and thus toward the exact nodes (within the basis set limit). The determinant selection is energy-based and it is thus not necessarily optimal for improving the nodal surface. With our approach, we explore the possibility to obtain sufficiently accurate DMC energies by retaining compact trial wave functions and by varying the nodal surface through orbital and CI coefficients optimization in the presence of a Jastrow factor. While this is also an energy-based optimization and not a direct optimization of the nodal surface, the wave function, and thus the nodal surface, has additional flexibility through the optimization of all parameters. Although the restriction to a compact active space will not allow for accurate total energies, the contributions of higher excitations to the nodal surface may not be important for the calculation of dissociation energies or other energy differences. In this paper, ZnO, FeO, FeH as well as CrS are reexamined. Truhlar and coworkers 23 performed DFT and coupled cluster (CC) calculations for a set of 20 transition metal dimers. They found a significant discrepancy to the experimental dissociation energies of ZnO and FeH for both methods. The authors of that work also state a prominent multi-reference character for those compounds. ZnO and FeH were also examined by Doblhoff-Dier et

4

ACS Paragon Plus Environment

Page 4 of 30

Page 5 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

al. 17 using DMC with Kohn-Sham (KS) initial orbitals. The significant deviation to the experimental dissociation energies again suggests that a single-reference approach is not able to correctly describe these systems. For FeH, several works 24–26 even suggest a breakdown of the Born-Oppenheimer approximation. The investigation of FeO is motivated by our previous work regarding the FeS system. The aim is to obtain as accurate results for FeO as were obtained for the isovalent FeS system. In 2011, one of the current authors presented the evaluation of transition metal sulfides with single-determinant DMC. 14 The largest deviations were found for FeS and CrS. Therefore, we reexamined CrS with MR-DMC.

Methods Since complete reviews on QMC methods as well as on the nature of the wave function ansatz used in this study are available 8,10,27 , only a brief overview shall be given in this section. The MR-DMC method employed in this work is thoroughly discussed in our previous publication. 18

Trial Wave Function The trial wave function used in the QMC calculations is of a Slater-Jastrow type,

|ΨT i = eU ·

X

ci |Φi i

(1)

i

with a Jastrow correlation function eU , the configuration interaction (CI) coefficients ci , and the configuration state functions (CSFs) |Φi i. The Jastrow function is totally symmetric with respect to electron permutations and is responsible for the inclusion of the dynamic electron correlation. Throughout this work, a Jastrow factor with electron-electron and electron-nucleus (two-particle terms) as well as with electron-electron-nucleus contributions (three-particle term), developed by Lüchow et al. 28 , is used. The anti-symmetric part of the trial wave function consists of CSFs which are linear combinations of products of spin-up and 5

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 30

spin-down Slater determinants |Dk↑ i and |Dk↓ i, with the coefficients di,k being determined by the spatial and spin symmetries of the state at hand.

|Φi i =

X

di,k |Dk↑ i |Dk↓ i

(2)

k

The determinants are constructed from molecular orbitals with each MO being expanded into a standard basis set. The trial wave function depends on the coordinates of the electrons as well as on the parameter vector p = {α, c, κ},

|ΨT i ≡ |ΨT (R, p)i

(3)

with α describing the Jastrow parameters, c being the CI coefficients and κ corresponding to the orbital rotation parameters. The latter are optimized alternatingly with respect to the VMC energy. There are different methods suitable for energy minimization. Throughout this work, the linear method 29 is used to optimize the Jastrow and CI parameters while the orbital parameters are optimized with the perturbative method (POPT) 29 . Effective core potentials (ECPs) are used for all calculations in order to include scalar relativistic effects and to decrease the computational demand of QMC by making larger time steps feasible.

Computational Approach The trial wave functions were generated with the Molpro 30 package. The initial wave functions were obtained from ab initio calculations, such as HF, KS-DFT (with the B3LYP 31,32 functional), and CASSCF. The active space for the latter included the 4s and 3d orbitals of the metals and the valence p orbitals of the main group elements (1s for H). The QMC calculations were performed with the program Amolqc 28,33,34 , developed in our group. A 69-term Jastrow correlation function (denoted as sm666 in ref. 28) with cusp-less three-particle terms was used for all the calculations.

6

ACS Paragon Plus Environment

Page 7 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

All calculations were performed using the ECPs with the triple-ζ basis sets of Burkatzki, Filippi and Dolg 35,36 , referred to in the following by BFD-VTZ. The non-local part of the pseudopotentials was localized on a spherical grid by using the trial wave function. 37,38 The dissociation energies were spin-orbit (SO) corrected, and the core-valence (CV) correlation contribution was added. The first-order SO corrections for the atoms were derived from experimental splittings. 39 The ones for the molecules were taken from literature. The CV correlation was estimated by means of multi-reference perturbation theory (MR-MP2) 40 , as implemented in GAMESS 41 . The core-valence basis set TK+NOSeC-V-QZP with all diffuse functions 42,43 was used for these calculations. The dissociation energy was calculated with and without correlating the core electrons to estimate that quantity. The same active space as for the QMC calculations was chosen. The importance of including the core-valence correlation contribution in order to obtain accurate transition metal properties was also noted in other studies. 44 The zero-point energy (ZPE) of CrS was determined with the electron structure modelling program Gaussian 45 at B3LYP/BFD-VTZ 35,36 level. The ZPEs of the other compounds were obtained from a Morse fit to the MR-DMC potential energy curves, computed with fully optimized guide functions.

Results and Discussion In this section, the single- and multi-determinant QMC calculations for different transition metal dimers will be discussed. The effect of the MO optimization on the VMC and DMC energies for both approaches is to be investigated in particular. Finally, the dissociation energies of the different compounds are evaluated. The absolute energies will only be discussed for the first system, the ones of the remaining compounds can be found in the supporting information. The DMC energies of the atomic species are given there as well. Elementary information about the dimers is summarized in Table 1. Further details are provided in the sections describing the different compounds. The CSF column in Table 1

7

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 30

represents the number of CSFs with non-zero coefficients, necessary to describe the ground state of the respective compound for a given CAS. Table 1: Ground state, active space, CSFs, equilibrium bond distance (Å), spin-orbit correction, core-valence correlation contribution, and zero point energy for all four compounds. The energy corrections are given in eV. Compound Ground state ZnO FeO FeH CrS a

1

Σ+ 46,47 ∆ 48,49 4 51,52 ∆ 5 53,54 Π 5

CAS

CSFs

re

SO

[16,9] [12,9] [9,7] [10,9]

10 184 30 670

1.709 1.623 1.567 2.0781 54

n/a -0.0558 50 -0.0477 23 -0.0118 54

CV

ZPE

0.0923 0.0462 0.126 0.0537 0.0675 0.114 0.140 0.0278a

calculated at B3LYP/BFD-VTZ level.

ZnO ZnO exhibits an electronic configuration of π 4 σ 2 δ 4 π 4 σ 2 . The [16,9]-CAS is built from the 4s and 3d orbitals of zinc, and from the 2p orbitals of oxygen. The equilibrium bond distance, see Table 1, is obtained from a potential curve at MR-DMC level. The curve was recorded for a small time step of τ = 0.0005 a.u. in order to reach an acceptance ratio of about 99%. The different wave function parameters were optimized alternatingly with respect to the VMC energies. The latter are illustrated in Table 2 together with the zero time step extrapolated DMC energies. The non-optimized parameters, such as the CI and the MO coefficients, are taken from the respective ab initio calculations for some optimization levels. Table 2: ZnO VMC and DMC energies in Eh at various optimization levels (Jas = Jastrow only), using different starting orbitals and BFD-VTZ/sm666. Ansatz

Orbitals Optimization level VMC energy DMC energy

Single det

HF B3LYP opt

Jas Jas Jas+MO

-242.8836(3) -242.8944(3) -242.9013(3)

-242.9931(5) -243.0022(5) -243.0065(6)

[16,9]-CAS

CAS CAS opt

Jas Jas+CI Jas+MO+CI

-242.8971(3) -242.9047(3) -242.9176(3)

-242.9950(5) -243.0023(6) -243.0111(5)

8

ACS Paragon Plus Environment

Page 9 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

A systematic lowering of the VMC energies can be observed from HF over KS (B3LYP) to CAS orbitals. The MO optimization leads to an improvement of the energies in both cases, with the change being more significant for the CAS orbitals. It is interesting to see that the coupling between dynamic and static correlation, which will only be included by optimizing the orbitals in the presence of the Jastrow correlation function, has a substantial impact on the energies. The use of KS orbitals in QMC, on the other hand, partly captures this effect. The fully optimized wave function thus captures the dynamic correlation explicitly through the Jastrow factor while the static one is included through the different configurations. The optimization of the molecular orbital parameters in the presence of a Jastrow correlation function leads to a substantial improvement of the nodal surface for the CAS ansatz. For the KS orbitals, the DMC energy is only slightly lowered by the optimization of the MO parameters which emphasizes that the nodal surface was already almost optimal for a single-determinant ansatz before the optimization. The optimization of all parameters (Jas, MO, and CI) is necessary to obtain DMC energies that are lower than the ones with the optimized KS orbitals. This reveals the influence of the dynamic correlation on the quality of the nodal surface. Without optimizing the MO and CI parameters of the CAS initial wave function, the Jastrow optimization does not change the nodal surface. The VMC as well as the DMC energies follow similar trends for all dimers and are thus only discussed for ZnO. Table 3 contains the dissociation energies of ZnO at different optimization levels. As for the atomic species, the DMC energies for the different starting orbitals (and optimization levels) only differed within the order of the statistical error. The core-valence correlation contribution of ZnO amounts to 0.0923 eV. It yields a substantial contribution to the dissociation energy, making it a non-negligible quantity. ZnO does not exhibit first-order spin-orbit coupling due to its totally symmetric ground state. The ZPE was determined by means of the potential energy curve and corresponds to 0.0462 eV.

9

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 30

Table 3: DMC dissociation energies of ZnO in eV at various optimization levels, using different starting orbitals and BFD-VTZ/sm666. The energies are CV and SO corrected. Ansatz

Orbitals Optimization level

D0

Single det

HF B3LYP opt

Jas Jas Jas+MO

1.20(2) 1.45(2) 1.57(2)

CAS

CAS CAS opt

Jas Jas+CI Jas+MO+CI

1.25(2) 1.45(2) 1.69(2)

Table 4: Bond dissociation energies in eV calculated and measured for ZnO. Investigators This work This work Clemmer et al. 55 Zhang et al. 56 Krogel et al. 57 Weaver et al. 58 Aoto et al. 44 Xu et al. 23

Method SR-DMC MR-DMC Mass Spectrometry from ∆Hf

De

D0

1.61(2) 1.57(2) 1.74(2) 1.69(2) 1.61(4) 1.65(4)

DMC 1.45(2) CASPT2 1.54 CCSD(T)(CV)/CBS-DK 1.55 CCSDT(2)Q /apTZ-DK(3) 1.36 DFT/B97-1-DK 1.30 DFT/M06-L-DK 1.35

1.45

The DMC ansatz with HF nodes fails to reproduce the dissociation energy of ZnO, see Table 4. At the Jas optimization level, the KS orbitals yield a significantly more accurate dissociation energy than the HF and CAS nodes. Without the MO optimization, the DFT and the CAS guide functions yield comparable results. The optimization of the orbital parameters in the presence of a Jastrow correlation function improves the dissociation energy substantially. In comparison to the experimental dissociation energy 55,56 , the single determinant dissociation energy is improved by more than 0.1 eV when optimizing the B3LYP orbitals together with the Jastrow function. The resulting dissociation energy is accurate to less than 0.1 eV. The fully optimized MR-DMC dissociation energy is slightly less accurate but still better than 0.1 eV. The significant improvement of the dissociation energy for the 10

ACS Paragon Plus Environment

Page 11 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

MR-DMC approach by optimizing not only the Jastrow parameters but also the CI and the MO parameters is visualized in Fig. 1. A single-reference treatment thus seems to be accurate enough for the ZnO system. This is in contrast to the findings by Xu et al. 23 based on different multi-reference diagnostics and the discrepancy they found in comparison with the experimental dissociation energies for CC as well as for DFT. Doblhoff-Dier et al. 17 also reported a large deviation to the experiment for their single-reference DMC dissociation energies (about 0.3 eV for ZnO), using KS determinants with different functionals. We believe however, that the discrepancy of our SR-DMC dissociation energy to this work is due to the orbital optimization and the core-valence correlation contribution to the dissociation energy. Krogel et al. 57 published single-determinant DMC potential energy curves of transition metal oxides. Their dissociation energy for ZnO is significantly lower than the one obtained in this work with SR-DMC. The CASPT2 method is not able to reproduce the experimental dissociation energy of Clemmer and coworkers. The CC dissociation energy of Aoto et al. 44 , extrapolated to the complete basis set, core-valence corrected and including scalar relativistic effects underestimates the experimental dissociation energy by about 0.1 eV. Xu et al. 23 studied the performance of CC compared to DFT calculations. Their CC and DFT dissociation energies are significantly lower than the one computed in this work.

11

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Jas Jas+CI Jas+MO+CI

0.8

|D0 − D0exp| [eV]

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 30

0.6

0.4

0.2

0.0

ZnO

FeO

FeH

CrS

Figure 1: Deviation to experimental dissociation energies at the various optimization levels for all four compounds.

FeO The ground state of FeO corresponds to an electronic configuration of σ 2 π 4 σ 2 δ 3 σ 1 π 2 . The potential energy curve at MR-DMC level is computed with a fixed time step of τ = 0.001 a.u. The spin-orbit correction of -0.0558 eV, see Table 1, is taken from ref. 50, the authors of which made use of an earlier experimental study. 59 The wave function is fully optimized with an sm666 Jastrow factor. The VMC and DMC energies of FeO are listed in the supporting information. They are not discussed since they show trends comparable to the ones that were already observed for the ZnO compound.

12

ACS Paragon Plus Environment

Page 13 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Table 5: DMC dissociation energies of FeO in eV at various optimization levels, using different starting orbitals and BFD-VTZ/sm666. The energies are CV and SO corrected. Ansatz

Orbitals Optimization level

D0

Single det

HF B3LYP opt

Jas Jas Jas+MO

2.88(2) 3.69(2) 3.83(2)

CAS

CAS CAS opt

Jas Jas+CI Jas+MO+CI

3.27(2) 3.76(2) 4.11(2)

The dissociation energies of FeO are listed in Table 5. Similarly to the ZnO compound, the optimization of the MO parameters improves the dissociation energy to a large extent. The approach used in this work is able to systematically improve the dissociation energy for the different wave function ansätze and optimization levels, see Figure 1. The singledeterminant as well as the CAS guide functions without MO optimization yield similar results, they all underestimate the dissociation energy substantially. Only the dissociation energy of the fully optimized multi-reference guide function is in good agreement with the experimental results, realized by means of various methods 60–62 (cf. Table 6). Table 6: Bond dissociation energies in eV calculated and measured for FeO. Investigators This work Chestakov et al. 60 Li et al. 61 Smoes and Drowart 62 Krogel et al. 57 Aoto et al. 44 Sakellaris et al. 63 Jensen et al. 64

Method MR-DMC

De

4.17(2) 4.11(2)

Photodissociation Collision-induced dissociation Mass spectrometry DMC CCSD(T) MRCI+DKH2+Q DFT/B3LYP DFT/BP86 DFT/PBE

D0

4.18(1) 4.18(1) 4.16(8) 4.25(1) 4.21 3.69 3.96 5.21 5.31

The DMC approach of Krogel et al. 57 yields a dissociation energy that is comparable to the one computed in this study. The CC dissociation energy of Aoto et al. 44 agrees well 13

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 30

with the experimental ones as well as with the MR-DMC value of this work, while the MRCI result of Sakellaris et al. 63 is significantly lower. As for the DFT approaches of Jensen et al. 64 , the different functionals are not able to yield consistent results, either severely underor overestimating the experimental dissociation energies.

FeH The ground state of FeH is described by the electron configuration σ 2 π 4 δ 2 σ 1 . The [9,7]-CAS, see Table 1, is constructed from the 4s and 3d orbitals of iron, and the 1s orbital of hydrogen. The MR-DMC potential energy curve is recorded at a fixed time step of τ = 0.001 a.u. Table 7 shows a significantly higher dissociation energy for the KS nodes than for the HF nodes at the Jastrow optimization level. The optimization of the molecular orbital parameters shows no improvement of the dissociation energy for the single determinant guide function. For the multi-reference approach, a systematic improvement of the dissociation energy can be observed for the different optimization levels, which is visualized by Figure 1. The MO optimization of the CAS guide function in the presence of a Jastrow correlation factor has a significant effect on the dissociation energy of FeH, it is increased by about 0.4 eV. Note that even if the anti-symmetric part of the multi-reference guide function is not optimized, a more accurate dissociation energy is obtained compared to the one from the fully optimized single-reference wave function, which speaks for FeH exhibiting a prominent multi-reference character. Table 7: DMC dissociation energies of FeH in eV at various optimization levels, using different starting orbitals and BFD-VTZ/sm666. The energies are CV and SO corrected. Ansatz

Orbitals Optimization level

D0

Single det

HF B3LYP opt

Jas Jas Jas+MO

0.81(2) 1.02(2) 1.02(2)

CAS

CAS CAS opt

Jas Jas+CI Jas+MO+CI

1.10(2) 1.37(2) 1.79(2)

14

ACS Paragon Plus Environment

Page 15 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

A good agreement is achieved when comparing the MR-DMC dissociation energy to the experimental one of Schultz and Armentrout 52 , see Table 8. The breakdown of the BO approximation, which was mentioned in several studies 24–26 , can be refuted by the accurate MR-DMC results. The DFT calculations with different functionals, performed by Jensen et al. 64 , fail to yield satisfactory results since they severely overestimate the dissociation energy. The CC results of Aoto et al. 44 and Cheng et al. 65 agree well with our dissociation energy. The focal point analysis (FPA) of DeYonker and Allen 66 yields a substantial deviation to the experimental dissociation energy. Xu and coworkers 23 confirmed the multi-reference character of FeH by different diagnostics. However, they obtained an accurate dissociation energy with CCSDT(2)Q including scalar relativistic effects, while their reported DFT dissociation energies are severely larger than the experimental value. Nonetheless, they argued that KS DFT yields overall comparable results to CC theory for the twenty transition metal compounds that they investigated. Table 8: Bond dissociation energies in eV calculated and measured for FeH. Investigators

Method

This work

MR-DMC

Schultz and Armentrout 52 Jensen et al. 64

Xu et al.

23

Aoto et al. 44 Cheng et al. 65 DeYonker and Allen 66

De

D0

1.90(2) 1.79(2)

Mass Spectrometry

1.63(8)

DFT/B3LYP DFT/BP86 DFT/PBE

2.10 2.41 2.30

CCSDT(2)Q /apTZ-DK(3) DFT/B97-1-DK DFT/M06-L-DK CCSD(T)(CV)/CBS-DK CCSD(T) FPA

1.78 2.00 2.17 1.95 1.99 1.86

CrS A slightly modified initial wave function compared to the usual CAS ansatz is chosen for the CrS system due to the inability to converge the MO parameters with QMC when starting 15

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 30

from CAS orbitals generated by a [10,9]-CASSCF calculation. The active orbitals of the CAS wave function are further relaxed by performing a RASSCF calculation with single and double excitations into a set of virtual orbitals. The original [10,9]-CAS corresponds to the RAS2, where all possible excitations are performed while a RAS3 with 11 virtual orbitals is created for the single and double excitations from the RAS2. The RAS1 remains empty. The RASSCF calculation will henceforth be referred to as [10,9;2,11]-RAS calculation. The aim of this approach is to obtain better initial orbitals that can then be further optimized in a QMC energy minimization calculation. The CAS (=RAS2) orbitals are optimized in the partial presence of dynamic correlation through excitations to the RAS3. The RAS2 orbitals are hence expected to be closer to the converged orbitals in a full MR-VMC optimization. The CAS for the QMC calculations is however built similarly to the one of the other compounds, namely from the 4s and the 3d orbitals of chromium, and from the 3p orbitals of sulfide. The ground state of CrS is described by the electron configuration σ 2 π 4 σ 1 δ 2 π 1 . The calculations for CrS are performed with the experimental bond length of 2.0781 Å. 54 Table 9: DMC dissociation energies of CrS in eV at various optimization levels, using different starting orbitals and BFD-VTZ/sm666. The energies are CV and SO corrected. Ansatz

Orbitals Optimization level

D0

Single det

HF B3LYP opt

Jas Jas Jas+MO

2.05(2) 2.77(2) 2.77(2)

CAS

CAS CAS RAS2 RAS2 opt

Jas Jas+CI Jas Jas+CI Jas+MO+CI

2.43(2) 2.70(2) 2.80(2) 3.04(2) 3.10(2)

The dissociation energies of CrS for the different approaches are listed in Table 9. Similarly to the other compounds, a systematic improvement of the dissociation energy can be observed for the different methods and optimization levels. The KS nodes appear ideal since the MO optimization does not improve the dissociation energy. The ansatz with CAS or16

ACS Paragon Plus Environment

Page 17 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

bitals yields lower dissociation energies than the one obtained with KS orbitals at the same optimization level. Relaxing the initial active orbitals through a RAS calculation has a substantial effect on the dissociation energy. Not only are the dissociation energies significantly improved when comparing them to the ones obtained with CAS initial orbitals but also the molecular orbital parameters could be successfully optimized. At a given optimization level, the dissociation energies for the different CAS nodes differ by about 0.3 eV. When further optimizing the orbitals initially taken from the RAS calculation, the dissociation energy can be improved by 0.05 eV. Figure 1 shows that the deviations between the experimental dissociation energies and the ones computed with MR-DMC can be systematically reduced by increasing the level of optimization. Table 10: Bond dissociation energies in eV calculated or measured for CrS. Investigators This work Drowart et al. 67

Method

D0

MR-DMC

3.10(2)

Mass Spectrometry 3.36(15)

Petz and Lüchow 14 Bauschlicher and Maitre 68 Liang and Andrews 69

DMC/PPII CCSD(T) DFT/BPW91

2.969(9) 2.89 3.33

Table 10 yields experimental and theoretical dissociation energies for CrS. The calculated dissociation energy is larger than both, the single-determinant DMC 14 and the CCSD(T) 68 values, but still smaller than the experimental D0 of Drowart et al. 67 . Assessing the accuracy of the MR-DMC result proves challenging due to the large experimental error bar. Our dissociation energy is about 0.1 eV below the lower bound of Drowart and coworkers. In order to estimate the accuracy of the obtained MR-DMC result, experimental data with smaller error bars are needed.

17

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Spectroscopic Constants The potential energy curves of ZnO, FeO, and FeH were computed at the fully optimized MR-DMC level and fitted to Morse functions from which spectroscopic constants, such as the equilibrium bond distance (minimum of the Morse curve), the harmonic frequency as well as the anharmonicity could be deduced. The evaluation of those constants allows an assessment of the employed method. Table 11 illustrates the obtained quantities and compares them for different methods. For ZnO, the equilibrium bond distance is in good agreement with the experimental bond length of Zack et al. 70 and it is slightly shorter than the CC and DFT ones. The MR-DMC bond distance of FeO is slightly larger than the one obtained from other theoretical methods and it agrees with the experiment. As for FeH, the equilibrium bond distance obtained from the Morse fit is similar to the one from other theoretical studies. The harmonic frequencies and the anharmonicities obtained from the Morse fit are in good agreement as well with the experimental as with the theoretical results for all three compounds.

18

ACS Paragon Plus Environment

Page 18 of 30

Page 19 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Table 11: Spectroscopic constants for the different transition metal compounds. The equilibrium bond distance is given in Å, the harmonic frequency and the anharmonicity in cm−1 . System Investigators This work ZnO

Zack et al. 70 Fancher et al. 46 Weaver et al. 58 Bauschlicher and Partridge 47 This work

FeO

Allen et al.a Drechsler et al. 49 Hendrickx and Anam 71 Sakellaris et al. 63 This work

FeH

Philips et al. 72 Dulick et al. 73 DeYonker and Allen 66 Jensen et al. 64

a

Method

re

ωe

ωe x e

MR-DMC

1.709

746(8)

4.4(1)

738 805(40)

4.88

Direct-absorption methods 1.7047(2) Photoelectron Spectrum CASPT2 CCSD(T) DFT/B3LYP

1.7 1.719 1.713

742 727.2 741

5.83

MR-DMC

1.623

866(79)

4.7(7)

882

4

1.619 anion-ZEKE CASPT2 MRCI RCCSD(T)

1.612 1.612 1.607

887 864 905

7.2 5.9

MR-DMC

1.567

1842(27)

38.9(9)

Near IR Spectrum CCSDT DFT/B3LYP

1826.86 31.96 1831.8(19) 34.9(9) 1.5660 1.57

1798.8

37.8

derived from Allen et al. 74

Conclusion The dissociation energies of ZnO, FeO, FeH, and CrS were determined through single- and multi-determinant DMC calculations. The Jastrow, CI, and MO parameters of the wave functions were both partially and fully optimized with respect to the energy. A systematic improvement of the dissociation energy could be observed for all compounds for the different ansätze. In the single determinant approach, optimizing the KS orbitals led for all four systems to either minor or no significant improvement of the nodal surface of the guide functions. For the multi-reference ansatz, on the other hand, the optimization of the molecular orbital parameters in the presence of a Jastrow correlation function is the key 19

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

contribution. A good agreement of the MR-DMC dissociation energy with the experimental ones was achieved for ZnO, FeO, and FeH. We found that the ZnO dissociation energy could be obtained within 0.1 eV already with a single-reference ansatz, but only after MO optimization. In addition, potential energy curves at MR-DMC level were recorded for these three compounds, which yielded equilibrium bond distances and spectroscopic constants that agree well with literature. As for CrS, the complex MO optimization was tackled by employing more accurate initial orbitals, generated by a RASSCF calculation. The calculated dissociation energy of CrS agrees well with other theoretical methods. Unfortunately, the error bar of the experimental dissociation energy is rather large which is why the accuracy of the obtained result is difficult to assess. Our results show that it is possible to obtain accurate dissociation energies and properties by compact wave functions generated from a small, physically motivated CAS.

Acknowledgement The authors thank the Jülich-Aachen research alliance (JARA) and the RWTH Aachen University for the granted computing time under project rwth0278. The authors thank Christina Zitlau for the performance of preliminary calculations.

Author Information Corresponding Author E-mail: [email protected] (Arne Lüchow)

Funding The authors declare no funding for this project.

20

ACS Paragon Plus Environment

Page 20 of 30

Page 21 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Supporting Information Available The supporting information contains the DMC energies and the first-order spin-orbit corrections of the atoms as well as the VMC and DMC energies for FeO, FeH, and CrS. This information is available free of charge via the Internet at http://pubs.acs.org

References (1) Stiefel, E. I. Transition Metal Sulfur Chemistry: Biological and Industrial Significance and Key Trends. In Transit. Met. Sulfur Chem.; American Chemical Society, 1996. (2) Xiao, S.; Li, X.; Sun, W.; Guan, B.; Wang, Y. General and facile synthesis of metal sulfide nanostructures: In situ microwave synthesis and application as binder-free cathode for Li-ion batteries. Chem. Eng. J. 2016, 306, 251–259. (3) Kunkes, E. L.; Studt, F.; Abild-Pedersen, F.; Schlögl, R.; Behrens, M. Hydrogenation of CO2 to methanol and CO on Cu/ZnO/Al2 O3 : Is there a common intermediate or not? J. Catal. 2015, 328, 43–48. (4) Eshghi, H.; Rahimizadeh, M.; Saberi, S. Fe(HSO4)3 as an inexpensive, eco-friendly, heterogeneous and reusable catalyst for acetal/ketal formation and their facile regeneration. Catal. Commun. 2008, 9, 2460–2466. (5) Zhao, Y.; Truhlar, D. G. A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions. J. Chem. Phys. 2006, 125, 194101. (6) Fuchs, M.; Niquet, Y.-M.; Gonze, X.; Burke, K. Describing static correlation in bond dissociation by Kohn–Sham density functional theory. J. Chem. Phys. 2005, 122, 094116.

21

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(7) Musia, M.; Perera, A.; Bartlett, R. J. Multireference coupled-cluster theory: The easy way. J. Chem. Phys. 2011, 134, 114108. (8) Lüchow, A. Quantum Monte Carlo methods. WIREs Comput. Mol. Sci. 2011, 1, 388– 402. (9) Lüchow, A.; Petz, R.; Schwarz, A. Electron Structure Quantum Monte Carlo. Zeitschrift für Phys. Chemie 2010, 224, 343–355. (10) Austin, B. M.; Zubarev, D. Y.; Lester, W. A. Quantum Monte Carlo and Related Approaches. Chem. Rev. 2012, 112, 263–288. (11) Wagner, L. K. Quantum Monte Carlo for Ab Initio calculations of energy-relevant materials. Int. J. Quantum Chem. 2014, 114, 94–101. (12) Foulkes, W. M. C.; Mitas, L.; Needs, R. J.; Rajagopal, G. Quantum Monte Carlo simulations of solids. Rev. Mod. Phys. 2001, 73, 33–83. (13) Wagner, L. K.; Mitas, L. Energetics and dipole moment of transition metal monoxides by quantum Monte Carlo. J. Chem. Phys. 2007, 126, 034105. (14) Petz, R.; Lüchow, A. Energetics of diatomic transition metal sulfides ScS to FeS with diffusion quantum Monte Carlo. ChemPhysChem 2011, 12, 2031–2034. (15) Diedrich, C.; Lüchow, A.; Grimme, S. Performance of diffusion Monte Carlo for the first dissociation energies of transition metal carbonyls. J. Chem. Phys. 2005, 122, 021101. (16) Horváthová, L.; Dubecký, M.; Mitas, L.; Štich, I. Quantum Monte Carlo Study of πBonded Transition Metal Organometallics: Neutral and Cationic Vanadium–Benzene and Cobalt–Benzene Half Sandwiches. J. Chem. Theory Comput. 2013, 9, 390–400. (17) Doblhoff-Dier, K.; Meyer, J.; Hoggan, P. E.; Kroes, G. J.; Wagner, L. K. Diffusion Monte Carlo for Accurate Dissociation Energies of 3d Transition Metal Containing Molecules. J. Chem. Theory Comput. 2016, 12, 2583–2597. 22

ACS Paragon Plus Environment

Page 22 of 30

Page 23 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

(18) Haghighi Mood, K.; Lüchow, A. Full Wave Function Optimization with Quantum Monte Carlo and Its Effect on the Dissociation Energy of FeS. J. Phys. Chem. A 2017, 121, 6165–6171. (19) Giner, E.; Scemama, A.; Caffarel, M. Using perturbatively selected configuration interaction in quantum Monte Carlo calculations. Can. J. Chem. 2013, 91, 879–885. (20) Giner, E.; Scemama, A.; Caffarel, M. Fixed-node diffusion Monte Carlo potential energy curve of the fluorine molecule F2 using selected configuration interaction trial wavefunctions. J. Chem. Phys. 2015, 142, 044115. (21) Scemama, A.; Applencourt, T.; Giner, E.; Caffarel, M. Accurate nonrelativistic groundstate energies of 3d transition metal atoms. J. Chem. Phys. 2014, 141, 244110. (22) Scemama, A.; Garniron, Y.; Caffarel, M.; Loos, P.-F. Deterministic construction of nodal surfaces within quantum Monte Carlo: the case of FeS. J. Chem. Theory Comput. 2018, 14, 1395–1402. (23) Xu, X.; Zhang, W.; Tang, M.; Truhlar, D. G. 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? J. Chem. Theory Comput. 2015, 11, 2036–2052. (24) Brown, J. M.; Körsgen, H.; Beaton, S. P.; Evenson, K. M. The rotational and finestructure spectrum of FeH, studied by far-infrared laser magnetic resonance. J. Chem. Phys. 2006, 124, 234309. (25) Harrison, J. J.; Brown, J. M. Measurement of the Magnetic Properties of FeH in Its X 4 ∆ and F 4 ∆ States from Sunspot Spectra. Astrophys. J. 2008, 686, 1426–1431. (26) Shulyak, D.; Reiners, A.; Wende, S.; Kochukhov, O.; Piskunov, N.; Seifahrt, A. Mod-

23

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

elling the molecular Zeeman-effect in M-dwarfs: methods and first results. Astron. Astrophys. 2010, 523, A37. (27) Hammond, B. L.; Lester, Jr., W. A.; Reynolds, P. J. Monte Carlo methods in ab initio quantum chemistry; World Scientific, 1994. (28) Lüchow, A.; Sturm, A.; Schulte, C.; Haghighi Mood, K. Generic expansion of the Jastrow correlation factor in polynomials satisfying symmetry and cusp conditions. J. Chem. Phys. 2015, 142, 084111. (29) Toulouse, J.; Umrigar, C. J. Optimization of quantum Monte Carlo wave functions by energy minimization. J. Chem. Phys. 2007, 126, 084102. (30) Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G.; Shamasundar, K. R.; Adler, T. B.; Amos, R. D.; Bernhardsson, A.; Berning, A.; Cooper, D. L.; And, M. J. O. D.; Dobbyn, A. J.; Eckert, F.; Goll, E.; Hampel, C.; Hesselmann, A.; Hetzer, G.; Hrenar, T.; Jansen, G.; Köoppl, C.; Liu, Y.; Lloyd, A. W.; Mata, R. A.; May, A. J.; McNicholas, S. J.; Meyer, W.; Mura, M. E.; Nicklass, A.; O’Neill, D. P.; Palmieri, P.; Peng, D.; Pflüger, K.; Pitzer, R.; Reiher, M.; Shiozaki, T.; Stoll, H.; Stone, A. J.; Tarroni, R.; Thorsteinsson, T.; Wang, M. MOLPRO, version 2015.1, a package of ab initio programs. (31) Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. (32) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. (33) Manten, S.; Lüchow, A. On the accuracy of the fixed-node diffusion quantum Monte Carlo method. J. Chem. Phys. 2001, 115, 5362–5366.

24

ACS Paragon Plus Environment

Page 24 of 30

Page 25 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

(34) Lüchow, A. Maxima of |Ψ|2 : A connection between quantum mechanics and Lewis structures. J. Comput. Chem. 2014, 35, 854–864. (35) Burkatzki, M.; Filippi, C.; Dolg, M. Energy-consistent pseudopotentials for quantum Monte Carlo calculations. J. Chem. Phys. 2007, 126, 234105. (36) Burkatzki, M.; Filippi, C.; Dolg, M. Energy-consistent small-core pseudopotentials for 3d-transition metals adapted to quantum Monte Carlo calculations. J. Chem. Phys. 2008, 129, 164115. (37) Fahy, S.; Wang, X. W.; Louie, S. G. Variational quantum Monte Carlo nonlocal pseudopotential approach to solids: Formulation and application to diamond, graphite, and silicon. Phys. Rev. B 1990, 42, 3503–3522. (38) Mitáš, L.; Shirley, E. L.; Ceperley, D. M. Nonlocal pseudopotentials and diffusion Monte Carlo. J. Chem. Phys. 1991, 95, 3467–3475. (39) Kramida, A.; Ralchenko, Y.; Reader, J.; NIST ASD Team, NIST Atomic Spectra Database (ver. 5.5) physics.nist.gov/asd. 2017. (40) Nakano, H. Quasidegenerate perturbation theory with multiconfigurational selfconsistent-field reference functions. J. Chem. Phys. 1993, 99, 7983–7992. (41) Schmidt, M. W.; Baldridge, K. K.; Boatz, J. A.; Elbert, S. T.; Gordon, M. S.; Jensen, J. H.; Koseki, S.; Matsunaga, N.; Nguyen, K. A.; Su, S.; Windus, T. L.; Dupuis, M.; Montgomery, J. A. General atomic and molecular electronic structure system. J. Comput. Chem. 1993, 14, 1347–1363. (42) Koga, T.; Tatewaki, H.; Matsuyama, H.; Satoh, Y. Contracted Gaussian-type basis functions revisited. III. Atoms K through Kr. Theor. Chem. Acc. 1999, 102, 105–111. (43) Noro, T.; Sekiya, M.; Koga, T.; Matsuyama, H. Valence and correlated basis sets for the first-row transition atoms from Sc to Zn. Theor. Chem. Acc. 2000, 104, 146–152. 25

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(44) Aoto, Y. A.; De Lima Batista, A. P.; Köhn, A.; De Oliveira-Filho, A. G. How to Arrive at Accurate Benchmark Values for Transition Metal Compounds: Computation or Experiment? J. Chem. Theory Comput. 2017, 13, 5291–5316. (45) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery , Jr., J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, Revision A02 ; Gaussian, Inc.: Wallingford CT, 2009. (46) Fancher, C. A.; de Clercq, H. L.; Thomas, O. C.; Robinson, D. W.; Bowen, K. H. Zinc oxide and its anion: A negative ion photoelectron spectroscopic study. J. Chem. Phys. 1998, 109, 8426–8429. (47) Bauschlicher, C. W.; Partridge, H. A comparison of ZnO and ZnO− . J. Chem. Phys. 1998, 109, 8430–8434. (48) Cheung, A.-C.; Lee, N.; Lyyra, A.; Merer, A.; Taylor, A. Spectroscopic properties of the 5 ∆i ground state of FeO. J. Mol. Spectrosc. 1982, 95, 213–225. (49) Drechsler, G.; Boesl, U.; Bäßmann, C.; Schlag, E. W. Mass selected anion-zero kinetic

26

ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

energy photoelectron spectroscopy (anion-ZEKE): Ground and low excited states of FeO. J. Chem. Phys. 1997, 107, 2284–2291. (50) Schultz, N. E.; Zhao, Y.; Truhlar, D. G. Density Functionals for Inorganometallic and Organometallic Chemistry. J. Chem. Phys. A 2005, 109, 11127–11143. (51) Stevens, A. E.; Feigerle, C. S.; Lineberger, W. C. Laser photoelectron spectroscopy of MnH− and FeH− : Electronic structures of the metal hydrides, identification of a lowspin excited state of MnH, and evidence for a low-spin ground state of FeH. J. Chem. Phys. 1983, 78, 5420–5431. (52) Schultz, R. H.; Armentrout, P. B. The gas-phase thermochemistry of FeH. J. Chem. Phys. 1991, 94, 2262–2268. (53) Shi, Q.; Ran, Q.; Tam, W.; Leung, J.-H.; Cheung, A.-C. Laser-induced fluorescence spectroscopy of CrS. Chem. Phys. Lett. 2001, 339, 154–160. (54) Pulliam, R. L.; Ziurys, L. M. The pure rotational spectrum of the CrS radical in its X 5

Πr state. J. Chem. Phys. 2010, 133, 174313.

(55) Clemmer, D. E.; Dalleska, N. F.; Armentrout, P. B. Reaction of Zn+ with NO2 . The gas-phase thermochemistry of ZnO. J. Chem. Phys. 1991, 95, 7263–7268. (56) Zhang, W.; Truhlar, D. G.; Tang, M. Tests of exchange-correlation functional approximations against reliable experimental data for average bond energies of 3d transition metal compounds. J. Chem. Theory Comput. 2013, 9, 3965–3977. (57) Krogel, J. T.; Santana, J. A.; Reboredo, F. A. Pseudopotentials for quantum Monte Carlo studies of transition metal oxides. Phys. Rev. B 2016, 93, 075143. (58) Weaver, M. N.; Merz, K. M.; Ma, D.; Kim, H. J.; Gagliardi, L. Calculation of Heats of Formation for Zn Complexes: Comparison of Density Functional Theory, Second Order

27

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Perturbation Theory, Coupled-Cluster and Complete Active Space Methods. J. Chem. Theory Comput. 2013, 9, 5277–5285. (59) Merer, A. J. Spectroscopy of the Diatomic 3d Transition Metal Oxides. Annu. Rev. Phys. Chem. 1989, 40, 407–438. (60) Chestakov, D. A.; Parker, D. H.; Baklanov, A. V. Iron monoxide photodissociation. J. Chem. Phys. 2005, 122, 084302. (61) Li, M.; Liu, S.-R.; Armentrout, P. B. Collision-induced dissociation studies of Fem O+ n: Bond energies in small iron oxide cluster cations, Fem O+ n (m=1–3, n=1–6). J. Chem. Phys. 2009, 131, 144310. (62) Smoes, S.; Drowart, J. Determination of the Dissociation Energies of Gaseous Iron Monoxide and Manganese Monoxide by the Mass Spectrometric Knudsen Cell Method. In Mod. High Temp. Sci.; Humana Press, 1984; pp 31–52. (63) Sakellaris, C. N.; Miliordos, E.; Mavridis, A. First principles study of the ground and excited states of FeO, FeO+ , and FeO− . J. Chem. Phys. 2011, 134, 234308. (64) Jensen, K. P.; Roos, B. O.; Ryde, U. Performance of density functionals for first row transition metal systems. J. Chem. Phys. 2007, 126, 014103. (65) Cheng, L.; Gauss, J.; Ruscic, B.; Armentrout, P. B.; Stanton, J. F. Bond Dissociation Energies for Diatomic Molecules Containing 3d Transition Metals: Benchmark ScalarRelativistic Coupled-Cluster Calculations for 20 Molecules. J. Chem. Theory Comput. 2017, 13, 1044–1056. (66) DeYonker, N. J.; Allen, W. D. Taming the low-lying electronic states of FeH. J. Chem. Phys. 2012, 137, 234303. (67) Drowart, J.; Pattoret, A.; Smoes, S. Mass Spectrometric Studies Of The Vaporization Of Refractory Compounds. In Proc. Brit. Ceram. Soc.; 1967; Vol. 8; pp 67–89. 28

ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

(68) Bauschlicher, C. W.; Maitre, P. Theoretical study of the first transition row oxides and sulfides. Theor. Chim. Acta 1995, 90, 189–203. (69) Liang, B.; Andrews, L. Infrared Spectra and Density Functional Theory Calculations of Group 6 Transition Metal Sulfides in Solid Argon. J. Phys. Chem. A 2002, 106, 6945–6951. (70) Zack, L. N.; Pulliam, R. L.; Ziurys, L. M. The pure rotational spectrum of ZnO in the X1 Σ+ and a3 Πi states. J. Mol. Spectrosc. 2009, 256, 186–191. (71) Hendrickx, M. F. A.; Anam, K. R. A New Proposal for the Ground State of the FeO− Cluster in the Gas Phase and for the Assignment of Its Photoelectron Spectra. J. Phys. Chem. A 2009, 113, 8746–8753. (72) Phillips, J. G.; Davis, S. P.; Lindgren, B.; Balfour, W. J. The near-infrared spectrum of the FeH molecule. Astrophys. J. Suppl. Ser. 1987, 65, 721–778. (73) Dulick, M.; Bauschlicher, Jr., C. W.; Burrows, A.; Sharp, C. M.; Ram, R. S.; Bernath, P. Line Intensities and Molecular Opacities of the FeH F 4 ∆i – X 4 ∆i Transition. Astrophys. J. 2003, 594, 651–663. (74) Allen, M.; Ziurys, L.; Brown, J. The millimeter-wave spectrum of FeO in its X5 ∆i state (v = 0): a study of all five spin components. Chem. Phys. Lett. 1996, 257, 130–136.

29

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

30

ACS Paragon Plus Environment

Page 30 of 30