计算化学公社

标题: 用gic做扫描6步后显示New curvilinear step failed, RedCar failed. [打印本页]

作者
Author: 隔壁老张sir 时间: 2023-3-29 19:03
标题: 用gic做扫描6步后显示New curvilinear step failed, RedCar failed.
本帖最后由隔壁老张sir 于 2023-3-29 19:11 编辑

输入文件：
%nprocshared=12
%mem=24GB
%chk=relax_bentTshape_cas12_dummyatom_dke_GIC.chk
#p opt=addgic casscf(12,12) maug-cc-pvtz guess(mix,always) int=ultrafine Guess=Permute SCF=maxcyc=500 nosymm

Title Card Required

0 1
N             -1    0.47544795 0.00000000 -0.27450000
N             -1 -0.47544795 0.00000000 0.27450000
N             0    0.00000000 0.00000000  -10.40000000
N             0    0.00000000 0.00000000 -9.60000000

XC1=XCntr(1,2)
YC1=YCntr(1,2)
ZC1=ZCntr(1,2)
XC2=XCntr(3,4)
YC2=YCntr(3,4)
ZC2=ZCntr(3,4)
F1F2(Frozen)=sqrt[(XC1-XC2)^2+(YC1-YC2)^2+(ZC1-ZC2)^2]*0.529177
D34(StepSize=0.05,NSteps=40)=R(3,4)

“casscf(12,12)/maug-cc-pvtz”是参考的一篇文献，“int=ultrafine Guess=Permute SCF=maxcyc=500”是为了收敛加的。

报错：
NOTE: GIC-type coordinates are in arbitrary units.
Lowest energy point so far.  Saving SCF results.
Iteration  1 RMS(Cart)=  0.01434970 RMS(Int)=  0.00802064
Iteration  2 RMS(Cart)=  0.00000225 RMS(Int)=  0.00802064
New curvilinear step failed, DQL= 3.39D-03 SP=-1.09D-07.
RedCar failed.
Error termination via Lnk1e in /home/zhangsir/g16/l103.exe at Wed Mar 29 17:00:31 2023.
Job cpu time:    1 days  8 hours 49 minutes 29.2 seconds.
Elapsed time:    0 days  2 hours 44 minutes 58.7 seconds.
File lengths (MBytes):  RWF= 1353 Int=    0 D2E=    0 Chk=    4 Scr=    1

请问各位老师，这是什么原因呢？有什么解决办法吗？

作者
Author: Deepast 时间: 2023-3-30 02:39
gic有时候就会这样，程序不知道该怎么优化
我不知道有什么好的解决方法，有时候甚至交换原子顺序就能跑
先续算看看能不能跑，不行就手动调下，再不行换别的扫描方法

作者
Author: Deepast 时间: 2023-3-30 02:49
最好把广义坐标也简化下，我不确定程序崩和这个有没有关系
比如冻结的话没必要写这么多
开方和单位转换那部分可以删

作者
Author: 隔壁老张sir 时间: 2023-3-30 08:34

Deepast 发表于 2023-3-30 02:39
gic有时候就会这样，程序不知道该怎么优化
我不知道有什么好的解决方法，有时候甚至交换原子顺序就能跑
...

您说的“先续算看看能不能跑，不行就手动调下”可以麻烦您说的具体一点吗？
我刚接触量化不久，是个小白，不太知道具体是要怎么做

作者
Author: 隔壁老张sir 时间: 2023-3-30 08:38

Deepast 发表于 2023-3-30 02:49
最好把广义坐标也简化下，我不确定程序崩和这个有没有关系
比如冻结的话没必要写这么多
开方和单位转换那 ...

如果不固定两个原子，那个分子就会旋转

作者
Author: sobereva 时间: 2023-3-30 09:27
不知道你想算什么，这关键词完全莫名其妙，看不出这计算有任何实际意义

GIC的算法并不稳健，尝试不同定义方式，或者根据实际情况考虑不用GIC的可能

作者
Author: 隔壁老张sir 时间: 2023-3-30 10:05

sobereva 发表于 2023-3-30 09:27
不知道你想算什么，这关键词完全莫名其妙，看不出这计算有任何实际意义

GIC的算法并不稳健，尝试不同定 ...

sob老师好。我是参考《Global Ab Initio Ground-State Potential Energy Surface of N4》这篇文献，尝试构建N4体系的势能面。我这里是选了文章中的一种N2+N2的构型，按文章中给的键长变化参数，尝试计算一下势能变化。
“casscf(12,12)/maug-cc-pvtz”是该文献中用的，用它的实际意义我也不太清楚；“int=ultrafine Guess=Permute SCF=maxcyc=500”是我为了收敛加的；过程中有共价键断裂，就又加了“guess(mix,always)”；为了防止扫描过程分子跳变，加的“nosymm”。

sob老师说的“尝试不同定义方式”是指什么呀？
我一开始也想着不用GIC，但是没有找到其他办法使扫描过程中两氮分子的中心距离保持不变，sob老师有什么好的建议吗？

作者
Author: Deepast 时间: 2023-3-30 11:44

隔壁老张sir 发表于 2023-3-30 10:05
sob老师好。我是参考《Global Ab Initio Ground-State Potential Energy Surface of N4》这篇文献，尝试 ...

准经典轨线的研究是吧，这个研究不是靠这种扫描做的
做这个的比较少，建议多问问导师和师兄师姐

作者
Author: 隔壁老张sir 时间: 2023-3-30 11:53

Deepast 发表于 2023-3-30 11:44
准经典轨线的研究是吧，这个研究不是靠这种扫描做的
做这个的比较少，建议多问问导师和师兄师姐

唉我是组里第一个做这个的人，导师也不懂这个方向（方向是大老板硬给我的）

作者
Author: Deepast 时间: 2023-3-30 13:25

隔壁老张sir 发表于 2023-3-30 11:53
唉我是组里第一个做这个的人，导师也不懂这个方向（方向是大老板硬给我的）

这个是要算一系列单点能，然后构建全维势能面
结构可以通过编程生成
算了单点后可以根据一些算法拟合势能面，这个一般也要自己写
你想练手的话建议开始把计算精度调低点
我懂的也不多，只能说到这

作者
Author: 隔壁老张sir 时间: 2023-3-31 10:07

Deepast 发表于 2023-3-30 13:25
这个是要算一系列单点能，然后构建全维势能面
结构可以通过编程生成
算了单点后可以根据一些算法拟合势 ...

嗯嗯十分感谢

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