C[n]环结构优化时达到C[24]以后卡在L101.exe不动了？

liyunlu

最近尝试重复文献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
计算时卡着不动：
输出文件截图如下：

liyunlu

只计算opt或者只计算freq都卡在L101不动，计算单点可以正常结束

zjxitcc

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

liyunlu

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不动。但只跑单点，上述都可以~还请老师们指点

biogon

你先用#p看到底是停在什么步骤上了

zjxitcc

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

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

sobereva

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

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

liyunlu

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