| 查看: 3008 | 回复: 9 | |||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||
beefly专家顾问 (职业作家)
地沟油冶炼专家
|
[交流]
Gaussian用冗余内坐标优化大分子结构时的问题已有7人参与
|
||||||||
|
出错信息： GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Eigenvalue 1920 is 9.98D-07 should be greater than 0.000001 Eigenvector: A602 A605 A526 A529 A590 1 -0.23218 0.21568 -0.17301 0.17222 -0.16746 A684 A687 D849 A603 D896 1 0.13759 -0.13443 -0.13232 0.11672 0.11547 NTrRot= -1 NTRed= 3567 NAtoms= 642 NSkip= 1647 IsLin=F Error in internal coordinate system. 问题分析： Gaussian默认使用冗余内坐标做结构优化，计算中检查使用的冗余坐标能否构成3N-L（L=6或5）个独立内坐标空间，也就是检查Wilson G矩阵是否有3N-L个非零特征值。由于产生冗余内坐标的算法本身存在缺陷，随着分子尺寸不断增加（通常几百个原子以上），有些重要的弱化学键（如氢键）会被丢掉，导致G矩阵的非零特征值个数小于3N-L，就会遇到上面的出错信息。在Gaussian09.d中，由于重新定义了冗余内坐标中的成键参数，丢失的弱键可能更多，因此程序报错的机会就更多一些。 解决方法： 1，自己写个小程序，给定原子类型后，找出键长在一定范围的所有的弱键。如果不包含在Gaussian产生的冗余内坐标中，就把它用opt=ModRedundant加到输入文件里。 2，更好的办法是改用直角坐标做优化，不仅不会遇到上述问题，还可以避免耗时的G矩阵对角化过程。 |
回复此楼» 收录本帖的淘帖专辑推荐
 我学习计算的一些帖子 | 高斯细节 | 经验 | Gaussian |
» 本帖已获得的红花（最新10朵）
枪下游魂» 猜你喜欢
超导转变温度NaN
已经有4人回复
-
华南师范大学先进光电子研究院电子纸团队诚聘2026级博士研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有134人回复
-
QE利用声子计算的一系列lambda和Tc如何选
已经有8人回复
-
理论计算合作
已经有2人回复
-
小木虫的论文辅导靠谱吗？有没有用过的同学？
已经有1人回复
-
QE计算声子ph.out的Γ点出现虚频（-30cm-1）vasp计算没有
已经有0人回复
-
40-50万/年，中国散裂中子源诚聘计算模拟方向博士后
已经有81人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
书籍求助：汽车市场营销理论与实务（电子版）——章小平
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
高斯大侠留步，看看这个高斯输入文件有何问题？
已经有8人回复- 用Gaussian算有机共轭大分子时，怎样构建分子的共轭键？
已经有14人回复
- 固定键长优化问题
已经有7人回复
- 高斯09寻找过渡态，急需帮助，十分感谢
已经有6人回复
- 【求助】TDDFT计算中基组的问题
已经有4人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +457 - QC强帖: 12
- 应助: 86 (初中生)
- 金币: 12735.4
- 散金: 27307
- 红花: 300
- 帖子: 4663
- 在线: 925.3小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
10楼2017-09-29 12:30:35
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
2楼2014-11-19 06:51:34
agent99
禁虫 (正式写手)
物化程序猿
- 应助: 75 (初中生)
- 贵宾: 0.001
- 金币: 3630.7
- 散金: 70
- 红花: 8
- 沙发: 1
- 帖子: 392
- 在线: 177.5小时
- 虫号: 1247746
- 注册: 2011-03-28
- 性别: GG
- 专业: 理论和计算化学
3楼2014-11-19 07:16:47
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
4楼2014-11-19 08:15:07