wzkchem5 发表于 2022-4-13 19:05 是的 确实容易忽视犯错 |
scf 发表于 2022-4-13 18:56 笔误，已纠正 |
granvia 发表于 2022-4-13 12:02 对，只不过因为咱们论坛上经常有人从老帖学习知识，我加那一句话主要是为了避免误导后人 |
wzkchem5 发表于 2022-4-13 18:46 有道理，确实应根据自旋量子数严格考虑。但具体到楼主所关心的计算，ta做的确是单参考态的计算吧 |
granvia 发表于 2022-4-13 11:41 因为自旋多重度是自旋量子数的2倍加1，只有对于单参考态波函数，它才等于alpha电子数减去beta电子数加1，对于多参考态波函数未必。所以不能认为氧原子的两个单电子要么自旋都向上，要么自旋都向下，有可能虽然一个自旋向上一个自旋向下，但是仍然是三重态，只不过这个态没办法用单参考态方法描述，所以那些只做过单参考态计算的人会下意识地认为不存在这样的情况。类似地，两个孤立的氧原子组成的总体系，即便自旋向上的电子数等于向下的电子数，也不必然代表是单重态。 |
| 参与人数Participants 2 | eV +10 | 收起 理由Reason |
|---|---|---|
|
| + 5 | 谢谢 |
|
| + 5 | 我很赞同 |
本帖最后由 granvia 于 2022-4-13 19:03 编辑 wzkchem5 发表于 2022-4-13 14:21 可是两个相距足够远的氧原子，可认为已是孤立原子了，每个孤立氧原子的基态自旋多重度都应该是3吧。因此总体系基态总的多重度怎可能是3呢？ |
granvia 发表于 2022-4-13 05:50 准确来说，仅当用单参考态方法直接算基态时是这样的。用多参考态方法，或者用单参考的SF-TDDFT以高自旋态为参考态算低自旋态时，自旋多重度是可以取3的 |
changxiaoya 发表于 2022-4-13 03:24 能量只有一半应该是能量读错了，或者输入文件写错了。 除了氢分子之类的电子数极其少的体系以外，自旋多重度即使设错了，也不可能导致绝对能量差一半这么多。 |
本帖最后由 granvia 于 2022-4-13 12:52 编辑 changxiaoya 发表于 2022-4-13 10:24 基态的自旋多重度只可能是1或5 |
wzkchem5 发表于 2022-4-12 20:23 我检查一下，感谢！  |
granvia 发表于 2022-4-12 20:29 请问一下，距离远的时候，两个单独的三重态氧原子的能量，总的自旋多重态应该是1 还是 3 呢？ 我尝试用片段组合的方法，氧原子分别为三重态，总的多重度设置为1和3 ， 能量都只有一半（O2为150 hart， 片段组合计算出来只有75 hart）且键长改变，键能几乎不变。 |
本帖最后由 changxiaoya 于 2022-4-13 10:27 编辑 scf 发表于 2022-4-13 05:48 没有用scan，是拆分成单点来计算的。我试一下限制性计算，谢谢！ 另外，感谢分享这个网站。 |
|
本帖最后由 scf 于 2022-4-13 06:02 编辑 用的scan吗？印象里高斯的scan用上一个位置的scf的波函数作为初始猜测，有的时候不如单独算单点的默认初猜得到的势能面看上去平滑。另外unrestricted和guess=mix和stable=opt放在一起算单重态算的能量, 有时候不如 restricted + 没有guess=mix + stable=opt更低。O2基态是三重态，觉得三重态的键能应该比单重态更大 https://www.ucl.ac.uk/~ucapsy0/diatomics.html |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-14 04:50 , Processed in 0.183495 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites