Moreover, even though there was agreement between the all-atom CHARMM36 aggregation behaviour and the Martini 3 coarse-grained simulations, we decided to further compare the final Martini 3 model with its predecessor Martini 2. To this end, we performed three further replica simulations of the Glcα(1→1)Glcα system using the Martini 2 force field, utilizing the standard parameters for trehalose. Similar to the Martini 3 simulations, a unimodal distribution of aggregate sizes was observed in all three replicates with a peak in aggregate sizes corresponding to disaccharide dimer formation. Notably, however, the largest aggregate observed in the Martini 2 simulations comprised 15 individual disaccharides, indicating that while both Martini force fields did not reproduce the bimodal distribution obtained in the GLYCAM-06j simulation, Martini 2 exhibited a greater tendency towards the formation of large disaccharide aggregates when compared to Martini 3. It has been reported previously that the nonbonded interactions between glycans are too attractive in Martini 2 . The solution of this issue was a key target during development of the Martini 3 force field . It is therefore reassuring that our results indicate that the disaccharides are readily soluble in both Martini models, and that the Martini 3 force field exhibits a reduced extent of carbohydrate aggregation compared to Martini 2.