The CHARMM36 carbohydrate parameter set did not adequately reproduce experimental thermodynamic data of carbohydrate interactions with water or proteins or carbohydrate self-association; thus, a new nonbonded parameter set for carbohydrates was developed. The parameters were developed to reproduce experimental Kirkwood-Buff integral values, defined by the Kirkwood-Buff theory of solutions, and applied to simulations of glycerol, sorbitol, glucose, sucrose, and trehalose. Compared to the CHARMM36 carbohydrate parameters, these new Kirkwood-Buff-based parameters reproduced accurately carbohydrate self-association and the trend of activity coefficient derivative changes with concentration. When using these parameters, preferential interaction coefficients calculated from simulations of these carbohydrates and the proteins lysozyme, bovine serum albumin, alpha-chymotrypsinogen A, and RNase A agreed well with the experimental data, whereas use of the CHARMM36 parameters indicated preferential inclusion of carbohydrates, in disagreement with the experiment. Thus, calculating preferential interaction coefficients from simulations requires using a force field that accurately reproduces trends in the thermodynamic properties of binary excipient-water solutions, and in particular the trend in the activity coefficient derivative. Finally, the carbohydrate-protein simulations using the new parameters indicated that the carbohydrate size was a major factor in the distribution of different carbohydrates around a protein surface.