Mitali Damle
We analyse the circumgalactic medium (CGM) surrounding the Milky Way-Andromeda system in three different highest resolution realizations from the HESTIA suite of cosmological constrained simulations of the Local Group. Each of the three realizations, along with an MW-M31 pair, also comprise of a host of satellite systems. The positions and kinematics of all the constituents of each realization depend both on the initial conditions as well as the subsequent hydrodynamic interactions among themselves.
An ensemble of five ions, two (HI and SiIII- hereafter, low ions) which trace the cold and cool-ionized gas and three (OVI, OVII and OVIII- hereafter, high ions) which primarily track the hot-diffuse gas, is modelled using the Cloudy photoionization code. Mollweide projection maps are generated for all the three realizations and for all five ions using the Healpy python package. Dependence of ion column densities as a function of impact parameter (out to 4.0 Rvir) is plotted.
Not only do we see cold gas clumps (as expected), but we also see extended diffuse, cool gas streams that stretch far beyond the virial radii of MW and M31. Low ions show an overall larger column density dispersion as opposed to the high ions which exhibit rather flat column density-impact parameter profiles. Our SiIII column density profiles match well with corresponding results from the EAGLE simulations but are found to be a factor of two higher than those seen in the Project AMIGA survey.
Simulations have been an excellent counterpart to the observations since the past few decades. Modern cosmological as well as idealized simulations not only provide with possible guiding directions to the future observational projects, but also act as robust testing grounds for existing and new theories. Cosmological simulations have been able to well-replicate the large-scale structure of the Universe, the Cosmic Web pattern and the morphology of galaxies. However, in the context of the LG, most of these simulations have so far had modest success in getting the local small-scale structure right. The HESTIA suite of constrained simulations goes a step further in this regard, wherein, they manage to get both the large-scale as well as small-scale cosmography of the LG correct.
Our aim here is to use three of the highest-resolution simulations from HESTIA in order to obtain a better understanding about the tenuous gas (more commonly called as the Circumgalactic medium or CGM) that surrounds each galaxy in our LG. The gas comprising this CGM is observed to be multiphase in nature i.e. the presence of a wide range of temperatures and densities gives rise to an ensemble of ions, all of which co-exist in the same physical space. A rigorous knowledge of the CGM, thus, necessitates a thorough ionization modeling implementation; in our case we use the Cloudy photoionization modeling code and model five ions- two of which trace the cold and cool-ionized CGM (HI and SiIII; Low ions) and three, which trace the hot, diffuse CGM (OVI, OVII and OVIII; High ions). Thereafter, we set the midpoint of MW-M31 as our frame of reference and generate All-sky projection maps (Mollweide maps) for all the five ions and three simulations. We also plot the 2D column densities for each ion as a function of the impact parameter (Rvir; a radius characterizing the physical extent of any galaxy), in order to track the changes in respective ion column densities as we move from the inner, denser regions to the outer, more diffuse regions of the galaxies.
We see that the clumpiness of the gas decreases as we go from the low ions to the high ions (i.e. the gas becomes increasingly uniform). Low ions also show a wider range of column density distributions (seen as thicker column density profiles) as compared to High ions (seen as thinner column density profiles). Many of the satellite galaxies, due to the gravitational influence of the bigger MW-M31 systems, show significant distortion and stripping features (as is evident in the 2D projection maps in Page 1). Finally, our results of the SiIII column densities for M31 are seen to match rather well with those from EAGLE cosmological simulations. However, our results are about a factor of two higher than the corresponding results from the Project AMIGA observational survey of the Andromeda galaxy.