高斯无法优化环C36

archer

初始结构经过xtb预优化，计算单点可以正常结束。加上opt关键词，无后续输出，卡在读取文件输入文件：
....
7//1,2,3,16;
1/14=-1,18=20,19=15,26=4/3(-5);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
-------
C18-C18
-------
Symbolic Z-matrix:
Charge =  0 Multiplicity = 1
C                     7.33219  -0.02159   0.00077

.....

Stardust0831

确实很奇怪，我这边可以重复出这一问题。精简关键词、换基组和泛函都没用。
甚至freq反而可以正常跑。

archer

Stardust0831 发表于 2025-5-8 16:02
确实很奇怪，我这边可以重复出这一问题。精简关键词、换基组和泛函都没用。
甚至freq反而可以正常跑。

不加opt，单算频率，可以跑，请问你用的也是G09吗。

Stardust0831

我用http://sobereva.com/610中提到的选项18生成了随机位移结构破坏了对称性，可以正常跑opt了。

Stardust0831

archer 发表于 2025-5-8 16:10
不加opt，单算频率，可以跑，请问你用的也是G09吗。

我是linux的g16的c版本

Stardust0831

我找到问题的最优解决了，使用opt=cartesian可以正常优化，估计是内坐标与笛卡尔坐标相互转换这部分的代码细节有问题

不好意思刚刚想修改原来的回帖的，不小心把这条消息当新的回帖发出来了，有点丑陋hhh

archer

Stardust0831 发表于 2025-5-8 16:18
我找到问题的最优解决了，使用opt=cartesian可以正常优化，估计是内坐标与笛卡尔坐标相互转换这部分的代码 ...

这种算是Gaussian的bug吗

wal

我是g16 B01 用主楼那个直接(除了titel行少个空行之外)就能正常跑

是某些版本特殊bug吗

Stardust0831

wal 发表于 2025-5-8 17:29
我是g16 B01 用主楼那个直接(除了titel行少个空行之外)就能正常跑

是某些版本特殊bug吗

好神奇。我和stw这边用高斯跑都发现会卡在读取坐标这一步不动。stw是g16的A.03；我是g16的C.01

Stardust0831

archer 发表于 2025-5-8 16:49
这种算是Gaussian的bug吗

这是我群友stw的一个猜测：

想到一个可能性，他这三个原子的夹角是169.99，所以gaussian构建二面角的时候，可能认为这三个原子在一条直线上，然后再以这三个原子为基扩展搜索二面角，结果发现下一个原子也“在一条直线上”所以就不断循环，一直转圈

archer

Stardust0831 发表于 2025-5-8 17:49
这是我群友stw的一个猜测：

可能是的，我手动修改几个原子坐标，就能跑起来了。

fatpig

有两种解决方法：
1. opt=cartesian
2. opt=z-matrix
当然，分子坐标需要写成对应模式...