The polarizable ion model (PIM)用CP2K如何实现？

duoconsole

在做熔盐的分子动力学相关的计算。文献当中使用了cp2k的Fist模块实现了极化离子力场（The polarizable ion model，PIM），力场的相关势函数见图片。
PIM力场主要包括四部分，库伦项，极化项、排斥项和色散项。我查阅cp2k手册，库伦项应该主要是在poisson部分设置，排斥项和色散项刚好是与cp2k中非键相互作用中的BMHFTD这个力场一样，但是目前不知道如何把极化项考虑进去，极化项涉及偶极矩的变化，在cp2k手册中也看到PIPOLE项可以迭代计算偶极矩，我应该如何设置CP2K的输入文件在得到偶极矩的情况下计算极化项，并把这个极化项考虑到相互作用势中呢？