关于团簇结构与有机分子间吸附能计算时的BSSE处理问题

卡开发发

喵星大佬 发表于 2021-8-19 19:12
我觉得出现这种问题应该先检查模型和SCF的问题

确实如此，有时候我们总认为其他人会自己多检查几遍才会到论坛上来问，不过实际遇到的低级错误是不少（比如POTCAR顺序写错这种事情就发生了n次）。

telunsu

卡开发发 发表于 2021-8-19 19:26
确实如此，有时候我们总认为其他人会自己多检查几遍才会到论坛上来问，不过实际遇到的低级错误是不少（比 ...

实在抱歉，引起大家的误会了，我只是拿这个楼主的公式举了个例子，我实际是在用DMOL3计算MIL-101与水、呋喃的相互作用
DMol3 BSSE Energy Calculations (Ha)
Total energy A(B) = -232.219935
Total energy B(A) = -1756.972062
Total energy A    = -232.205937
Total energy B    = -1756.970660
BSSE correction energy = -0.015400
Total energy BSSE corrected = -1989.452953
Counterpoise corrected interaction energy = -0.276356

还是想问下： 1，是否可以Total energy BSSE corrected = -1989.452953代替E_1+2和E_BSSE 这两项，计算总能量？

因为我有时发现Total energy BSSE corrected =E_AB -BSSE correction energy,而有时相差又有点大,体系含Cr，spin为2的时候能量最低，选pw91 是参照了别人之前在JPCC的工作
2，如果不能代替，那么Total energy BSSE corrected  代表什么？十分感谢各位老师的回答！

# Task parameters
Calculate                     optimize
Opt_energy_convergence        1.0000e-005
Opt_gradient_convergence      2.0000e-003 A
Opt_displacement_convergence  5.0000e-003 A
Opt_iterations                800
Opt_max_displacement          0.3000 A

Symmetry                      off
Max_memory                    204800
File_usage                    smart
Scf_density_convergence       1.000000e-005
Scf_charge_mixing             2.000000e-001
Scf_spin_mixing               5.000000e-001
Scf_diis                      6 pulay
Scf_iterations                300

# Electronic parameters
Spin_polarization             unrestricted
Spin                          2
Charge                        0
Basis                         dnp
basis_version                 basfile_v4.4
Pseudopotential               dspp
Functional                    gga(p91)
Aux_density                   octupole
Dftd                          OBS
Integration_grid              medium
Occupation                    thermal 0.0020
Cutoff_Global                 4.4000 angstrom

卡开发发

telunsu 发表于 2021-8-19 19:43
实在抱歉，引起大家的误会了，我只是拿这个楼主的公式举了个例子，我实际是在用DMOL3计算MIL-101与水、呋 ...

按照前面给的BSSE correction energy去做应该肯定没问题，这个值和上面的A(B)和A，B(A)和B之间的关系算出来是没问题的，能确保这个就可了。Total energy BSSE corrected他是咋定义的我不太清楚。