VBPCM: A Valence Bond Method that Incorporates a Polarizable Continuum Model

Valence Bond Theory-Its Birth, Struggles with Molecular Orbital Theory, Its Present State and Future Prospects

This essay describes the successive births of valence bond (VB) theory during 1916–1931. The alternative molecular orbital (MO) theory was born in the late 1920s. The presence of two seemingly different descriptions of molecules by the two theories led to struggles between the main proponents, Linus Pauling and Robert Mulliken, and their supporters. Until the 1950s, VB theory was dominant, and then it was eclipsed by MO theory. The struggles will be discussed, as well as the new dawn of VB theory, and its future.

Ab initio valence bond theory: A brief history, recent developments, and near future

This Perspective presents a survey of several issues in ab initio valence bond (VB) theory with a primary focus on recent advances made by the Xiamen VB group, including a brief review of the earlier history of the ab initio VB methods, in-depth discussion of algorithms for nonorthogonal orbital optimization in the VB self-consistent field method and VB methods incorporating dynamic electron correlation, along with a concise overview of VB methods for complex systems and VB models for chemical bonding and reactivity, and an outlook of opportunities and challenges for the near future of the VB theory.

Nonorthogonal orbital based n-body reduced density matrices and their applications to valence bond theory. III. Second-order perturbation theory using valence bond self-consistent field function as reference

Using the formulas and techniques developed in Papers I and II of this series, the recently developed second-order perturbation theory based on a valence bond self-consistent field reference function (VBPT2) has been extended by using the internally contracted correction wave function. This ansatz strongly reduces the size of the interaction space compared to the uncontracted wave function and thus improves the capability of the VBPT2 method dramatically. Test calculations show that internally contracted VBPT2 using only a small number of reference valence bond functions, can give results as accuracy as the VBPT2 method and other more sophisticated methods such as full configuration interaction and multireference configuration interaction.

How solvent influences the anomeric effect: roles of hyperconjugative versus steric interactions on the conformational preference

The block-localized wave function (BLW) method, which can derive optimal electron-localized state with intramolecular electron delocalization completely deactivated, has been combined with the polarizable continuum model (PCM) to probe the variation of the anomeric effect in solution. Currently both the hyperconjugation and electrostatic models have been called to interpret the anomeric effect in carbohydrate molecules. Here we employed the BLW-PCM scheme to analyze the energy differences between α and β anomers of substituted tetrahydropyran C5OH9Y (Y = F, Cl, OH, NH2, and CH3) and tetrahydrothiopyran C5SH9Y (Y = F, Cl, OH, and CH3) in solvents including chloroform, acetone, and water. In accord with literature, our computations show that for anomeric systems the conformational preference is reduced in solution and the magnitude of reduction increases as the solvent polarity increases. Significantly, on one hand the solute-solvent interaction diminishes the intramolecular electron delocalization in β anomers more than in α anomers, thus destabilizing β anomers relatively. But on the other hand, it reduces the steric effect in β anomers much more than α anomers and thus stabilizes β anomers relatively more, leading to the overall reduction of the anomeric effect in anomeric systems in solutions.

Perspective on Diabatic Models of Chemical Reactivity as Illustrated by the Gas-Phase S(N)2 Reaction of Acetate Ion with 1,2-Dichloroethane

Diabatic models are widely employed for studying chemical reactivity in condensed phases and enzymes, but there has been little discussion of the pros and cons of various diabatic representations for this purpose. Here we discuss and contrast six different schemes for computing diabatic potentials for a charge rearrangement reaction. They include (i) the variational diabatic configurations (VDC) constructed by variationally optimizing individual valence bond structures and (ii) the consistent diabatic configurations (CDC) obtained by variationally optimizing the ground-state adiabatic energy, both in the nonorthogonal molecular orbital valence bond (MOVB) method, along with the orthogonalized (iii) VDC-MOVB and (iv) CDC-MOVB models. In addition, we consider (v) the fourfold way (based on diabatic molecular orbitals and configuration uniformity), and (vi) empirical valence bond (EVB) theory. To make the considerations concrete, we calculate diabatic electronic states and diabatic potential energies along the reaction path that connects the reactant and the product ion-molecule complexes of the gas-phase bimolecular nucleophilic substitution (S(N)2) reaction of 1,2-dichloethane (DCE) with acetate ion, which is a model reaction corresponding to the reaction catalyzed by haloalkane dehalogenase. We utilize ab initio block-localized molecular orbital theory to construct the MOVB diabatic states and ab initio multi-configuration quasidegenerate perturbation theory to construct the fourfold-way diabatic states; the latter are calculated at reaction path geometries obtained with the M06-2X density functional. The EVB diabatic states are computed with parameters taken from the literature. The MOVB and fourfold-way adiabatic and diabatic potential energy profiles along the reaction path are in qualitative but not quantitative agreement with each other. In order to validate that these wave-function-based diabatic states are qualitatively correct, we show that the reaction energy and barrier for the adiabatic ground state, obtained with these methods, agree reasonably well with the results of high-level calculations using the composite G3SX and G3SX(MP3) methods and the BMC-CCSD multi-coefficient correlation method. However, a comparison of the EVB gas-phase adiabatic ground-state reaction path with those obtained from MOVB and with the fourfold way reveals that the EVB reaction path geometries show a systematic shift towards the products region, and that the EVB lowest-energy path has a much lower barrier. The free energies of solvation and activation energy in water reported from dynamical calculations based on EVB also imply a low activation barrier in the gas phase. In addition, calculations of the free energy of solvation using the recently proposed SM8 continuum solvation model with CM4M partial atomic charges lead to an activation barrier in reasonable agreement with experiment only when the geometries and the gas-phase barrier are those obtained from electronic structure calculations, i.e., methods i-v. These comparisons show the danger of basing the diabatic states on molecular mechanics without the explicit calculation of electronic wave functions. Furthermore, comparison of schemes i-v with one another shows that significantly different quantitative results can be obtained by using different methods for extracting diabatic states from wave function calculations, and it is important for each user to justify the choice of diabatization method in the context of its intended use.

On the construction of diabatic and adiabatic potential energy surfaces based on ab initio valence bond theory

A theoretical model is presented for deriving effective diabatic states based on ab initio valence bond self-consistent field (VBSCF) theory by reducing the multiconfigurational VB Hamiltonian into an effective two-state model. We describe two computational approaches for the optimization of the effective diabatic configurations, resulting in two ways of interpreting such effective diabatic states. In the variational diabatic configuration (VDC) method, the energies of the diabatic states are variationally minimized. In the consistent diabatic configuration (CDC) method, both the configuration coefficients and orbital coefficients are simultaneously optimized to minimize the adiabatic ground-state energy in VBSCF calculations. In addition, we describe a mixed molecular orbital and valence bond (MOVB) approach to construct the CDC diabatic and adiabatic states for a chemical reaction. Note that the VDC-MOVB method has been described previously. Employing the symmetric S(N)2 reaction between NH(3) and CH(3)NH(3)(+) as a test system, we found that the results from ab initio VBSCF and from ab initio MOVB calculations using the same basis set are in good agreement, suggesting that the computationally efficient MOVB method is a reasonable model for VB simulations of condensed phase reactions. The results indicate that CDC and VDC diabatic states converge, respectively, to covalent and ionic states as the molecular geometries are distorted from the minimum of the respective diabatic state along the reaction coordinate. Furthermore, the resonance energy that stabilizes the energy of crossing between the two diabatic states, resulting in the transition state of the adiabatic ground-state reaction, has a strong dependence on the overlap integral between the two diabatic states and is a function of both the exchange integral and the total diabatic ground-state energy.

On the interpretation of valence bond wavefunctions

Valence bond wavefunctions are naturally geared to the chemist's idea of chemical bonding. In a structure one may distinguish different electron pair bonds and possible radical character. A structure may correspond to a covalent bond, where all electrons are equally divided over the atoms, or may describe an excess charge on a discrete part of the molecule, which indicates ionic bond character. From the weights of the structure in a variationally optimised multi-structure valence bond wavefunction one may derive the importance of the different bonding types. The individual structures could be considered to represent the different physical situations. We explore this concept for simple diatomic molecules and for polyatomics, and we discuss the relation to Lewis structures. We show that the assumption of individual properties for the individual structures leads to inconsistencies.

Identity SN2 Reactions X−+CH3X→XCH3+X− (X=F, Cl, Br, and I) in Vacuum and in Aqueous Solution: A Valence Bond Study

The recently developed (L. Song, W. Wu, Q. Zhang, S. Shaik, J. Phys. Chem. A 2004, 108, 6017) valence bond method coupled with a polarized continuum model (VBPCM) has been applied to the identity SN2 reaction of halides in the gas phase and in aqueous solution. The barriers computed at the level of the breathing orbital VB method (P. C. Hiberty, J. P. Flament, E. Noizet, Chem. Phys. Lett. 1992, 189, 259), BOVB and VBPCM//BOVB, are comparable to CCSD(T) and CCSD(T)//PCM results and to experimentally derived barriers in solution (W. J. Albery, M. M. Kreevoy, Adv. Phys. Org. Chem. 1978, 16, 85). The reactivity parameters needed to apply the valence bond state correlation diagram (VBSCD) method (S. Shaik, J. Am. Chem. Soc. 1984, 106, 1227), were also determined by VB calculations. It has been shown that the reactivity parameters along with their semiempirical derivations provide a satisfactory qualitative and quantitative account of the barriers.

Valence bond theory for chemical dynamics

This essay provides a perspective on several issues in valence bond theory: the physical significance of semilocal bonding orbitals, the capability of valence bond concepts to explain systems with multireferences character, the use of valence bond theory to provide analytical representations of potential energy surfaces for chemical dynamics by the method of semiempirical valence bond potential energy surfaces (an early example of specific reaction parameters), by multiconfiguration molecular mechanics, by the combined valence bond-molecular mechanics method, and by the use of valence bond states as coupled diabatic states for describing electronically nonadiabatic processes (photochemistry). The essay includes both ab initio and semiempirical approaches.

A survey of recent developments in ab initio valence bond theory

Starting from the 1980s and onwards, Valence Bond theory has been enjoying renaissance that is characterized by the development of a growing number of ab initio methods, and by many applications to chemical reactivity and to the central paradigms of chemistry. Owing the increase of computational power of modern computers and to significant advances in the methodology, valence bond theory begins to offer a sound and attractive alternative to Molecular Orbital theory. This review aims at summarizing the most important developments of ab initio valence bond methods during the last two or three decades, and is primarily devoted to a description of what the various methods can actually achieve within their specific scopes and limitations. Key available softwares are surveyed.

Hybrid ab initio VB/MM Method − A Valence Bond Ride through Classical Landscapes

The development of a new hybrid (QM/MM) method, where the QM part is treated by ab initio valence bond (VB) theory is presented. This VB/MM method has the advantages of empirical VB (EVB) methodology but does not rely on empirical parameterization for the quantum part. The method implements embedding of the quantum region of each diabatic state separately, by treating the electrostatic interactions between QM and MM regions classically. Additionally, it assumes that changes of the off diagonal matrix element due to different environments are such that the overall resonance integral does not change. These assumptions are discussed in detail and the validity of the method is shown to be successful in three different bond dissociation processes, where bond dissociation as well as solvation energies compare very well with the experimental data.

XMVB: a program for ab initio nonorthogonal valence bond computations

An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available.