用Multiwfn通过ETS-NOCV分析时复合物片段的电荷和自旋多重度设置问题

dahai7719 · 发表于 Post on 2023-5-16 21:20:12

本帖最后由 dahai7719 于 2023-5-16 23:45 编辑

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有重复校正的问题？谢谢

sobereva · 发表于 Post on 2023-5-16 22:21:58

先整理好帖子内容，别留一堆未解析的格式标签，严重妨碍阅读

dahai7719 · 发表于 Post on 2023-5-16 23:35:58

本帖最后由 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有重复校正的问题？谢谢

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[波函数分析求助] 用Multiwfn通过ETS-NOCV分析时复合物片段的电荷和自旋多重度设置问题