A Stellar Feedback Origin for Neutral Hydrogen in High-Redshift Quasar-Mass Halos

Observations of quasar pairs reveal that quasar host halos at z ~ 2 have large covering fractions of cool dense gas (≳ 60% for Lyman limit systems within a projected virial radius). Most simulations have so far failed to explain these large observed covering fractions. We analyze a new set of 15 simulated massive halos with explicit stellar feedback from the FIRE project, covering the halo mass range M_h ≈ 2 x 10^(12) - 10^(13) M_☉ at z = 2. This extends our previous analysis of the circum-galactic medium of high-redshift galaxies to more massive halos. Feedback from active galactic nuclei (AGN) is not included in these simulations. We find covering fractions consistent with those observed around z ~ 2 quasars. The large HI covering fractions arise from star formation-driven galactic winds, including winds from low-mass satellite galaxies that interact with the cosmological infalling filaments in which they are typically embedded. The simulated covering fractions increase with both halo mass and redshift over the ranges covered, as well as with resolution. Our simulations predict that galaxies occupying dark matter halos of mass similar to quasars but without a luminous AGN should have Lyman limit system covering fractions comparable to quasars. This prediction can be tested by measuring covering fractions transverse to sub-millimeter galaxies or to more quiescent galaxies selected based on their high stellar mass.