ik Domin

Zori 1.0: A parallel quantum Monte Carlo electronic structure package

The Zori 1.0 package for electronic structure computations is described. Zori performs variational and diffusion Monte Carlo computations as well as correlated wave function optimization. This article presents an overview of the implemented methods and code capabilities.

A Diffusion Monte Carlo Study of the O−H Bond Dissociation of Phenol†

The homolytic O-H bond dissociation energy (BDE) of phenol was determined from diffusion Monte Carlo (DMC) calculations using single determinant trial wave functions. DMC gives an O-H BDE of 87.0 +/- 0.3 kcal/mol when restricted Hartree-Fock orbitals are used and a BDE of 87.5 +/- 0.3 kcal/mol with restricted B3LYP Kohn-Sham orbitals. These results are in good agreement with the extrapolated B3P86 results of Costa Cabral and Canuto (88.3 kcal/mol), the recommended experimental value of Borges dos Santos and Martinho Simões (88.7 +/- 0.5 kcal/mol), and the G3 (88.2 kcal/mol), CBS-APNO (88.2 kcal/mol), CBS-QB3 (87.1 kcal/mol) results of Mulder.

Multicenter Bonding in Ditetracyanoethylene Dianion: A Simple Aromatic Picture in Terms of Three-Electron Bonds.

The nature of the multicenter, long bond in ditetracyanoethylene dianion complex [TCNE]2(2-) is elucidated using high level ab initio Valence Bond (VB) theory coupled with Quantum Monte Carlo (QMC) methods. This dimer is the prototype of the general family of pancake-bonded dimers with large interplanar separations. Quantitative results obtained with a compact wave function in terms of only six VB structures match the reference CCSD(T) bonding energies. Analysis of the VB wave function shows that the weights of the VB structures are not compatible with a covalent bond between the π* orbitals of the fragments. On the other hand, these weights are consistent with a simple picture in terms of two resonating bonding schemes, one displaying a pair of interfragment three-electron σ bonds and the other displaying intrafragment three-electron π bonds. This simple picture explains at once (1) the long interfragment bond length, which is independent of the countercations but typical of three-electron (3-e) CC σ bonds, (2) the interfragment orbital overlaps which are very close to the theoretical optimal overlap of 1/6 for a 3-e σ bond, and (3) the unusual importance of dynamic correlation, which is precisely the main bonding component of 3-e bonds. Moreover, it is shown that the [TCNE]2(2-) system is topologically equivalent to the square C4H4(2-) dianion, a well-established aromatic system. To better understand the role of the cyano substituents, the unsubstituted diethylenic Na(+)2[C2H4]2(2-) complex is studied and shown to be only metastable and topologically equivalent to a rectangular C4H4(2-) dianion, devoid of aromaticity.

Quantum Monte Carlo study of first-row atoms using transcorrelated variational Monte Carlo trial functions

The effect of using the transcorrelated variational Monte Carlo (TC-VMC) approach to construct a trial function for fixed node diffusion Monte Carlo (DMC) energy calculations has been investigated for the first-row atoms, Li to Ne. The computed energies are compared with fixed node DMC energies obtained using trial functions constructed from Hartree-Fock and density functional levels of theory. Despite major VMC energy improvement with TC-VMC trial functions, no improvement in DMC energy was observed using these trial functions for the first-row atoms studied. The implications of these results on the nodes of the trial wave functions are discussed.

Breathing Orbital Valence Bond Method in Diffusion Monte Carlo: C-H Bond Dissociation of Acetylene †

This study explores the use of breathing orbital valence bond (BOVB) trial wave functions for diffusion Monte Carlo (DMC). The approach is applied to the computation of the carbon-hydrogen (C-H) bond dissociation energy (BDE) of acetylene. DMC with BOVB trial wave functions yields a C-H BDE of 132.4 +/- 0.9 kcal/mol, which is in excellent accord with the recommended experimental value of 132.8 +/- 0.7 kcal/mol. These values are to be compared with DMC results obtained with single determinant trial wave functions, using Hartree-Fock orbitals (137.5 +/- 0.5 kcal/mol) and local spin density (LDA) Kohn-Sham orbitals (135.6 +/- 0.5 kcal/mol).

Isomer Energy Differences for the C4H3 and C4H5 Isomers Using Diffusion Monte Carlo

A new diffusion Monte Carlo study is performed on the isomers of C4H3 and C4H5 emulating the methodology of a previous study (Int. J. Chem. Kinet. 2001, 33, 808). Using the same trial wave function form of the previous study, substantially different isomerization energies were found owing to the use of larger walker populations in the present work. The energy differences between the E and i isomers of C4H3 were found to be 10.5 +/- 0.5 kcal/mol and for C4H5, 9.7 +/- 0.6 kcal/mol. These results are in reasonable accord with recent MRCI and CCSD(T) findings.

Quantitative characteristics of qualitative localized bonding patterns.

Patterns of localized and delocalized chemical bonding obtained using the recently proposed adaptive natural density partitioning (AdNDP) provide a qualitative description of electronic structure but miss any quantitative information. Descriptors, such as the electron localization function (ELF), provide quantitative characteristics of bonding and can enhance the usefulness of qualitative patterns. In the present study, we used ELF and a related construct, charge-density-weighted ELF (ELF rho), to characterize localized and delocalized bonding in a variety of systems. It is demonstrated that ELF rho yields a more detailed description than ELF when used to analyze bonding in aromatic, conflicting aromatic, and antiaromatic systems. Both canonical molecular orbitals (CMOs) and localized multicenter two-electron (nc-2e) bonds obtained in the latter case by AdNDP localization are used to calculate ELF rho.

Influence of Water on the Oxidation of Dimethyl Sulfide by the ·OH Radical

Oxidative stress of sulfur-containing biological molecules in aqueous environments may lead to the formation of adduct intermediates that are too short-lived to be experimentally detectable. In this study we have modeled the simplest of such oxidative reactions: the attack of dimethyl sulfide (DMS) by a hydroxyl radical (·OH) to form a radical adduct, whose subsequent heterolytic dissociation leads to a radical cation (DMS+) that is important for further reactions. We have modeled the aqueous environment with a limited number of discrete water molecules, selected after an original multistep procedure, and further embedded in a polarizable continuum model, to observe the impact of the water configuration on the heterolytic dissociation of the radical adduct. Molecular dynamics and quantum chemical methods (DFT, MP2, and CCSD) were used to elucidate the lowest energy structures resulting from the ·OH attack on DMS. Subsequent high level ab initio valence bond (BOVB) calculations revealed the possibility for the occurrence of subsequent heterolytic dissociation.

Some Recent Developments of Quantum Monte Carlo to Molecular Systems.

A brief review of recent advances in quantum Monte Carlo for the electronic structure of molecules at the University of California, Berkeley, is given based on an invited talk presented at the Workshop on Nuclear and Mesoscopic Physics.

Recent Developments in Quantum Monte Carlo: Methods and Applications

The quantum Monte Carlo method in the diffusion Monte Carlo form has become recognized for its capability of describing the electronic structure of atomic, molecular and condensed matter systems to high accuracy. This talk will briefly outline the method with emphasis on recent developments connected with trial function construction, linear scaling, and applications to selected systems.