MPWB1K calculations of stepwise encapsulations: Lix@C60

Computations are carried out on Li@C60, Li2@C60 and Li3@C60 with the newly introduced MPWB1K density functional. The computed equilibrium thermodynamics shows that Li2@C60 is at least by two orders of magnitude less populated than Li@C60 when the saturated Li vapor is considered while for Li3@C60 it is at least four orders of magnitude compared to Li@C60. The evaluations correspond to the available experimental findings.