| 查看: 3831 | 回复: 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朵）
» 猜你喜欢
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有193人回复
-
美国杰克逊州立大学周剑歌教授的化学计算与AI实验室诚聘博士后研究员
已经有4人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有3人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Tc数据拟合
已经有1人回复
-
韩国釜山大学储能电池与功能材料实验室诚招博后（电池/电化学/器件制备/MOF合成）
已经有187人回复
-
祈福！
已经有79人回复
-
散金求好运
已经有96人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯大侠留步，看看这个高斯输入文件有何问题？
已经有8人回复
- 用Gaussian算有机共轭大分子时，怎样构建分子的共轭键？
已经有14人回复
- 固定键长优化问题
已经有7人回复
- 高斯09寻找过渡态，急需帮助，十分感谢
已经有6人回复
- 【求助】TDDFT计算中基组的问题
已经有4人回复
枪下游魂
木虫 (著名写手)
- 应助: 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
abdoman
木虫 (著名写手)
- 应助: 56 (初中生)
- 金币: 5369.5
- 散金: 918
- 红花: 3
- 沙发: 1
- 帖子: 1321
- 在线: 451.5小时
- 虫号: 301521
- 注册: 2006-12-01
- 性别: GG
- 专业: 理论和计算化学
5楼2014-11-19 15:41:01
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13331
- 散金: 27307
- 红花: 301
- 帖子: 4683
- 在线: 928.6小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
6楼2014-11-20 01:23:53
7楼2015-07-08 09:14:22
送红花一朵 |
本帖内容被屏蔽 |
8楼2016-02-02 12:30:45
9楼2017-09-29 11:41:06
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13331
- 散金: 27307
- 红花: 301
- 帖子: 4683
- 在线: 928.6小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
10楼2017-09-29 12:30:35
