Carbodiimides are important crosslinkers in organic synthesis and are used in the isocyanate industry as modifier additives. Therefore, the understanding of their formation is of high importance. In this work, we present a theoretical B3LYP/6-31G(d) and SMD solvent model and experimental investigation of the formation of diphenylcarbodiimide (CDI) from phenyl isocyanate using a phosphorus-based catalyst (MPPO) in ortho-dichlorobenzene (ODCB) solvent. Kinetic experiments were based on the volumetric quantitation of CO2 evolved, at different temperatures between 40 and 80 °C. Based on DFT calculations, we managed to construct a more detailed reaction mechanism compared to previous studies which is supported by experimental results. DFT calculations revealed that the mechanism is composed of two main parts, and the rate determining step of the first part, controlling the CO2 formation, is the first transition state with a 52.9 kJ mol-1 enthalpy barrier. The experimental activation energy was obtained from the Arrhenius plot (ln k vs. 1/T) using the observed second-order kinetics, and the obtained 55.8 ± 2.1 kJ mol-1 was in excellent agreement with the computational one, validating the complete mechanism, giving a better understanding of carbodiimide production from isocyanates.
Keywords: DFT; activation energy; carbodiimide; gas volumetry; isocyanate; organocatalysis; reaction kinetics.