Block-Localized Wavefunction (BLW) Based Two-State Approach for Charge Transfers between Phenyl Rings.

The block-localized wave function (BLW) method is the simplest and most efficient variant of ab initio valence bond (VB) theory which defines electron-localized resonance states following the conventional VB concepts. Here, a BLW-based two-state approach is proposed to probe the charge/hole transfer reactions within the Marcus-Hush model. With this approach, both the electronic coupling and reorganization energies can be derived at the ab initio level.

On the Interfragment Exchange in the X-Pol Method.

The inclusion of exchange repulsion terms in the explicit polarization (X-Pol) model is examined by antisymmetrizing the X-Pol Hartree-product wave function; this yields XPol with full eXchange, called X-Pol-X. When the monomers are treated by Hartree-Fock theory, this calculation can be accomplished by using the formalism of block-localized wave functions (BLW) that has been used in a variety of applications. In this case the block-localized structure in the X-Pol-X wave function allows for decomposition of the full Fock matrix of a dimension of M blocks into M smaller Fock matrices.

A Non-Orthogonal Block-Localized Effective Hamiltonian Approach for Chemical and Enzymatic Reactions.

The effective Hamiltonian-molecular orbital and valence bond (EH-MOVB) method based on non-orthogonal block-localized fragment orbitals has been implemented into the program CHARMM for molecular dynamics simulations of chemical and enzymatic reactions, making use of semiempirical quantum mechanical models. Building upon ab initio MOVB theory, we make use of two parameters in the EH-MOVB method to fit the barrier height and the relative energy between the reactant and product state for a given chemical reaction to be in agreement with experiment or high-level ab initio or density functional results. Consequently, the EH-MOVB method provides a highly accurate and computationally efficient QM/MM model for dynamics simulation of chemical reactions in solution.

Block-Localized Density Functional Theory (BLDFT), Diabatic Coupling, and Their Use in Valence Bond Theory for Representing Reactive Potential Energy Surfaces.

A multistate density functional theory in the framework of the valence bond model is described. The method is based on a block-localized density functional theory (BLDFT) for the construction of valence-bond-like diabatic electronic states and is suitable for the study of electron transfer reactions and for the representation of reactive potential energy surfaces. The method is equivalent to a valence bond theory with the treatment of the localized configurations by using density functional theory (VBDFT).

An Effective Hamiltonian Molecular Orbital-Valence Bond (MOVB) Approach for Chemical Reactions Applied to the Nucleophilic Substitution Reaction of Hydrosulfide Ion and Chloromethane.

An effective Hamiltonian mixed molecular orbital and valence bond (EH-MOVB) method is described to obtain an accurate potential energy surface for chemical reactions. Building upon previous results on the construction of diabatic and adiabatic potential surfaces using ab initio MOVB theory, we introduce a diabatic-coupling scaling factor to uniformly scale the ab initio off-diagonal matrix element H(12) such that the computed energy of reaction from the EH-MOVB method is in agreement with the target value. The scaling factor is very close to unity, resulting in minimal alteration of the potential energy surface of the original MOVB model.

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.

Explicit polarization (X-Pol) potential using ab initio molecular orbital theory and density functional theory.

The explicit polarization (X-Pol) method has been examined using ab initio molecular orbital theory and density functional theory. The X-Pol potential was designed to provide a novel theoretical framework for developing next-generation force fields for biomolecular simulations. Importantly, the X-Pol potential is a general method, which can be employed with any level of electronic structure theory.

Molecular dynamics simulations of the detoxification of paraoxon catalyzed by phosphotriesterase.

Combined QM(PM3)/MM molecular dynamics simulations together with QM(DFT)/MM optimizations for key configurations have been performed to elucidate the enzymatic catalysis mechanism on the detoxification of paraoxon by phosphotriesterase (PTE). In the simulations, the PM3 parameters for the phosphorous atom were reoptimized. The equilibrated configuration of the enzyme/substrate complex showed that paraoxon can strongly bind to the more solvent-exposed metal ion Zn(beta), but the free energy profile along the binding path demonstrated that the binding is thermodynamically unfavorable.

Incorporation of a QM/MM buffer zone in the variational double self-consistent field method.

The explicit polarization (X-Pol) potential is an electronic-structure-based polarization force field, designed for molecular dynamics simulations and modeling of biopolymers. In this approach, molecular polarization and charge transfer effects are explicitly treated by a combined quantum mechanical and molecular mechanical (QM/MM) scheme, and the wave function of the entire system is variationally optimized by a double self-consistent field (DSCF) method. In the present article, we introduce a QM buffer zone for a smooth transition from a QM region to an MM region.

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.

An efficient algorithm for energy gradients and orbital optimization in valence bond theory.

An efficient algorithm for energy gradients in valence bond theory with nonorthogonal orbitals is presented. A general Hartree-Fock-like expression for the Hamiltonian matrix element between valence bond (VB) determinants is derived by introducing a transition density matrix. Analytical expressions for the energy gradients with respect to the orbital coefficients are obtained explicitly, whose scaling for computational cost is m(4), where m is the number of basis functions, and is thus approximately the same as in HF method.

The variational explicit polarization potential and analytical first derivative of energy: Towards a next generation force field.

A previous article proposed an electronic structure-based polarizable potential, called the explicit polarization (X-POL) potential, to treat many-body polarization and charge delocalization effects in polypeptides. Here, we present a variational version of the X-POL potential, in which the wave function of the entire molecular system is variationally optimized to yield the minimum total electronic energy. This allows the calculation of analytic gradients, a necessity for efficient molecular dynamics simulations.

Heterolytic bond dissociation in water: why is it so easy for C4H9Cl but not for C3H9SiCl?

The recently developed (Song, L.; Wu, W.; Zhang, Q.

Block-localized wavefunction (BLW) method at the density functional theory (DFT) level.

The block-localized wavefunction (BLW) approach is an ab initio valence bond (VB) method incorporating the efficiency of molecular orbital (MO) theory. It can generate the wavefunction for a resonance structure or diabatic state self-consistently by partitioning the overall electrons and primitive orbitals into several subgroups and expanding each block-localized molecular orbital in only one subspace. Although block-localized molecular orbitals in the same subspace are constrained to be orthogonal (a feature of MO theory), orbitals between different subspaces are generally nonorthogonal (a feature of VB theory).

The physical origin of large covalent-ionic resonance energies in some two-electron bonds.

This study uses valence bond (VB) theory to analyze in detail the previously established finding that alongside the two classical bond families of covalent and ionic bonds, which describe the electron-pair bond, there exists a distinct class of charge-shift bonds (CS-bonds) in which the fluctuation of the electron pair density plays a dominant role. Such bonds are characterized by weak binding, or even a repulsive, covalent component, and by a large covalent-ionic resonance energy RE(cs) that is responsible for the major part, or even for the totality, of the bonding energy. In the present work, the nature of CS-bonding and its fundamental mechanisms are analyzed in detail by means of a VB study of some typical homonuclear bonds (H-H, H3C-CH3, H2N-NH2, HO-OH, F-F, and Cl-Cl), ranging from classical-covalent to fully charge-shift bonds.

A valence bond study of the dioxygen molecule.

The dioxygen molecule has been the subject of valence bond (VB) studies since 1930s, as it was considered as the first "failure" of VB theory. The object of this article is to provide an unambiguous VB interpretation for the nature of chemical bonding of the molecule by means of modern VB computational methods, VBSCF, BOVB, and VBCI. It is shown that though the VBSCF method can not provide quantitative accuracy for the strongly electronegative and electron-delocalized molecule because of the lack of dynamic correlation, it still gives a correct qualitative analysis for wave function of the molecule and provides intuitive insights into chemical bonding.

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.

Steric strain versus hyperconjugative stabilization in ethane congeners.

Rotation barriers in the group IVB ethane congeners H(3)X-YH(3) (X, Y = C, Si, Ge, Sn, Pb) have been systematically studied and deciphered using the ab initio valence bond theory in terms of the steric strain and hyperconjugation effect. Our results show that in all cases the rotation barriers are dominated by the steric repulsion whereas the hyperconjugative interaction between the X-H bond orbitals and the vicinal Y-H antibond orbitals (and vice versa) plays a secondary role, although indeed the hyperconjugation effect favors staggered structures. By the independent estimations of the hyperconjugative and steric interactions in the process of rotations, we found that the structural effect which mainly refers to the central X-Y bond relaxation makes a small contribution to the rotational barriers.

Barriers of hydrogen abstraction vs halogen exchange: an experimental manifestation of charge-shift bonding.

This paper shows that the differences between the barriers of the halogen exchange reactions, in the H + XH systems, and the hydrogen abstraction reactions, in the X + HX systems (X = F, Cl, Br), measure the covalent-ionic resonance energies of the corresponding X-H bonds. These processes are investigated using CCSD(T) calculations as well as the breathing-orbital valence bond (BOVB) method. Thus, the VB analysis shows that (i) at the level of covalent structures the barriers are the same for the two series and (ii) the higher barriers for halogen exchange processes originate solely from the less efficient mixing of the ionic structures into the respective covalent structures.

Origins of Rotational Barriers in Hydrogen Peroxide and Hydrazine.

Compared with their isoelectronic system ethane, both hydrogen peroxide and hydrazine exhibit a double well torsional energy curve where skew conformers are favored over trans conformers and cis conformers are energy-maximum states. Clearly, the involvement of the lone oxygen and nitrogen pairs, or more specifically, the enhanced stabilizing n→σ* negative hyperconjugation effect and destabilizing repulsion among lone pairs, complicates the conformational analysis. In this work, the modern ab initio valence bond (VB) method is employed to quantitatively investigate the torsional energy curves of hydrogen peroxide and hydrazine in terms of hyperconjugative stabilization, steric repulsion, and structural and electronic relaxations.

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.

Valence bond calculations of hydrogen transfer reactions: a general predictive pattern derived from theory.

Hydrogen abstraction reactions of the type X(*) + H-H' --> X-H + H'(*) (X = F, Cl, Br, I) are studied by ab initio valence bond methods and the VB state correlation diagram (VBSCD) model. The reaction barriers and VB parameters of the VBSCD are computed by using the breathing orbital valence bond and valence bond configuration interaction methods. The combination of the VBSCD model and semiempirical VB theory leads to analytical expressions for the barriers and other VB quantities that match the ab initio VB calculations fairly well.

The ground and excited states of polyenyl radicals C2n-1H2n + 1 (n = 2-13): a valence bond study.

The semiempirical valence bond (VB) method, VBDFT(s), is applied to the ground states and the covalent excited states of polyenyl radicals C2n - 1H2n + 1 (n = 2-13). The method uses a single scalable parameter with a value that carries over from the study of the covalent excited states of polyenes (W. Wu, D.