RESUMO
Recently, thermoset vitrimer polymers have shown significant promise for structural applications because of their ability to be reshaped and remolded due to their covalent adaptive network (CAN). In these vitrimers, the transesterification reaction is responsible for the CAN, where the efficiency of the reaction is controlled either by organic or by organometallic catalysts. Understanding the mechanism of the transesterification reaction in the bulk phase using direct experimental techniques is extremely difficult due to the highly cross-linked complex structure of thermosetting vitrimers. Therefore, we use solution-phase experiments to investigate the catalytic efficiency and to guide density functional theory (DFT) simulations of the transesterification reaction mechanism with catalysts triazabicyclodecene (TBD), zinc acetate (Zn(OAc)2), 1-methylimidazole (1-MI), and dibutyltin oxide (DBTO). The estimated catalytic efficiency from the detailed DFT reaction path calculations follows the order TBD â³ DBTO â³ Zn(OAc)2 > 1-MI, which agrees with the experimental results. In addition to reaction path modeling, the mechanism and the relative rates of the transesterification reaction are analyzed with the assistance of Fukui indices as a measure of electrophilicity and nucleophilicity of atomic sites and with partial charges. It was found that the sum of the nucleophilicity index of the base and the electrophilicity index of the acid of the bifunctional catalysts correlates with the SN2 transition state and tetrahedral intermediate energies, which are related to the barrier of the rate-limiting step. This correlation provides a hypothesis for computational prescreening of potentially better catalysts that have an index in a range of values. These results provide a basis for understanding an important part of the mechanism of transesterification in vitrimer systems and may assist with designing new catalysts.