In addition to RDF analysis, the locally averaged Steinhardt order parameters, 〈𝑞 @ 〉 and 〈𝑞 B 〉, are calculated for each water molecule, to further address the possibility of ice crystallization . The Steinhardt order parameters are local bond order parameters based on spherical harmonics that have been frequently used in the literature to monitor ice formation . In the locally averaged form, 𝑞 @ and 𝑞 B of a given molecule is averaged over all neighbors within the first RDF shell, to yield 〈𝑞 @ 〉 and 〈𝑞 B 〉. These order parameters have been demonstrated to yield good separation between ice Ic, ice Ih, and water 〈𝑞 @ 〉 vs. 〈𝑞 B 〉 distributions at lower temperatures, allowing for the identification of ice molecules based on the 〈𝑞 @ 〉 and 〈𝑞 B 〉 values of an individual molecule . In this work, the locally averaged Steinhardt order parameters were calculated using the position of oxygen atoms from all 10 simulation trajectories at the conclusion of the simulations.