Winnky 发表于 2024-10-20 19:09 你对自旋密度的理解大错特错 仔细看此文 谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算 http://sobereva.com/353(http://bbs.keinsci.com/thread-4422-1-1.html) |
sobereva 发表于 2017-1-2 12:06 卢老师,请问我可以理解成自旋密度+1吗? 因为我希望在做完开壳层计算之后确定选用的多重态是否合理 谢谢! |
wzkchem5 发表于 2024-3-23 00:57 感谢!看来是只能暴力枚举了 |
liaojianhan 发表于 2024-3-22 15:37 首先不难自动由总核电荷数确定自旋多重度的奇偶。 其次,如果需要确定自旋多重度的物种不是多得离谱,比如100个左右,那么一个一个根据化学常识确定也并不慢。 如果要跑几万个分子,那么首先考虑一下,这些分子是不是从某个已有的数据集批量下载下来的?如果是,在原数据集里往往可以找到所有分子的自旋多重度列表。 如果既要跑大量分子,又没有现成的自旋多重度可查,那就只能批量用比较便宜的方法试所有可能的自旋多重度,自动挑能量最低的了。但最好人工抽查几个分子,确认判断结果正确,尤其是要考虑到: (1)该做对称性破缺的时候有没有做对称性破缺?不仅包括开壳层单重态,即使二重态、三重态等等也可能存在很多反铁磁耦合的解,比如你研究一个体系,含一个Fe(II)和一个Fe(III),可能真实的自旋多重度是2,对应两者都是高自旋且是反铁磁耦合,结果你的自旋多重度2的计算收敛到了两者都是低自旋的态,导致高估自旋多重度2的能量,误把自旋多重度判断成10(两者都是高自旋,但铁磁耦合)。 (2)你试的最大的自旋多重度够不够大?仍然考虑以上的例子,假如你只试2、4、6、8,以为这就够了,不试10,那结果也不可靠。 |
zjxitcc 发表于 2024-3-22 21:46 我主要是因为要批量跑,然后每个产物都不一样,甚至可能自旋多重度也不一样。想寻找一个自动化一点的方法 |
本帖最后由 zjxitcc 于 2024-3-22 21:47 编辑 liaojianhan 发表于 2024-3-22 21:43 自旋多重度和电荷没写对,则计算完全错误,不是你想研究的东西(这时候你可能还没想清楚你想计算具体什么物种,需要的不是计算,而是停下来想一想“我要研究什么、我想算什么”)。 |
想请教一下各位大佬,对于单电子比较多的过渡金属体系,难以确认多重度和单电子数,我如果输入的多重度和单电子数不太合理的话,对于几何优任务,过渡态搜索任务是否会造成很严重的影响呢?(使用的方法有半经验方法和DFT方法) |
wzkchem5 发表于 2022-6-27 20:10 嗯嗯,谢谢 |
wzkchem5 发表于 2022-6-29 14:06 谢谢! |
tmdzd 发表于 2022-6-28 21:39 确实不一样。普通TDDFT只有后者,要得到前者需要做spin-flip TDDFT。 实际涉及的激发是哪个,就用哪个激发对应的轨道 |
本帖最后由 tmdzd 于 2022-6-29 05:21 编辑 wzkchem5 发表于 2022-6-29 03:46 感谢回复! 还有一个问题是 9 alpha→10 beta和9 beta → 10 beta的能量差应该不一样吧,毕竟9 alpha和9 beta的能量应该不一样, 如果做一个tddft计算会出现前者这种激发吗? 如果要画出MO来看电子转移,要用9a还是9b呢 |
tmdzd 发表于 2022-6-28 20:19 有可能,但因为交换alpha、beta电子不改变体系的总能量,所以一般总是把比较多的那种自旋定义为alpha。所以9个alpha、10个beta电子的体系,虽然存在,但一般是当作10个alpha、9个beta算的 |
sobereva 发表于 2017-1-2 12:06 老师想问一下有可能出现beta比alpha电子多的情况吗 比如一个开壳层体系,10个alpha 9个beta 然后激发的时候有可能出现9 alpha→10 beta的情况吗 还是一定是9 beta → 10 beta? |
wzkchem5 发表于 2022-6-27 20:22 感谢老师的解答 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 02:21 , Processed in 0.180595 second(s), 26 queries , Gzip On.