I added a new pair style to LAMMPS for simulation Hydrophobic interactions.
I simulated a system containing 5 polymers (including 5 beads). In the absence of electrostatic interaction (coul/debye) our results are correct (Fig1). When I turned on the coul/debye interaction, our simulation results dramatically changed even though we didn’t have any charge in our simulation box (Fig 2).
All conditions in Fig 1 & 2 are the same ( no charge-hydrophobicity-temperature) but in Fig 2 I turned on the coul/debye pair style.