计算化学公社

 找回密码 Forget password
 注册 Register
Views: 7596|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[Molclus] molclus结合Gaussian计算报错

[复制链接 Copy URL]

51

帖子

0

威望

135

eV
积分
186

Level 3 能力者

老师们好,最近做molclus结合Gaussian做DFT计算时,总是出现这种报错。不知道如何处理,请老师们帮助。

图片5.png (20.01 KB, 下载次数 Times of downloads: 37)

图片5.png

settings.ini

2.66 KB, 下载次数 Times of downloads: 5

template.gjf

79 Bytes, 下载次数 Times of downloads: 13

template2.gjf

103 Bytes, 下载次数 Times of downloads: 7

traj.xyz

29.76 KB, 下载次数 Times of downloads: 11

51

帖子

0

威望

135

eV
积分
186

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2020-4-29 09:55:19 | 只看该作者 Only view this author
sobereva 发表于 2020-4-27 01:39
这样不靠谱。
分子力场精度很低,其相对能量为0的那个,未必经过DFT进一步优化后一定能收敛到全局最小点 ...

知道了,感谢您

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

7#
发表于 Post on 2020-4-27 01:39:57 | 只看该作者 Only view this author
三三33 发表于 2020-4-26 10:37
老师,您好,我现在改画了小分子,并用molclus结合Gaussian做分子力场计算,团簇数目为20,已经计算结束 ...

这样不靠谱。
分子力场精度很低,其相对能量为0的那个,未必经过DFT进一步优化后一定能收敛到全局最小点
用更高精度方法refine的时候应当把低精度方法算出来的能量较低的一批都拿出来继续做
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

51

帖子

0

威望

135

eV
积分
186

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2020-4-26 10:37:50 | 只看该作者 Only view this author
sobereva 发表于 2020-4-23 19:33
槽点太多了,计算方式太野蛮了

这种复合物优化直接用B3LYP等于往枪口上撞,完全是自寻死路。仔细看下文 ...

老师,您好,我现在改画了小分子,并用molclus结合Gaussian做分子力场计算,团簇数目为20,已经计算结束。如果想让molclus调用量子化学程序进一步优化,可以只把能量为0的那个拿出来进行优化吗

51

帖子

0

威望

135

eV
积分
186

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2020-4-24 09:02:28 | 只看该作者 Only view this author
sobereva 发表于 2020-4-23 19:33
槽点太多了,计算方式太野蛮了

这种复合物优化直接用B3LYP等于往枪口上撞,完全是自寻死路。仔细看下文 ...

好的,谢谢老师了

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

4#
发表于 Post on 2020-4-23 19:33:29 | 只看该作者 Only view this author
槽点太多了,计算方式太野蛮了

这种复合物优化直接用B3LYP等于往枪口上撞,完全是自寻死路。仔细看下文补充弱相互作用最基本的常识
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413
DFT-D色散校正的使用
http://sobereva.com/210
乱谈DFT-D
http://sobereva.com/83
大体系弱相互作用计算的解决之道
http://sobereva.com/214

这么巨大体系在B3LYP/6-31G*下连极好的服务器都算不动,还指望Windows版Gaussian能算得动?这属于没有Gaussian基本计算常识。非要直接用这么大模型做优化,必须让molclus调用MOPAC或者xtb做。(调用Gaussian做考虑色散校正的半经验也可以,但效率低)

当前traj.xyz里就一个结构,拿molclus做毫无意义。如果你是想考察小分子在DNA什么地方什么方式结合比较稳定,起码得产生20个构型,而且必须在genmer里恰当设置参数,让小分子只出现在DNA的靠中间区域,避免有限长度的DNA模型的边界效应导致虚假的结果(你当前的这一帧里小分子太靠边了)。

我不晓得你的体系的细节,自旋多重度为2有极大嫌疑不合理。而且你的体系里DNA的骨架上的磷酸基都是质子解离状态,整体带巨量负电荷,总电荷怎么可能是0?而且整体带显著电荷的体系优化时必须带着隐式溶剂模型。

你的小分子上带的真的是硅元素么?

没事不要在opt里写cartesian,通常只会阻碍收敛


问题实在太多...感觉你的基础太欠缺了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

164

帖子

0

威望

1991

eV
积分
2155

Level 5 (御坂)

3#
发表于 Post on 2020-4-23 14:09:40 | 只看该作者 Only view this author
看看高斯输出文件有无报错

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

2#
发表于 Post on 2020-4-23 13:52:17 | 只看该作者 Only view this author
700多个原子用DFT优化还是用Gaussian你确定算的动么?
这里没有理由用opt=Cartesian,用了只会降低算法的效率。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-25 23:49 , Processed in 0.192425 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list