Numerical simulations of the polymerization reaction of deuterated diacetylene 2,4‐hexadiynylene bis(p‐toluenesulfonate) crystal

In this paper, we developed two types of programs in order to simulate the polymerization reaction of a fully deuterated crystal of diacetylene 2,4‐hexadiynylene bis(p‐toluenesulfonate) (pTS‐D). The first simulation is based on a modification of Baughman's model, a classical model for simulating the polymerization of diacetylene crystals. The agreement between the simulated and experimental results concerning the reaction kinetics is satisfactory. With this simulation algorithm, we take into account the experimental observation that the polymerization of pTS‐H and pTS‐D crystals is really a random process of formation of polymer chains along the crystallographic axis b . The second simulation is based on the Monte Carlo method, which permits not only to simulate the kinetics of the reaction, but also the chain‐length distribution in the hydrogenated and deuterated compounds. These two types of simulations were already developed for the hydrogenated crystal of diacetylene, named pTS‐H. Two main modifications are applied in the case of pTS‐D for taking into account experimental results: in the first the rate constants of chain‐terminating microscopic processes are different in pTS‐H and pTS‐D which must be considered. The second modification concerns the evolution of the lattice deformation during the course of polymerization. The experimental variation of the b parameter as a function of polymer content X in pTS‐D is different from that in pTS‐H; this result is important to consider when calculating the activation energy of the initiation and propagation microscopic processes.