Supplemental Material for Molecular simulations of heterogeneous ice nucleation I: Controlling ice nucleation through surface hydrophilicity Stephen J. Cox,1, 2 Shawn M. Kathmann,3 Ben Slater,1 and Angelos Michaelides1, 2, a) 1)Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, U.K. 2)London Centre for Nanotechnology, 17{19 Gordon Street, London WC1H 0AH, U.K. 3)Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States (Dated: 22 April 2015) This Supplemental Material contains further details of the simulation methods and the data fitting procedure used. A movie and still images of a nucleation event at the NP with Eads=∆Hvap ≈ 0:9 are also provided. a)Electronic mail:
[email protected] 1 I. FURTHER SIMULATIONS DETAILS To construct the NP, the ASE package1 was used: 6 atomic layers were used in the f1,0,0g family of directions; 9 layers in the f1,1,0g family; and 5 layers in the f1,1,1g family, except along the (1¯; 1¯; 1)¯ direction where no layers were used. As the equations of motion for the atoms in the NP were not integrated (i.e. they were fixed) no interaction potential was defined between them, and the adsorption energy Eads of a water monomer to the NP was therefore simply defined as the total energy after geometry optimization of a single water molecule at the center of the (111) face. The velocity Verlet algorithm was used to propagate the equations of motion of the water molecules, using a 10 fs time step.