As computing power and force field accuracy increase, large multicomponent and multiphasic systems are increasingly coming into focus. Here it is even more crucial to carefully consider finite size effects, related to the low copy number of each of the different components in the simulation system. Furthermore, given the sluggish dynamics in the dense phases and the slow exchange across the phase boundaries, accelerating the exploration of phase space via suitable enhanced sampling methods will be important, not only for atomistic but also for CG simulations.