计算有机铁电材料的偶极子，cp2k的电场设置问题

jasondeng

想请教各位大佬，最近我在计算有机铁电材料的偶极子的时候，想通过外加高斯电场，来观察变化，但是我一加电场  &EFIELD
   ENVELOP GAUSSIAN
   INTENSITY 1E5
   PHASE 0
   POLARISATION: 0.0 1.0 0.0
   WAVELENGTH 3E5
   &GAUSSIAN_ENV
    SIGMA 1.06E1
    T0 2.5E2
   &END GAUSSIAN_ENV
  &END EFIELD
就计算不收敛，我意识到应该是&EFIELD只能作用于非周期性体系的原因，因为berry相位法只能计算周期性结构，但是&PERIODIC_EFEILD只能设置恒定电场，无法解决我需要设置外加高斯电场的需求，请问大佬们，这种问题应该怎么解决啊，求助。还想问一下，如果是计算体系的自发极化，我是只需要计算该原胞的偶极子，开非周期；还是要考虑周期性啊。