本帖最后由 dahai7719 于 2023-5-16 23:48 编辑 Sobereva老师您好: 不好意思,当时编辑的时候没有乱,我又重新整理了: 1:我对这个复合物分子(figure1 图)用Gaussian做了优化:opt=(gdiis,maxstep=3,notrust,recalc=5) freq b3lyp/6-311g(d,p) em=gd3bj] 2:然后用M06-2X方法,结合您在509贴子中给def2以ma-方式加弥散函数的方法基组ma-TZVPP做了单点,得到fch文件后,做了IRI分析(结果见figure2), 发现里面有很强的H键和弱的范德华作用,然后做了拓扑分析,找了所有的临界点,然后按照您帖子522中基于键临界点的电子密度估计内氢键键能 的方法计算了 figure3中(3,-1)关键点H键的键能(见figure4),发现有-19kcal/mol的H键键能,判断应该是有电子转移,我想继续做ETS-NOCV分析,把两边的NH3OH作为一个片段,中间C2N4O2做一个片段,再做这两个片段的单点时电荷和自选多重度拿不准怎么设,是按照自由基呢,还是按照正负离子把NH3OH 按照2个正电荷、C2N4O2按2个负电荷来设置比较合理? 按照文献中(figure1)有明显的正负部分,和自己计算的H键键能(figure4中有-19kcal/mol的H键键能) 也应该有电荷的转移, 是不是应该按照正-负两个电荷分别给两个片段来设置,谢谢您在百忙中能给与指点! 3:另外还有一个问题不太明白,如果用上面M06-2X方法结合给TZVPP加了弥散和D3校正,如果再加上Counterpoise计算单点后,再用522中基于键临界点的电子密度估计内氢键键能的方法估算H键(拟合公式中已经考虑的BSSE的问题),会不会和单点计算中的Counterpoise有重复校正的问题?谢谢 |
先整理好帖子内容,别留一堆未解析的格式标签,严重妨碍阅读 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-27 17:52 , Processed in 0.231355 second(s), 26 queries , Gzip On.