计算化学公社

标题: C[n]环结构优化时达到C[24]以后卡在L101.exe不动了？ [打印本页]

作者
Author: liyunlu 时间: 2021-6-2 17:32
标题: C[n]环结构优化时达到C[24]以后卡在L101.exe不动了？
本帖最后由 liyunlu 于 2021-6-2 17:53 编辑

最近尝试重复文献C[n]环结构优化( Aromaticity of Even-Number Cyclo[n]carbons (n=6-100), J. Phys. Chem. A 2020, 124, 51, 10849–10855. )。在n=6-22时，结构优化和频率计算都能正常进行。但n=24以后的结构优化和频率计算都不能进行，都卡在L101.exe不懂，不论是windows的G09D01，还是Linux的G16A03都是这样。但只进行单点计算可以正常计算。
这里贴出计算出入文件：
%chk=C26-opt.chk
# opt freq wb97xd/def2tzvp out=wfn
Title Card Required
0 1
C    -2.590324000       4.650570000       -0.000305000
C    -1.468481000       5.117110000       -0.000364000
C       2.354235000       4.773576000       -0.000197000
C       3.375091000       4.114686000       -0.000211000
C       4.301706000       3.131710000       -0.000100000
C       4.899518000       2.073916000       -0.000074000
C       5.020394000       -1.762606000       0.000268000
C       5.301920000       -0.441372000       0.000116000
C       4.490326000       -2.855883000       0.000242000
C       2.649878000       -4.616857000       0.000291000
C       1.402577000       -5.135657000       0.000314000
C       0.201633000       -5.319988000       0.000262000
C    -1.143820000       -5.198997000       0.000242000
C    -2.292671000       -4.803521000       0.000182000
C    -3.427568000       -4.070812000       0.000075000
C    -4.261073000       -3.186747000       0.000107000
C    -4.925725000       -2.010694000       -0.000022000
C    -5.042705000       1.697935000       -0.000127000
C    -5.295805000       0.509575000       -0.000154000
C       3.627295000       -3.895160000       0.000330000
C    -3.677105000       3.848178000       -0.000235000
C    -4.453371000       2.913503000       -0.000201000
C       1.076897000       5.213193000       -0.000256000
C    -0.133223000       5.321956000       -0.000283000
C       5.263650000       0.773035000       0.000091000
C    -5.253245000       -0.840646000       0.000009000
C26-opt.wfn
计算时卡着不动：
输出文件截图如下：

作者
Author: liyunlu 时间: 2021-6-2 17:40
只计算opt或者只计算freq都卡在L101不动，计算单点可以正常结束

作者
Author: zjxitcc 时间: 2021-6-2 17:44
用#p，而非#，可以看到更多输出、计算过程，当前输出太少。没必要尝试Win版Gaussian，直接在linux下用g16算即可，接下来也无需再贴Win下的任何内容（因为没用）。你还得交待你的linux下所用的核数和内存，你不可能真的只用1个核在算啊，那当然是动不了了，体系和基组都不小。

作者
Author: liyunlu 时间: 2021-6-2 17:51

zjxitcc 发表于 2021-6-2 17:44
用#p，而非#，可以看到更多输出、计算过程，当前输出太少。没必要尝试Win版Gaussian，直接在linux下用g16算 ...

做过的尝试如下：
（1）换#p和#，
（2）换版本G09D01(Windows 4c,1.5G), G16a03(36c, 180G),
（3）换泛函m062x-wb97xd.
（4）单独计算opt，单独计算freq。
上述尝试都只停在L101.exe不动。但只跑单点，上述都可以~还请老师们指点

作者
Author: biogon 时间: 2021-6-2 19:51
你先用#p看到底是停在什么步骤上了

作者
Author: zjxitcc 时间: 2021-6-2 20:08

liyunlu 发表于 2021-6-2 17:51
做过的尝试如下：
（1）换#p和#，
（2）换版本G09D01(Windows 4c,1.5G), G16a03(36c, 180G),

不用列很多种情况。直接给出#p下，linux下某个“卡住”的文件末尾截图。

作者
Author: sobereva 时间: 2021-6-3 07:10
没必要自己优化。我在这篇博文介绍的文章的补充材料里都已经提供好了<=C30的优化完的坐标了
揭示各种新奇的碳环体系的振动特征
http://sobereva.com/578（http://bbs.keinsci.com/thread-20919-1-1.html）

这和Gaussian对于某些高阶点群处理上的bug有关，之前我在上面这篇文章里提到的研究中算大环的时候就已经遇到了

作者
Author: liyunlu 时间: 2021-6-3 09:35
感谢各位老师的解答，这个的确与sob老师说的，可能是高阶对称性的问题，使用opt=z-matrix可以解决这个问题，就是让优化时使用内坐标进行计算（输入文件的坐标时内坐标还是笛卡尔坐标）

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