All simulations use all-atom representations and molecular dynamics (MD) simulations were performed using LAMMPS 40 with a 0.025 ps timestep under NVE dynamics. The ReaxFF 41 reactive force field was used and partial atomic charges were calculated using the charge equilibration scheme at every step in the simulations with a tolerance of 1x10 -6 .