We perform a suite of simulations with realistic gravity and thermal balance in shells to quantify the role of the ratio of cooling time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) and the amplitude of density perturbations ($\delta \rho/\rho$) in the production of multiphase gas in the circumgalactic medium (CGM). Previous idealized simulations, focussing on small amplitude perturbations in the intracluster medium (ICM), found that cold gas can condense out of the hot ICM in global thermal balance when the background $t_{\rm cool}/t_{\rm ff} \lesssim 10$. Recent observations suggest the presence of cold gas even when the background profiles have somewhat large values of ${t_{\rm cool}/t_{\rm ff}}$. This partly motivates a better understanding of additional factors such as large density perturbations that can enhance the propensity for cooling and condensation even when the background ${t_{\rm cool}/t_{\rm ff}}$ is high. Such large density contrasts can be seeded by galaxy wakes or dense cosmological filaments. From our simulations, we introduce a condensation curve in the $(\delta \rho/\rho)$ - min$(t_{\rm cool}/t_{\rm ff})$ space, that defines the threshold for condensation of multiphase gas in the CGM. We show that this condensation curve corresponds to ${(t_{\rm cool}/t_{\rm ff})}_{\rm blob} \lesssim 10$ applied to the overdense blob instead of the background for which $t_{\rm cool}/t_{\rm ff}$ can be higher. We also study the modification in the condensation curve by varying entropy stratification. Steeper (positive) entropy gradients shift the condensation curve to higher amplitudes of perturbations (i.e., make condensation difficult). A constant entropy core, applicable to the CGM in smaller halos, shows condensation over a larger range of radii as compared to the steeper entropy profiles in the ICM.

We study an idealized 1D model for the evolution of hot gas in dark matter halos for redshifts $z=[0,6]$. We introduce a numerical setup incorporating cosmological accretion of gas, along with the growth of the halo, based on the Van den Bosch model for the average growth of halos as a function of cosmic time. We evolve one-dimensional Lagrangian shells with radiative cooling of the gas and heating due to feedback from the gas cooling and moving in toward the center. A simple Bondi accretion model on to a central black hole is used to include feedback heating. The setup captures some of the key characteristics of spherically symmetric accretion onto the halos: formation of virial shocks slightly outside $r_{200}$ and long-term thermal balance in the form of cooling and heating cycles. The gas density outside our initial halos at $z=6$ is constrained by requiring that the baryon fraction within the virial radius for non-radiative evolution be equal to the universal value at almost all times. The total mass in the cold phase (taken to be $\sim 10^4$ K) within $40$ kpc is tracked as a function of the halo mass and redshift. We compare the evolution of the cold gas mass to the observed stellar-mass versus halo mass relations, following which, we can constrain the feedback energy required for different halo masses and redshifts. We also compare and match the hot gas density and temperature profiles for our most massive halo to those of clusters observed upto redshift $2$. Our model is thus an improvement over the semi-analytic models in which isothermal condition and $\rho \propto r^{-2}$ are assumed.

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g., spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in halos critically depends on the ratio of the cooling time to the free-fall time ($t_{cool}/t_{ff}$). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focussed on {\em in situ} condensation of multiphase gas but spherical simulations studied condensation {\em anywhere} in the box. Unlike previous claims, our nonlinear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas condensing depends on the shape of the gravitational potential well; gas has more time to condense if gravitational acceleration decreases toward the center. In our idealized simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive $t_{cool}/t_{ff}$ far beyond the critical value close to 10. Triggered condensation is very prominent in plane-parallel simulations, in which a large amount of cold gas condenses out.