The specific number of glycine and water molecules in the simulation box for each polymorph can be found in supplementary table 1 in the supplementary material. To create an initial fluid, 2700 glycine molecules were solvated with SPC/E water using the gromacs/5.0 36 solvate command at a concentration of 0.34 g/ml. The fluid was then equilibrated by first performing NVT MD for 6 ns, and then NPT MD for an additional 6 ns at 380 K to prevent any crystallization from occurring. The box length of the equilibrated fluid was 9.8 nm. A single a-, b-, or 𝛾-glycine spherical nanonanocrystal, with a 2.3 nm radius, was then embedded in the equilibrated fluid using the gromacs/5.0 solvate command. Nanocrystals were sliced out of their respective bulk crystal structures obtained from the Cambridge Structural Database . The nanocrystals had an initial net zero dipole-alignment with the x-axis, as discussed below in section 3c. The number of solid molecules in the nanocrystals cut from their respective bulk structures can be obtained in supplementary table . The energy of the nanocrystal embedded in fluid system was then minimized in gromacs/5.0, using steepest descent minimization, with an energy minimization tolerance of 100.0 kJ/mol/nm. Production runs were then launched from this energy minimized configuration. Production simulations were performed for 4.2 ns.