电荷平衡问题：生物磷脂双分子层在非水溶剂条件下的电荷平衡问题

ysrhhhh

目前，想通过LAMMPS中的NEMD方法来计算生物界面和不同材料之间的界面热阻值，因此构建了磷脂双分子层，但是构建的磷脂双分子层的电量是非中性的，因此需要添加离子进行电荷平衡；但是我的研究目的是计算生物双分子层和金属/液氮，两种不同的接触材料之间的界面热阻，因此无法像水溶液一样加入na、cl离子进行电荷平衡。看到一些资料可以使用dummy ions（虚拟离子）的形式来实现电荷平衡：1.带有电荷（±1e），用于电荷补偿；
2.无质量或质量极小（如0.0001 amu），不影响体系动量；
3.Lennard-Jones 参数极小，不参与结构排斥；
4.不与系统内原子发生排斥/吸引相互作用（通过 pair_coeff 设置）
但是在目前的研究中好像没有找到相关的论文，不知道这种方法的可行性，所以想请大家帮忙看是否合理或者存在别的方法实现系统的电荷平衡