Clustering of MgII absorption line systems around massive galaxies: an important constraint on feedback processes in galaxy formation

We use the latest version of the metal line absorption catalogue of Zhu & M\'enard (2013) to study the clustering of MgII absorbers around massive galaxies (~10^11.5 M_sun), quasars and radio-loud AGN with redshifts between 0.4 and 0.75. Clustering is evaluated in two dimensions, by binning absorbers both in projected radius and in velocity separation. Excess MgII is detected around massive galaxies out to R_p=20 Mpc. At projected radii less than 800 kpc, the excess extends out to velocity separations of 10,000 km/s. The extent of the high velocity tail within this radius is independent of the mean stellar age of the galaxy and whether or not it harbours an active galactic nucleus. We interpret our results using the publicly available Illustris and Millennium simulations. Models where the MgII absorbers trace the dark matter particle or subhalo distributions do not fit the data. They overpredict the clustering on small scales and do not reproduce the excess high velocity separation MgII absorbers seen within the virial radius of the halo. The Illustris simulations which include thermal, but not mechanical feedback from AGN, also do not provide an adequate fit to the properties of the cool halo gas within the virial radius. We propose that the large velocity separation MgII absorbers trace gas that has been pushed out of the dark matter halos, possibly by multiple episodes of AGN-driven mechanical feedback acting over long timescales.Comment: 10 pages, 11 figures, accepted in MNRA

Zooming in on accretion - II. Cold Circumgalactic Gas Simulated with a super-Lagrangian Refinement Scheme

In this study we explore the complex multi-phase gas of the circumgalactic medium (CGM) surrounding galaxies. We propose and implement a novel, super-Lagrangian 'CGM zoom' scheme in the moving-mesh code AREPO, which focuses more resolution into the CGM and intentionally lowers resolution in the dense ISM. We run two cosmological simulations of the same galaxy halo, once with a simple 'no feedback' model, and separately with a more comprehensive physical model including galactic-scale outflows as in the Illustris simulation. Our chosen halo has a total mass of ~10^12 Msun at z ~ 2, and we achieve a median gas mass (spatial) resolution of ~2,200 solar masses (~95 parsecs) in the CGM, six-hundred (fourteen) times better than in the Illustris-1 simulation, a higher spatial resolution than any cosmological simulation at this mass scale to date. We explore the primary channel(s) of cold-phase CGM gas production in this regime. We find that winds substantially enhance the amount of cold gas in the halo, also evidenced in the covering fractions of HI and the equivalent widths of MgII out to large radii, in better agreement with observations than the case without galactic winds. Using a tracer particle analysis to follow the thermodynamic history of gas, we demonstrate how the majority of this cold, dense gas arises due to rapid cooling of the wind material interacting with the hot halo, and how large amounts of cold, ~10^4 K gas can be produced and persist in galactic halos with Tvir ~ 10^6 K. At the resolutions presently considered, the quantitative properties of the CGM we explore are not appreciably affected by the refinement scheme.Comment: MNRAS submitted, comments welcome. High-res version at http://www.mpa-garching.mpg.de/~dnelson/papers/Suresh19_zooming2.pd

Zooming in on the circumgalactic medium: resolving small-scale gas structure with the GIBLE cosmological simulations

We introduce Project GIBLE (Gas Is Better resoLved around galaxiEs), a suite of cosmological zoom-in simulations where gas in the circumgalactic medium (CGM) is preferentially simulated at ultra-high numerical resolution. Our initial sample consists of eight galaxies, all selected as Milky Way-like galaxies at $z=0$ from the TNG50 simulation. Using the same galaxy formation model as IllustrisTNG, and the moving-mesh code AREPO, we re-simulate each of these eight galaxies maintaining a resolution equivalent to TNG50-2 ($m_{\rm{gas}}$ $\sim$ $8 \times 10^5 {\rm M}_{\odot}$). However, we use our super-Lagrangian refinement scheme to more finely resolve gas in the CGM around these galaxies. Our highest resolution runs achieve 512 times better mass resolution ($\sim$ $10^3 {\rm M}_{\odot}$). This corresponds to a median spatial resolution of $\sim$ $75$ pc at $0.15~R_{\rm{200,c}}$, which coarsens with increasing distance to $\sim$ $700$ pc at the virial radius. We make predictions for the covering fractions of several observational tracers of multi-phase CGM gas: HI, MgII, CIV and OVII. We then study the impact of improved resolution on small scale structure. While the abundance of the smallest cold, dense gas clouds continues to increase with improving resolution, the number of massive clouds is well converged. We conclude by quantifying small scale structure with the velocity structure function and the auto-correlation function of the density field, assessing their resolution dependence. The GIBLE cosmological hydrodynamical simulations enable us to improve resolution in a computationally efficient manner, thereby achieving numerical convergence of a subset of key CGM gas properties and observables.Comment: Submitted to MNRAS. Comments welcom

The cosmic web in Lyman-alpha emission

We develop a comprehensive theoretical model for Lyman-alpha emission, from the scale of individual Lyman-alpha emitters (LAEs) to Lyman-alpha halos (LAHs), Lyman-alpha blobs (LABs), and Lyman-alpha filaments (LAFs) of the diffuse cosmic web itself. To do so, we post-process the high-resolution TNG50 cosmological magnetohydrodynamical simulation with a Monte Carlo radiative transfer method to capture the resonant scattering process of Lyman-alpha photons. We build an emission model incorporating recombinations and collisions in diffuse gas, including radiative effects from nearby AGN, as well as emission sourced by stellar populations. Our treatment includes a physically motivated dust model, which we empirically calibrate to the observed LAE luminosity function. We then focus on the observability, and physical origin, of the $z=2$ Lyman-alpha cosmic web, studying the dominant emission mechanisms and spatial origins. We find that diffuse Lyman-alpha filaments are, in fact, illuminated by photons which originate, not from the intergalactic medium itself, but from within galaxies and their gaseous halos. In our model, this emission is primarily sourced by intermediate mass halos ($10^{10} - 10^{11}\,$M$_{\odot}$), principally due to collisional excitations in their circumgalactic media as well as central, young stellar populations. Observationally, we make predictions for the abundance, area, linear size, and embedded halo/emitter populations within filaments. Adopting an isophotal surface brightness threshold of $10^{-20}\,$erg$\,$s$^{-1}\,$cm$^{-2}\,$arcsec$^{-2}$, we predict a volume abundance of Lyman-alpha filaments of ${\sim}10^{-3}$ cMpc$^{-3}\,$ for lengths above $400\,$pkpc. Given sufficiently large survey footprints, detection of the Lyman-alpha cosmic web is within reach of modern integral field spectrographs, including MUSE, VIRUS, and KCWI.Comment: Submitted to MNRA

scida: scalable analysis for scientific big data

scida is a Python package for reading and analyzing large scientific data sets with support for various cosmological and galaxy formation simulations out-of-the-box. Data access is provided through a hierarchical dictionary-like data structure after a simple load() function. Using the dask library for scalable, parallel and out-of-core computation, all computation requests from a user session are first collected in a task graph. Arbitrary custom analysis, as well as all available dask (array) operations, can be performed. The subsequent computation is executed only upon request, on a target resource (e.g. a HPC cluster).Comment: recommended for acceptance in the Journal of Open Source Software; open-source development at https://github.com/cbyrohl/scid

Following the flow: tracer particles in astrophysical fluid simulations

We present two numerical schemes for passive tracer particles in the hydrodynamical moving-mesh code AREPO, and compare their performance for various problems, from simple setups to cosmological simulations. The purpose of tracer particles is to allow the flow to be followed in a Lagrangian way, tracing the evolution of the fluid with time, and allowing the thermodynamical history of individual fluid parcels to be recorded. We find that the commonly-used `velocity field tracers', which are advected using the fluid velocity field, do not in general follow the mass flow correctly, and explain why this is the case. This method can result in orders-of-magnitude biases in simulations of driven turbulence and in cosmological simulations, rendering the velocity field tracers inappropriate for following these flows. We then discuss a novel implementation of `Monte Carlo tracers', which are moved along with fluid cells, and are exchanged probabilistically between them following the mass flux. This method reproduces the mass distribution of the fluid correctly. The main limitation of this approach is that it is more diffusive than the fluid itself. Nonetheless, we show that this novel approach is more reliable than what has been employed previously and demonstrate that it is appropriate for following hydrodynamical flows in mesh-based codes. The Monte Carlo tracers can also naturally be transferred between fluid cells and other types of particles, such as stellar particles, so that the mass flow in cosmological simulations can be followed in its entirety.Comment: Accepted for publication in MNRAS, minor updates to match accepted version. 19 pages, 14 figure

Zooming in on accretion - I. The structure of halo gas

We study the properties of gas in and around 10^12 solar mass halos at z=2 using a suite of high-resolution cosmological hydrodynamic 'zoom' simulations. We quantify the thermal and dynamical structure of these gaseous reservoirs in terms of their mean radial distributions and angular variability along different sightlines. With each halo simulated at three levels of increasing resolution, the highest reaching a baryon mass resolution of ~10,000 solar masses, we study the interaction of filamentary inflow and the quasi-static hot halo atmosphere. We highlight the discrepancy between the spatial resolution available in the halo gas as opposed to within the galaxy itself, and find that stream morphologies become increasingly complex at higher resolution, with large coherent flows revealing density and temperature structure at progressively smaller scales. Moreover, multiple gas components co-exist at the same radius within the halo, making radially averaged analyses misleading. This is particularly true where the hot, quasi-static, high entropy halo atmosphere interacts with cold, rapidly inflowing, low entropy accretion. We investigate the process of gas virialization and identify different regimes for the heating of gas as it accretes from the intergalactic medium. Haloes at this mass have a well-defined virial shock, associated with a sharp jump in temperature and entropy at ~1.25 r_vir. The presence, radius, and radial width of this boundary feature, however, vary not only from halo to halo, but also as a function of angular direction, covering roughly ~85% of the 4pi sphere. Our findings are relevant for the proper interpretation of observations pertaining to the circumgalactic medium, including evidence for large amounts of cold gas surrounding massive haloes at intermediate redshifts.Comment: High-res PDF and simulation movies available at http://www.cfa.harvard.edu/~dnelson/#research (MNRAS submitted, comments welcome

Moving mesh cosmology: tracing cosmological gas accretion

We investigate the nature of gas accretion onto haloes and galaxies at z=2 using cosmological hydrodynamic simulations run with the moving mesh code AREPO. Implementing a Monte Carlo tracer particle scheme to determine the origin and thermodynamic history of accreting gas, we make quantitative comparisons to an otherwise identical simulation run with the smoothed particle hydrodynamics (SPH) code GADGET-3. Contrasting these two numerical approaches, we find significant physical differences in the thermodynamic history of accreted gas in haloes above 10^10.5 solar masses. In agreement with previous work, GADGET simulations show a cold fraction near unity for galaxies forming in massive haloes, implying that only a small percentage of accreted gas heats to an appreciable fraction of the virial temperature during accretion. The same galaxies in AREPO show a much lower cold fraction, <20% in haloes above 10^11 solar masses. This results from a hot gas accretion rate which, at this same halo mass, is an order of magnitude larger than with GADGET, while the cold accretion rate is also lower. These discrepancies increase for more massive systems, and we explain both as due to numerical inaccuracies in the standard formulation of SPH. We also observe that the relatively sharp transition from cold to hot mode dominated accretion, at a halo mass of ~10^11, is a consequence of comparing past gas temperatures to a constant threshold value independent of virial temperature. Examining the spatial distribution of accreting gas, we find that gas filaments in GADGET tend to remain collimated and flow coherently to small radii, or artificially fragment and form a large number of purely numerical "blobs". Similar gas streams in AREPO show increased heating and disruption at 0.25-0.5 virial radii and contribute to the hot gas accretion rate in a manner distinct from classical cooling flows.Comment: 21 pages, 12 figures. MNRAS accepted (in press). High-resolution images can be found at http://www.cfa.harvard.edu/itc/research/movingmeshcosmology

The Circumgalactic Medium of Milky Way-like Galaxies in the TNG50 Simulation -- II: Cold, Dense Gas Clouds and High-Velocity Cloud Analogs

We use the TNG50 simulation of the IllustrisTNG project to study cold, dense clouds of gas in the circumgalactic media (CGM) of Milky Way-like galaxies. We find that their CGM is typically filled with of order one hundred (thousand) reasonably (marginally) resolved clouds, possible analogs of high-velocity clouds (HVCs). There is a large variation in cloud abundance from galaxy to galaxy, and the physical properties of clouds that we explore -- mass, size, metallicity, pressure, and kinematics -- are also diverse. We quantify the distributions of cloud properties and cloud-background contrasts, providing cosmological inputs for idealized simulations. Clouds characteristically have sub-solar metallicities, diverse shapes, small overdensities ($\chi = n_{\rm cold} / n_{\rm hot} \lesssim 10$), are mostly inflowing, and have sub-virial rotation. At TNG50 resolution, resolved clouds have median masses of $\sim 10^6\,\rm{M_\odot}$ and sizes of $\sim 10$ kpc. Larger clouds are well converged numerically, while the abundance of the smallest clouds increases with resolution, as expected. In TNG50 MW-like haloes, clouds are slightly (severely) under-pressurised relative to their surroundings with respect to total (thermal) pressure, implying that magnetic fields may be important. Clouds are not distributed uniformly throughout the CGM, but are clustered around other clouds, often near baryon-rich satellite galaxies. This suggests that at least some clouds originate from satellites, via direct ram-pressure stripping or otherwise. Finally, we compare with observations of intermediate and high velocity clouds from the real Milky Way halo. TNG50 shows a similar cloud velocity distribution as observations, and predicts a significant population of currently difficult-to-detect low velocity clouds.Comment: Accepted for publication (MNRAS). Part of a set of papers based on TNG50 MW/M31-like galaxies. Additional visuals and data products at www.tng-project.org/ramesh23