Recently, it has been shown that the reaction mechanism in self-healing diphenyl dichalcogenide-based polymers involves the formation of sulfenyl and selenyl radicals. These radicals are able to attack a neighbouring dichalcogenide bond via a three-membered transition state, leading to the interchange of chalcogen atoms. Hence, the chain mobility is crucial for the exchange reaction to take place. In this work, molecular dynamics simulations have been performed in a set of disulfide- and diselenide-based materials to analyze the effect of the molecular structure in the chain mobility. First of all, a validation of the computational protocol has been carried out, and different simulation parameters like initial guess, length of the molecular chains, size of the simulation box and simulation time, have been evaluated. This protocol has been used to study the chain mobility and also the self-healing capacity, which depends on the probability to generate radicals ( ρ ), the barrier of the exchange reaction ( Δ G ) and the mobility of the chains ( ω ). The first two parameters have been obtained in previous quantum chemical calculations on the systems under study in this work. After analyzing the self-healing capacity, it is concluded that aromatic diselenides (PD-SeSe) are the best candidates among those studied to show self-healing, due to lower reaction barriers and larger ω values.
Keywords: molecular dynamics; non-covalent interactions; self-healing polymers.