This data is a supplement to the publication https://doi.org/10.1016/j.fuel.2018.12.047. The data includes Turbomole input and output files. The calculations are performed using DFT/TPSSh-D3/TZVP method and Turbomole version 7.2. The equilibrium structures for the reactions of polyaromatics are named as following: CXHY (S1) + O -> CXHYO (S2) -> CXHY-1O (S3) + H (i) O addition and H abstraction CXHY-1O (S3) [-> S4 -> S5 ] -> CX-1HY-1 (S6) + CO (ii) Single, two, or three step CO elimination The transition state structures are named according to the naming of the corresponding reactant and product. For example, the transition state connecting the structure S3 to S5 is named as T35. Under some of the transition state directories, intrinsic reaction coordinate calculation output can be found under the directories named as "IRC". The LibreOffice Calc spreadsheet "SUPPINFO.ods" includes the activation and reaction energies to the reaction steps.