从能量和信息理论视角理解单取代烷烃的异构化

Transcription

1 物理化学学报 cta Phys. -Chim. Sin. 2018, 34 (3), [rticle] doi: /PKU.WHXB 从能量和信息理论视角理解单取代烷烃的异构化 钟爱国 *, 李嵘嵘, 洪琴, 张杰, 陈丹 台州学院医药化工学院, 浙江台州 摘要 : 对直链烷烃和支链烷烃的相对稳定性统一的解释仍然没有定论, 并且一直在进行着 以单取代的烷烃体系 CnH2n+1 R (n = 3, 4, 5, 6;R = OH, OCH3, NH2, NO2, F, Cl, CN, CHO) 为例, 本文对支链效应的有效性和本质进行了研究 与传统的基于轨道的描述不同的是, 本文采用了密度泛函理论的总能量和基于新能量分配方案的能量分量 [ 见 Liu, S. B. J. Chem. Phys. 2007, 126, ] 新型能量分解方法计算结果表明, 静电效应和立体效应等对支链效应的存在都起着重要作用, 但是它们均不能单独用来解释支链效应的本质 用双变量 ( 静电势和空间位阻 ) 组合, 发现单取代烷烃衍生物的异构化反应主要影响因子是静电势作用, 空间位阻效应的影响是次要的 此外还发现了香农熵差与 Fisher 信息差之间的线性关系, 未能发现总能量差或者分能量差值和 Fisher 信息或者 Shannon 熵之间的关系 这与前人发现是一致的 关键词 : 分子稳定性 ; 空间位阻 ; 密度泛函理论 ; 能量分项中图分类号 :O641 Understanding the Isomerization of Monosubstituted lkanes from Energetic and Information-Theoretic Perspectives ZHONG iguo *, LI Rongrong, HONG Qin, ZHNG Jie, CHEN Dan College of Medicine and Chemical Engineering, Taizhou University, Taizhou , Zhejiang Province, P. R. China. bstract: unified explanation for the relative stability of linear and branched alkanes is still lacking, and research in this direction is ongoing. Unlike the conventional orbital-based description, we have employed the density functional theory (DFT) based on the total energy and its components as well as a newly proposed energy partitioning scheme [Liu, S. B. J. Chem. Phys. 2007, 126, ]. Taking monosubstituted hydrocarbons CnH2n+1 R (n = 3, 4, 5, 6; R = OH, OCH3, NH2, NO2, F, Cl, CN, CHO) as examples, we have investigated the molecular origin and nature within the framework of the DFT. Our computational results revealed that no such a single energy component that dictates the transformation of the mono-substituted alkane derivatives. By the binary fit of the electrostatic potential and steric hindrance from the new energy decomposition scheme, we unraveled that isomerization is mainly governed by the electrostatic potential, while the steric effect has a minor influence. Moreover, we established a linear relationship between the Shannon entropy difference and the Fisher information difference, which is accordance with previous findings. Further, there was no linear relationship between Ee and the Fisher information or Shannon entropy. Key Words: Molecular stability; Steric hindrance; Density functional theory; Energy component 1 Introduction It is widely accepted both by the experimental measurements 1 and high-level ab initio computations 2 that branched alkane hydrocarbons are generally Received: July 19, 2017; Revised: ugust 16, 2017; ccepted: ugust 24, 2017; Published online: ugust 30, * Corresponding author. zhongaiguo@tzc.edu.cn. Tel.: The project was supported by the Natural Science Foundation of Zhejiang Province, China (LY15B030001, LY15B060001) and Zhejiang Province Undergraduate Student Innovation Project, China (2017R430014). 浙江省自然科学基金 (LY15B030001; LY15B060001) 和浙江省大学生科技创新 (2017R430014) 资助项目 Editorial office of cta Physico-Chimica Sinica

2 304 cta Physico-Chimica Sinica Vol.34 thermodynamically more stable than its straight-chain counterpart. Nevertheless, this phenomenon is somewhat counterintuitive at first glance. In chemistry, steric effects are a consequence of the space required to accommodate the atoms and groups within a molecule, and are often thought to be predominated by repulsive forces arising from overlapping electron densities. In this sense, more branched alkanes normally engender more steric destabilization energy than the normal ones. However, we have just touched upon one side of the coin. Simultaneously, there exist some attracting interactions between alky groups, stemming from electron correlation. Thus, inclusion of attractive interactions such as London dispersion forces is indispensable to understand electronic and geometrical structures. Finally, these two forces coexist that makes such a phenomenon a hard nut to crack. To provide a unified explanation for this abnormal phenomenon, a surge of interest has been aroused. Schleyer et al. pointed out that intramolecular 1,3-alkyl-alkyl interactions make positive contributions to the molecular stability of hydrocarbons 3,4. Nonetheless, this protobranching or prototypical branching model was severely criticized by Gronert. He 5 argued that it cannot be rationalized by a through-space model of dispersive London forces. Pitzer and Catalano 6 have also provided an explanation for the relative stability of branched alkanes based on dispersion interactions. But this argument does not hold because the Hartree-Fock (HF) approach intrinsically has no electron correlation but is able to predict the correct stability trend 3. In this regard, the overall stabilization effect does not originate from a simple physical picture. In the literature, there are some other orbital-based intramolecular stabilization explanations. Ma and Inagaki 7 proposed an orbital phase rational and found that σ cross conjugation between two C H bonds and one C C bond in an antiperiplanar conformation in branched alkanes results in greater orbital stabilization than straight-chain σ orbital conjugation 8. Kemnitz et al. employed natural bond orbital (NBO) analysis 9 and valence bond theory to argue that 10 germinal carbon-carbon σ to σ* second-order donor-acceptor orbital interactions account for the branching stability. However, this hyperconjugation model has an inherent drawback since this perturbative type of energy analysis is referenced to a hypothetical localized state that does not analyze/partition the total energy of each molecule. The objective of this paper is two-fold. First, we aim to find out the physical differences of mono-substituted linear and branched alkanes alkane derivatives (see Fig.1 for details). It is intuitive that the isomerization from a linear to a branched alkane involves differences in steric energy. Like many concepts in chemistry, actually there is no exact definition of steric effect. Liu has proposed that the Weizäcker kinetic energy can be employed as a novel definition of steric effect in DFT 11,12. This density-based quantification of the steric effect differs from its conventional wavefunction-based counterpart. Fig.1 Schematic representations of mono-substituted alkane derivatives under study.

3 No.3 doi: /PKU.WHXB The latter definition is resulted only from orthogonality and exchange antisymmetry of the wavefunction (the Fermi hole) 8. This density-based steric definition is quite in accordance with the early definition by Weisskopf who ascribed the steric repulsion to kinetic energy pressure 13. In this work, we employ two energy decomposition schemes in Kohn-Sham density functional theory (KS-DFT) to understand the total physical differences between branched and linear alkanes. Isolation of a specific intramolecular interaction is beyond the scope of this paper. Second, other two density-based key quantities from information theory, Shannon entropy and Fisher information are taken into consideration at molecular level to gain insights into the transformation of mono-substituted alkane derivatives. The questions we will address are (i) whether the branched effect is always valid and (ii) what effect(s) dominantly contribute(s) to the validity of this branched effect. This paper is mainly laid out as follows. In section 2, we outline the total energy decomposition schemes, followed by computational details. fter the results and discussion in section 3, a brief summary is given in section 4. 2 Theory 2.1 n outline of the two energy partition schemes In conventional KS-DFT, the total energy of a system, E[ρ], can be cast in the form of three independent contributions Ts[ρ], Ee[ρ], and Exc[ρ] (Eq.(1)) 14, E[ρ] = Ts[ρ] + Ee[ρ] + Exc[ρ] (1) where the first term Ts[ρ] signifies the non-interacting one-electron kinetic energy, the second term Ee[ρ] is the electrostatic potential, as shown in Eq.(2), a combination of the nuclear to electron attraction stabilization, Vne[ρ], the classical inter-electron Coulombic repulsion, J[ρ], and the nuclear-nuclear repulsion, Vnn, and the last term Exc[ρ] is a sum of the exchange and correlation energy components. Ee[ρ] = Vne[ρ] + J[ρ] + Vnn (2) Unlike the definition in Eq.(1), Liu 12 proposed a novel energy partition scheme (Eq.(3)), to unambiguously quantify the omnipresent steric repulsion, as defined by the Weizäcker kinetic energy, Tw[ρ] with and E[ρ] = Es[ρ] + Ee[ρ] + Eq[ρ] (3) Es[ρ] TW[ρ] = (1/8) ρ(r) 2 /ρ(r)dr (4) Eq[ρ] = Exc[ρ] + Ts[ρ] Es[ρ] (5) where Es and Eq are the steric effect and fermionic quantum effects, respectively. The electrostatic potential Ee is exactly the same as that defined in Eq.(2), ρ(r) is the total electron density distribution, and ρ(r) is the first-order derivative. It is worth mentioning that the Weizäcker kinetic functional is the exact kinetic functional for one-electron atoms and two-electron Hartree-Fock atoms. It has become an essential ingredient in orbital-free DFT 15. This novel density-based definition of steric effect is totally different from its wavefunction-based counterpart, resulted from the Pauli Exclusion Principle that solely accounts for the electrons possessing the same spin. In this work, these two energy partition schemes will be utilized to analyze the isomerization energy between the linear and branched hydrocarbon derivatives, and its components to find out which component is the predominant term that governs the transformation of an alkane. The new scheme has recently been applied to a couple of systems 16 25, among which are bond rotation barriers, cis-effect, anomeric effect, SN2 reaction barriers, water clusters, beryllium bonding interactions, and so on. key point from these different investigations is that the electrostatic potential is the leading term while other contributions from the steric repulsion and quantum exchange-correlation interactions play minor but nontrivial roles. It is worthwhile to mention that energy partitioning scheme with Weizäcker functional (in our case Es[ρ]) as the dominant term was shown elsewhere as well Shannon entropy and Fisher information In information theory, Shannon entropy SS is defined as the following 27 : (6) SS[ ] ( r)ln ( r)d r ss( r)dr where SS(r) is the Shannon entropy density and ρ(r) is the total electron density of a molecular system, satisfying the following condition: ( r)dr N (7) with N as the total number of electrons in the system. In addition, Fisher information, IF, is defined as follows 28 : 2 ( r ) IF[ ]= d r if( )d ( ) r r (8) r where if(r) is the Fisher information density and ρ(r) is the density gradient. To calculate the atomic values of Shannon entropy in a molecule, Eq.(11) can be recast as follows: S [ ] S [ ] (9) S where S S S [ ] ( r)ln ( r)dr (10) with Ω the atomic basin of atom. Bader s zero-flux partition condition is employed to partition atoms in molecules 29. The same is true for Fisher information, I [ ] I [ ] (11) F F

4 306 cta Physico-Chimica Sinica Vol.34 with 2 ( r ) I F [ ] dr (12) ( r ) 2.3 Computational details ll mono-substituted hydrocarbons CnH2n+1 R (n = 3, 4, 5, 6; R = OH, OCH3, NH2, NO2, F, Cl, CN, CHO) (see Fig.1 for details) were fully optimized at the M06-2X/6-311+G(d,p) 28 level of theory. No symmetry constraint was imposed to allow for full variational degrees of freedom. ll local minima on the hyper-potential surface were further verified by vibrational frequency calculations (e.g., no imaginary frequency). DFT total energy evaluations and steric energy decomposition calculations were carried out by employing the M06-2X functional with Dunning s augmented correlation-consistent valence polarized triple-ζ basis set aug-cc-pvtz (denoted as avtz) in NWChem 6.6, 30 a freely accessible computational chemistry package. Furthermore, a relatively larger basis set aug-cc-pvqz (avqz), was taken into consideration to guarantee the rationalization of our results. M06-2X is superior to B3LYP 31,32 in many aspects, especially in the description of ubiquitous dispersion forces, irrespective of its over-parameterization. The tight self-consistent field (SCF) convergence criteria and ultrafine integration grids were employed throughout to rule out numerical problems. Suffice to note that the M06-2X functional gives isomerization energies very close to those predicted by the gold standard CCSD(T) in quantum chemistry and experiment. 2 ll MP2 and CCSD(T) computations were executed with the frozen core (FC) technique to reduce the computational cost without compromising the accuracy too much. One more point to mention is that zero-point corrections were not considered in our computations. The linear alkane derivatives are taken as reference so that the isomerization energy is negative in all cases. Unless otherwise stated, all energy differences and its components are in unit of kj mol 1. dditionally, we employed a multifunctional wavefunction analyzer, Multiwfn 33,34 to calculate the Shannon entropy and Fisher information at molecular level by utilizing the M06-2X/aVTZ molecular wavefunction as an input file in molden format. To eliminate the numeric instabilities and guarantee the rationalization of our results, we adopted 150 points in radial and 1454 points on spherical surface of Becke s numerical quadrature. 3 Results and discussion 3.1 Validation Table 1 exhibits the isomerization energy of a serials of mono-substituted alkanes C3H7-R with both electron-withdrawing and electron-donating groups R = OH, OCH3, NH2, NO2, F, Cl, CN, and CHO. Here we have chosen two typical density functionals M06-2X and B3LYP. It is well-established that B3LYP inherently fails to depict the dispersion forces which have a far-reaching effect on the molecular stability. In this sense, for comparison we have adopted some B3LYP variants with a portion of empirical dispersion, such as Grimme s density-independent, atomic pair wise corrections, D3BJ (B3LYP_D3BJ), 35 the exchange-hole dipole moment (B3LYP_XDM) 36 theory proposed by Becke and Johnson, and the dispersion-core-potential (B3LYP_DCP) 37 model, which contains both local and semilocal terms. ll the structures both reactants and products were optimized by using each method with Pople s standard basis set G(d,p) while all post-hf optimizations were carried out at the MP2/6-311+G(d,p) level of theory. The avtz basis set was employed to obtain more accurate energy differences. n even larger basis set avqz was employed to confirm the rationalization of the avtz results. The CCSD(T)/aVTZ results were taken as a reference to evaluate the overall performances of HF, MP2 and DFT B3LYP and M06-2X. Two statistical parameters, mean signed error (MSE) and mean absolute error (ME) were adopted. quick inspection of Table 1 can give rise to some points in order: (i) all the methods under consideration can give a qualitative prediction of the isomerization energies (negative in Table 1 Total energy difference (in kj mol 1 ) of the isomerization reaction of C3 serials a. R M06-2X M06-2X b B3LYP B3LYP_D3BJ B3LYP_XDM B3LYP_DCP HF MP2 CCSD(T) OH OCH NH NO F Cl n.a. c CN CHO ME d MSE e a Dunning s avtz basis set is adopted to obtain the total energy differences and the CCSD(T) data are taken as reference. b Dunning's avqz basis set is used here. c n.a. means not available due to the lack of dispersion-corrected potential for Cl. d ME is mean absolute error. e MSE is mean signed error.

5 No.3 doi: /PKU.WHXB sign), indicating that electron correlation is not the only factor that is responsible for the relative stability of mono-substituted alkanes since HF intrinsically has no electron correlation; (ii) compared with CCSD(T), MP2 and B3LYP_DCP overestimate the isomerization energies while the rest approaches adopted [M06-2X, B3LYP, B3LYP_D3BJ, B3LYP_XDM, and HF] underestimate the data; (iii) dispersion-corrected B3LYP variants improve the energy differences to some extent, with mean absolute errors (MEs) decreasing from 3.4 kj mol 1 of B3LYP to 2.2 kj mol 1 of B3LYP_D3BJ, 2.0 kj mol 1 of B3LYP_XDM, and 2.3 kj mol 1 of B3LYP_DCP, respectively; however, they are still far from satisfactory when compared with CCSD(T) and entail more substantial improvement; (iv) compared with the CCSD(T) data, M06-2X has the best performance if combined with the avtz basis set. In a nutshell, we have screened a rational and cost-efficient combination of method and basis set M06-2X/aVTZ and it will be employed for further energy decomposition analysis. 3.2 Total energy decompositions Table 2 shows the total energy difference ( Etot) via the conventional and newly proposed energy partition schemes and its components: kinetic energy ( Ts), steric hindrance ( Es), exchange ( Ex), correlation ( Ec), electrostatic potential ( Ee) and quantum effect ( Eq), which are obtained at the M06-2X/aVTZ level of theory, with the geometries completely optimized at the M06-2X/6-311+G(d,p) level of theory. From Table 2, one can readily see that the total energy difference is negative in sign, indicating that the branched alkane derivatives are energetically favorable. mong all the energy components, steric potential ( Es) and exchangecorrelation potential (exceptions do exist in C3 serials) possess negative values, meaning that they make positive contributions to the total energy difference, while they are largely compensated by the quantum effect ( Eq) and kinetic energy ( Ts), respectively. n intriguing phenomenon has been discovered that for the C3 serials the exchange Table 2 Total energy difference ( Etot) and its components with kinetic energy ( Ts), steric hindrance ( Es), exchange ( Ex), correlation ( Ec), electrostatic potential ( Ee) and quantum effect ( Eq), evaluated at the M06-2X/aVTZ level a. Reaction Ts Es Ex Ec Exc Ee Eq Etot C3-serials C4-serials to be continued

6 308 cta Physico-Chimica Sinica Vol.34 continued Table 2 Reactions Ts Es Ex Ec Exc Ee Eq Etot C5-serials to be continued

7 No.3 doi: /PKU.WHXB continued Table 2 Reaction Ts Es Ex Ec Exc Ee Eq Etot C6-serials to be continued

8 310 cta Physico-Chimica Sinica Vol.34 continued Table 2 Reaction Ts Es Ex Ec Exc Ee Eq Etot to be continued

9 No.3 doi: /PKU.WHXB continued Table 2 Reaction Ts Es Ex Ec Exc Ee Eq Etot a The conventional KS-DFT total energy ΔE = Ts + Exc + Ee and the newly proposed total energy partition scheme ΔE = Es + Eq + Ee. The total energy difference and its components are in kj mol 1. potential ( Ex) is larger than the correlation ( Ec) in absolute values and the results are reverse for these two energy components in the rest serials. It is worthwhile to point out that only the total effect of exchange and correlation potential has a physical meaning. Meanwhile, we find that there exists a strong linear relationship between the quantum effect and steric hindrance, with the correlation coefficient R 2 = 0.99 (Fig.2a), serving as a confirmation of no systematic errors of our computational results. Similar trends have been observed in previous studies Furthermore, we aim to find out if there exists a single energy component that dictates the isomerization of alkane derivatives. To make that happen, in Fig.2(b, c), we have plotted the electrostatic potential ( Ee) and the kinetic energy ( Ts) vs the total energy difference ( Etot), respectively. The corresponding correlation coefficient R 2 is only 0.65 and 0.71 while those for the remaining energy components are less than 0.5 and ignored. In a previous work for alkanes without Fig.2 Correlations between (a) steric potential and quantum effect (y = 0.98 x , R 2 = 0.99); (b) total energy difference and electrostatic potential (y = 0.40x 8.61, R 2 = 0.64), and (c) total energy difference and kinetic energy (y = 0.68 x 0.84, R 2 = 0.71), the y axis spans from (b) to (c). ll data are obtained at the M06-2X/aVTZ level of theory with all the geometries completely optimized at the M06-2X/6-311+G(d,p) level of theory.

10 312 cta Physico-Chimica Sinica Vol.34 substituents, Ee can be linearly related to Etot, but that's a special case 21. Up to now, we have revealed that no such a single energy component that governs the isomerization of alkane derivatives, which is reminiscent of the resultant force in nature for such a phenomenon. Furthermore, we want to look into the relative importance of energy components contributing the total energy, thus rendering the so-called branching effect. To this end, we have employed a two-variable strategy. In Fig.3, we plotted the total energy difference vs the two components: the electrostatic potential ( Ee) and steric hindrance ( Es), with R 2 = Moreover, we can find out the relative significance of the two energy components via the coefficients. It is clearly shown that the electrostatic potential plays a predominant role while the steric effect has some minor effect in the transformation of mono-substituted alkane derivatives, with the coefficients 0.33 and 0.04, respectively. In other words, though steric hindrance is engendered due to the branching of alkanes, this portion is overwhelmingly compensated by the electrostatic potential. 3.3 Shannon entropy and Fisher information Having touched upon one side of the coin, the total energy and its components, next we will switch our gear from an information-theoretic viewpoint by analyzing the M06-2X/aVTZ wavefunction. Here in this work, suffice to note that our Shannon entropy and Fisher information data are based on the electron probability density, rather than the shape function. From the information theory point of view, Shannon entropy measures the spatial delocalization of the electronic density, and Fisher information measures its sharpness or concentration. They sound to be vastly different measurements but apparently in electronic systems, these two density-based quantities are not absolutely independent. They are strongly correlated to each other, as theoretically proven earlier 38 and numerically shown 39,40. We have plotted the Shannon entropy difference and Fisher information difference in Fig.4 with the Fig.4 Linear correlation between Shannon entropy difference and Fisher information difference. correlation coefficient R 2 = Once again, we have verified that Shannon entropy and Fisher information are interrelated as a windfall of this work. The linear correlations of these two quantities can be used an evidence of weak interactions as previously shown in beryllium bonds 24, though it is not observed here. We have also plotted the Shannon entropy difference vs the electrostatic potential energy difference ΔEe or Fisher information difference vs total energy difference ΔEtot with the correlation coefficient R 2 < That the total energy difference ΔEtot does not strongly correlate any of these information-theoretic quantities suggests that the branched alkanes effect is of the complex nature. 4 Concluding remarks Isomerization of mono-substituted alkanes from energetic and information-theoretic perspectives was investigated within the KS-DFT framework. To this end, the KS-DFT total energy and a newly proposed energy partitioning scheme were employed as well as Shannon entropy and Fisher information from information theory. We have found that no such a single energy component dictates the transformation of the mono-substituted alkane derivatives. Molecular stability is governed by a resultant force. From the binary fit, we have unraveled that the electrostatic potential and the steric repulsion are responsible for the relative stability of linear and branched alkanes. Moreover, from the regression coefficient, we have shown that the electrostatic potential has played a predominant role while the steric effect has a trivial role. Shannon entropy and Fisher information at molecular level are linearly correlated, which is in line with our previous work. References (1) Olah, G..; Molnur,. Hydrocarbon Chemistry; Wiley: New York, Fig.3 binary fit of total energy difference ( Etot) with electrostatic potential ( Ee) and steric hindrance ( Es) (2) Miao, J.; Hua, S.; Li, S. Chem. Phys. Lett. 2012, 541, 7.

11 No.3 doi: /PKU.WHXB doi: /j.cplett (3) Wodrich, M. D.; Wannere, C. S.; Mo, Y.; Jarowski, P. D.; Houk, K. N.; Schleyer, P. V. R. Chem. Eur. J. 2007, 13 (27), doi: /chem (4) llen, T. L. J. Chem. Phys. 1959, 31, doi: / (5) Gronert, S. Chem. Eur. J. 2009, 15 (21), doi: /chem (6) Pitzer, K. S.; Catalano, E. J. m. Chem. Soc. 1956, 78, doi: /ja01600a006 (7) Ma, J.; Inagaki, S. J. m. Chem. Soc. 2001, 123, doi: /ja003067v (8) Ito, K. J. m. Chem. Soc. 1953, 75, doi: /ja01106a046 (9) Badenhoop, J. K.; Weinhold, F. J. Chem. Phys. 1997, 107 (14), doi: / (10) Kemnitz, C. R.; Mackey, J. L.; Loewen, M. J.; Hargrove, J. L.; Lewis, J. L.; Hawkins, W. E.; Nielsen,. F. Chem. Eur. J. 2010, 16 (23), doi: /chem (11) Weizsäcker, C. F. V. Z. Phys. 1935, 96, 431. doi: /BF (12) Liu, S. B. J. Chem. Phys. 2007, 126, doi: / (13) Weisskopf, V. F. Science 1975, 187, 605. doi: /science (14) Parr, R. G.; Yang, W. Density Functional Theory of toms and Molecules; Oxford University Press: Oxford, (15) Chakraborty, D.; Kar, S.; Chattaraj, P. K. Phys. Chem. Chem. Phys. 2015, 17, doi: /C5CP00995B (16) Liu, S. B.; Govind, N. J. Phys. Chem. 2008, 112 (29), doi: /jp800376a (17) Liu, S. B.; Govind, N; Pedersen, L.G. J. Chem. Phys. 2008, 129 (9), doi: / (18) Uğur, İ.; Vleeschouwer, F. D.; Tüzün, N.; viyente, V.; Geerlings, P.; Liu, S. B.; yers, P. W.; DeProft, F. J. Phys. Chem. 2009, 113 (30), doi: /jp903371b (19) Torrent-Sucarrat, M.; Liu, S. B.; DeProft, F. J. Phys. Chem. 2009, 113 (15), doi: /jp (20) Liu, S. B.; Hu, H.; Pedersen, L. G. J. Phys. Chem. 2010, 114, doi: /jp101329f (21) Tsirelson, V. G.; Stash,. I.; Liu, S. B. J. Chem. Phys. 2010, 133 (11), doi: / (22) Ess, D. H.; Liu, S. B.; DeProft, F. J. Phys. Chem. 2010, 114, doi: /jp108577g Liu, S. B. cta Phys. -Chim. Sin. 2013, 29 (1), 43. [ 赵东波, 荣春英, 苏曼, 苏文, 尹笃林, 刘述斌. 物理化学学报, 2013, 29 (1), 43.] doi: /PKU.WHXB (24) Wang, Y. J.; Zhao, D. B.; Rong, C. Y.; Liu, S. B. cta Phys. -Chim. Sin. 2013, 29 (10), [ 王友娟, 赵东波, 荣春英, 刘述斌. 物理化学学报, 2013, 29 (10), 2173.] doi: /PKU.WHXB (25) Zhong,. G.; Chen, D.; Li, R. R. Chem. Phys. Lett. 2015, 633, 265. doi: /j.cplett (26) Deb, B. M; Chattaraj, P. K. Phys. Rev. 1989, 39, doi: /PhysRev (27) Shannon, C. E. Bell Syst. Tech. J. 1948, 27, 379. doi: /bltj issue-3 (28) Fisher, R.. Proc. Camb. Philos. Soc. 1925, 22, 700. doi: /S (29) Bader, R. F. W. toms in Molecules: Quantum Theory; Oxford University Press: Oxford, (30) Zhao, Y.; Truhlar, D. G. Theor. Chem. cc. 2008, 120, 215. doi: /s x (31) Valiev, M.; Bylaska, E. J.; Govind, N.; Kowalski, K.; Straatsma,T. P.; van Dam, H. J. J.; Wang, D.; Nieplocha, J.; pra, E.; Windus, T. L.; et al. Comput. Phys. Commun. 2010, 181, doi: /j.cpc (32) Becke,. D. J. Chem. Phys. 1993, 98, doi: / (33) Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B 1988, 37, 785. doi: /PhysRevB (34) Lu, T.; Chen, F. J. Comput. Chem. 2012, 33 (5), 580. doi: /jcc (35) Lu, T. Multiwfn, Version 3.3; Multifunctional Wavefunction nalyzer, 2016, (accessed Sep 18, 2016). (36) Grimme, S.; Ehrlich, S.; Goerigk, L. J. Comput. Chem. 2011, 32, doi: /jcc (37) Becke,. D.; Johnson, E. R. J. Chem. Phys. 2005, 122, doi: /acs.chemrev.5b00533 (38) Heßelmann,. J. Chem. Phys. 2009, 130 (8), doi: / (39) Liu, S. B. J. Chem. Phys. 2007, 126, doi: / (40) Rong, C. Y.; Lu, T.; Liu, S. B. J. Chem. Phys. 2014, 140 (2), doi: / (23) Zhao, D. B.; Rong, C. Y.; Jenkins, S.; Kirk, R. S.; Yin, D. L.;