求助TDDFT中不同级别XC kernel的求法

枫沨

我想想大家请教一个问题，通过libxc直接得到的都是能量对rho以及sigma的各种组合到处。
以GGA泛函为例（下面给出的式子），在计算TDDFT中的xc kernel的时候 f_xc(r,r')是如何由libxc的结果组合成真正GGA级别的f_xc的数值的（而不是停留在LDA级别）。由具体的算法指导与例子吗，提前谢谢大家！
$$
f_{xc}(r,r')=\frac{\partial^2 G}{\partial \rho \partial \rho} - \nabla_j \frac{\partial^2 G}{\partial \nabla_j \rho \partial \rho} - \nabla_j \nabla_k \frac{\partial^2 G}{\partial \nabla_j \rho \partial \nabla_k \rho}
$$

hebrewsnabla

https://github.com/pyscf/pyscf/issues/1539 这个也是你问的？

可以阅读一下scf._response_functions里面RKS的get_response，进一步地他会调用ni.cache_xc_kernel

枫沨

hebrewsnabla 发表于 2022-12-14 16:28
https://github.com/pyscf/pyscf/issues/1539 这个也是你问的？

可以阅读一下scf._response_functions里 ...

是的，我以为大家都过节回家了。所以就在这里又问了一下。谢谢解答，我自己先去试试！