For gaining insights into interactions in periodic systems, an analysis is developed based on the fragment molecular orbital method combined with periodic boundary conditions. The adsorption energy is decomposed into guest and surface polarization and deformation energy, guest−surface and guest−guest interactions, and the vibrational free energy. The analysis is applied to the adsorption of guest molecules to Ih (001) ice surface. The cooperativity effects result in a non-linear change in the adsorption energy with coverage due to many-body effects. The role of dispersion is found to be dominant for guests with long hydrophobic tails. A rule is proposed relating the length of the alkyl tail with the formation of the guest layer. The computed binding enthalpies are in good agreement with experimental values. For high coverage, adsorbed molecules can form an ordered layer known as self-assembled monolayer.