The computations in this study resulted in a total of eighteen fully optimized energy minima and thirteen associated intervening fully optimized first-order saddle points that connect the paths of the reactions of interest; these are all referred to as "stationary points". The fully optimized minima represent reactants, products, and intermediates, and the first-order saddle points represent transition states. These stationary points are in four distinct systems each with their own elemental compositions and assigned basis sets. Relative energy changes were calculated between stationary points along a reaction path within each system only. The complete set of configurations are listed in the supplement. Note that these discovered pathways determine the current best upper bounds to activation energies because competing pathways would need to have reaction rates that are equal to or greater than these mechanisms.