Dissociation and energy transfer in high-energy collisions of O2 playimportant roles in simulating thermal energy content and heat flux in flows around hypersonicvehicles. Furthermore, atomic oxygen reactions on the vehicle surface are an importantcontributor to heat shield erosion. Molecular dynamics modeling is needed to better understandthe relevant rate processes. Because it is necessary to model the gas flows in high-temperatureshock waves, electronically excited states of O2 and O can be populated, and molecular dynamicssimulations should include collisions of electronically excited species and electronicallynonadiabatic collisions. This requires potential energy surfaces and state couplings for many energetically accessible electronic states. Here we report a systematic strategy to calculate such surfaces and couplings. We have applied this method to the fourteen lowest-energy potential energy surfaces in the 3A′ manifold of O3, and we report a neural-network fit to diabatic potential energy matrix (DPEM). We illustrate the use of the resulting DPEM by carrying out semiclassical dynamics calculations of cross sections for excitation of O2 inat two collision energies; these dynamics calculations are carried out by the 3A′ collisions with Ocurvature-driven coherent switching with decay of mixing method.