计算化学公社

标题: 过渡态使用直角坐标出错 [打印本页]

作者
Author: xujian 时间: 2016-9-28 10:28
标题: 过渡态使用直角坐标出错
刚开始做量化计算的时候，找过渡态使用的内坐标，一直也没在意，有一次保存成直角坐标，刚提交任务就出错，请问老师，用直角坐标找过渡态是不是要添加什么关键词啊，烦请格外老师不吝赐教，谢谢！

作者
Author: librakitty 时间: 2016-9-28 11:58
这个和你的初始坐标输入没有关系，guassian 优化默认就用的内坐标，倒是你想用直角坐标优化得输入opt=cartesian。你最好检查一下报错信息再考虑到底是啥问题。

作者
Author: xujian 时间: 2016-9-28 12:21

librakitty 发表于 2016-9-28 11:58
这个和你的初始坐标输入没有关系，guassian 优化默认就用的内坐标，倒是你想用直角坐标优化得输入opt=carte ...

好的。谢谢

作者
Author: xujian 时间: 2016-9-28 12:31

librakitty 发表于 2016-9-28 11:58
这个和你的初始坐标输入没有关系，guassian 优化默认就用的内坐标，倒是你想用直角坐标优化得输入opt=carte ...

我关键词是opt=(calcfc,ts,ef,maxcycle=400,cartesian)，用的直角坐标，然后直接报错，Error: segmentation violation

作者
Author: liyuanhe211 时间: 2016-9-28 12:45

xujian 发表于 2016-9-28 12:31
我关键词是opt=(calcfc,ts,ef,maxcycle=400,cartesian)，用的直角坐标，然后直接报错，Error: segmentati ...

压缩后上传输出，segmentation violation毛线信息都没有。
还用eigenvalue following，要么是根本不明白它什么意思，要么是用了啥毫无意义的水平来优化。
maxcycle垃圾关键词
不用cartesian不出错就别用

作者
Author: xujian 时间: 2016-9-28 17:06

liyuanhe211 发表于 2016-9-28 12:45
压缩后上传输出，segmentation violation毛线信息都没有。
还用eigenvalue following，要么是根本不明白 ...

谢谢老师，附件如下，麻烦您查看解答

(, 下载次数 Times of downloads: 4)

作者
Author: xujian 时间: 2016-9-28 17:52

liyuanhe211 发表于 2016-9-28 12:45
压缩后上传输出，segmentation violation毛线信息都没有。
还用eigenvalue following，要么是根本不明白 ...

高斯说明书上说。对于所有不能使用解析梯度的计算方法，几何优化计算的默认算法是本征值跟踪算(Opt=EF)，EF是什么意思我确实不理解，不知道老师能否解答一下

作者
Author: sobereva 时间: 2016-9-28 19:08

xujian 发表于 2016-9-28 17:52
高斯说明书上说。对于所有不能使用解析梯度的计算方法，几何优化计算的默认算法是本征值跟踪算(Opt=EF)， ...

和默认的优化方法大同小异
没必要手动去指定之

作者
Author: liyuanhe211 时间: 2016-9-28 20:56

xujian 发表于 2016-9-28 17:06
谢谢老师，附件如下，麻烦您查看解答

你这就是本来用不着EF非要用，然后瞎用EF还不会用。不要写任何自己不理解的关键词

(, 下载次数 Times of downloads: 65)

(, 下载次数 Times of downloads: 51)

作者
Author: xujian 时间: 2016-9-29 09:33

liyuanhe211 发表于 2016-9-28 20:56
你这就是本来用不着EF非要用，然后瞎用EF还不会用。不要写任何自己不理解的关键词

谢谢老师，我把ef去掉就可以运算了，还有请教一下，说明书中对MaxCycle=N的解释是：指定优化计算最大迭代次数为N。默认是（默认的计算程序）所用冗余内坐标数目的两倍，最大为20，或是（其它计算程序）做优化计算的变量数目的两倍。它说最大值是20，为什么我写400不会出错？

作者
Author: 绫波レイ 时间: 2016-9-29 09:37

xujian 发表于 2016-9-29 09:33
谢谢老师，我把ef去掉就可以运算了，还有请教一下，说明书中对MaxCycle=N的解释是：指定优化计算最大迭代 ...

最大迭代次数一般对于优化没什么帮助的，具体可以看http://sobereva.com/164

作者
Author: KiritsuguPapa 时间: 2016-9-29 11:52
做过渡态要学学此文：　　http://sobereva.com/44

作者
Author: xujian 时间: 2016-9-29 12:44

KiritsuguPapa 发表于 2016-9-29 11:52
做过渡态要学学此文：　　http://sobereva.com/44

谢谢，请问过渡态出现 Error termination request processed by link 9999重新保存之后依然是这个问题，是不是跟基态处理方法不同

作者
Author: KiritsuguPapa 时间: 2016-9-30 11:01

xujian 发表于 2016-9-29 12:44
谢谢，请问过渡态出现 Error termination request processed by link 9999重新保存之后依然是这个问题， ...

仔细看11楼发的那个地址的文章，优化不收敛的解决办法。肯定不是直接保存最后结构重跑就行的，这样和用maxcyc没有区别。

作者
Author: xujian 时间: 2016-9-30 11:23

KiritsuguPapa 发表于 2016-9-30 11:01
仔细看11楼发的那个地址的文章，优化不收敛的解决办法。肯定不是直接保存最后结构重跑就行的，这样和用ma ...

优化不收敛显示的不应该是l502么。如果跟maxcyc没有区别，那为什么在基态出现这种情况的时候这么处理一般都是可行的，如果您知道的话，是否可以详细解答一下，谢谢！

作者
Author: sobereva 时间: 2016-9-30 11:57

xujian 发表于 2016-9-30 11:23
优化不收敛显示的不应该是l502么。如果跟maxcyc没有区别，那为什么在基态出现这种情况的时候这么处理一般 ...

你都没搞清楚SCF和几何优化的基本过程。
SCF和opt都是迭代过程，opt的每一轮循环都要将SCF迭代至收敛。opt算是外循环，SCF算是内循环。
二者都有maxcyc控制迭代次数上限，502报错是SCF不收敛，9999报错是opt不收敛
SCF不收敛看
解决SCF不收敛问题的方法
http://sobereva.com/61
opt不收敛看
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

作者
Author: xujian 时间: 2016-9-30 12:18

sobereva 发表于 2016-9-30 11:57
你都没搞清楚SCF和几何优化的基本过程。
SCF和opt都是迭代过程，opt的每一轮循环都要将SCF迭代至收敛。o ...

非常感谢sob老师！

作者
Author: KiritsuguPapa 时间: 2016-9-30 12:44

xujian 发表于 2016-9-30 11:23
优化不收敛显示的不应该是l502么。如果跟maxcyc没有区别，那为什么在基态出现这种情况的时候这么处理一般 ...

502是SCF不收敛，和优化不收敛不是一个东西
没有“一般都是可行的”这一说，偶尔有那也是巧合，是刚好能优化完只是步数不够，和maxcyc一样的；如果最后振荡了，你再取最后结构继续用原方法优化，那也还是震荡。
一般取最后结构，调整方法（比如calcfc,gdiis)继续优化。
11楼的那个文末写了优化不收敛的处理方式和尝试顺序。

作者
Author: KiritsuguPapa 时间: 2016-9-30 12:44

xujian 发表于 2016-9-30 11:23
优化不收敛显示的不应该是l502么。如果跟maxcyc没有区别，那为什么在基态出现这种情况的时候这么处理一般 ...

那个文章真的是够详细了，你先看了再说

作者
Author: xujian 时间: 2016-9-30 17:01

KiritsuguPapa 发表于 2016-9-30 12:44
那个文章真的是够详细了，你先看了再说

好的，非常感谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)