$\texttt{slick}$: Modeling a Universe of Molecular Line Luminosities in Hydrodynamical Simulations

We present {\sc slick} (the Scalable Line Intensity Computation Kit), a software package that calculates realistic CO, [\ion{C}{1}], and [\ion{C}{2}] luminosities for clouds and galaxies formed in hydrodynamic simulations. Built on the radiative transfer code {\sc despotic}, {\sc slick} computes the thermal, radiative, and statistical equilibrium in concentric zones of model clouds, based on their physical properties and individual environments. We validate our results applying {\sc slick} to the high-resolution run of the {\sc Simba} simulations, testing the derived luminosities against empirical and theoretical/analytic relations. To simulate the line emission from a universe of emitting clouds, we have incorporated random forest machine learning (ML) methods into our approach, allowing us to predict cosmologically evolving properties of CO, [\ion{C}{1}] and [\ion{C}{2}] emission from galaxies such as luminosity functions. We tested this model in 100,000 gas particles, and 2,500 galaxies, reaching an average accuracy of $\sim$99.8\% for all lines. Finally, we present the first model light cones created with realistic and ML-predicted CO, [\ion{C}{1}], and [\ion{C}{2}] luminosities in cosmological hydrodynamical simulations, from $z=0$ to $z=10$.