Viscoelasticity near the gel-point: a molecular dynamics study

We report on extensive molecular dynamics simulations on systems of soft spheres of functionality f, i.e. particles that are capable of bonding irreversibly with a maximum of f other particles. These bonds are randomly distributed throughout the system and imposed with probability p. At a critical concentration of bonds, p_c approximately equal to 0.2488 for f=6, a gel is formed and the shear viscosity \eta diverges according to \eta ~ (p_c-p)^{-s}. We find s is approximately 0.7 in agreement with some experiments and with a recent theoretical prediction based on Rouse dynamics of phantom chains. The diffusion constant decreases as the gel point is approached but does not display a well-defined power law.