计算化学公社

标题: 求助CP2K激发态性质计算考虑SOC时的赝势设置 [打印本页]

作者
Author: 许旭辉 时间: 2025-6-16 20:33
标题: 求助CP2K激发态性质计算考虑SOC时的赝势设置
在用cp2k计算考虑SOC的TDDFT时，报错“Spin-orbit coupling selected, but GTH potential without SOC parameters provided” ，看来以往的帖子显示只是warning，但是这里直接就是停了。具体输入文件如下。
还有一个问题就是，考虑soc后不是应该是旋量形式吗，为什么输出还是显示一个轨道占据数是2啊？希望大佬指教。

作者
Author: f294540290 时间: 2025-6-16 21:33
考虑SOC时不能带有赝势~ 老老实实用全电子基组来算

作者
Author: 许旭辉 时间: 2025-6-16 21:55

f294540290 发表于 2025-6-16 21:33
考虑SOC时不能带有赝势~ 老老实实用全电子基组来算

谢谢您，那用全电子的话，会是轨道单占据形式吗？

作者
Author: 许旭辉 时间: 2025-6-16 22:06
本帖最后由许旭辉于 2025-6-16 22:09 编辑

f294540290 发表于 2025-6-16 21:33
考虑SOC时不能带有赝势~ 老老实实用全电子基组来算

大佬，没有看到Mo的全电子基，能用GTH_SOC_POTENTIALS吗

作者
Author: sobereva 时间: 2025-6-17 06:52

f294540290 发表于 2025-6-16 21:33
考虑SOC时不能带有赝势~ 老老实实用全电子基组来算

用GTH赝势基组也可以做SOC-TDDFT，此时赝势文件要用带旋轨势的GTH_SOC_POTENTIALS

作者
Author: sobereva 时间: 2025-6-17 06:57

考虑soc后不是应该是旋量形式吗，为什么输出还是显示一个轨道占据数是2啊？

SOC-TDDFT仅仅是对TDDFT算的激发态考虑旋轨耦合，原理和下文说的一样。参考态计算依然是标量的。
使用Gaussian+PySOC在TDDFT下计算旋轨耦合矩阵元
http://sobereva.com/411（http://bbs.keinsci.com/thread-9442-1-1.html）
使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱
http://sobereva.com/462（http://bbs.keinsci.com/thread-12114-1-1.html）

作者
Author: 许旭辉 时间: 2025-6-17 09:24

sobereva 发表于 2025-6-17 06:57
SOC-TDDFT仅仅是对TDDFT算的激发态考虑旋轨耦合，原理和下文说的一样。参考态计算依然是标量的。
使用 ...

好的好的，谢谢老师

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