Multiscale simulation on product distribution from pyrolysis of styrene-butadiene rubber

Abstract

Pyrolysis of styrene-butadiene rubber receives renewed attention due to its application in tackling the waste tire disposal problem while allowing energy recovery. The density functional theory calculation (DFT) and ReaxFF molecular dynamics simulation (MD) are adopted to study the pyrolysis process with the variation of temperature and pressure. The bond dissociation energies of intramonomer and intermonomer bonds in trimers with different linking methods are calculated by DFT, where the bond with low energy tends to break during the pyrolysis process. The following MD simulation shows the pyrolysis product distribution of chain segments in styrene-butadiene rubber, where bond breaking positions in MD agree well with corresponding results in DFT and experiment. The next nearest neighbor bonds (single bonds) connected with double bond or benzene usually have lower dissociation energies than other single bonds and prone to break during the pyrolysis process. And thus, the intermonomer bonds tend to break at relatively low temperatures (around 650 K in experiment) prior to intramonomer bonds, which result in the emergence of monomers. With the temperature increase, intramonomer bonds are broken and thus large fragments are further pyrolyzed into small ones (e.g., C2 and C). Besides, the pressure strongly influences the product distribution, where high pressures promote the occurrence of secondary reactions.