In Figure , we present the radial distribution function of the hydroxyl oxygens to the water oxygens for the β anomers (see Figure for the α anomers plot). In Table , we show the position of the first peak, position of the first minima and coordination number at the first minima, for β -Glc, as reference. The six monosaccharides exhibit a very similar hydration pattern with water oxygens, most likely, 2.85 Å away from the hydroxyl oxygens (Table ). Looking at the radial distribution functions from our MD simulations performed without polarization (dotted lines in Figure ), we see shorter and broader peaks. These results indicate that polarization enhances the structure of proximal waters, but does not result in an expanded hydration shell as reported elsewhere. Rather, we compute an hydration shell of 3.5 Å for β -Glc, compared to 4.05 Å without polarization (Table ), and 3.8 Å as predicted by the polarizable Drude force field. In our simulations, we see a tighter, more ordered, hydration layer with higher coordination number, corresponding to a significantly denser hydration shell than has been previously observed with other force fields. In Table , we list the hydration number of each monomer, computed with and without mutual polarization, as well as the number of acceptor and donor hydrogen bonds (H-bonds) to the hydroxyls and nitrogen groups. We find that Glc is surrounded by 25 ± 8 water molecules, which is consistent with the 21 hydration number found with dielectric spectroscopy. Further, we observe a relatively large standard deviation of the hydration number for all monosaccharides, regardless of the anomeric form. This suggests a dynamic hydration layer where the water molecules exchange position on a timescale To further test the performance of these new AMOEBA force field parameters for monosaccharides, we ran MD simulations (150 ns in the NPT ensemble) of multiple β -Glc in water. Following previous work, we simulated six different concentrations, all below the concentration at which β -Glc aggregates: 0.25, 0.5, 1, 2, 3 and 4 mol/kg. In each case, we computed the radial distribution function of the center-ofmass (COM) of β -Glc, as well as their mean-square displacement. We fitted the mean-square displacement as a function of time with a linear function (Figure ) to extract the diffusion coefficient for each concentration: