| 查看: 2435 | 回复: 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朵）
» 猜你喜欢
华算科技技术工程师招聘（生物模拟方向）
已经有0人回复
-
求助，stern层的电势降如何计算？
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有192人回复
-
祈福国自然基金通过！
已经有12人回复
-
云母片的计算构型可以从哪里获取？
已经有0人回复
-
Biomaterials期刊返稿时间求助
已经有11人回复
-
深圳湾实验室杨奕课题组招聘博士后、助理研究员及软件工程师
已经有18人回复
-
山东省基金：作为项目负责人承担面上项目数量累计不得超过2项
已经有29人回复
-
再次祈福国自然基金高中！
已经有20人回复
-
fluent带有化学反应的模拟
已经有0人回复
-
求《物理化学学习进阶与考研攻略》pdf文件
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯大侠留步，看看这个高斯输入文件有何问题？
已经有8人回复
- 用Gaussian算有机共轭大分子时，怎样构建分子的共轭键？
已经有14人回复
- 关于Gaussian Scan和opt=modredundant的问题，求大神帮忙
已经有6人回复
- 固定键长优化问题
已经有7人回复
- 量子化学计算中帮助几何优化收敛的常用方法
已经有38人回复
- 高斯09寻找过渡态，急需帮助，十分感谢
已经有6人回复
- 【求助】TDDFT计算中基组的问题
已经有4人回复
- 【coolrainbow个人文集】Gaussian机器人-真的存在黑箱吗-理论FAQ
已经有56人回复
- 【sobereva个人文集】转载一篇原创精品文章：过渡态、反应路径的计算方法及相关问题
已经有42人回复
枪下游魂
木虫 (著名写手)
- 应助: 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
5楼2014-11-19 15:41:01
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +453 - QC强帖: 12
- 应助: 85 (初中生)
- 金币: 10868.3
- 散金: 27307
- 红花: 296
- 帖子: 4588
- 在线: 916.8小时
- 虫号: 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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +453 - QC强帖: 12
- 应助: 85 (初中生)
- 金币: 10868.3
- 散金: 27307
- 红花: 296
- 帖子: 4588
- 在线: 916.8小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
10楼2017-09-29 12:30:35
10