用CP2K计算偶极矩，输入文件如何写

刘由牛

请问各位大佬，如果想利用CP2K跑AIMD，计算偶极矩随时间的变化，输入文件该如何写？

洁然不同

DFT PRINT部分有个moments选项设置ON即可，默认计算到二阶，即电偶极矩和电四极矩张量同时输出，这是体系的总偶极矩。还可以在 DFT LOCALIZE中用Wannier函数计算，PRINT打开TOTAL_DIPOLE 显示总偶极矩（矢量和）。USE_HISTORY在MD中设置成TRUE,减少计算量。