==== Front ChemistryOpenChemistryOpen10.1002/(ISSN)2191-1363OPENChemistryOpen2191-1363John Wiley and Sons Inc. Hoboken 3009412410.1002/open.201800051OPEN201800051Full PaperFull PapersIntrinsic Dynamic Nature of Neutral Hydrogen Bonds Elucidated with QTAIM Dual Functional Analysis: Role of the Compliance Force Constants and QTAIM‐DFA Parameters in Stability Nishide Taro 1 Hayashi Satoko Dr.hayashi3@sys.wakayama-u.ac.jp 1 Nakanishi Waro Prof.http://orcid.org/0000-0003-2163-7134nakanisi@sys.wakayama-u.ac.jp 1 1 Wakayama University Faculty of Systems Engineering Wakayama Japan 06 6 2018 8 2018 7 8 10.1002/open.v7.8565 575 05 4 2018 © 2018 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.This is an open access article under the terms of the http://creativecommons.org/licenses/by-nc-nd/4.0/ License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.Abstract The dynamic and static nature of various neutral hydrogen bonds (nHBs) is elucidated with quantum theory of atoms‐in‐molecules dual functional analysis (QTAIM‐DFA). The perturbed structures generated by using the coordinates derived from the compliance force constants (Cij) of internal vibrations are employed for QTAIM‐DFA. The method is called CIV. The dynamic nature of CIV is described as the “intrinsic dynamic nature”, as the coordinates are invariant to the choice of the coordinate system. nHBs are, for example, predicted to be van der Waals (H2Se−✶−HSeH; ✶=bond critical point), t‐HBnc (typical‐HBs with no covalency: HI−✶−HI), t‐HBwc (t‐HBs with covalency: H2C=O−✶−HI), CT‐MC [molecular complex formation through charge transfer (CT): H2C=O−✶−HF], and CT‐TBP (trigonal bipyramidal adduct formation through CT: H3N−✶−HI) in nature. The results with CIV were the same as those with POM in the calculation errors, for which the perturbed structures were generated by partial optimization, and the interaction distances in question were fixed suitably in POM. The highly excellent applicability of CIV for QTAIM‐DFA was demonstrated for the various nHBs, as well as for the standard interactions previously reported. The stability of the HBs, evaluated by ΔE, is well correlated with Cij (ΔE×Cij=constant value of −165.64), and the QTAIM parameters, although a few deviations were detected. ab initio calculationsatoms-in-molecules dual functional analysis (QTAIM-DFA)compliance force constantsdynamicshydrogen bonds source-schema-version-number2.0component-idopen201800051cover-dateAugust 2018details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_NLMPMC version:version=5.4.9 mode:remove_FC converted:19.09.2018 T. Nishide, S. Hayashi, W. Nakanishi, ChemistryOpen 2018, 7, 565. ==== Body 1 Introduction Hydrogen bonds (HBs) are fundamentally important because of their ability to form molecular associations, which stabilizes a system in terms of energy; the direction of the interacting three atoms in B⋅⋅⋅H−Y (see below) is controlled through the formation of a HB that is almost a linear asymmetric σ bond (3c–4e: three‐center, four‐electron bond).1, 2, 3, 4, 5, 6 Weak HBs can be considered to be van der Waals (vdW) interactions, whereas strong HBs tend to be more covalent (Cov) in nature. The formation of HBs plays a crucial role in all fields of chemical and biological sciences. HBs control various chemical processes depending on their strength. It is imperative to clarify the nature of HBs for better understanding of chemical processes controlled by HBs.7, 8, 9, 10 We previously reported the dynamic and static nature of HBs in the neutral and charged forms by applying the quantum theory of atoms‐in‐molecules dual functional analysis (QTAIM‐DFA).11, 12, 13, 14, 15, 16 Perturbed structures were employed for QTAIM‐DFA to clarify the dynamic behavior of the interactions. The perturbed structures were generated by using the normal coordinates of the (best‐fitted) internal vibrations and/or by partial optimization, for which the interaction distances in question were suitably fixed. The methods are called NIV14, 15, 16 and POM,11, 12, 13 respectively. Neutral HBs (nHBs) are predicted to be vdW to CT‐TBP [trigonal bipyramidal adduct formation through charge transfer (CT)] in nature, whereas charged HBs are typically Cov in nature.17 The QTAIM approach, introduced by Bader,18, 19 enables us to analyze the nature of chemical bonds and interactions.20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 Many QTAIM investigations have been reported, mainly by theoretical researchers, and experimental scientists have recently used QTAIM to explain their results by considering chemical bonding and interactions.20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39 A bond critical point (BCP, ✶)40 is an important concept in QTAIM, for which ρ(r) (charge density) reaches a minimum along the interatomic (bond) path, whereas it is a maximum on the interatomic surface separating the atomic basins. The ρ(r) at the BCP is described by ρ b(r c) and so are other QTAIM functions, such as the total electron energy densities H b(r c), potential energy densities V b(r c), and kinetic energy densities G b(r c) at the BCPs. A chemical bond or interaction between A and B is denoted by A−B, which corresponds to a bond path (BP) in QTAIM. We will use A−✶−B for a BP, in which the asterisk emphasizes the existence of a BCP in A−B.18, 19, 40 Equations (1), (2), (2′) represent the relations between G b(r c), V b(r c), H b(r c), and ∇2 ρ b(r c); H b(r c) must be negative if ∇2 ρ b(r c) <0, as confirmed by Equation (2) with negative V b(r c) at all BCPs. (1) Hb(rc)=Gb(rc)+Vb(rc) (2) (ℏ28m)∇2ρb(rc)=Hb(rc)-Vb(rc)2 (2′) =Gb(rc)+Vb(rc)2 QTAIM‐DFA was recently formulated on the basis of QTAIM.41, 42, 43, 44, 45, 46, 47, 48 The H b(r c) parameters are plotted versus H b(r c)−V b(r c)/2 [=(ħ2/8 m)∇2 ρ b(r c)] at the BCPs in QTAIM‐DFA. In our treatment, data from the perturbed structures around the fully optimized structures are employed, in addition to those from the fully optimized ones. Data from the fully optimized structures are analyzed by using the polar coordinate (R, θ) representation, which corresponds to the static natures of the interactions.12, 13, 14, 15, 16 Each interaction plot, which contains data from both the perturbed structures and the fully optimized ones, includes a specific curve that provides important information about the interaction. This plot is expressed by (θ p, κ p), for which θ p corresponds to the tangent line of the plot and κ p is the curvature. The dynamic nature of the interactions has been proposed on the basis of (θ p, κ p).12, 13, 14, 15, 16 Namely, the signs of the first and second derivatives of H b(r c) and H b(r c)−V b(r c)/2 [=(ħ2/8 m)∇2 ρ b(r c)] are also employed to characterize the interactions in QTAIM‐DFA, whereas the signs of H b(r c) and H b(r c)−V b(r c)/2 are employed for classification in QTAIM. We call (R, θ) and (θ p, κ p) the QTAIM‐DFA parameters, which are drawn in Figure 2, exemplified by H3N−✶−HI (26: C 3v). QTAIM‐DFA is applied to standard interactions, and rough criteria that distinguish the interaction in question from others are obtained. QTAIM‐DFA and the criteria are explained in the Supporting Information by using Schemes S1 and S2, Figures S1, and Equations (S1)–(S7). The basic concept of the QTAIM approach is also discussed. QTAIM‐DFA has excellent potential for evaluating, classifying, characterizing, and understanding weak to strong interactions according to a unified form.12, 13, 14, 15, 16 However, the predicted dynamic nature would somewhat change depending on the perturbed structures employed. Very recently, we proposed a new method to generate the perturbed structures, other than NIV and POM, for QTAIM‐DFA.49 The method employs the coordinates corresponding to the compliance force constants (Cij) for the internal vibrations.50, 51, 52, 53, 54 We call the method CIV.49 The compliance force constants (Cij)50, 51, 52, 53, 55 are defined as the partial second derivatives of the potential energy due to an external force, as shown in Equation (3), for which i and j refer to internal coordinates, and the force constants fi and fj correspond to i and j, respectively. The value in Equation (3) corresponds to a lower numerical value (i) of a compliance constant representing a stronger bond (j), that is, the compliance constants measure the flexibility (or compliance) of a particular bond. The applications of CIV to the closed‐shell (CS) interactions are substantially more effective than those to the shared‐shell (SS) interactions in QTAIM‐DFA.49 (3) Cij=∂2E∂fi∂fj The very high applicability of CIV is demonstrated to generate the perturbed structures for QTAIM‐DFA. The dynamic nature of the interactions based on the perturbed structures with CIV is described as the “intrinsic dynamic nature of interactions”, as the coordinates corresponding to Cij are invariant to the choice of the coordinate system. The results with CIV are the same as those with POM in terms of the calculation errors. However, CIV has been applied only to the typical interactions of a limited number of HBs, and the default in NIV seems large for HBs.49 The establishment of QTAIM‐DFA on the firm basis of employing the perturbed structures with CIV for the wide range of nHBs is another purpose of this paper. The neutral HBs in the species examined in this work are denoted by B−✶−HY (1–29), containing HI adducts, which are found in Table 1. Table 1 QTAIM functions and QTAIM‐DFA parameters evaluated for the neutral hydrogen bonds (nHBs) in 1–29 by applying the QTAIM dual functional analysis by employing the perturbed structures generated with CIV, NIV, and POM.[a,b] Species (X−✶−Y) c∇2 ρ b(r c)[c] [au] H b(r c) [au] k b(r c)[d] R [au] θ [°] Cij [e] [mDyn−1 Å3] θ p:CIV [°] κ p:CIV [au−1] H2Se−✶−HSeH (1) 0.0026 0.0006 −0.858 0.0027 76.0 [f] 23.4 88.1 194 H2S−✶−HSH (2) 0.0032 0.0008 −0.861 0.0033 76.3 19.2 91.8 229 H3N−✶−HNH2 (3) 0.0059 0.0016 −0.844 0.0062 74.9 12.1 87.5 188 H2O−✶−HOH (4) 0.0106 0.0005 −0.976 0.0107 87.3 6.4 123.7 159 H3N−✶−HOH (5) 0.0094 −0.0020 −1.096 0.0096 102.1 [f] 5.5 157.3 87.3 HI−✶−HI (6) 0.0034 0.0004 −0.945 0.0034 84.1 13.8 102.6 304 HBr−✶−HBr (7) 0.0038 0.0010 −0.853 0.0039 75.6 17.9 91.4 269 HCl−✶−HCl (8) 0.0049 0.0015 −0.828 0.0052 73.6 16.1 95.0 294 HF−✶−HF (9) 0.0125 −0.0002 −1.007 0.0125 90.8 5.9 128.2 107 H2Se−✶−HI (10) 0.0040 0.0001 −0.986 0.0040 88.5 12.7 126.5 [f] 464 H2Se−✶−HBr (11) 0.0040 0.0002 −0.978 0.0040 87.6 13.0 130.1 [f] 488 H2Se−✶−HCl (12) 0.0044 0.0001 −0.989 0.0044 88.7 11.2 137.3 [f] 431 H2Se−✶−HF (13) 0.0051 −0.0013 −1.113 0.0053 104.3 [f] 7.3 164.5 146 H2S−✶−HI (14) 0.0043 0.0001 −0.991 0.0043 89.0 13.4 124.5 334 H2S−✶−HBr (15) 0.0047 −0.0001 −1.010 0.0047 91.1 12.0 133.9 309 H2S−✶−HCl (16) 0.0051 −0.0002 −1.024 0.0051 92.8 10.3 140.5 269 H2S−✶−HF (17) 0.0061 −0.0020 −1.143 0.0064 108.5 [f] 6.6 165.1 120 H2O−✶−HI (18) 0.0091 0.0009 −0.949 0.0091 84.5 10.1 113.6 217 H2O−✶−HBr (19) 0.0103 −0.0006 −1.028 0.0103 93.2 8.2 138.6 182 H2O−✶−HCl (20) 0.0112 −0.0018 −1.072 0.0114 98.9 6.4 149.9 116 H2O−✶−HF (21) 0.0131 −0.0089 −1.252 0.0158 124.0 3.4 166.1 6.9 H2C=O−✶−HI (22) 0.0102 −0.0009 −1.044 0.0103 95.3 9.7 139.7 216 H2C=O−✶−HBr (23) 0.0108 −0.0022 −1.093 0.0111 101.6 [f] 8.2 154.4 138 H2C=O−✶−HCl (24) 0.0115 −0.0032 −1.122 0.0119 105.9 [f] 6.6 160.4 92.0 H2C=O−✶−HF (25) 0.0127 −0.0099 −1.279 0.0161 127.8 3.5 170.1 6.7 H3N−✶−HI (26) 0.0050 −0.0268 −1.728 0.0272 169.4 19.8 194.1 4.2 H3N−✶−HBr (27) 0.0069 −0.0189 −1.579 0.0201 160.0 7.9 190.3 6.4 H3N−✶−HCl (28) 0.0080 −0.0155 −1.492 0.0174 152.7 5.5 186.9 9.3 H3N−✶−HF (29) 0.0085 −0.0195 −1.533 0.0213 156.4 2.8 182.0 2.8 Species (X−✶−Y) θ p:POM [°] κ p:POM [au] ν [g] [cm−1] k f [h] [mDyn Å−1] θ p:NIV [°] κ p:NIV [au] ΔE [i] [kJ mol−1] Predicted nature H2Se−✶−HSeH (1) 88.0 202 41.8 0.016 88.3 196 −7.6 p‐CS/vdW H2S−✶−HSH (2) 91.8 232 69.1 0.009 93.3 263 −8.7 p‐CS/t‐HBnc H3N−✶−HNH2 (3) 87.5 151 141.2 0.036 86.6 159 −13.8 p‐CS/vdW H2O−✶−HOH (4) 123.8 159 188.1 0.043 116.7 158 −22.2 p‐CS/t‐HBnc H3N−✶−HOH (5) 157.5 88.1 200.2 0.050 158.6 83.0 −28.2 r‐CS/CT‐MC HI−✶−HI (6) 102.7 309 43.5 0.024 102.5 296 −12.9 p‐CS/t‐HBnc HBr−✶−HBr (7) 91.4 269 48.8 0.028 91.2 259 −8.3 p‐CS/t‐HBnc HCl−✶−HCl (8) 95.0 295 76.4 0.021 94.8 267 −10.0 p‐CS/t‐HBnc HF−✶−HF (9) 128.3 109 166.9 0.081 128.5 103 −20.7 r‐CS/t‐HBwc H2Se−✶−HI (10) 126.5 [f] 464 52.5 0.031 126.4 [f] 454 −14.5 p‐CS/t‐HBnc H2Se−✶−HBr (11) 130.0 498 57.9 0.044 129.9 [f] 480 −13.9 p‐CS/t‐HBnc H2Se−✶−HCl (12) 137.4 438 79.3 0.057 137.1 423 −15.5 p‐CS/t‐HBnc H2Se−✶−HF (13) 164.5 151 123.0 0.101 163.9 144 −21.3 r‐CS/CT‐MC H2S−✶−HI (14) 124.3 340 68.2 0.017 125.3 [f] 325 −13.9 p‐CS/t‐HBnc H2S−✶−HBr (15) 133.9 317 77.5 0.028 134.3 301 −14.6 r‐CS/t‐HBwc H2S−✶−HCl (16) 140.6 274 98.2 0.042 140.7 260 −16.6 r‐CS/t‐HBwc H2S−✶−HF (17) 165.1 121 145.7 0.096 164.6 117 −23.2 r‐CS/CT‐MC H2O−✶−HI (18) 112.9 212 97.6 0.013 122.9 227 −18.1 p‐CS/t‐HBnc H2O−✶−HBr (19) 138.1 186 119.6 0.034 140.7 168 −20.7 r‐CS/t‐HBwc H2O−✶−HCl (20) 149.9 120 150.1 0.048 152.0 [j] 104 −24.7 r‐CS/t‐HBwc [k] H2O−✶−HF (21) 166.1 8.5 229.9 0.079 167.6 7.1 −38.4 r‐CS/CT‐MC H2C=O−✶−HI (22) 139.8 202 141.9 0.049 138.3 219 −21.5 r‐CS/t‐HBwc H2C=O−✶−HBr (23) 154.5 135 152.0 0.070 152.8 140 −22.4 r‐CS/CT‐MC H2C=O−✶−HCl (24) 160.4 91.7 176.0 0.115 158.6 91.8 −25.9 r‐CS/CT‐MC H2C=O−✶−HF (25) 170.0 8.0 246.7 0.267 168.3 5.7 −36.3 r‐CS/CT‐MC H3N−✶−HI (26) 194.2 5.3 100.7 0.025 193.9 4.2 −30.5 r‐CS/CT‐TBP H3N−✶−HBr (27) 190.3 8.1 148.1 0.059 189.8 6.3 −33.7 r‐CS/CT‐TBP H3N−✶−HCl (28) 186.9 11.7 186.8 0.105 186.2 9.3 −38.0 r‐CS/CT‐TBP H3N−✶−HF (29) 181.9 5.4 227.0 0.241 180.6 1.8 −54.8 r‐CS/CT‐TBP [a] The functions and parameters were evaluated at the BCPs of the nHBs in the fully optimized structures. [b] With MP2/6–311++G(3df,3pd), except for I, for which calculations were performed with (7433111/743111/7411/2+1s1p1d1f) from the Sapporo Basis Set Factory, which is called MP2/BSS‐A. [c] c∇2 ρ b(r c)=H b(r c)−V b(r c)/2, for which c=ħ2/8 m. [d] k b(r c)=V b(r c)/G b(r c). [e] Defined in Equation (3). [f] Minor values that do not satisfy the characterization from the major ones are shown in italics. [g] Internal vibrational frequency corresponding to the interaction. [h] Force constant corresponding to the frequency. [i] From the components. [j] The nature of r‐CS/CT‐MC was predicted with NIV. [k] On the borderline area between r‐CS/t‐HBwc and r‐CS/CT‐MC if evaluated with CIV and POM. Wiley‐VCH Verlag GmbH & Co. KGaAHerein, we present the results of investigations on the “intrinsic dynamic nature of nHBs”, together with the static nature in B−✶−HY (1–29). To elucidate the nature, QTAIM‐DFA is applied to B−✶−HY (1–29) by employing the perturbed structures generated with CIV. Each HB interaction in 1–29 is classified and characterized by employing the criteria obtained on the basis of the standard interactions as a reference. The applicability of CIV to QTAIM‐DFA is also established in the nHBs of 1–29, for which the QTAIM‐DFA parameters elucidated by using CIV are compared with those elucidated by using NIV and POM. As a result, a firm basis for QTAIM‐DFA by employing the perturbed structures generated with CIV is established over the wide range of the nHBs in 1–29. The stability of 1–29 is discussed by examining the relations between the stability and the Cij, (R, θ), and (θ p, κ p) values. Computational Details The 6‐311++G(3df,3pd) basis sets of the Gaussian 09 programs56 were employed for the calculations of 1–29, together with the basis set of the 7433111/743111/7411/2+1s1p1d1f type for I, as implemented in the Sapporo Basis Set Factory.57 The basis set system is called BSS‐A. All calculations were performed under nonrelativistic conditions. The Møller–Plesset second‐order energy correlation (MP2) level58, 59, 60 was applied to the calculations (MP2/BSS‐A). The optimized structures were confirmed by frequency analysis. To obtain the perturbed structures with POM, the optimized structures were further (partially) optimized by employing Z matrices. The distances in question in the perturbed structures (r) were fixed to satisfy Equation (4), in which r o is the distance in the fully optimized structure with a o of Bohr radius (0.52918 Å). The method to generate the perturbed structures with NIV is explained by Equation (5). A k‐th perturbed structure in question (S kw) is generated by the addition of the normal coordinates of the k‐th internal vibration (N k) to the standard orientation of a fully optimized structure (S o) in the matrix representation. The coefficient fkw in Equation (5) controls the structural difference between S kw and S o: fkw is determined to satisfy Equation (4) for r. The selected motion must be most effectively localized on the interaction in question among the zero‐point internal vibrations in NIV. Equations (4) and (6) are similarly applied to generate the perturbed structures with CIV, for which Ci is the coordinates corresponding to Cij in Equation (3).50, 51, 52, 53, 55 The N k and Ci values of five digits are used to predict S kw and S iw, respectively: (4) r=ro+wao[w=(0),±0.05,and±0.1;ao=0.52918Å] (5) Skw=So+fkw·Nk (6) Siw=So+fiw·Ci (7) y=co+c1x+c2x2+c3x3(Rc2:squareofcorrelationcoefficient) QTAIM functions are calculated by using the same basis sets at the MP2 level as in the optimizations (MP2/BSS‐A), unless otherwise noted, and are analyzed with the AIM2000 program.18, 61 H b(r c) are plotted versus H b(r c)−V b(r c)/2 for the data of five data points of w=0, ±0.05, and ±0.1 in Equation (4) in QTAIM‐DFA. Each plot is analyzed by using a regression curve of the cubic function, shown in Equation (7), for which (x, y)=[H b(r c)−V b(r c)/2, H b(r c)] (R c 2>0.99999 per usual).16 2 Results and Discussion 2.1 Optimized Structures of Neutral Hydrogen‐Bonded Species and Stability Neutral HB species, B−HY [1–29: B=H2Se, H2S, H2O, H2CO, H3N, and HX (=HI, HBr, HCl, and HF); HY=HSeH, HSH, HOH, HNH2, and HX], were optimized with MP2/BSS‐A, although some were optimized in a previous study.17 The optimized B−H distances [r o(B, H)] are collected in Table S1, together with the sum of the vdW radii62 [Δr=r o(B, H)−Σr vdW(B, H)]. The symmetries are also given in Table S1. The energies for 1–29 on the energy surface (E) and the relative energies from the components (ΔE) [=E(HB)−E(components)] are collected in Table S1. The ΔE values of 1–29 are also shown in Table 1 for convenience of discussion. The ΔE values are plotted versus Δr for 1–29. The plot is shown in Figure S2, and the correlation is given in the figure.63 The ΔE values of B−✶−HX (B=H2Se, H2S, H2C=O, and H3N; HX=HF, HCl, HBr, and HI) are plotted versus those of H2O−✶−HX. The plot is shown in Figure S3, which also contains the plot of ΔE(H2O−✶−HX) versus ΔE(H2O−✶−HX) for convenience of comparison. The correlations are very good (Table 2, entries 1–4). The results show that the ΔE values of B−✶−HX (B=H2Se, H2S, H2C=O, H2O, and H3N) are well correlated with each other if the ΔE values of common HX are compared, although the ΔE value of H2Se−✶−HI seems somewhat smaller (more stable) than that predicted from the correlation for H2Se−✶−HX. The magnitudes of ΔE become larger in the order H2Se≤H2S≪H2O≤H2C=O≪H3N, although ΔE(H2O−✶−HF)<ΔE(H2C=O−✶−HF). The relations between ΔE in B−✶−HX are also confirmed in this work for HX=HI in addition to HX=HF, HCl, and HBr, although the E values are all evaluated under nonrelativistic conditions. Table 2 Correlations in 1–29, evaluated with NIV, POM, and CIV, under the MP2/BSS‐A conditions.[a] Entry Correlation a b R c 2 n 1 ΔE H2Se-HX vs. ΔE H2O-HX 0.371 −6.83 0.960 4 2 ΔE H2S-HX vs. ΔE H2O-HX 0.468 −5.13 0.997 4 3 ΔE H2CO-HX vs. ΔE H2O-HX 0.751 −7.39 0.996 4 4 ΔE H3N-HX vs. ΔE H2O-HX 1.200 −8.72 0.9997 4 5 θ p:NIV vs. θ p:CIV 0.988 1.71 0.994 29 6 θ p:NIV vs. θ p:CIV 0.992 1.02 0.999 27[b] 7 θ p:POM vs. θ p:CIV 1.001 −0.15 0.99997 29 8 κ p:NIV vs. κ p:CIV 0.980 −0.31 0.994 29 9 κ p:POM vs. κ p:CIV 1.009 −0.42 0.998 29 10 ΔE vs. R −2012.0 −3.83 0.866 27[c] 11 ΔE vs. θ −0.479 25.70 0.891 27[c] 12 θ p:CIV vs. θ 2.390 −86.95 0.957 23[d] 13 ΔE vs. θ p:CIV −0.314 19.69 0.971 8[e] 14 ΔE vs. θ p:CIV −0.219 14.01 0.957 8[f] 15 ΔE vs. θ p:CIV −0.155 −0.05 0.838 6[g] 16 ΔE vs. θ p:CIV 1.994 −414.66 0.898 4[h] [a] Analyzed by assuming the linear correlation y=ax+b (R c 2: square of correlation coefficient). [b] For 1–29, except for 4 and 18. [c] For 1–29, except for 26 and 27. [d] For 1–29, except for 21, 24, and 25–29. [e] For 1–9, except for 3. [f] For 10–17. [g] For 18–25, except for 21 and 25. [h] For 24–29. Wiley‐VCH Verlag GmbH & Co. KGaAAfter clarifying the basic behavior in Δr and ΔE, molecular graphs with contour plots of ρ(r) are examined before detailed discussion of the nature of the nHBs in 1–29. 2.2 Molecular Graphs with Contour Plots for B−✶−HX Figure 1 illustrates molecular graphs for B−✶−HI (B=H2Se, H2S, H2O, H2C=O, H3N, and HI) containing the contour plots of ρ(r). All of the BCPs expected for B−✶−HI are clearly detected. They seem to be well located at three‐dimensional saddle points of ρ(r). The molecular graphs of 1–29, other than B−✶−HI, were similarly drawn, and although they are not shown, they are very close to those of B−✶−HI. Figure 1 Molecular graphs, with contour plots of ρ(r) for a) HI−✶−HI (6), b) H2Se−✶−HI (10), c) H2S−✶−HI (14), d) H2O−✶−HI (18), e) H2C=O−✶−HI (22), and f) H3N−✶−HI (26). 2.3 Survey of the B−✶−HY Interactions The HB interactions seem straight for B−✶−HX on the basis of the BPs, as shown in Figure 1. To examine the linearity of the BPs, further, the lengths of the BPs (r BP) in question are collected in Table S2 for 1–29, together with the corresponding straight‐line distances (R SL). The differences between them (Δr BP=r BP−R SL) are less than 0.04 Å. Consequently, the BPs for all B−✶−HY of 1–29 can be described by straight lines. QTAIM functions were calculated for B−✶−HY (1–29) at the BCPs. Table 1 collects the H b(r c)−V b(r c)/2 [=(ħ2/8 m)∇2 ρ b(r c)] and H b(r c) values, and the ρ b(r c) values are given in Table S3, whereas some were reported previously.17 The H b(r c) values are plotted versus H b(r c)−V b(r c)/2 for the data shown in Table 1 together with those from the perturbed structures generated with CIV. Figure 2 shows the plots. The plots appear in the region of H b(r c)−V b(r c)/2>0; therefore, the HBs are all classified by CS interactions. The CS interactions will be further classified by the signs of H b(r c). They are specifically called pure CS (p‐CS) interactions if they appear in the region of H b(r c)>0, whereas they will be regular CS (r‐CS) interactions for H b(r c)<0. The behavior of the nHBs in 1–29 will be discussed in detail after evaluations of the QTAIM‐DFA parameters (see Table 1). Figure 2 Plots of H b(r c) versus H b(r c)−V b(r c)/2 for 1–29, for which data from the perturbed structures generated with CIV were employed, in addition to the data from the optimized structures. Definitions of (R, θ) and (θ p, κ p) are illustrated, as exemplified by H3N−✶−HI (26: C 3v). 2.4 QTAIM‐DFA Parameters of (R, θ) and (θ p, κ p) for Neutral HBs in 1–29 Evaluated with POM, NIV, and CIV The QTAIM‐DFA parameters of (R, θ) and (θ p, κ p) were obtained by analyzing the plots of H b(r c) versus H b(r c)−V b(r c)/2 according to Equations (S3)–(S6). The (θ p, κ p) values evaluated by employing the perturbed structures generated with CIV, POM, and NIV are denoted by (θ p:CIV, κ p:CIV), (θ p:POM, κ p:POM), and (θ p:NIV, κ p:NIV), respectively. The (θ p:CIV, κ p:CIV) values can be obtained if the plots shown in Figure 2 are analyzed. Table 1 collects the QTAIM‐DFA parameters for 1–29. Table 1 contains the Cij values for the nHBs in 1–29 together with the frequencies correlated to the NIVs employed to generate the perturbed structures and the force constants (k f). 2.5 Behavior of θ p:CIV, θ p:POM, and θ p:NIV Together with That of κ p:CIV, κ p:POM, and κ p:NIV Figure 3 a shows the plot of θ p:NIV versus θ p:CIV, which gives very good correlation. The correlation is shown in entry 5 of Table 2 (see also Figure 3 a). The magnitudes of the differences between θ p:NIV and θ p:CIV (Δθ p:NIV−CIV=θ p:NIV−θ p:CIV) are less than 2° for most of the interactions. The magnitudes of Δθ p:NIV−CIV are larger than 2.0° for H2O−✶−HOH (Δθ p:NIV−CIV=−7.0°),49 H2O−✶−HI (9.3°), H2O−✶−HBr (2.1°), and H2O−✶−HCl (2.1°). Large deviations are detected for H2O−✶−HX (HX=HOH and HI). The selected internal vibration for H2O−✶−HX could not be located effectively on O−✶−H by mixing with some other vibrational modes in the same symmetry,49 although the selected mode is the best fit for the O−✶−H interaction. The correlation for the plot is much improved (Table 2, entry 6; see also Figure 3 a) if the data for H2O−✶−HX (X=HOH and HI) are omitted from the correlation. On the other hand, excellent correlation is obtained if θ p:POM is plotted versus θ p:CIV, as shown in Figure 3 b (for the correlation, also see Table 2, entry 7). The magnitudes of Δθ p:POM−CIV are equal to or less than 0.1° for all HB adducts examined, except for H3N−✶−HOH (Δθ p:POM−CIV=0.2°), H2S−✶−HI (−0.2°), H2O−✶−HI (−0.7°), and H2O−✶−HBr (−0.5°). The results must be a reflection of the fact that the perturbed structures generated with POM and CIV are very similar.49 The results demonstrate the excellent applicability of CIV to generate the perturbed structures also for the nHB species in QTAIM‐DFA. Figure 3 Plots of a) θ p:NIV versus θ p:CIV and b) θ p:POM versus θ p:CIV. Figure 4 a, b shows the plots of κ p:NIV versus κ p:CIV and κ p:POM versus κ p:CIV. The correlations are given in entries 8 and 9 of Table 2 (see also Figure 4 a, b, respectively). The correlations seem very good, although substantial deviations are observed in the plots. The magnitudes of Δκ p:NIV−CIV are larger than 10 au−1 for H2S−✶−HSH (Δκ p:NIV−CIV=34 au−1), H3N−✶−HNH2 (−29 au−1), HBr−✶−HBr (−10 au−1), HCl−✶−HCl (−27 au−1), H2O−✶−HI (10 au−1), H2O−✶−HBr (−14 au−1), and H2O−✶−HCl (−12 au−1), together with magnitudes of 5 to 10 au−1 for HI−✶−HI (−8.7 au−1), H2Se−✶−HI (−9.3 au−1), H2Se−✶−HBr (−8.5 au−1), H2Se−✶−HCl (−7.6 au−1), H2S−✶−HI (−9.2 au−1), H2S−✶− HBr (−8.7 au−1), and H2S−✶−HCl (−9.1 au−1). In the case of Δκ p:POM−CIV, the magnitudes are less than 5 au−1 for most cases. The values are larger than 10 au−1 for H3N−✶−HNH2 (Δκ p:POM−CIV=−37 au−1), H2Se−✶−HBr (10 au−1), and H2C=O−✶−HI (−14 au−1), together with magnitudes of 5 to 10 au−1 for H2Se−✶−HSeH (8.0 au−1), H2Se−✶−HCl (7.1 au−1), H2S−✶−HI (6.2 au−1), and H2S−✶−HBr (8.0 au−1). The magnitudes for Δκ p:NIV−CIV seem very large at first glance. However, the very large values of κ p would be responsible for the large magnitudes of Δκ p as a whole. The magnitudes of Δκ p:POM−CIV seem to be much improved relative to the case of Δκ p:NIV−CIV; however, there are some severe deviations, such as H3N−✶−HNH2 (−37 au−1). Figure 4 Plots of a) κ p:NIV versus κ p:CIV and b) κ p:POM versus κ p:CIV. The correlation of θ p:POM versus θ p:CIV is much better than that of κ p:POM versus κ p:CIV (see Table 2, entries 7 and 9). This observation seems curious at first glance, as the same perturbed structures are employed to evaluate θ p and κ p in QTAIM‐DFA. The differences may be mainly attributable to the much more complex route to evaluate κ p (=[d2 y/dx 2]/[1+(dy/dx)2]3/2) relative to the case of θ p (=90°−tan−1 (dy/dx)). The θ p and κ p values are evaluated by using the common regression curve shown in Equation (7), as pointed out in a previous paper.49 The small differences in the QTAIM functions based on the perturbed structures generated with CIV and POM will be magnified in the second derivatives of the regression curves used to evaluate κ p. As discussed above, the θ p:POM values can be recognized to be the same as the θ p:CIV values in terms of the calculation errors as a whole, although the Δθ p:POM−CIV values of −0.7° for H2O−✶−HI and −0.5° for H2O−✶−HBr seem slightly larger than the calculation errors. There must be a reason for these results, but we did not examine this point further. Larger magnitudes of Δκ p:POM−CIV are usually detected if κ p is very large. However, the results will not damage the excellent reliability in the characterization of the nHBs, as the κ p values are not used to characterize the interactions. Namely, the excellent applicability of CIV to generate the perturbed structures for QTAIM‐DFA is also well established for the various nHBs, as discussed above. 2.6 Nature of Neutral HBs Evaluated with the (θ, θ p) Values The nature of the wide range of HBs for neutral forms 1–29 is now classified and characterized on the basis of the θ and θ p:CIV values obtained in this work. It is instructive to survey the criteria shown in Scheme S2 before a detailed discussion. The θ values classify interactions, whereas the θ p values predict the character of these interactions. The criteria tell us that 45°<θ<180° [00] for the pure CS interactions (p‐CS) and 90°<θ<180° [H b(r c)<0] for the regular CS interactions (r‐CS). In the p‐CS region of 45°<θ<90°, the character of interactions will be of the vdW type for 45°<θ p<90° (45°<θ<75°), whereas it will be of the typical hydrogen bond (t‐HB) type with no covalency (t‐HBnc) for 90°<θ p<125° (75°<θ<90°), for which θ=75° and θ p=125° are tentatively given for θ p=90° and θ=90°, respectively. The CT interaction will appear in the r‐CS region of 90°<θ<180°. The t‐HB interactions with covalency (t‐HBwc) appear in the range of 125°<θ p<150° (90°<θ<115°), for which (θ, θ p)=(115°, 150°) are tentatively given as the borderline between the t‐HBwc and CT‐MC (interactions in the molecular complex formation through CT) natures. The borderline in the interactions between CT‐MC and CT‐TBP is defined by (θ, θ p)=(150, 180°), for which θ=150° is tentatively given corresponding to θ p=180°. As a result, CT‐MC and CT‐TBP will appear in the ranges of 150°<θ p<180° (115°<θ<150°) and 180°<θ p<190° (150°<θ<180°), respectively. Namely, the (θ, θ p) values of (75°, 90°), (90°, 125°), (115°, 150°), (150°, 180°), and (180°,190°) correspond to the borderlines between the interactions for vdW/t‐HBnc, t‐HBnc/t‐HBwc, t‐HBwc/CT‐MC, CT‐MC/CT‐TBP, and CT‐TBP/Cov‐w (weak covalent bonds), respectively. The basic parameters of θ and θ p, described in bold, are superior to the tentatively given parameters in the classification and characterization of the interactions. The classical chemical bonds of the SS (180°<θ) will be strong if R>0.15 au (Cov‐s: strong covalent bonds), whereas they will be weak for R<0.15 au (Cov‐w), although SS interactions are not detected in the nHBs studied in this work. The (θ, θ p:CIV) values of H2Se−✶−HSeH (1) are (76.0°, 88.1°), and therefore, it is classified by the p‐CS interaction and is characterized by its vdW nature (p‐CS/vdW). The θ p:CIV value of 88.1° should be superior to θ=76.0° (>75°) to predict the nature. The HB interaction in H3N−✶−HNH2 (3) is also predicted to be p‐CS/vdW in nature with (θ, θ p:CIV)=(74.9°, 87.5°) for the interaction. However, the HB interactions in 1 and 3 would be close to the borderline area between p‐CS/vdW and p‐CS/t‐HBnc judging from the (θ, θ p:CIV) values. HB interactions other than these two were similarly classified and characterized. The nHB interactions are predicted to have the p‐CS/t‐HBnc nature for H2S−✶−HSH (2), H2O−✶−HOH (4), HX−✶−HX [6 (HX=HI), 7 (HBr), and 8 (HCl)], H2Se−✶− HX [10 (HX=HI), 11 (HBr), and 12 (HCl)], H2S−✶−HI (14), and H2O−✶−HI (18). The r‐CS/t‐HBwc nature is predicted for HF−✶−HF (9), H2S−✶−HX [15 (HX=HBr) and 16 (HCl)], H2O−✶−HX [19 (HX=HBr) and 20 (HCl)], and H2C=O−✶−HI (22). On the other hand, the r‐CS/CT‐MC nature is predicted for H3N−✶−HOH (5), H2Se−✶−HF (13), H2S−✶−HF (17), H2O−✶−HF (21), and H2C=O−✶−HX [23 (HX=Br), 24 (HCl), and 25 (HF)], whereas the r‐CS/CT‐TBP nature is predicted for H3N−✶−HX [26 (HX=HI), 27 (HBr), 28 (HCl), and 29 (HF)]. The results are summarized in Table 1. The superior values of θ or θ p can be employed to predict the nature if either θ or θ p:CIV does not satisfy the categories to determine the nature. The characterization based on POM is the same as that based on CIV, and the characterization based on NIV is equal to that based on CIV and POM, except for 20 (H2O−✶−HCl). The nature of r‐CS/CT‐MC is predicted for 20 with NIV, whereas it is just borderline between r‐CS/t‐HBwc and r‐CS/CT‐MC if evaluated with CIV and POM. The results show that the HB interactions can also be characterized satisfactorily by employing θ p:NIV for most cases, irrespective of the substantial differences between θ p:NIV and θ p:CIV in some cases. The predicted nature for B−✶−HX is summarized in Table 3, exemplified by the formation of B−✶−HX from B (=H2Se, H2S, H2O, H2C=O, and H3N) and HX (=HI, HBr, HCl, and HF). It enables us to visualize the roles of B and HX in the formation of B−✶−HX. The HB interactions are predicted to be stronger in the order shown in Equations (8), (9). The results shown in Table 3 and Equations (8) and (9) can be essentially explained on the basis of the results shown in Figure S3, although there are some differences, namely, the order shown in Equation (9) holds for B=H2Se, H2S, H2O, and H2C=O in BH−✶−HX, but it is reversed for B=H3N. The indirect B⋅⋅⋅(H)−X soft–soft interactions may affect the (θ, θ p:CIV) values in H2Se−✶−HI and H2S−✶−HI. (8) B=H2Se0; see Table 2, entry 16), contrary to the cases of G (B) and G (C) with negative correlations (a<0; see Table 2, entries 13–15). As shown in Figure 7, ΔE correlates rather well with θ p as a whole, with a few deviations. The correlation of ΔE versus θ p should be a reflection of the correlation of ΔE versus θ through the correlation of θ p versus θ. Figure 7 Plot of ΔE versus θ p for 1–29. Black dots for 1–9 belong to G (A), red triangles for 10–17 to G (B), blue squares for 18–25 to G (C), and green diamonds for 26–29 to G (D), although a few deviations are also included. The plot of ΔE versus κ p for 1–29 is shown in Figure S6. A linear correlation is not detected between κ p and ΔE, as expected. The plot shows a characteristic shape. The κ p values are less than 10 au−1 if ΔE<−30 kJ mol−1, 80 au−1<κ p<300 au−1 for −30 kJ mol−1<ΔE<−7 kJ mol−1, and especially 270 au−1<κ p<490 au−1 if −17 kJ mol−1<ΔE<−10 kJ mol−1, except for κ p=188 au−1 with ΔE=−13.8 kJ mol−1 for 6 (HI−✶−HI). Very large values of κ p are observed around the borderline area for HBs between the vdW and t‐HBnc types. The stability of the HBs in 1–29, evaluated by ΔE, is well explained on the basis of Cij, R, θ, and θ p. For the plots of ΔE versus Cij, the magnitudes of ΔE for 26 and 27 seem to be overestimated relative to those expected on the basis of the correlations for 1–25, 28, and 29. On the other hand, the magnitudes of ΔE for 26 and 27 would be underestimated relative to those expected from the correlations of ΔE versus R and θ for 1–25, 28, and 29. As shown in Figure 7, the plot for ΔE versus θ p could be recognized as a correlation as a whole, with deviation for G (D) (H3N−✶−HX) from the whole correlation for G (A)–G (C), if the correlation constants for the a values for the groups are compared. 3 Conclusions Hydrogen bonds (HBs) are fundamentally important in all fields of chemical and biological sciences. Therefore, HBs have been variously investigated. However, it has been difficult to characterize the nature of HBs spread over the range of van der Waals (vdW) type for pure closed‐shell (CS) interactions to the covalent type of shared‐shell (SS) interactions. In this paper, HBs of the neutral form were characterized by applying quantum theory of atoms‐in‐molecules dual functional analysis (QTAIM‐DFA) by employing the perturbed structures generated with CIV. The neutral hydrogen bond (nHB) interactions were characterized on the basis of the static and dynamic behavior predicted with QTAIM‐DFA. The static nature arises from the data of the fully optimized structures, whereas the dynamic nature originates from the data of the perturbed structures around the fully optimized structures. The dynamic nature of the interactions could be described as the “intrinsic dynamic nature of interactions” if the perturbed structures were generated with CIV, as the coordinates corresponding to the compliance force constants (Cij), used in CIV, are invariant to the choice of the coordinate system. The method was applied to nHBs and the interactions were characterized. Some of them are as follows: nHBs in H2Se−✶−HSeH and H3N−✶−HNH2 were characterized by the p‐CS (pure CS)/vdW nature. The p‐CS/t‐HBnc (typical hydrogen bond with no covalency) nature was predicted for H2S−✶−HSH, H2O−✶−HOH, and HX−✶−HX (HX=HI, HBr, and HCl), whereas the r‐CS (regular CS)/t‐HBwc (typical‐HB interactions with covalency) nature was predicted for HF−✶−HF, H2S−✶−HX (HX=HBr and HCl), and H2O−✶−HX (HX=HBr and HCl). On the other hand, HBs in H2C=O−✶−HX (HX=HBr, HCl, and HF) were predicted to have the r‐CS/CT‐MC (interactions in the molecular complex formation through CT) nature, whereas the r‐CS/CT‐TBP (trigonal bipyramidal adduct formation through CT) nature was predicted for H3N−✶−HX (HX=HI, HBr, HCl, and HF). Characterization based on POM was the same as that based on CIV, and characterization based on NIV was equal to that based on CIV and POM, except for 20 (H2O−✶−HCl). The r‐CS/CT‐MC nature was predicted for 20 with NIV, whereas it was borderline between r‐CS/t‐HBwc and r‐CS/CT‐MC if evaluated with CIV and POM. The highly excellent applicability of CIV is well demonstrated in QTAIM‐DFA by applying the method to nHBs, in addition to the standard interactions in a previous paper. Relations between ΔE and Cij or the QTAIM‐DFA parameters were examined, together with the reasons. A relation of ΔE×Cij=−165.64 was found for 1–25, 28, and 29. Namely, ΔE could be well described by the inverse nature of Cij. Similarly, the R, θ, and θ p values correlated linearly well with ΔE. The results showed that the values became larger as the stability of the HBs, described by ΔE, increased in the region examined, although there were a few deviations. Conflict of interest The authors declare no conflict of interest. Supporting information As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors. Supplementary Click here for additional data file. Acknowledgements This work was partially supported by the Ministry of Education, Culture, Sports, Science and Technology of Japan by a Grant‐in‐Aid for Scientific Research (No. 17K05785). ==== Refs 1 L. Pauling , The Nature of the Chemical Bond , Cornell University Press , Ithaca, NY , 1960 . 2 G. C. Pimentel , A. L. McClellan , The Hydrogen Bond , W. H. Freeman , San Francisco, CA , 1960 . 3 P. Schuster , G. Zundel , C. Sandorfy , The Hydrogen Bond, Recent Developments in Theory and Experiments , North-Holland Publishing Company , Amsterdam , 1976 . 4 G. A. Jeffrey , An Introduction to Hydrogen Bonding , Oxford University Press , New York , 1997 . 5 S. Scheiner , Hydrogen Bonding, A Theoretical Perspective , Oxford University Press , Oxford , 1997 . 6 G. R. Desiraju , T. Steiner , The Weak Hydrogen Bond in Structural Chemistry and Biology , International Union of Crystallography Monographs on Crystallography , Oxford University Press, New York , 1999 . 7 Hydrogen Bonding: New Insights (Challenges and Advances in Computational Chemistry and Physics) (Ed.: S. J. Grabowski ), Springer , The Netherlands , 2006 , Vol. 3 . 8 “Intramolecular Hydrogen Bonds. Methodologies and Strategies for Their Strength Evaluation”: G. Buemi in Hydrogen Bonding: New Insights (Challenges and Advances in Computational Chemistry and Physics) (Eds.: S. J. Grabowski ), Springer , New York , 2006 , Vol. 3 , Ch. 2. 9 K.-L. Han , G.-J. Zhao , Hydrogen Bonding and Transfer in the Excited State , Wiley , Chichester, UK , 2010 . 10 See the references cited in ref. [41]. 11 W. Nakanishi , S. Hayashi , K. Narahara , J. Phys. Chem. A 2008 , 112 , 13593 –13599 .19053566 12 W. Nakanishi , S. Hayashi , K. Narahara , J. Phys. Chem. A 2009 , 113 , 10050 –10057 .19621871 13 W. Nakanishi , S. Hayashi , Curr. Org. Chem. 2010 , 14 , 181 –197 . 14 W. Nakanishi , S. Hayashi , J. Phys. Chem. A 2010 , 114 , 7423 –7430 .20540587 15 W. Nakanishi , S. Hayashi , Bull. Chem. Soc. Jpn. 2012 , 85 , 1293 –1305 . 16 W. Nakanishi , S. Hayashi , J. Phys. Chem. A 2013 , 117 , 1795 –1803 .23347251 17 S. Hayashi , K. Matsuiwa , M. Kitamoto , W. Nakanishi , J. Phys. Chem. A 2013 , 117 , 1804 –1816 .23347280 18 Atoms in Molecules: A Quantum Theory (Ed.: R. F. W. Bader ), Oxford University Press , Oxford , 1990 . 19 “An Introduction to the Quantum Theory of Atoms in Molecules”: C. F. Matta , R. J. Boyd in The Quantum Theory of Atoms in Molecules: From Solid State to DNA and Drug Design (Eds.: C. F. Matta, R. J. Boyd ), Wiley-VCH , Weinheim, Germany , 2007 , Ch. 1. 20 R. F. W. Bader , T. S. Slee , D. Cremer , E. Kraka , J. Am. Chem. Soc. 1983 , 105 , 5061 –5068 . 21 R. F. W. Bader , Chem. Rev. 1991 , 91 , 893 –926 . 22 R. F. W. Bader , J. Phys. Chem. A 1998 , 102 , 7314 –7323 . 23 R. F. W. Bader , Acc. Chem. Res. 1985 , 18 , 9 –15 . 24 T. H. Tang , R. F. W. Bader , P. MacDougall , Inorg. Chem. 1985 , 24 , 2047 –2053 . 25 F. Biegler-König , R. F. W. Bader , T. H. Tang , J. Comput. Chem. 1982 , 3 , 317 –328 . 26 F. Biegler-König , J. Schönbohm , D. Bayles , J. Comput. Chem. 2001 , 22 , 545 –559 . 27 F. Biegler-König , J. Schönbohm , J. Comput. Chem. 2002 , 23 , 1489 –1494 .12370951 28 J. A. Dobado , H. Martinez-Garcia , J. Molina , M. R. Sundberg , J. Am. Chem. Soc. 2000 , 122 , 1144 –1149 . 29 J. Molina , J. A. Dobado , Theor. Chem. Acc. 2001 , 105 , 328 –337 . 30 S. K. Ignatov , N. H. Rees , B. R. Tyrrell , S. R. Dubberley , A. G. Razuvaev , P. Mountford , G. I. Nikonov , Chem. Eur. J. 2004 , 10 , 4991 –4999 .15372659 31 W. Nakanishi , T. Nakamoto , S. Hayashi , T. Sasamori , N. Tokitoh , Chem. Eur. J. 2007 , 13 , 255 –268 .17066493 32 R. J. Boyd , S. C. Choi , Chem. Phys. Lett. 1986 , 129 , 62 –65 . 33 M. T. Carroll , R. F. W. Bader , Mol. Phys. 1988 , 65 , 695 –722 . 34 S. J. Grabowski , J. Phys. Chem. A 2001 , 105 , 10739 –10746 . 35 E. Espinosa , I. Alkorta , J. Elguero , E. Molins , J. Chem. Phys. 2002 , 117 , 5529 –5542 . 36 I. Rozas , I. Alkorta , J. Elguero , J. Am. Chem. Soc. 2000 , 122 , 11154 –11161 . 37 M. Domagała , S. Grabowski , K. Urbaniak , G. Mloston , J. Phys. Chem. A 2003 , 107 , 2730 –2736 . 38 S. Grabowski , W. A. Sokalski , J. Leszczynski , J. Phys. Chem. A 2005 , 109 , 4331 –4341 .16833763 39 M. Domagala , S. Grabowski , J. Phys. Chem. A 2005 , 109 , 5683 –5688 .16833901 40 Dots are usually employed to show BCPs in molecular graphs. Therefore, A−⋅−B would be more suitable to describe the BP with a BCP. Nevertheless, A−✶−B is employed to emphasize the existence of a BCP on the BP in question in our case. 41 W. Nakanishi , S. Hayashi , M. B. Pitak , M. B. Hursthouse , S. J. Coles , J. Phys. Chem. A 2011 , 115 , 11775 –11787 .21902277 42 Y. Sugibayashi , S. Hayashi , W. Nakanishi , Phys. Chem. Chem. Phys. 2015 , 17 , 28879 –28891 .26451525 43 S. Hayashi , Y. Sugibayashi , W. Nakanishi , Phys. Chem. Chem. Phys. 2016 , 18 , 9948 –9960 .26818845 44 Y. Sugibayashi , S. Hayashi , W. Nakanishi , ChemPhysChem 2016 , 17 , 1804 –1816 . 45 S. Hayashi , Y. Sugibayashi , W. Nakanishi , RSC Adv. 2016 , 6 , 49651 –49660 . 46 S. Hayashi , Y. Sugibayashi , W. Nakanishi , RSC Adv. 2017 , 7 , 31858 –31865 . 47 S. Hayashi , K. Nagata , S. Otsuki , W. Nakanishi , J. Phys. Chem. A 2017 , 121 , 2482 –2496 .28257204 48 Y. Tsubomoto , S. Hayashi , W. Nakanishi , L. K. Mapp , S. J. Coles , RSC Adv. 2018 , 8 , 9651 –9660 . 49 W. Nakanishi , S. Hayashi , Int. J. Quantum Chem. 2018 , 118 , e25590 . 50 K. Brandhorst , J. Grunenberg , J. Chem. Phys. 2010 , 132 , 184101 –184107 . 51 K. Brandhorst , J. Grunenberg , Chem. Soc. Rev. 2008 , 37 , 1558 –1567 .18648681 52 J. Grunenberg , J. Am. Chem. Soc. 2004 , 126 , 16310 –16311 .15600318 53 J. Grunenberg , Angew. Chem. Int. Ed. 2001 , 40 , 4027 –4029 ; Angew. Chem. 2001 , 113 , 4150 –4153 . 54 J. Böhnke , H. Braunschweig , P. Constantinidis , T. Dellermann , W. C. Ewing , I. Fischer , K. Hammond , F. Hupp , J. Mies , H.-C. Schmitt , A. Vargas , J. Am. Chem. Soc. 2015 , 137 , 1766 –1769 .25625817 55 The Cij values and the coordinates corresponding to Cij were calculated by using the Compliance 3.0.2 program released by Grunenberg and Brandhorst. 56 Gaussian 09, Revision D.01, M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, Ö. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, D. J. Fox, Gaussian, Inc., Wallingford, CT, 2009. 57 T. Noro , M. Sekiya , T. Koga , Theor. Chem. Acc. 2012 , 131 , 1124 –1128 . 58 C. Møller , M. S. Plesset , Phys. Rev. 1934 , 46 , 618 –622 . 59 J. Gauss , J. Chem. Phys. 1993 , 99 , 3629 –3643 . 60 J. Gauss , Ber. Bunsenges, Phys. Chem. 1995 , 99 , 1001 –1008 . 61 F. Biegler-König , J. Comput. Chem. 2000 , 21 , 1040 –1048 . 62 A. Bondi , J. Phys. Chem. 1964 , 68 , 441 –451 . 63 The correlation is ΔE=1.21+32.41×Δr (R c 2=0.876), if the data for H3N−✶−HF are omitted. The results show that the HB adducts become more stable as the distances are shortened, although H3N−✶−HF is much more stable than that expected from the correlation.