Full six-dimensional nonadiabatic quantum dynamics calculation for the energy pooling reaction O2(aΔ1)+O2(aΔ1)→O2(bΣ1)+O2(XΣ3)

Six new potential energy surfaces of four singlet states and two triplet states for the title oxygen molecule reaction along with the spin-orbit coupling among them have been constructed from the complete active space second-order perturbation theory with a 6-311+G(d) basis. Accurate integral cross sections are calculated with a full six-dimensional nonadiabatic time-dependent quantum wave packet method. The thermal rate constant based on the integral cross sections agrees well with the result of the experimental measurements, and the intersystem crossing effects are also discussed in this electronic energy-transfer process.