Empirical force fields for computer simulations of carbohydrates are often implicitly assumed to be valid also at temperatures different from room temperature for which they were optimited: Herein, the temperature dependence of the hydroxymethyl group rotamer populations in short oligogaccharides is invegtigated using Molecular dynamics simulations and NMR spectroscopy. Two oligosaccharides, methyl beta-cellobioside and beta-cellotetraose were simulated using three different carbohydrate force fields (CHARMM C35, GLYCAM06, and GROMOS 56A(carbo)) in combination with different water models (SPC, SPC/E, and TIP3P) using replica exchange molecular dynamics simulations. For comparison, hydroxymethyl group rotamer populations were investigated for methyl beta-cellobioside and cellopentaose based- on measured NMR (3)J(H5,H6) coupling constants, in the latter case by using a chemical shift selective NMR-filter. Molecular dynamics simulations in combination with NMR spectroscopy show that the temperature dependence of the hydroxymethyl rotamer population in these short cellooligomers, in the range 263-344 K, generally becomes exaggerated in simulations when compared to experimental data, but also that it is dependent on simulation conditions, and most notably properties of the water model.