计算化学公社

标题: 高斯无法优化环C36 [打印本页]

作者
Author: archer 时间: 2025-5-8 15:20
标题: 高斯无法优化环C36
初始结构经过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

.....

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

作者
Author: archer 时间: 2025-5-8 16:10

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

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

作者
Author: Stardust0831 时间: 2025-5-8 16:12
我用http://sobereva.com/610中提到的选项18生成了随机位移结构破坏了对称性，可以正常跑opt了。

作者
Author: Stardust0831 时间: 2025-5-8 16:14

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

我是linux的g16的c版本

作者
Author: Stardust0831 时间: 2025-5-8 16:18
本帖最后由 Stardust0831 于 2025-5-8 16:21 编辑

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

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

作者
Author: archer 时间: 2025-5-8 16:49

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

这种算是Gaussian的bug吗

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

(, 下载次数 Times of downloads: 12)
是某些版本特殊bug吗

作者
Author: Stardust0831 时间: 2025-5-8 17:47

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

是某些版本特殊bug吗

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

作者
Author: Stardust0831 时间: 2025-5-8 17:49

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

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

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

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

作者
Author: archer 时间: 2025-5-8 18:31

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

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

作者
Author: fatpig 时间: 2025-5-9 09:16
有两种解决方法：
1. opt=cartesian
2. opt=z-matrix
当然，分子坐标需要写成对应模式...

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