wzkchem5 发表于 2021-2-8 13:15 好的，谢谢老师的指点 |
xiaos 发表于 2021-2-7 15:55 注意金属离子必须用显式溶剂化。比如假设你是计算水里的吸附，又假设Pb2+是和x个水配位（应该是6个，你自己查文献确认一下Pb2+是不是配6个水），吸附以后还剩y个水配位，那么计算方法应该是E(abs)=E(吸附后的复合物)+(x-y)E(H2O)-E([Pb(OH2)x]2+)-E(fragment)。文献的方法仅适用于没有使用显式溶剂化的情况，也可能是文献默认使用者知道在使用显式溶剂化的情况下应该怎么改。如果不清楚y是多少，应该计算Gibbs自由能，找到配多少个水最稳定。 |
sobereva 发表于 2021-2-7 21:57 好的 谢谢sob老师 |
xiaos 发表于 2021-2-7 15:55 仔细看此文 谈谈该从Gaussian输出文件中的什么地方读电子能量 http://sobereva.com/488（http://bbs.keinsci.com/thread-13450-1-1.html） 铅离子就照常算就完了，没有任何特殊性 |
wzkchem5 发表于 2021-2-7 16:08 好的好的，非常感谢您的回复 |
xiaos 发表于 2021-2-7 15:55 再者，如果你用的基组不是特别大，比如小于def2TZVP，那么一般还需要计算BSSE |
xiaos 发表于 2021-2-7 15:55 那个不叫HF能量。虽然在输出文件里写成HF=xxx的，但是如果你是用DFT算的，那就叫DFT能量，不叫HF能量。 只有当你用的是HF或DFT，且用的不是双杂化泛函的时候，才可以这样做。另外建议优化之后加大基组（有时还可以改用更好的泛函）算个单点能，用这个单点能作差，不要直接用结构优化最后一步的能量作差。 |
本帖最后由 xiaos 于 2021-2-7 16:00 编辑 sobereva 发表于 2020-4-21 18:09 sob老师您好，请问计算化学吸附能的时候，分别优化吸附前和吸附后的两分子结构，是否可以取分子的HF能量求差获得？在文献中也找到一种吸附能的计算方法：E(abs)=E(total)-E(fragment)-E(atom)；E(total)和E(fragment)分别为吸附和未吸附的能量，E(atom)为原子的能量。但在我的体系里，吸附的是铅离子，所以想请问sob老师，这两种方法哪个更合理？如果使用第二种方法，那铅离子的能量如何计算？ |
sobereva 发表于 2020-4-23 06:54 谢谢社长，根据博文的指导，ESP图在Multiwfn中已经画出来了，现在正在研究怎么弄到VMD中。 |
rabbitkiller 发表于 2020-4-22 17:35 让Multiwfn调用cubegen节约时间 Multiwfn现已可以调用cubegen使静电势分析耗时有飞跃式的下降！ http://sobereva.com/435（http://bbs.keinsci.com/thread-10682-1-1.html） 对ORCA用户大幅节约Multiwfn做分子表面静电势分析耗时的方法 http://bbs.keinsci.com/thread-16499-1-1.html 不要用QZ级别算静电势，毫无必要，白浪费时间。 def2-SVP糙了些 |
sobereva 发表于 2020-4-21 18:09 请教社长，我看了绘制两个分子的范德华表面穿透图的方法（http://sobereva.com/443），想按照您介绍的方法画图。 我是用ORCA的B3LYP/def2-SVP做的优化（70多个碳原子），然后B3LYP/def2-QZVPP计算的单点能。用Multiwfn绘制表面静电势图用的就是QZVPP计算出来的波函数，开始绘图的时候Multiwfn就停在1%进度不动了，估计是数据太大，处理慢。 所以我想请问下，可否用B3LYP/def2-SVP计算单点能的波函数绘制静电势图？ 谢谢！ |
sobereva 发表于 2020-4-21 18:09 明白了，谢谢社长~ |
|
这个体系势能面极小点应该非常多，氮原子也不是决定吸附的唯一因素，所以“蓝色氮原子上的吸附能”这种说法有些太牵强。推荐先按照此文绘制一下石墨烯部分的表面静电势，估计一下吸附特征 使用Multiwfn+VMD快速地绘制静电势着色的分子范德华表面图和分子间穿透图（含视频演示） http://sobereva.com/443（http://bbs.keinsci.com/thread-11080-1-1.html） 可。 合适 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-25 05:33 , Processed in 0.169832 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites