LIMFAST. I. A Seminumerical Tool for Line Intensity Mapping

We present LIMFAST, a seminumerical code for simulating high-redshift galaxy formation and cosmic reionization as revealed by multitracer line intensity mapping (LIM) signals. LIMFAST builds upon and extends the 21cmFAST code widely used for 21 cm cosmology by implementing state-of-the-art models of galaxy formation and evolution. The metagalactic radiation background, including the production of various star formation lines, together with the 21 cm line signal tracing the neutral intergalactic medium (IGM), is self-consistently described by photoionization modeling and stellar population synthesis coupled to the galaxy formation model. We introduce basic structure and functionalities of the code, and demonstrate its validity and capabilities by showing broad agreements between the predicted and observed evolution of cosmic star formation, IGM neutral fraction, and metal enrichment. We also present the LIM signals of 21 cm, Lyα, Hα, Hβ, [O ii], and [O iii] lines simulated by LIMFAST, and compare them with results from the literature. We elaborate on how several major aspects of our modeling framework, including models of star formation, chemical enrichment, and photoionization, may impact different LIM observables and thus become testable once applied to observational data. LIMFAST aims at being an efficient and resourceful tool for intensity mapping studies in general, exploring a wide range of scenarios of galaxy evolution and reionization and frequencies over which useful cosmological signals can be measured.