Acknowledgements: L.M. acknowledges funding support from the DAE through the NISER project RNI 4011. L.M. appreciates discussions with Sukrit Ranjan regarding common convergence issues in solvers used in photochemical kinetics codes. Thanks are also due to Jeehyun Yang for valuable discussions on the use of the Reaction Mechanism Generator (RMG). L.M. further thanks Karen Willacy for insightful discussions related to hydrocarbon chemistry. P.G. and D.D. thank Dwaipayan Dubey for helpful discussions on developing an improved version of the Dijkstra algorithm for rate analysis, and they also thank Namrata Rani for assistance with Gaussian 16. P.G. also thanks S. Maitrey for discussions on the pylsodes solver, which was developed for the PEGASIS astrochemical code and adapted for use in this work. S.M. expresses sincere gratitude to Mr. Nishil Mehta for sharing his experience and expertize in disequilibrium chemistry in the early stages of the project. L.M. gratefully acknowledges support from Breakthrough Listen at the University of Oxford through a sub-award to NISER under Agreement R82799/CN013, provided as part of the global Breakthrough Listen collaboration funded by the Breakthrough Prize Foundation. PBR received funding for this work from the Leverhulme Centrefor Life in the Universe Joint Collaborations Research Project GrantG112026, Project KKZA/237. We thank the anonymous referee for constructive comments that helped improve the manuscript.
Ghosh et al. (Sun,) studied this question.