块定域波函数方法及其应用

提出了块定域波函数方法以定量分析分子内的电子定域现象或分子间的电荷传递效应。对于一个假想的严格定域的分子,我们通过将全部的电子和基轨道配分成几个子空间来构造其相应的波函数。其中每一个分子轨道只对某一个子空间展开,各子空间内的分子轨道相互正交,但不同子空间内的分子轨道间是非正交的。Hartree-Fock波函数和块定域波函数之间的能量之差即为分子内的电子定域能或分子间的电荷传递能。我们应用块定域波函数方法讨论了丁二烯分子中的旋转势垒。     

The generalized block-localized wavefunction method: a case study on the conformational preference and C-O rotational barrier of formic acid

A Lewis structure corresponding to the most stable electron-localized state is often used as a reference for the measure of electron delocalization effect in the valence bond (VB) theory. As the simplest variant of ab initio VB theory, the generalized block-localized wavefunction (BLW) method defines the wavefunction for an electron-localized state with block-localized orbitals without the orthogonalization constraint on different blocks. The validity of the method can be critically examined with experimental evidences. Here the BLW method has been applied to the investigation of the roles of both the π conjugation and σ hyperconjugation effects in the conformational preference of formic acid for the trans (Z) conformer over the cis (E) conformer. On one hand, our computations showed that the deactivation of the π conjugation or σ hyperconjugation has little impact on the Z-E energy gap, thus neither is decisive and instead the local dipole-dipole electrostatic interaction between the carbonyl and hydroxyl groups is the key factor determining the Z-E energy gap. On the other hand, the present study supported the conventional view that π conjugation is largely responsible for the C-O rotation barrier in formic acid, though the existence of hyperconjugative interactions in the perpendicular structure lowers the barrier considerably.     

Energy decomposition analysis based on a block-localized wavefunction and multistate density functional theory

An interaction energy decomposition analysis method based on the block-localized wavefunction (BLW-ED) approach is described. The first main feature of the BLW-ED method is that it combines concepts of valence bond and molecular orbital theories such that the intermediate and physically intuitive electron-localized states are variationally optimized by self-consistent field calculations. Furthermore, the block-localization scheme can be used both in wave function theory and in density functional theory, providing a useful tool to gain insights on intermolecular interactions that would otherwise be difficult to obtain using the delocalized Kohn-Sham DFT. These features allow broad applications of the BLW method to energy decomposition (BLW-ED) analysis for intermolecular interactions. In this perspective, we outline theoretical aspects of the BLW-ED method, and illustrate its applications in hydrogen-bonding and π-cation intermolecular interactions as well as metal-carbonyl complexes. Future prospects on the development of a multistate density functional theory (MSDFT) are presented, making use of block-localized electronic states as the basis configurations.     

Strength of Intramolecular Hydrogen Bonds

The concept of resonance-assisted hydrogen bonds(RAHBs)highlights the synergistic interplay between theπ-resonance and hydrogen bonding interactions.This concept has been well-accepted in academia and is widely used in practice.However,it has been argued that the seemingly enhanced intramolecular hydrogen bonding(IMHB)in unsaturated compounds may simply be a result of the constraints imposed by theσ-skeleton framework.Thus,it is crucial to estimate the strength of IMHBs.In this work,we used two approaches to probe the resonance effect and estimate the strength of the IMHBs in the two exemplary cases of the enol forms of acetylacetone and o-hydroxyacetophenone.One approach is the block-localized wavefunction(BLW)method,which is a variant of the ab initio valence bond(VB)theory.Using this approach,it is possible to derive the geometries and energetics with resonance shut down.The other approach is Edmiston’s truncated localized molecular orbital(TLMO)technique,which monitors the energy changes by removing the delocalization tails from localized molecular orbitals.The integrated BLW and TLMO studies confirmed that the hydrogen bonding in these two molecules is indeed enhanced byπ-resonance,and that this enhancement is not a result ofσconstraints.     

Two-state model based on the block-localized wave function method

The block-localized wave function (BLW) method is a variant of ab initio valence bond method but retains the efficiency of molecular orbital methods. It can derive the wave function for a diabatic (resonance) state self-consistently and is available at the Hartree-Fock (HF) and density functional theory (DFT) levels. In this work we present a two-state model based on the BLW method. Although numerous empirical and semiempirical two-state models, such as the Marcus-Hush two-state model, have been proposed to describe a chemical reaction process, the advantage of this BLW-based two-state model is that no empirical parameter is required. Important quantities such as the electronic coupling energy, structural weights of two diabatic states, and excitation energy can be uniquely derived from the energies of two diabatic states and the adiabatic state at the same HF or DFT level. Two simple examples of formamide and thioformamide in the gas phase and aqueous solution were presented and discussed. The solvation of formamide and thioformamide was studied with the combined ab initio quantum mechanical and molecular mechanical Monte Carlo simulations, together with the BLW-DFT calculations and analyses. Due to the favorable solute-solvent electrostatic interaction, the contribution of the ionic resonance structure to the ground state of formamide and thioformamide significantly increases, and for thioformamide the ionic form is even more stable than the covalent form. Thus, thioformamide in aqueous solution is essentially ionic rather than covalent. Although our two-state model in general underestimates the electronic excitation energies, it can predict relative solvatochromic shifts well. For instance, the intense pi-->pi* transition for formamide upon solvation undergoes a redshift of 0.3 eV, compared with the experimental data (0.40-0.5 eV).     

Hydrogen and Halogen Bonding in Homogeneous External Electric Fields: Modulating the Bond Strengths

Recent years have witnessed various fascinating phenomena arising from the interactions of noncovalent bonds with homogeneous external electric fields (EEFs). Here we performed a computational study to interpret the sensitivity of intrinsic bond strengths to EEFs in terms of steric effect and orbital interactions. The block-localized wavefunction (BLW) method, which combines the advantages of both ab initio valence bond (VB) theory and molecular orbital (MO) theory, and the subsequent energy decomposition (BLW-ED) approach were adopted. The sensitivity was monitored and analyzed using the induced energy term, which is the variation in each energy component along the EEF strength. Systems with single or multiple hydrogen (H) or halogen (X) bond(s) were also examined. It was found that the X-bond strength change to EEFs mainly stems from the covalency change, while generally the steric effect rules the response of H-bonds to EEFs. Furthermore, X-bonds are more sensitive to EEFs, with the key difference between H- and X-bonds lying in the charge transfer interaction. Since phenylboronic acid has been experimentally used as a smart linker in EEFs, switchable sensitivity was scrutinized with the example of the phenylboronic acid dimer, which exhibits two conformations with either antiparallel or parallel H-bonds, thereby, opposite or consistent responses to EEFs. Among the studied systems, the quadruple X-bonds in molecular capsules exhibit remarkable sensitivity, with its interaction energy increased by −95.2 kJ mol⁻¹ at the EEF strength 0.005 a.u.     

How to properly compute the resonance energy within the ab initio valence bond theory: a response to the ZHJVL paper

Valence bond (VB) theory describes a conjugated system by a set of electron-localized Lewis resonance structures. VB assumes that the magnitude of the intramolecular electron delocalization can be measured in terms of a resonance energy (RE), taken to be the energy difference between the real conjugated system (delocalized) and the corresponding most stable virtual resonance structure (localized). Proper RE estimates within VB theory require both delocalized and localized states to be defined at the same theoretical level, and the definition of the localized state to closely correspond to the intuitive picture of the corresponding VB structure. In contrast, the VB-delocal and VB-local computational approaches adopted by Zielinski, et al. [preceding paper in this issue] used definitions for either the delocalized or the localized states which, in our view, depart from the intuitive chemical picture. Consequently, their RE estimates are much lower than seemingly appropriate experimental evaluations with which they strongly disagree. Very large basis sets approaching completeness blur the boundaries among resonance structures and result in “basis set artifact” problems within any variant of VB theory. However, block-localized wavefunction (BLW) computations with mid-size basis sets not only exhibit insignificant variations with theoretical levels, but the resulting RE estimates also are justified by comparisons with those employing experimental data and MO computations. We stress that RE differs from the aromatic stabilization energy (ASE). The RE measures the total stabilization of an aromatic system, whereas ASE measures only the part of the RE that exceeds that of appropriate conjugated (but non-aromatic) reference molecules.     

Natural Bond Orbital Study on the Nature of Binding in the H2 Molecule

We use the natural bond orbital (NBO) method to decompose a MO wavefunction into the intuitive valence bond (VB) structures. At least two natural orbital type MO are required to describe the essential binding of the H₂ molecule at all inter nuclear distances. At first the MO wavefunction is transformed into an unrestricted Hartree-Fock wave-function consisted of non-orthogonal localized orbitals u' and v', and then the NBO method is used to decompose u' and v' into the physical meaningful orthogonal localized orbitals. Our results show that the orbitals u' and v' are decomposed into an atomic and an overlap parts. The latter part gives rise to the conventional ionic structure in the VB picture.     

Can QTAIM topological parameters be a measure of hydrogen bonding strength?

The block-localized wave function (BLW) method, which is the simplest variant of ab initio valence bond (VB) theory, together with the quantum theory of atoms in molecules (QTAIM) approach, have been used to probe the intramolecular hydrogen bonding interactions in a series of representative systems of resonance-assisted hydrogen bonds (RAHBs). RAHB is characteristic of the cooperativity between the π-electron delocalization and hydrogen bonding interactions and is identified in many biological systems. While the deactivation of the π resonance in these RAHB systems by the use of the BLW method is expected to considerably weaken the hydrogen bonding strength, little change on the topological properties of electron densities at hydrogen bond critical points (HBCPs) is observed. This raises a question of whether the QTAIM topological parameters can be an effective measure of hydrogen bond strength.     

Minimal Active Space for Diradicals Using Multistate Density Functional Theory

This work explores the electronic structure as well as the reactivity of singlet diradicals, making use of multistate density functional theory (MSDFT). In particular, we show that a minimal active space of two electrons in two orbitals is adequate to treat the relative energies of the singlet and triplet adiabatic ground state as well as the first singlet excited state in many cases. This is plausible because dynamic correlation is included in the first place in the optimization of orbitals in each determinant state via block-localized Kohn–Sham density functional theory. In addition, molecular fragment, i.e., block-localized Kohn–Sham orbitals, are optimized separately for each determinant, providing a variational diabatic representation of valence bond-like states, which are subsequently used in nonorthogonal state interactions (NOSIs). The computational procedure and its performance are illustrated on some prototypical diradical species. It is shown that NOSI calculations in MSDFT can be used to model bond dissociation and hydrogen-atom transfer reactions, employing a minimal number of configuration state functions as the basis states. For p- and s-types of diradicals, the closed-shell diradicals are found to be more reactive than the open-shell ones due to a larger diabatic coupling with the final product state. Such a diabatic representation may be useful to define reaction coordinates for electron transfer, proton transfer and coupled electron and proton transfer reactions in condensed-phase simulations.     

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.     

Adjacent Lone Pair (ALP) Effect: A Computational Approach for Its Origin

The adjacent lone pair (ALP) effect is an experimental phenomenon in certain nitrogenous heterocyclic systems exhibiting the preference of the products with lone pairs separated over other isomers with lone pairs adjacent. A theoretical elucidation of the ALP effect requires the decomposition of intramolecular energy terms and the isolation of lone pair–lone pair interactions. Here we used the block‐localized wavefunction (BLW) method within the ab initio valence bond (VB) theory to derive the strictly localized orbitals which are used to accommodate one‐atom centered lone pairs and two‐atom centered σ or π bonds. As such, interactions among electron pairs can be directly derived. Two‐electron integrals between adjacent lone pairs do not support the view that the lone pair–lone pair repulsion is responsible for the ALP effect. Instead, the disabling of π conjugation greatly diminishes the ALP effect, indicating that the reduction of π conjugation in deprotonated forms with two σ lone pairs adjacent is one of the major causes for the ALP effect. Further electrostatic potential analysis and intramolecular energy decomposition confirm that the other key factor is the favorable electrostatic attraction within the isomers with lone pairs separated.     

价键理论新进展

概要介绍了现代价键理论的几个主要方法,并讨论了它们各自的特点及其发展现状,并重点介绍了键表方法的基本理论、计算程序及一些应用。     

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.     

Is cyclobutadiene really highly destabilized by antiaromaticity?

The high energy of cyclobutadiene (CBD) is not due primarily to "anti-aromaticity," but rather to angle strain, torsional strain, and Pauli repulsion between the parallel CC bonds. Estimations including block-localized wavefunction (BLW) computations conclude that the enormous ring strain (ca. 60 kcal mol(-1)) far exceeds its antiaromatic destabilization (only 16.5 kcal mol(-1)).     

六氢吡啶和氨形成的氢键团簇C5H10NH(NH3)n(n=1～3)结构的从头算研究

在RHF/6-31G(d)水平下,对C5H10NH(NH3)n(n=1～3)氢键团簇的平衡构型进行了从头算研究,优化得到各种可能的平衡构型.C5H10NH(NH3)为线型氢键结构,而C5H10NH(NH3)2为三元环结构,C5H10NH(NH3)3为四元环结构.在MP2/6-31G(d)//B3LYP/6-31G(d)水平下,对最稳定构型C5H10NH(NH3)n(Ⅰ)(n=1～3)的分子轨道进行布居分析,并且对相应的占据轨道进行指认.C5H10NH(NH3)n(Ⅰ)(n=1～3)垂直电离势的计算结果表明,形成氢键团簇后,分子的垂直电离势降低.     

Lone-pair orbital interactions in polythiaadamantanes

The electronic structures of a series of polythiaadamantanes from thiaadamantane through 2,4,6,8,9,10-hexathiaadamantane (HTA) have been analyzed using density functional theory calculations in conjunction with Hückel and natural bond orbital analysis. The effects of multiple sulfur p-type lone-pair orbital interactions on ionization potentials, hole mobilities, and electronic coupling have been determined. An overall increase in the average energy of the lone-pair orbitals as the number of sulfur atoms increases is predicted, with the exact positioning of the HOMO depending on specific lone-pair interactions. Separation of through-bond (TB) and through-space (TS) interactions between intramolecular sulfur atoms has been performed using localized molecular orbitals and model systems based on interacting hydrogen sulfide molecules. TB interations were found to reduce orbital splitting, while TS interactions were found to increase orbital splitting. TS interactions were more or less constant from one polythiaadamantane to the next, and the contributions of TB effects to individual orbital energies vary depending on the relative orientation of sulfur atoms as determined by the sigma molecular framework. Electronic coupling between intermolecular sulfur lone-pair orbitals was determined by investigating unique dimer pairs observed in the crystal structure of HTA. Electronic coupling is not as strong as expected given the short intermolecular S-S distances observed in the crystal structure. In general, B3LYP/6-31G(d) and B3LYP/6-311+G(d,p) give very similar orbital energies and splittings.     

Molecular orbital studies of dinitrogen tetroxide and related molecules

Hartree-Fock LCAO MO calculations for N₂O₄ have been performed in a basis of symmetry orbitals formed from a minimal Slater basis set. Effects of rounding and truncation errors were minimized by the use of the symmetry basis, which also allowed the order or tilling of molecular orbitals to be specified independently of orbital energies. Convergence difficulties were overcome by combined use of the conjugate gradients method and Roothaan's iterative procedure; the method of steepest descents was less effective than either of these. Multicentre ‘non-NDDO’ two-electron integrals were evaluated by the gaussian expansion technique. The wavefunction obtained for the lowest state is NN antibonding, largely as a result of the filling of the 6b_1u antibonding sigma orbital in preference to the 6a_g bonding sigma orbital. There is only a small amount of NN pi-bonding. A bond energy analysis shows that the lowest state is markedly stabilized by NNO three-centre interactions.     

Characterization of pentacarbonyl(4-methylpyridine)chromium(0) complex using density functional theory (DFT) and Hartree-Fock (HF) computational methods

The optimized molecular structure, atomic charges, vibrational frequencies, thermodynamic properties, nuclear magnetic resonance (NMR) and ultraviolet-visible (UV-Vis) spectral data of pentacarbonyl(4-methylpyridine)chromium(0) complex have been investigated by performing ab initio Hartree-Fock (HF) and density functional theory, B3LYP, B3PW91 and BE1PBE methods with 6-311G, 6-311+G(3d,3p) and 6-31G(d,p) basis set. The calculated NMR data at 6-311G basis set, vibrational frequencies at 6-311+G(3d,3p) basis set and the optimized geometric bond lengths and bond angles at 6-31G(d,p) basis set are in good agreement with the corresponding experimental values. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) have been simulated. In addition, the transition state and energy band gap and infrared intensities have also been reported.