学到了 |
实际跑了一下,genmer.ini里rmax从你原先的15设小到7,把完全没必要写的rmin去掉了。得到的结构如下所示,可见得到的初始结构就已经很合理了,Li被DMSO充分包着![]() genmer产生的10个结构: ![]() 然后用molclus调用xtb 6.4.0在GFN1-xTB下批量优化,跑出来的结构都很合理,Li与氧或硫显著结合。我为了省时间就跑了10个结构,其中能量最低的如下,Li+与三个氧和两个硫结合(其它结构也有与四个氧、或四个氧一个硫结合的)。之后再用DFT refine一下就可以了 ![]() isostat产生的文件在此: ![]() 这个体系不是非得跑AIMD。用genmer就完全可以很好地搞定 |
你拿UFF=QEq算这个根本不靠谱,建议改用PM7或者GFN-xTB做预优化 另外,genmer.ini里建议把rmax设小点,使得genmer产生的结构中Li周围总是有DMSO包围 |
这种建议跑AIMD |
喵星大佬 发表于 2021-3-14 16:25 可是算出来就是这样呀。是100个团簇还不够多吗 |
显然不合理。。。。这样的话你的Li一边不就没有显式溶剂了么 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-14 13:12 , Processed in 2.826250 second(s), 30 queries , Gzip On.