Periodic local vibrational modes were calculated with the rev-vdW-DF2 density functional to quantify the intrinsic strength of the X-I⋯OA-type halogen bonding (X = I or Cl; OA: carbonyl, ether and N-oxide groups) in 32 model systems originating from 20 molecular crystals. We found that the halogen bonding between the donor dihalogen X-I and the wide collection of acceptor molecules OA features considerable variations of the local stretching force constants (0.1–0.8 mdyn/Å) for I⋯O halogen bonds, demonstrating its powerful tunability in bond strength. Strong correlations between bond length and local stretching force constant were observed in crystals for both the donor X-I bonds and I⋯O halogen bonds, extending for the first time the generalized Badger’s rule to crystals. It is demonstrated that the halogen atom X controlling the electrostatic attraction between the σ -hole on atom I and the acceptor atom O dominates the intrinsic strength of I⋯O halogen bonds. Different oxygen-containing acceptor molecules OA and even subtle changes induced by substituents can tweak the n → σ ∗ (X-I) charge transfer character, which is the second important factor determining the I⋯O bond strength. In addition, the presence of the second halogen bond with atom X of the donor X-I bond in crystals can substantially weaken the target I⋯O halogen bond. In summary, this study performing the in situ measurement of halogen bonding strength in crystalline structures demonstrates the vast potential of the periodic local vibrational mode theory for characterizing and understanding non-covalent interactions in materials.

The first combined quantum chemical energy-electron density description of bond breaking/forming events using URVA and QTAIM.

The analysis of chemical bonding in crystal structures and surfaces is an important research topic in theoretical chemistry. In this work, we present a PyMOL plugin, named LModeA-nano, as implementation of the local vibrational mode theory for periodic systems (Tao et al. J. Chem. Theory Comput. 2019, 15, 1761) assessing bond strength in terms of local stretching force constants in extended systems of one, two, and three dimensions. LModeA-nano can also analyze chemical bonds in isolated molecular systems thus enabling a head-to-head comparison of bond strength across systems with different dimensions in periodicity (0−3D). The new code is interfaced to the output generated by various solid-state modeling packages including VASP, CP2K, Quantum ESPRESSO, CASTEP, and CRYSTAL. LModeA-nano is crossplatform, open-source and freely available on GitHub: https://github.com/smutao/ LModeA-nano.

Local stretching force constants derived from periodic local vibrational modes at the vdW-DF2 density functional level have been employed to quantify the intrinsic hydrogen bond strength of 16 ice polymorphs, ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI, XII, XIII, XIV, XV, XVII, and XIX, that are stable under ambient to elevated pressures. Based on this characterization on 1820 hydrogen bonds, relationships between local stretching force constants and structural parameters such as hydrogen bond length and angle were identified. Moreover, different bond strength distributions, from uniform to inhomogeneous, were observed for the 16 ices and could be explained in relation to different local structural elements within ices, that is, rings, that consist of different hydrogen bond types. In addition, criteria for the classification of hydrogen bonds as strong, intermediate, and weak were introduced. The latter was used to explore a different dimension of the water–ice phase diagram. These findings will provide important guidelines for assessing the credibility of new ice structures.

The unified reaction valley approach combined with the local vibrational mode and ring puckering analysis is applied to investigate the hydrogen evolution from water in the presence of small hydrides such as BH3, metal hydrides as AlH3, and their derivatives. We studied a series of reactions involving BH3, AlH3, B2H6, Al2H6, and AlH3BH3 with one- and two-water molecules, considering multiple reaction paths. In addition, the influence of the aqueous medium was examined. A general reaction mechanism was identified for most of the reactions. Those that deviate could be associated with unusually high reaction barriers with no hydrogen release. The charge transfer along the reaction path suggests that a viable hydrogen release is achieved when the catalyst adopts the role of a charge donor during the chemical processes. The puckering analysis showed that twistboat and boat forms are the predominant configurations in the case of an intermediate six-membered ring formation, which influences the activation barrier. The local mode analysis was used as a tool to detect the H–H bond formation as well as to probe catalyst regenerability. Based on the correlation between the activation energy and the change in the charge separation for cleaving O–H and B­(Al)–H bonds, two promising subsets of reactions could be identified along with prescriptions for lowering the reaction barrier individually with electron-donating/withdrawing substituents.