Fixed-node diffusion Monte Carlo shows promise for modeling reaction thermochemistry of hydrocarbon-based radicals.

Abstract

Reliable thermodynamic and kinetic properties of free radical polymerization reactions are essential for synthesizing both primary polymeric materials and specialty polymers. The computational generation of these data from quantum chemistry requires a time-efficient method capable of capturing the essential physics. One such method, fixed-node diffusion Monte Carlo (FN-DMC) (using single Slater-Jastrow trial wavefunctions), has demonstrated the capability to recover 90%-95% of missing dynamic correlation energy for typical systems. In this study, methyl radical addition to ethylene serves as a simple model to test FN-DMC's ability to calculate enthalpies of reaction and activation energies with different time steps, antisymmetric trial wavefunctions, basis set sizes, and effective core potentials. The FN-DMC computational protocol thus defined for methyl radical addition to ethylene is subsequently benchmarked against Weizmann-1 and experimental reaction enthalpies from Lin et al.'s test set of 21 radical addition and 28 hydrogen abstraction enthalpies. Our findings reveal that FN-DMC consistently generates reaction enthalpies with chemical accuracy, exhibiting mean absolute deviation of 3.5(7) and 1.4(8) kJ/mol from the Weizmann-1 reference for radical addition and hydrogen abstraction reactions, respectively. Given its favorable computational scaling and high degree of parallelizability, we, therefore, recommend more comprehensive testing of FN-DMC with effective core potentials to address more extensive and intricate polymerization reactions and reactions with other radicals.