The Monte Carlo simulation started by defining the initial chain length distribution of the PEG molecules. The chain lengths, ni for i = 1, 2, 3, …, m (where m = number of chains with ni monomeric units), were assumed to follow a normal distribution with a mean chain length 𝜆 (representing the average degree of polymerization 𝑋 ! """" ), and a standard deviation s.