求助cp2k计算极化子和束缚激子的方法

Shana

各位老师好。由于课题需要，我需要计算极化子和束缚激子，也找到一些使用cp2k计算极化子和激子的文献（Chemistry of Materials 2020 32 (19), 8393-8400；[color=var(--sd-ui-anchor-colour)][size=1em]https://doi.org/10.1016/j.apsusc.2015.10.164）。但是对于如何使用cp2k计算极化子还不是很清楚，因此向各位有经验的老师请教。
（1）我原先是使用vasp计算极化子，在vasp中使用通过FERWE和FERDO，限制spin-up和spin-down 的轨道占据，进而指定电子占据到特定的轨道，那么在cp2k中有类似的指定轨道占据方法吗？

（2）我原先使用vasp计算极化子时，需要预先畸变结构，才能使顺利弛豫出特定的极化子构型，在cp2k中也需要预先畸变结构吗？

（3）我在论坛中看到有坛友使用LSD和UKS计算极化子，也看到有CDFT方法计算，这两种方法等价吗？哪一种更为方便一点呢？

（４）在vasp计算中，会生成PROCAR文件，里面记录了每个K-S轨道的组成，在ｃｐ２ｋ中能生成类似的文件吗？方便检查各个K－S轨道的组成．

谢谢各位老师和同学的指点．

leeru

一般来说，预先畸变本身带有主观性，不算是真正意义上的第一性原理。算极化子最好是epw，算激子最好用yambo或者berkeleygw。如果非要用cp2k算极化子，可以使用cdft+pbe0(hse06)，但依然不能消除电子的自相互作用SIC。当然也可以在PBE。

leeru

cp2k在pbe下加入SIC校正，结合UKS和预先畸变也凑合。只是SIC那几个参数要调，而且超胞体积比较大，不然不太好收敛。

Shana

leeru 发表于 2024-9-1 09:07
cp2k在pbe下加入SIC校正，结合UKS和预先畸变也凑合。只是SIC那几个参数要调，而且超胞体积比较大，不然不太 ...

感谢回复。我去试试

Shana

leeru 发表于 2024-9-1 09:05
一般来说，预先畸变本身带有主观性，不算是真正意义上的第一性原理。算极化子最好是epw，算激子最好用yambo ...

谢谢您的回复。目前只是用ΔSCF近似一下，可以请教您一下cdft怎么使用吗？

leeru

Shana 发表于 2024-9-4 20:08
谢谢您的回复。目前只是用ΔSCF近似一下，可以请教您一下cdft怎么使用吗？

https://sites.google.com/view/theorchem/lecture/tutorial参考这个教程，这个是用cp2k的cdft研究TiO2中的极化子。

Shana

leeru 发表于 2024-9-4 21:25
https://sites.google.com/view/theorchem/lecture/tutorial参考这个教程，这个是用cp2k的cdft研究TiO2中 ...

感谢您的指导