AIMD中如何通过外加电场(势场)模拟电场下双电层中离子对电极表面的吸附

12345321

各位老师好，我想通过AIMD模拟电场作用下，充电时阳离子吸附在负极表面的过程，研究阳离子吸附时的构象。
体系是胆碱阳离子，三氟甲烷磺酸根离子，锌离子，水，和切出来的Zn-slab，400~500原子(其中120原子被冻结)，Cutoff 500Ry,ref cutoff 55,64核服务器,PBE+D3(BJ)/DZVP-MOLOPT-SR-GTH，开smearing。
想请教一下各位老师该如何通过外加电场模拟充电过程中双电层中电场分布呢，还是说直接用卢老师在CP2K培训班中提到的静态外电场：


&EFIELD
   INTENSITY 0.02
   POLARISATION 0 0 1
   &END EFIELD
或是加外势：
&EXTERNAL_POTENTIAL
  FUNCTION  (A/B)*Z
  PARAMETERS A 5
VALUES  …………
  &END EXTERNAL_POTENTIAL
双电层中电势绝对值与离子到极板位置关系如图所示

求各位老师解答，谢谢