The simulation is performed with N P T ensemble and the temperature of 333 K is maintained using the Nosé-Hoover thermostat (coupling constant 0.5 ps) and a pressure of 1 atm is maintained using the Parrinello-Rahman barostat(coupling constant 2 ps). Electrostatic interaction is treated using particle-mesh Ewald summation with a Coulomb cutoff 1.6 nm same as van der Waals cutoff. All hydrogen bonds are constrained by the LINCS algo-rithm. We leave out the first 20 ns trajectory for equilibration (for details see supplementary information) and after equilibration a run length of 1µs is used for data analysis.