Federico Marinacci

Galactic fountains and the rotation of disc-galaxy coronae

In galaxies like the Milky Way, cold (˜104 K) gas ejected from the disc by stellar activity (the so-called galactic-fountain gas) is expected to interact with the virial-temperature (˜106 K) gas of the corona. The associated transfer of momentum between cold and hot gas has important consequences for the dynamics of both gas phases. We quantify the effects of such an interaction using hydrodynamical simulations of cold clouds travelling through a hot medium at different relative velocities. Our main finding is that there is a velocity threshold between clouds and corona, of about 75 km s-1, below which the hot gas ceases to absorb momentum from the cold clouds. It follows that in a disc galaxy like the Milky Way a static corona would be rapidly accelerated; the corona is expected to rotate and to lag, in the inner regions, by ˜80-120 km s-1 with respect to the cold disc. We also show how the existence of this velocity threshold can explain the observed kinematics of the cold extraplanar gas.

Simulating galaxy formation with black hole driven thermal and kinetic feedback

The inefficiency of star formation in massive elliptical galaxies is widely believed to be caused by the interactions of an active galactic nucleus (AGN) with the surrounding gas. Achieving a sufficiently rapid reddening of moderately massive galaxies without expelling too many baryons has however proven difficult for hydrodynamical simulations of galaxy formation, prompting us to explore a new model for the accretion and feedback effects of supermassive black holes. For high-accretion rates relative to the Eddington limit, we assume that a fraction of the accreted rest mass energy heats the surrounding gas thermally, similar to the ‘quasar mode’ in previous work. For low-accretion rates, we invoke a new, pure kinetic feedback model that imparts momentum to the surrounding gas in a stochastic manner. These two modes of feedback are motivated both by theoretical conjectures for the existence of different types of accretion flows as well as recent observational evidence for the importance of kinetic AGN winds in quenching galaxies. We find that a large fraction of the injected kinetic energy in this mode thermalizes via shocks in the surrounding gas, thereby providing a distributed heating channel. In cosmological simulations, the resulting model produces red, non-star-forming massive elliptical galaxies, and achieves realistic gas fractions, black hole growth histories and thermodynamic profiles in large haloes.

Unveiling the corona of the Milky Way via ram-pressure stripping of dwarf satellites

The spatial segregation between dwarf spheroidal (dSph) and dwarf irregular galaxies in the Local Group has long been regarded as evidence of an interaction with their host galaxies. In this paper, we assume that ram-pressure stripping is the dominant mechanism that removed gas from the dSphs and we use this to derive a lower bound on the density of the corona of the Milky Way at large distances (R similar to 50-90 kpc) from the Galactic Centre. At the same time, we derive an upper bound by demanding that the interstellar medium of the dSphs is in pressure equilibrium with the hot corona. We consider two dwarfs (Sextans and Carina) with well-determined orbits and star formation histories. Our approach introduces several novel features: (i) we use the measured star formation histories of the dwarfs to derive the time at which they last lost their gas and (via a modified version of the Kennicutt-Schmidt relation) their internal gas density at that time; (ii) we use a large suite of 2D hydrodynamical simulations to model the gas stripping; and (iii) we include supernova feedback tied to the gas content. Despite having very different orbits and star formation histories, we find results for the two dSphs that are in excellent agreement with one another. We derive an average particle density of the corona of the Milky Way at R = 50-90 kpc in the range n(cor) = 1.3-3.6 x 10(-4) cm(-3). Including additional constraints from X-ray emission limits and pulsar dispersion measurements (that strengthen our upper bound), we derive Galactic coronal density profiles. Extrapolating these to large radii, we estimate the fraction of baryons (missing baryons) that can exist within the virial radius of the Milky Way. For an isothermal corona (T-cor = 1.8 x 10(6) K), this is small - just 10-20 per cent of the expected missing baryon fraction, assuming a virial mass of 1-2 x 10(12) M-circle dot. Only a hot (T-cor = 3 x 10(6) K) and adiabatic corona can contain all of the Galaxys missing baryons. Models for the Milky Way must explain why its corona is in a hot adiabatic thermal state; or why a large fraction of its baryons lie beyond the virial radius.

Halo mass and assembly history exposed in the faint outskirts: the stellar and dark matter haloes of Illustris galaxies

We use the Illustris Simulations to gain insight into the build-up of the outer, low-surface brightness regions which surround galaxies. We characterize the stellar haloes by means of the logarithmic slope of the spherically-averaged stellar density profiles, alphaSTARS at z=0, and we relate these slopes to the properties of the underlying Dark-Matter (DM) haloes, their central galaxies, and their assembly histories. We analyze a sample of ~5,000 galaxies resolved with more than 5x10^4 particles each, and spanning a variety of morphologies and halo masses (3x10^11 < Mvir < 10^14 Msun). We find a strong trend between stellar halo slope and total halo mass, where more massive objects have shallower stellar haloes than the less massive ones (-5.5 \pm 0.5 < alphaSTARS <-3.5 \pm 0.2 in the studied mass range). At fixed halo mass, we show that disk-like, blue, young, and more massive galaxies are surrounded by significantly steeper stellar haloes than elliptical, red, older, and less massive galaxies. Overall, the stellar density profiles fall off much more steeply than the underlying DM, and no clear trend holds between stellar slope and DM halo concentration. However, DM haloes which formed more recently, or which accreted larger fractions of stellar mass from infalling satellites, exhibit shallower stellar haloes than their older analogs with similar masses, by up to Delta(alphaSTARS) ~ 0.5-0.7. Our findings, combined with the most recent measurements of the strikingly different stellar power-law indexes for M31 and the Milky Way, appear to favour a massive M31, and a Milky Way characterized by a much quieter accretion history over the past 10 Gyrs than its companion.

Vertical disc heating in Milky Way-sized galaxies in a cosmological context

Vertically extended, high velocity dispersion stellar distributions appear to be a ubiquitous feature of disc galaxies, and both internal and external mechanisms have been proposed to be the major driver of their formation. However, it is unclear to what extent each mechanism can generate such a distribution, which is likely to depend on the assembly history of the galaxy. To this end, we perform 16 high-resolution cosmological-zoom simulations of Milky Way-sized galaxies using the state-of-the-art cosmological magnetohydrodynamical code AREPO, and analyse the evolution of the vertical kinematics of the stellar disc in connection with various heating mechanisms. We find that the bar is the dominant heating mechanism in most cases, whereas spiral arms, radial migration and adiabatic heating from mid-plane density growth are all subdominant. The strongest source, though less prevalent than bars, originates from external perturbations from satellites/subhaloes of masses log10(M/M⊙) ≳ 10. However, in many simulations the orbits of newborn star particles become cooler with time, such that they dominate the shape of the age–velocity dispersion relation and overall vertical disc structure unless a strong external perturbation takes place.

Simulating galaxy formation with the IllustrisTNG model

We introduce an updated physical model to simulate the formation and evolution of galaxies in cosmological, large-scale gravity+magnetohydrodynamical simulations with the moving mesh code AREPO. The overall framework builds upon the successes of the Illustris galaxy formation model, and includes prescriptions for star formation, stellar evolution, chemical enrichment, primordial and metal-line cooling of the gas, stellar feedback with galactic outflows, and black hole formation, growth and multi-mode feedback. In this paper we give a comprehensive description of the physical and numerical advances which form the core of the IllustrisTNG (The Next Generation) framework. We focus on the revised implementation of the galactic winds, of which we modify the directionality, velocity, thermal content, and energy scalings, and explore its effects on the galaxy population. As described in earlier works, the model also includes a new black hole driven kinetic feedback at low accretion rates, magnetohydrodynamics, and improvements to the numerical scheme. Using a suite of (25 Mpc

The Auriga Project: the properties and formation mechanisms of disc galaxies across cosmic time

h^{-1}

Baryonic impact on the dark matter distribution in Milky Way-sized galaxies and their satellites

)

MAGNETIC FIELDS IN COSMOLOGICAL SIMULATIONS OF DISK GALAXIES

^3

Diffuse gas properties and stellar metallicities in cosmological simulations of disc galaxy formation

cosmological boxes we assess the outcome of the new model at our fiducial resolution. The presence of a self-consistently amplified magnetic field is shown to have an important impact on the stellar content of