Both proper, red-shifting and improper, blue-shifting hydrogen bonds have been well-recognized with enormous experimental and computational studies. The current consensus is that there is no difference in nature between these two kinds of hydrogen bonds, where the electrostatic interaction dominates. Since most if not all the computational studies are based on molecular orbital theory, it would be interesting to gain insight into the hydrogen bonds with modern valence bond (VB) theory. In this work, we performed ab initio VBSCF computations on a series of hydrogen-bonding systems, where the sole hydrogen bond donor CF3H interacts with ten hydrogen bond acceptors Y (═NH2CH3, NH3, NH2Cl, OH(-), H2O, CH3OH, (CH3)2O, F(-), HF, or CH3F). This series includes four red-shifting and six blue-shifting hydrogen bonds. Consistent with existing findings in literature, VB-based energy decomposition analyses show that electrostatic interaction plays the dominating role and polarization plays the secondary role in all these hydrogen-bonding systems, and the charge transfer interaction, which denotes the hyperconjugation effect, contributes only slightly to the total interaction energy. As VB theory describes any real chemical bond in terms of pure covalent and ionic structures, our fragment interaction analysis reveals that with the approaching of a hydrogen bond acceptor Y, the covalent state of the F3C-H bond tends to blue-shift, due to the strong repulsion between the hydrogen atom and Y. In contrast, the ionic state F3C(-) H(+) leads to the red-shifting of the C-H vibrational frequency, owing to the attraction between the proton and Y. Thus, the relative weights of the covalent and ionic structures essentially determine the direction of frequency change. Indeed, we find the correlation between the structural weights and vibrational frequency changes.

A novel valence bond based automatic diabatisation method by compression, called valence bond based compression approach for dibatization (VBCAD), is presented in this letter. It is a "blackbox" type method that provides an automatic diabatisation from a classical valence bond (VB) perspective. In VBCAD, a model space projection is performed by an eigenvalue decomposition algorithm followed by dimensional reduction based on a sequence of Householder transformations. Our diabaticity criterion is implemented in a way that maximizes the diversity of VB structure weights between different diabatic states. Owing to the rigorous Householder transformations employed in this entire procedure, the invariance of the target eigensubspace is preserved. This is illustrated on two prototypical examples.

In this work, a valence bond type multireference density functional theory (MRDFT) method, called the Hamiltonian matrix correction based density functional valence bond method (hc-DFVB), is presented. In hc-DFVB, the static electronic correlation is considered by the valence bond self-consistent field (VBSCF) strategy, while the dynamic correlation energy is taken into account by Kohn-Sham density functional theory (KS-DFT). Different from our previous version of DFVB (J. Chem. Theory Comput. 2012, 8, 1608), hc-DFVB corrects the dynamic correlation energy with a Hamiltonian correction matrix, improving the functional adaptability and computational accuracy. The method was tested for various physical and chemical properties, including spectroscopic constants, bond dissociation energies, reaction barriers, and singlet-triplet gaps. The accuracy of hc-DFVB matches that of KS-DFT and high level molecular orbital (MO) methods quite well. Furthermore, hc-DFVB keeps the advantages of VB methods, which are able to provide clear interpretations and chemical insights with compact wave functions.

A new energy decomposition analysis (EDA) scheme based on valence bond (VB) wave function, called VB-EDA, is presented. In VB-EDA, the total interaction energy is decomposed into frozen, charge transfer, polarization and dynamic correlation terms based on valence bond calculations. The frozen term is the energy variation of the unrelaxed VB wave function according to the change of an interaction distance. The charge transfer term is the contribution of the additional VB structures while the polarization term is due to the relaxation of VB orbitals. Dynamic correlation term is computed by post-VBSCF methods. Different from other existing VB based EDA schemes, which were used to analyze noncovalent interactions for some specific complexes, the newly developed VB-EDA is designed for the general use. Using VB-EDA, the bonding nature of cation-π interactions in a series of cation-π complexes (cations = Li, Na, K, Mg, and Ca; π systems = ethylene and benzene) is explored. Furthermore, a new covalency index, which demonstrates the covalency of cation-π interactions, is presented based on the VB-EDA results. The VB-EDA analysis reveals that the cation-π interactions in the Li, Na, and K complexes belong to the typical ionic bonds while the Mg and Ca complexes have the relatively large covalent characteristics. However, only the CH-Mg complex can be regarded as a covalent bonding complex while the other complexes belong to the typical ionic complexes. Thereupon, it must be careful in the cognition for the covalency of intermolecular interaction. Large nonelectrostatic interaction component does not always correspond to a covalent bond.