Energy conversion in cometary atmospheres-Hybrid modeling of 67P/Churyumov–Gerasimenko