“…In addition, the application of Buckingham‐type potential often leads to an unphysical phenomenon, in which the potential energy becomes negative infinity when r is too small. According to previous studie,
24,25,37–39 a correction term for
is therefore adopted and described as below:
where
,
, and
are parameters chosen to match the first and second derivatives of
and
at
. As for the long‐range Coulomb interaction, an effective partial charge was set to each species, with the charge of O equal to −1.2 e. The calculations were sped up by the damped shifted force method
40 with a cutoff of 8.0 Å and damping parameter of 0.25 Å −1 .…”