Thermal degradation of unstrained single polymer chain: non-linear effects at work

We examine the thermally induced fracture of an unstrained polymer chain of discrete segments coupled by an anharmonic potential by means of molecular dynamics simulation with a Langevin thermostat. Cases of both under- and over-damped dynamics are investigated, and a comparison with recent studies of bond scission in model polymers with harmonic interactions is performed. We find that the polymer degradation changes qualitatively between the inertial regime and that of heavily damped dynamics. The role of bond healing (recombination) is also studied and probability distributions for the recombination times and overstretched bond lengths are obtained. Our extensive simulations reveal many properties of the scission dynamics in agreement with the notion of random breakdown of independent bonds, e.g., the mean time of chain rupture, <τ> follows an Arrhenian behavior with temperature T, and depends on the number of bonds N in the polymer as <τ> ∝ N(-1). In contrast, the rupture rates of the individual bonds along the polymer backbone indicate clearly the presence of self-induced inhomogeneity resulting from the interplay of thermal noise and nonlinearity. Eventually we examine the fragmentation kinetics during thermolysis. We demonstrate that both the probability distribution function of fragment sizes as well as the mean length of fragments at subsequent times t characterize degradation as predominantly a first order reaction.