A spin–free approach for valence bond theory and its applications

The effect of immediate environment on bond strength of different bond types-A valence bond study

The ability to design catalysis largely depends on our understanding of the electrostatic effect of the surrounding on the bonds participating in the reaction. Here, we used a simplistic model of point charges (PCs) to determine a set of rules guiding how to construct PC-bond arrangement that can strengthen or weaken different chemical bonds. Using valence bond theory to calculate the in situ bond energies, we show that the effect of the PC mainly depends on the bond's dipole moment irrespective of its type (being covalent or charge shift). That is, polar bonds are getting stronger or weaker depending on the sign and location of the PC, whereas non- or weakly polar bonds become stronger or weaker depending only on the location of the PC and to a smaller extent compared with polar bonds. We also show that for polar bonds, the maximal bond strengthening and weakening effect can be achieved when the PC is placed along the bond axis, as close as possible to the more and less polarizable atom/fragment, respectively. Finally, due to the stabilizing effects of polarizability, we show that, overall, it is easier to cause bond strengthening compared with bond weakening. Particularly, for polar bonds, bond strengthening is larger than bond weakening obtained by an oppositely signed PC. These rules should be useful in the future design of catalysis in, e.g., enzyme active sites.

Assessing the Accuracy of Tailored Coupled Cluster Methods Corrected by Electronic Wave Functions of Polynomial Cost

Tailored coupled cluster theory represents a computationally inexpensive way to describe static and dynamical electron correlation effects. In this work, we scrutinize the performance of various coupled cluster methods tailored by electronic wave functions of polynomial cost. Specifically, we focus on frozen-pair coupled cluster (fpCC) methods, which are tailored by pair-coupled cluster doubles (pCCD), and coupled cluster theory tailored by matrix product state wave functions optimized by the density matrix renormalization group (DMRG) algorithm. As test system, we selected a set of various small- and medium-sized molecules containing diatomics (N₂, F₂, C₂, CN⁺, CO, BN, BO⁺, and Cr₂) and molecules (ammonia, ethylene, cyclobutadiene, benzene, hydrogen chains, rings, and cuboids) for which the conventional single-reference coupled cluster singles and doubles (CCSD) method is not able to produce accurate results for spectroscopic constants, potential energy surfaces, and barrier heights. Most importantly, DMRG-tailored and pCCD-tailored approaches yield similar errors in spectroscopic constants and potential energy surfaces compared to accurate theoretical and/or experimental reference data. Although fpCC methods provide a reliable description for the dissociation pathway of molecules featuring single and quadruple bonds, they fail in the description of triple or hextuple bond-breaking processes or avoided crossing regions.

Extended Mulliken-Hush Method with Applications to the Theoretical Study of Electron Transfer

A novel adiabatic-to-diabatic (ATD) transformation strategy, namely, the extended Mulliken-Hush (XMH) method, is proposed to evaluate diabatic properties including electronic couplings, potential energy surfaces, and their crossings. The XMH method is developed by adopting our recently proposed ATD transformation formula of a general vectorial physical observable, in which a useful ATD transformation is further determined by using an auxiliary dipole between localized frontier orbitals as a simple approximation of the diabatic transition dipole. The XMH method is simple and practical that provides a flexible way to construct diabatic states. To some extent, it can be regarded as an extension of the generalized Mulliken-Hush (GMH) method since the latter takes a stronger approximation, in which the diabatic transition dipole is assumed to be vanishing. Test calculations on the HeH₂⁺ system show that the electronic couplings predicted by the XMH method are closer to the ones calculated by the valence bond block-diagonalization approach than the GMH ones since the XMH method takes into account both the magnitude and direction of the diabatic transition dipole, which is consistent with the properties of this molecule. In the study of electron transfer in the two kinds of donor-bridge-acceptor systems, the XMH method maintains the simplicity of the GMH method and gives reasonable results even when the latter fails, wherein the diabatic transition dipole is nearly perpendicular to the difference of the initial and final adiabatic dipoles. More importantly, the XMH method can be easily combined with high-level electronic structure methods, in which the properties of the ground and excited states may be more accurately calculated, and hence, one may expect that further development of the XMH method would result in a general computational model for studying electron transfer reactions.

High-level coupled-cluster energetics by Monte Carlo sampling and moment expansions: Further details and comparisons

We recently proposed a novel approach to converging electronic energies equivalent to high-level coupled-cluster (CC) computations by combining the deterministic CC(P;Q) formalism with the stochastic configuration interaction (CI) and CC Quantum Monte Carlo (QMC) propagations. This article extends our initial study [J. E. Deustua, J. Shen, and P. Piecuch, Phys. Rev. Lett. 119, 223003 (2017)], which focused on recovering the energies obtained with the CC method with singles, doubles, and triples (CCSDT) using the information extracted from full CI QMC and CCSDT-MC, to the CIQMC approaches truncated at triples and quadruples. It also reports our first semi-stochastic CC(P;Q) calculations aimed at converging the energies that correspond to the CC method with singles, doubles, triples, and quadruples (CCSDTQ). The ability of the semi-stochastic CC(P;Q) formalism to recover the CCSDT and CCSDTQ energies, even when electronic quasi-degeneracies and triply and quadruply excited clusters become substantial, is illustrated by a few numerical examples, including the F-F bond breaking in F₂, the automerization of cyclobutadiene, and the double dissociation of the water molecule.

Covalent Excited States of Polyenes C2nH2n+2 (n = 2-8) and Polyenyl Radicals C2n-1H2n+1 (n = 2-8): An Ab Initio Valence Bond Study

The ab initio valence bond (VB) methods, VBSCF and VBCI, are applied to the ground states and the covalent excited states of polyenes C2nH2n+2 (n = 2-8) and polyenyl radicals C2n-1H2n+1 (n = 2-8). The excitation energy gap was computed at the ab initio VB level, which is in good agreement with the semiempirical VB method, VBDFT(s), and the experimental values as well as with the molecular orbital theory based methods, CASPT3 and MRCI. The ab initio VB wave functions of systems are also in very good agreement with those of the VBDFT(s) method, even though the former is based on the ab initio VB scheme while the latter is a semiempirical Hückel type method, in which no orbital optimization procedure is performed. The computational results show that the ab initio VB method is capable now of providing numerical accuracy not only for bond forming and breaking processes, as shown in the past, but also for excitation energies, as shown here. In addition, the computational results validate the efficiency of the VBDFT(s) method, which is a simple VB model with less computational effort but which provides intuitive insights into the excited states of conjugated molecules.

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.

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.

Combining active-space coupled-cluster methods with moment energy corrections via the CC(P;Q) methodology, with benchmark calculations for biradical transition states

We have recently suggested the CC(P;Q) methodology that can correct energies obtained in the active-space coupled-cluster (CC) or equation-of-motion (EOM) CC calculations, which recover much of the nondynamical and some dynamical electron correlation effects, for the higher-order, mostly dynamical, correlations missing in the active-space CC/EOMCC considerations. It is shown that one can greatly improve the description of biradical transition states, both in terms of the resulting energy barriers and total energies, by combining the CC approach with singles, doubles, and active-space triples, termed CCSDt, with the CC(P;Q)-style correction due to missing triple excitations defining the CC(t;3) approximation.

Valence bond perturbation theory. A valence bond method that incorporates perturbation theory

A post-VBSCF method, called valence bond second-order perturbation theory (VBPT2), is developed in this paper and is shown to be (i) economical and (ii) at par with more sophisticated VB and MO-based methods. The VBPT2 method starts with VBSCF using a minimal structure set. Subsequently, the Møller-Plesset (MP) partition of the zeroth-order Hamiltonian is obtained by introducing a generalized Fock matrix constructed from the VBSCF density matrix. The first-order wave function is expressed in terms of singly and doubly excited VB structures, which are generated by replacing occupied orbitals by virtual orbitals, the latter being defined as orthogonal to the occupied orbitals. The VBPT2 method retains the simplicity of a VB presentation by condensing contributions from the excited structures into the minimal number of fundamental structures that are involved in the VBSCF calculation. The method is tested by calculating the bond energies of H(2), F(2), N(2), O(2), the barrier of identity hydrogen abstraction reaction, the atomization energy and a potential energy curve for the water molecule and the structural weights and covalent-ionic resonance energy of F(2). It is shown that the VBPT2 method gives results in good agreement with those of the VBCI method and molecular-orbital based methods such as MRPT and MRCI at the same truncation levels. However, the computational effort is greatly reduced, compared to that of VBCI. Future potential directions for the development of the VBPT2 method are outlined.

Topology of electron charge density for chemical bonds from valence bond theory: a probe of bonding types

To characterize the nature of bonding we derive the topological properties of the electron charge density of a variety of bonds based on ab initio valence bond methods. The electron density and its associated Laplacian are partitioned into covalent, ionic, and resonance components in the valence bond spirit. The analysis provides a density-based signature of bonding types and reveals, along with the classical covalent and ionic bonds, the existence of two-electron bonds in which most of the bonding arises from the covalent-ionic resonance energy, so-called charge-shift bonds. As expected, the covalent component of the Laplacian at the bond critical point is found to be largely negative for classical covalent bonds. In contrast, for charge-shift bonds, the covalent part of the Laplacian is small or positive, in agreement with the weakly attractive or repulsive character of the covalent interaction in these bonds. On the other hand, the resonance component of the Laplacian is always negative or nearly zero, and it increases in absolute value with the charge-shift character of the bond, in agreement with the decrease of kinetic energy associated with covalent-ionic mixing. A new interpretation of the topology of the total density at the bond critical point is proposed to characterize covalent, ionic, and charge-shift bonding from the density point of view.

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 BLW method is particularly useful in the quantification of the electron delocalization (resonance) effect within a molecule and the charge-transfer effect between molecules. In this paper, we extend the BLW method to the density functional theory (DFT) level and implement the BLW-DFT method to the quantum mechanical software GAMESS. Test applications to the pi conjugation in the planar allyl radical and ions with the basis sets of 6-31G(d), 6-31+G(d), 6-311+G(d,p), and cc-pVTZ show that the basis set dependency is insignificant. In addition, the BLW-DFT method can also be used to elucidate the nature of intermolecular interactions. Examples of pi-cation interactions and solute-solvent interactions will be presented and discussed. By expressing each diabatic state with one BLW, the BLW method can be further used to study chemical reactions and electron-transfer processes whose potential energy surfaces are typically described by two or more diabatic states.

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. An accurate quantitative description for the molecule requires higher levels of VB methods that incorporate dynamic correlation. The potential energy curves of the molecule are computed at the various VB levels. It is shown that there exists a small hump in the PECs of VBSCF for the ground state, as found in previous studies. However, higher levels of VB methods dissolve the hump. The BOVB and VBCI methods reproduce the dissociation energies and other physical properties of the ground state and the two lowest excited states in very good agreement with experiment and with sophisticated MO based methods, such as the MRCI method.

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.

A practical valence bond method: a configuration interaction method approach with perturbation theoretic facility

The previously developed valence bond configuration interaction (VBCI) method (Wu, W.; Song, L.; Cao, Z.; Zhang, Q.; Shaik, S., J. Phys. Chem. A, 2002, 105, 2721) that borrows the general CI philosophy of the MO theory, is further extended in this article, and its methodological features are improved, resulting in three accurate and cost-effective procedures: (a) the effect of quadruplet excitation is incorporated using the Davidson correction, such that the new procedure reduces size consistency problems, with due improvement in the quality of the computational results. (b) A cost-effective procedure, named VBCI(D, S), is introduced. It includes doubly excited structures for active electrons and singly excited structures for inactive pairs. The computational results of VBCI(D, S) match those of VBCISD with much less computational effort than VBCISD. (c) Finally, a second-order perturbation theory is utilized as a means of configuration selection, and lead to considerable reduction of the computational cost, with little or no loss in accuracy. Applications of the new procedures to bond energies and barriers of chemical reactions are presented and discussed.

Charge Transfer in the Electron Donor−Acceptor Complex BH3NH3

As a simple yet strongly binding electron donor-acceptor (EDA) complex, BH(3)NH(3) serves as a good example to study the electron pair donor-acceptor complexes. We employed both the ab initio valence bond (VB) and block-localized wave function (BLW) methods to explore the electron transfer from NH(3) to BH(3). Conventionally, EDA complexes have been described by two diabatic states: one neutral state and one ionic charge-transferred state. Ab initio VB self-consistent field (VBSCF) computations generate the energy profiles of the two diabatic states together with the adiabatic (ground) state. Our calculations evidently demonstrated that the electron transfer between NH(3) and BH(3) falls in the abnormal regime where the reorganization energy is less than the exoergicity of the reaction. The nature of the NH(3)-BH(3) interaction is probed by an energy decomposition scheme based on the BLW method. We found that the variation of the charge-transfer energy with the donor-acceptor distance is insensitive to the computation levels and basis sets, but the estimation of the amount of electron transferred heavily depends on the population analysis procedures. The recent resurgence of interest in the nature of the rotation barrier in ethane prompted us to analyze the conformational change of BH(3)NH(3), which is an isoelectronic system with ethane. We found that the preference of the staggered structure over the eclipsed structure of BH(3)NH(3) is dominated by the Pauli exchange repulsion.