The accuracy of diffusion quantum Monte Carlo simulations in the determination of molecular equilibrium structures

For a test set of 17 first-row small molecules, the equilibrium structures are calculated with Ornstein-Uhlenbeck diffusion quantum Monte Carlo simulations guiding by trial wave functions constructed from floating spherical Gaussian orbitals and spherical Gaussian geminals. To measure performance of the Monte Carlo calculations, the mean deviation, the mean absolute deviation, the maximum absolute deviation, and the standard deviation of Monte Carlo calculated equilibrium structures with respect to empirical equilibrium structures are given. This approach is found to yield results having a uniformly high quality, being consistent with empirical equilibrium structures and surpassing calculated values from the coupled cluster model with single, double, and noniterative triple excitations [CCSD(T)] with the basis sets of cc-pCVQZ and cc-pVQZ. The nonrelativistic equilibrium atomization energies are also presented to assess performance of the calculated methods. The mean absolute deviations regarding experimental atomization energy are 0.16 and 0.21 kcal/mol for the Monte Carlo and CCSD(T)/cc-pCV(56)Z calculations, respectively.     

Prediction of the thermophysical properties of pure neon, pure argon, and the binary mixtures neon-argon and argon-krypton by Monte Carlo simulation using ab initio potentials

Gibbs ensemble Monte Carlo simulations were used to test the ability of intermolecular pair potentials derived ab initio from quantum mechanical principles, enhanced by Axilrod-Teller triple-dipole interactions, to predict the vapor-liquid phase equilibria of pure neon, pure argon, and the binary mixtures neon-argon and argon-krypton. The interaction potentials for Ne-Ne, Ar-Ar, Kr-Kr, and Ne-Ar were taken from literature; for Ar-Kr a different potential has been developed. In all cases the quantum mechanical calculations had been carried out with the coupled-cluster approach [CCSD(T) level of theory] and with correlation consistent basis sets; furthermore an extrapolation scheme had been applied to obtain the basis set limit of the interaction energies. The ab initio pair potentials as well as the thermodynamic data based on them are found to be in excellent agreement with experimental data; the only exception is neon. It is shown, however, that in this case the deviations can be quantitatively explained by quantum effects. The interaction potentials that have been developed permit quantitative predictions of high-pressure phase equilibria of noble-gas mixtures.     

A study of electron affinities using the initiator approach to full configuration interaction quantum Monte Carlo

For the atoms with Z ≤ 11, energies obtained using the "initiator" extension to full configuration interaction quantum Monte Carlo (i-FCIQMC) come to within statistical errors of the FCIQMC results. As these FCIQMC values have been shown to converge onto FCI results, the i-FCIQMC method allows similar accuracy to be achieved while significantly reducing the scaling with the size of the Slater determinant space. The i-FCIQMC electron affinities of the Z ≤ 11 atoms in the aug-cc-pVXZ basis sets are presented here. In every case, values are obtained to well within chemical accuracy [the mean absolute deviation (MAD) from the relativistically corrected experimental values is 0.41 mE(h)], and significantly improve on coupled cluster with singles, doubles and perturbative triples [CCSD(T)] results. Since the only remaining source of error is basis set incompleteness, we have investigated using CCSD(T)-F12 contributions to correct the i-FCIQMC results. By doing so, much faster convergence with respect to basis set size may be achieved for both the electron affinities and the FCIQMC ionization potentials presented in a previous paper. With this F12 correction, the MAD can be further reduced to 0.13 mE(h) for the electron affinities and 0.31 mE(h) for the ionization potentials.     

Ornstein-Uhlenbeck diffusion quantum Monte Carlo study on the bond lengths and harmonic frequencies of some first-row diatomic molecules

This article accesses the performance of the Ornstein-Uhlenbeck diffusion quantum Monte Carlo with regard to the calculation of molecular geometries and harmonic frequencies of H2, LiH, HF, Li2, LiF, CO, N2, and F2 molecules. A comparison of the results for the eight first-row diatomic molecules from experiments, CCSD(T)/6-311G(3df,3pd) and CCSD(T)/cc-pV5Z levels of theory as well as our work is given. The results presented show that quantum Monte Carlo is becoming powerful tools for ab initio electronic structure calculations.     

Quantum Monte Carlo for the x-ray absorption spectrum of pyrrole at the nitrogen K-edge

Fixed-node diffusion Monte Carlo (FNDMC) is used to simulate the x-ray absorption spectrum of a gas-phase pyrrole molecule at the nitrogen K-edge. Trial wave functions for core-excited states are constructed from ground-state Kohn-Sham determinants substituted with singly occupied natural orbitals from configuration interaction with single excitations calculations of the five lowest valence-excited triplet states. The FNDMC ionization potential (IP) is found to lie within 0.3 eV of the experimental value of 406.1 ± 0.1 eV. The transition energies to anti-bonding virtual orbitals match the experimental spectrum after alignment of IP values and agree with the existing assignments.     

Breaking the carbon dimer: the challenges of multiple bond dissociation with full configuration interaction quantum Monte Carlo methods

The full configuration interaction quantum Monte Carlo (FCIQMC) method, as well as its "initiator" extension (i-FCIQMC), is used to tackle the complex electronic structure of the carbon dimer across the entire dissociation reaction coordinate, as a prototypical example of a strongly correlated molecular system. Various basis sets of increasing size up to the large cc-pVQZ are used, spanning a fully accessible N-electron basis of over 10(12) Slater determinants, and the accuracy of the method is demonstrated in each basis set. Convergence to the FCI limit is achieved in the largest basis with only O[10(7)] walkers within random errorbars of a few tenths of a millihartree across the binding curve, and extensive comparisons to FCI, CCSD(T), MRCI, and CEEIS results are made where possible. A detailed exposition of the convergence properties of the FCIQMC methods is provided, considering convergence with elapsed imaginary time, number of walkers and size of the basis. Various symmetries which can be incorporated into the stochastic dynamic, beyond the standard abelian point group symmetry and spin polarisation are also described. These can have significant benefit to the computational effort of the calculations, as well as the ability to converge to various excited states. The results presented demonstrate a new benchmark accuracy in basis-set energies for systems of this size, significantly improving on previous state of the art estimates.     

Predicting atomic dopant solvation in helium clusters: the MgHe(n) case

We present a quantum Monte Carlo study of the solvation and spectroscopic properties of the Mg-doped helium clusters MgHe(n) with n=2-50. Three high-level [MP4, CCSD(T), and CCSDT] MgHe interaction potentials have been used to study the sensitivity of the dopant location on the shape of the pair interaction. Despite the similar MgHe well depth, the pair distribution functions obtained in the diffusion Monte Carlo simulations markedly differ for the three pair potentials, therefore indicating different solubility properties for Mg in He(n). Moreover, we found interesting size effects for the behavior of the Mg impurity. As a sensitive probe of the solvation properties, the Mg excitation spectra have been simulated for various cluster sizes and compared with the available experimental results. The interaction between the excited 1P Mg atom and the He moiety has been approximated using the diatomics-in-molecules method and the two excited 1pi and 1sigma MgHe potentials. The shape of the simulated MgHe50 spectra shows a substantial dependency on the location of the Mg impurity, and hence on the MgHe pair interaction employed. To unravel the dependency of the solvation behavior on the shape of the computed potentials, exact density-functional theory has been adapted to the case of doped He(n) and various energy distributions have been computed. The results indicate the shape of the repulsive part of the MgHe potential as an important cause of the different behaviors.     

The lithium–thiophene interaction: a critical study using highly correlated electronic structure approaches of quantum chemistry

The fundamental multicentric interaction of a lithium atom with a single thiophene ring is addressed. A systematic study of the interaction energy (IE) and geometry for the Li–T charge-transfer complex is done at the MP2 and CCSD(T) levels using increasingly large basis sets up to aug-cc-pVQZ (AVQZ). Basis set superposition errors (BSSE) are evaluated and shown to have a major impact on the value of the IE. The Fixed-Node Diffusion Monte Carlo (FN-DMC) method is used as an alternative basis-set-free approach to obtain what is likely to be the most accurate estimate of the IE obtained so far. While counterpoise-corrected MP2/AVQZ and CCSD(T)/AVTZ interaction energies are found to be ?3.8 and ?7.5 kcal/mol, the FN-DMC method yields +1.3 ± 1.7 kcal/mol. The slow convergence of the ab initio IE (and some key structural parameters) with respect to basis set quality and the discrepancy with the FN-DMC result is discussed. A visualization of the electron pairing using the electron pair localization function (EPLF) for the Li-doped versus undoped thiophene is also presented.     

The structures,binding energies and vibrational frequencies of Ca3 and Ca4 — An application of the CCSD(T) method

Summary The Ca₃ and Ca₄ metallic clusters have been investigated using state-of-the-artab initio quantum mechanical methods. Large atomic natural orbital basis sets have been used in conjunction with the singles and doubles coupled-cluster (CCSD) method, a coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), and the multireference configuration interaction (MRCI) method. The equilibrium geometries, binding energies and harmonic vibrational frequencies have been determined with each of the methods so that the accuracy of the coupled-cluster methods may be assessed. Since the CCSD(T) method reproduces the MRCI results very well, cubic and quartic force fields of Ca₃ and Ca₄ have been determined using this approach and used to evaluate the fundamental vibrational frequencies. The infrared intensities of both thee mode of Ca₃ and thet ₂ mode of Ca₄ are found to be small. The results obtained in this study are compared and contrasted with those from our earlier studies on small Be and Mg clusters.Dedicated to Prof. Klaus Ruedenberg on the occasion of his 70th birthday     

Auxiliary-field quantum Monte Carlo calculations of molecular systems with a Gaussian basis

We extend the recently introduced phaseless auxiliary-field quantum Monte Carlo (QMC) approach to any single-particle basis and apply it to molecular systems with Gaussian basis sets. QMC methods in general scale favorably with the system size as a low power. A QMC approach with auxiliary fields, in principle, allows an exact solution of the Schrodinger equation in the chosen basis. However, the well-known sign/phase problem causes the statistical noise to increase exponentially. The phaseless method controls this problem by constraining the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. In the present calculations, the trial wave function is a single Slater determinant from a Hartree-Fock calculation. The calculated all-electron total energies show typical systematic errors of no more than a few millihartrees compared to exact results. At equilibrium geometries in the molecules we studied, this accuracy is roughly comparable to that of coupled cluster with single and double excitations and with noniterative triples [CCSD(T)]. For stretched bonds in H(2)O, our method exhibits a better overall accuracy and a more uniform behavior than CCSD(T).     

Auxiliary-field quantum Monte Carlo study of first- and second-row post-d elements

A series of calculations for the first- and second-row post-d elements (Ga-Br and In-I) are presented using the phaseless auxiliary-field quantum Monte Carlo (AF QMC) method. This method is formulated in a Hilbert space defined by any chosen one-particle basis and maps the many-body problem into a linear combination of independent-particle solutions with external auxiliary fields. The phase/sign problem is handled approximately by the phaseless formalism using a trial wave function, which in our calculations was chosen to be the Hartree-Fock solution. We used the consistent correlated basis sets of Peterson et al. [J. Chem. Phys. 119, 11099 (2003); 119, 11113 (2003)], which employ a small-core relativistic pseudopotential. The AF QMC results are compared with experiment and with those from density functional (generalized gradient approximation and B3LYP) and CCSD(T) calculations. The AF QMC total energies agree with CCSD(T) to within a few millihartrees across the systems and over several basis sets. The calculated atomic electron affinities, ionization energies, and spectroscopic properties of dimers are, at large basis sets, in excellent agreement with experiment.     

Path integral ground state study of finite-size systems: application to small (parahydrogen)N (N=2-20) clusters

We use the path integral ground state method to study the energetic and structural properties of small para-H2 clusters of sizes ranging from 2 to 20 molecules. A fourth order formula is used to approximate the short imaginary-time propagator and two interaction potentials are considered. Our results are compared to those of exact basis set calculations and other quantum Monte Carlo methods when available. We find that for all cluster sizes considered, our results show a lower ground state energy than literature values obtained by diffusion Monte Carlo and variational Monte Carlo. For the dimer and trimer, ground state energies are in good agreement with exact results obtained using the discrete variable representation. Structural properties are found to be insensitive to the choice of interaction potential. We explore the use of Pekeris coordinates to analyze the importance of linear arrangement in trimers and for trimers within clusters of larger size.     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

Benzene-pyridine interactions predicted by the effective fragment potential method

The accurate representation of nitrogen-containing heterocycles is essential for modeling biological systems. In this study, the general effective fragment potential (EFP2) method is used to model dimers of benzene and pyridine, complexes for which high-level theoretical data -including large basis spin-component-scaled second-order perturbation theory (SCS-MP2), symmetry-adapted perturbation theory (SAPT), and coupled cluster with singles, doubles, and perturbative triples (CCSD(T))-are available. An extensive comparison of potential energy curves and components of the interaction energy is presented for sandwich, T-shaped, parallel displaced, and hydrogen-bonded structures of these dimers. EFP2 and CCSD(T) potential energy curves for the sandwich, T-shaped, and hydrogen-bonded dimers have an average root-mean-square deviation (RMSD) of 0.49 kcal/mol; EFP2 and SCS-MP2 curves for the parallel displaced dimers have an average RMSD of 0.52 kcal/mol. Additionally, results are presented from an EFP2 Monte Carlo/simulated annealing (MC/SA) computation to sample the potential energy surface of the benzene-pyridine and pyridine dimers.     

Vanadium oxide compounds with quantum Monte Carlo

Calculations with the diffusion quantum Monte Carlo method are presented for vanadium oxide molecules VO0/+0(n) with n = 1-4 and for V2O5. Atomization and ionization energies are calculated as well as oxygen abstraction energies. The fixed-node approximation is compared for guide functions with orbitals from B3LYP and BP86 calculations and higher accuracy was obtained with the latter orbitals. Additionally, all-electron and pseudopotential calculations are compared for the oxygen atom. The overall accuracy is found to be comparable to CCSD(T) calculations where experimental data is available.     

On the accuracy of explicitly correlated coupled-cluster interaction energies--have orbital results been beaten yet?

The basis set convergence of weak interaction energies for dimers of noble gases helium through krypton is studied for six variants of the explicitly correlated, frozen geminal coupled-cluster singles, doubles, and noniterative triples [CCSD(T)-F12] approach: the CCSD(T)-F12a, CCSD(T)-F12b, and CCSD(T)(F12*) methods with scaled and unscaled triples. These dimers were chosen because CCSD(T) complete-basis-set (CBS) limit benchmarks are available for them to a particularly high precision. The dependence of interaction energies on the auxiliary basis sets has been investigated and it was found that the default resolution-of-identity sets cc-pVXZ/JKFIT are far from adequate in this case. Overall, employing the explicitly correlated approach clearly speeds up the basis set convergence of CCSD(T) interaction energies, however, quite surprisingly, the improvement is not as large as the one achieved by a simple addition of bond functions to the orbital basis set. Bond functions substantially improve the CCSD(T)-F12 interaction energies as well. For small and moderate bases with bond functions, the accuracy delivered by the CCSD(T)-F12 approach cannot be matched by conventional CCSD(T). However, the latter method in the largest available bases still delivers the CBS limit to a better precision than CCSD(T)-F12 in the largest bases available for that approach. Our calculations suggest that the primary reason for the limited accuracy of the large-basis CCSD(T)-F12 treatment are the approximations made at the CCSD-F12 level and the non-explicitly correlated treatment of triples. In contrast, the explicitly correlated second-order Mo?ller-Plesset perturbation theory (MP2-F12) approach is able to pinpoint the complete-basis-set limit MP2 interaction energies of rare gas dimers to a better precision than conventional MP2. Finally, we report and analyze an unexpected failure of the CCSD(T)-F12 method to deliver the core-core and core-valence correlation corrections to interaction energies consistently and accurately.     

多维微正则蒙特卡罗方法计算异氰酸S0→T1系间窜跃速率

推导了通用的坐标空间多维蒙特卡罗微正则非绝热速率计算公式,并应用于计算异氰酸基5（S0）到第一激发三态（T1）的系间窜跃（ISC）速率．在CCSD水平上构建了一,二,三维势能面用于蒙特卡罗抽样．计算所得S0→T1 ISC速率与实验符合较好,因此可以预期多维蒙特卡罗微正则非绝热速率理论将成为计算多原子分子非绝热反应速率的有效手段之一．     

Quantum Monte Carlo calculations of the dissociation energies of three-electron hemibonded radical cationic dimers

We report variational and diffusion quantum Monte Carlo (VMC and DMC) calculations of the dissociation energies of the three-electron hemibonded radical cationic dimers of He, NH3, H2O, HF, and Ne. These systems are particularly difficult for standard density-functional methods such as the local-density approximation and the generalized gradient approximation. We have performed both all-electron (AE) and pseudopotential (PP) calculations using Slater-Jastrow wave functions with Hartree-Fock single-particle orbitals. Our results are in good agreement with coupled-cluster CCSD(T) calculations. We have also studied the relative stability of the hemibonded and hydrogen-bonded water radical dimer isomers. Our calculations indicate that the latter isomer is more stable, in agreement with post-Hartree-Fock methods. The excellent agreement between our AE and PP results demonstrates the high quality of the PPs used within our VMC and DMC calculations.