where 𝑟 6 denotes the position vector of the last atom in the polymer chain. The angular brackets ⟨...⟩ indicate the mean value over the polymer molecules and simulations timeframes. According to the Flory-Huggins theory, the gyration radius Rg of a polymer chain in a dilute solution scales as 𝑅 ' ∝ 𝑁 ? , where N is the number of monomers per chain (degree of polymerization) while v reflects the dependence of the polymer size on the solvent quality. In a good solvent (v = 3/5), the polymer swells due to the energetically more favorable monomer-solvent interaction. In a bad solvent (v = 1/3), the polymer collapses as the monomer-monomer attraction prevails over monomer-solvent interactions.