关于使用ORCA 6.0.0的GOAT模块进行多核计算时只有1个核心在运作的问题

QJW · 发表于 Post on yesterday 15:33

wzkchem5 发表于 2026-3-8 14:53
GOAT（包括GOAT-ENTROPY等所有衍生关键词）的输出已经是去重过后的结果了。后续用更高理论级别优化以后， ...

老师,gaot输出了三个看起来一样的结果

QJW · 发表于 Post on yesterday 15:34

wzkchem5 发表于 2026-3-8 14:53
GOAT（包括GOAT-ENTROPY等所有衍生关键词）的输出已经是去重过后的结果了。后续用更高理论级别优化以后， ...

请问这个该怎么整啊

wzkchem5 · 发表于 Post on yesterday 21:31

QJW 发表于 2026-3-8 15:33
老师,gaot输出了三个看起来一样的结果

什么叫看起来一样？确定不是对映异构体或非对映异构体？

Daniel_Arndt · 发表于 Post on 12 hour ago

wzkchem5 发表于 2026-3-8 21:31
什么叫看起来一样？确定不是对映异构体或非对映异构体？

可能类似我之前遇到的情况。我研究的体系里有三苯基膦配体，有时候两个结构的不同仅仅是三苯基膦转了约120度的差别，然后GOAT就会认为两个结构不一样，但只要用可视化软件看一下，就知道是一样的东西了。

wzkchem5 · 发表于 Post on 11 hour ago

Daniel_Arndt 发表于 2026-3-9 10:16
可能类似我之前遇到的情况。我研究的体系里有三苯基膦配体，有时候两个结构的不同仅仅是三苯基膦转了约12 ...

就是两个分子可以重合是吗？VMD等软件算出来的RMSD接近0？
如果RMSD很小，但也离0有一定距离（比如RMSD 0.1埃），那就是GOAT的结构优化标准太松了，或GOAT的RMSD标准太严了。可以用tightopt优化收敛后，再参考我之前给的链接去重

Daniel_Arndt · 发表于 Post on 7 hour ago

wzkchem5 发表于 2026-3-9 11:45
就是两个分子可以重合是吗？VMD等软件算出来的RMSD接近0？
如果RMSD很小，但也离0有一定距离（比如RMSD ...

实际情况是这样的，我用GaussView查看的时候，如果不看“Labels”（原子序号），两个构象没区别；如果看“Labels”（原子序号），才能看出来二者之间的不同就是三苯基膦绕着金属跟膦之间的键转了个120度（我当时算的是三苯基膦跟一个过渡金属试剂配位）。这两个构象如果直接用VMD算RMSD的话，肯定不接近0；但如果把其中一个构象中三个苯环的原子顺序调一下再算RMSD的话，结果就很接近0。

至于可能的解决办法，一个是让用户在GOAT计算的输入文件中设置fragment，要用到“%frag”、“ConstrainFragments”、“ConnectFragments”，把三苯基膦设置成一个不允许其变动的fragment，或者是仅允许氢原子变动的fragment（“optimizeHydrogens true”）；另一个办法就对开发者而言很麻烦了，更改计算RMSD的算法，仿照 https://doi.org/10.1186/s13321-020-00455-2 ，可能会导致很多代码要变动。

wzkchem5 · 发表于 Post on 3 hour ago

Daniel_Arndt 发表于 2026-3-9 14:54
实际情况是这样的，我用GaussView查看的时候，如果不看“Labels”（原子序号），两个构象没区别；如果看 ...

GOAT的RMSD算法是考虑了原子序号的变化的。考虑了原子序号变化后RMSD接近0，但不考虑原子序号变化时RMSD接近0的两个构象，叫做rotamer。参见https://www.faccts.de/docs/orca/ ... the-delta-s-rm-conf
所以如果两个rotamer被判断成了两个不同的conformer，应该还是RMSD阈值太严，或结构优化收敛得不够严

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] 关于使用ORCA 6.0.0的GOAT模块进行多核计算时只有1个核心在运作的问题

浏览过的版块