计算化学公社

标题: 想对比分子同位素替换的影响，替后换电子能量一致，可我计算的吉布斯自由能也一致 [打印本页]

作者
Author: cdw32156 时间: 2025-2-11 12:53
标题: 想对比分子同位素替换的影响，替后换电子能量一致，可我计算的吉布斯自由能也一致
想问一下， (, 下载次数 Times of downloads: 4) 我用这个分子，两个氘原子在不同位置进行替换。我只对未替换的原分子进行了优化，后面直接修改的原子质量，然后分别进行频率分析，计算得到的吉布斯自由能都是一样的。不懂是为什么，是因为替换的氘原子数量一致吗

作者
Author: dnlx 时间: 2025-2-11 22:27
1. 只替换一个氘、和替换两个氘，得到的自由能是否相同？如果相同，说明输入文件有问题。
2.替换两个氘的同位素异构体的频率是否有区别？如果没区别，说明输入文件有问题。如果有区别，说明吉布斯能的差异被抵消了。另外你是不是用quasiRRHO计算的振动熵。

作者
Author: sobereva 时间: 2025-2-12 05:34
几乎一定是设置同位素的做法不当、犯了低级错误
拿振动分析的具体输入、输出文件说事

作者
Author: cdw32156 时间: 2025-2-15 21:23

sobereva 发表于 2025-2-12 05:34
几乎一定是设置同位素的做法不当、犯了低级错误
拿振动分析的具体输入、输出文件说事

这是振动分析的具体输入、输出文件。找不到是什么原因

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-2-15 21:47

cdw32156 发表于 2025-2-15 21:23
这是振动分析的具体输入、输出文件。找不到是什么原因

20核才给4GB内存太离谱了，加弥散、pop=full和iop(9/40=4)的意义也不明。而且为啥几何优化结果只取结构而不从chk里读取波函数呢？

作者
Author: cdw32156 时间: 2025-2-15 22:02

Uus/pMeC6H4-/キ发表于 2025-2-15 21:47
20核才给4GB内存太离谱了，加弥散、pop=full和iop(9/40=4)的意义也不明。而且为啥几何优化结果只取结构而 ...

您好，给的内存小是因为这是一个单分子。加弥散是因为后续计算需要极化率。但是这些关键词对我的计算来说，除了所需时间变长，对结果本身是不是没有什么影响。我的问题是：同位素在不同位置替后换电子能量一致，可我计算的吉布斯自由能也一致，这是什么原因

作者
Author: sobereva 时间: 2025-2-16 00:56

cdw32156 发表于 2025-2-15 21:23
这是振动分析的具体输入、输出文件。找不到是什么原因

热力学校正量并不一样
CRdui.log
Zero-point vibrational energy    711227.0 (Joules/Mol)
                              169.98733 (Kcal/Mol)

CRlin.log
Zero-point vibrational energy    711222.3 (Joules/Mol)
                              169.98622 (Kcal/Mol)

保留小数位数太少时看不出差异

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