was imposed to adjacent triplets of monomers, where k ba is the bond-angle bending constant, ϕ the bond angle, and ϕ 0 = π the equilibrium bond angle. The value of k ba = 1ϵ was used as a fitting parameter to match the radius of gyration of the AA PEG molecules (Fig. ). The system was thermalised to a temperature of T = 298.15 K by coupling the particles to a heat bath using the Langevin thermostat with friction coefficient γ = 3r θ πησ, where r θ is a numerical prefactor. The value of η = 1 mPa • s was chosen to mimic the interaction with an implicit solvent corresponding to bulk water, and the coefficient r θ was chosen so that R 1 (f ) spectra show the same dependence in the frequency f as the one extracted from CoM calculations (Fig. ): for low PEG concentration c PEG < 1, r θ was taken equal to 1.5, and for c PEG ≥ 1, r θ was taken equal to 5. Finally the stochastic equations of motion were integrated using the Velocity-Verlet method . The scripts are available in the DaRUS repository .