We use the FIRE-2 cosmological simulations to study the formation of a virial temperature, quasi-static gas phase in the circumgalactic medium (CGM) at redshifts 0<z<5, and how the formation of this virialized phase affects the evolution of galactic discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average sub-virial temperatures (T<<T_vir), large pressure fluctuations and supersonic inflow/outflow velocities, to virial temperatures (T~T_vir), uniform pressures and subsonic velocities. This transition occurs when the outer CGM (~0.5 R_vir) is already subsonic and has a temperature ~T_vir, indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We demonstrate that virialization of the inner CGM coincides with abrupt changes in the properties of the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disc, star formation transitions from `bursty' to `steady,' and stellar-driven galaxy-scale outflows are suppressed. Our results hence suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic discs, rather than with the quenching of star formation as often assumed.

Star formation is observed to be strongly correlated to dense regions of molecular gas. Although the exact nature of the link between star formation and molecular hydrogen is still unclear, some have suggested that shielding of dense gas by dust grains is the key factor enabling the presence of both. We present a sub-grid model for use in galaxy formation simulations in which star formation is linked explicitly to local dust shielding. We developed and tested our shielding and star formation models using smoothed particle hydrodynamic simulations of solar and sub-solar metallicity isolated Milky Way-mass disk galaxies. We compared our dust shielding-based star formation model to two other star formation recipes that used gas temperature and H$_2$ fraction as star formation criteria. We further followed the evolution of a dwarf galaxy within a cosmological context using both the shielding and H$_2$-based star formation models. We find that the shielding-based model allows for star formation at higher temperatures and lower densities than a model in which star formation is tied directly to H$_2$ abundance, as requiring H$_2$ formation leads the gas to undergo additional gravitational collapse before star formation. However, the resulting galaxies are very similar for both the shielding and H$_2$-based star formation models, and both models reproduce the resolved Kennicutt-Schmidt law. Therefore, both star formation models appear viable in the context of galaxy formation simulations.