zjxitcc 发表于 2025-6-12 16:25 老师您好,我又想了一下,如果是开壳层单重态的话,我直接将原子的自旋布局乘以-1,这会导致电子自旋翻转,电子自旋翻转需要能量啊,这怎么解释呢? 感谢老师的耐心解答! |
zjxitcc 发表于 2025-6-12 20:57 好的好的!感谢老师 |
直接定义片段初猜就完了。仔细看: 谈谈片段组合波函数与自旋极化单重态 http://sobereva.com/82 |
本帖最后由 zjxitcc 于 2025-6-12 20:59 编辑 AuAuAu 发表于 2025-6-12 20:50 都适用(这个其实你自己花几秒思考就知道了)。足够。 |
zjxitcc 发表于 2025-6-12 16:25 好的,谢谢老师,之后会注意用词! 那请问老师(2)中的“但对于单重态体系,自旋布居的符号不是固定的,可以将所有原子的自旋布居乘以-1(因为乘-1前后,总和仍然等于0),无需再次计算,也无需换用计算软件,无需FlipSpin”是指开壳层单重态和闭壳层单重态均适用吗 以及(3)中检验波函数稳定性”#p ub3lyp/def2svp guess=mix stable=opt empiricaldispersion=gd3bj “是否就足够了 |
本帖最后由 zjxitcc 于 2025-6-13 08:11 编辑 (1)自旋密度是三维空间的函数,在三维空间中每一点处都有值。你所展示的数据是每个原子上有一个值,没有原子的地方没有值,这叫自旋布居,不叫自旋密度,也不可以简称为spin。注意:大量文献和程序输出搞错了这一基本概念,或者明知故犯,造成大量理论计算化学基础不扎实的科研人员用词错误。 (2)自旋布居是可以有正、有负、等于0的。如果体系是闭壳层单重态,每个原子的自旋布居都是0;如果体系是开壳层单重态,每个原子的自旋布居有正有负。如果全为正 或 全为负,说明此单重态计算是错误的(即用户某一步操作错了)。但对于单重态体系,自旋布居的符号不是固定的,可以将所有原子的自旋布居乘以-1(因为乘-1前后,总和仍然等于0),无需再次计算,也无需换用计算软件,无需FlipSpin,如果您看一下自旋布居的计算过程,便知这是合理操作。 但对于二重态等 不是单重态的计算,不可以将所有原子的自旋布居乘以-1。 (3)guess=mix仅仅是构造对称破缺初猜,不保证一定收敛到对称破缺UKS解。一定要记得检验波函数稳定性,确保收敛到稳定波函数。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-13 03:13 , Processed in 0.172493 second(s), 25 queries , Gzip On.