Atomic-Scale Mechanistic Insights into the Ring-Opening Polymerization of Elemental Sulfur

A detailed understanding of the composition and polymerization mechanism of elemental sulfur remains a decades long unresolved question for modern chemistry. However, the dynamic nature of molten sulfur significantly complicates its accurate characterization. To overcome this challenge, we performed the first comprehensive molecular dynamics (MD) simulations using a ReaxFF reactive force field specifically parameterized to capture the complex ring-opening polymerization dynamics of elemental sulfur. Rigorous development of the force field parameters, trained against extensive quantum mechanical datasets, was key to enabling these large-scale (>10 000 atoms) reactive MD simulations at polymerization-relevant temperatures. This study provides the first detailed atomic-level description of liquid sulfur, elucidating temperature-dependent molecular composition and offering unprecedented insights into sulfur polymerization mechanisms. These are the first simulations to reveal the formation of remarkably large macrocyclic sulfur rings at the onset of polymerization-a discovery that challenges longstanding mechanistic misconceptions, thus reshaping our understanding of sulfur polymerization. Our findings highlight the power of molecular dynamics in exploring complex polymerization processes, with broad impact in dynamic covalent chemistry and covalent adaptable polymers.