Showing 1–50 of 68 results for author: Scannapieco, C

Cosmological gas accretion history onto the stellar discs of Milky Way-like galaxies in the Auriga simulations -- (II) The inside-out growth of discs

Authors: Federico G. Iza, Sebastián E. Nuza, Cecilia Scannapieco, Robert J. J. Grand, Facundo A. Gómez, Volker Springel, Rüdiger Pakmor, Federico Marinacci, Francesca Fragkoudi

Abstract: We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar di… ▽ More We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar distributions of around 70% of the simulated galaxies exhibit an ``inside-out'' pattern, with older (younger) stellar populations preferentially located in the inner (outer) disc regions. In all cases, we find a very tight correlation between the infall, outflow and net accretion rates, as well as between these three quantities and the star formation rate. This is because the amount of gas which is ultimately available for star formation in each radial ring depends not only on the infall rates, but also on the amount of gas leaving the disc in outflows, which directly relates to the local star formation level. Therefore, any of these rates can be used to identify galaxies with inside-out growth. For these galaxies, the correlation between the dominant times of accretion/star formation and disc radius is well fitted by a linear function. We also find that, when averaged over galaxies with formation histories similar to the Milky Way, the simulated accretion rates show a similar evolution (both temporally- and radially-integrated) to the usual accretion prescriptions used in chemical evolution models, although some major differences arise at early times and in the inner disc regions. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 19 pages, 13 figures

Unveiling the origins of galactic bars: insights from barred and unbarred galaxies

Authors: Paula D. López, Cecilia Scannapieco, Sofía A. Cora, Ignacio D. Gargiulo

Abstract: A significant fraction of local galaxies exhibit stellar bars, non-axisymmetric structures composed of stars, gas, and dust. Identifying key differences between the properties of barred and unbarred galaxies can uncover clues about the conditions for triggering bar formation. We explore the early stages of bar formation in a small sample of disc barred galaxies extracted from the TNG50 cosmologica… ▽ More A significant fraction of local galaxies exhibit stellar bars, non-axisymmetric structures composed of stars, gas, and dust. Identifying key differences between the properties of barred and unbarred galaxies can uncover clues about the conditions for triggering bar formation. We explore the early stages of bar formation in a small sample of disc barred galaxies extracted from the TNG50 cosmological simulation, and compare their properties to those of unbarred galaxies. According to our results, the most important difference between barred and unbarred galaxies is that the former have systematically higher fractions of stellar to dark matter mass in their inner regions, from very early stages and prior to the formation of the bars. They harbour high initial gas content, fostering increased star formation rates and leading to a central mass concentration that grows faster over time compared to unbarred galaxies. Examining the evolution of the halo spin within 10 ckpc reveals that barred galaxies have higher angular momentum transfer from the disc to the halo. Curiously, both barred and unbarred galaxies share similar initial low values of the halo spin, consistent with those proposed in the literature for bar formation. Furthermore, we evaluate existing stability criteria to capture the complexity of the process, and investigate the effects of mergers, flybys, and environment as possible drivers of bar formation. We find no clear link between mergers and disc instabilities resulting in the formation of bars, even though some of the simulated barred galaxies might have been influenced by these events. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 20 pages, 20 figures, 2 tables

Journal ref: MNRAS, Volume 529, Issue 2, April 2024, Pages from 979 to 998

A model for molecular hydrogen-dependent star formation in simulations of galaxy evolution

Authors: Ezequiel Lozano, Cecilia Scannapieco, Sebastian E. Nuza

Abstract: Star formation, together with the associated chemical and energy feedback, is one of the most important processes in galaxy evolution. The star formation activity in galaxies defines and affects many of their fundamental properties, such as stellar mass, morphology and chemical enrichment levels. Simple models for star formation in cosmological hydrodynamical simulations have shown to be successfu… ▽ More Star formation, together with the associated chemical and energy feedback, is one of the most important processes in galaxy evolution. The star formation activity in galaxies defines and affects many of their fundamental properties, such as stellar mass, morphology and chemical enrichment levels. Simple models for star formation in cosmological hydrodynamical simulations have shown to be successful in reproducing the star formation rate (SFR) levels and shapes of different types of galaxies. However, with the advent of high-resolution simulations and more detailed observations, more sophisticated star formation models are needed; in particular, to better understand the relation between star formation and the amount of gas in the atomic and molecular phases. In this work, we apply a novel star formation model, recently developed to work in the context of hydrodynamical simulations, to the study of the SFR in Milky Way-mass galaxies. The new implementation describes the formation of molecular hydrogen from atomic material, considering also possible dependencies with the chemical abundance of the gas. This allows to implement various star formation models, where the SFR of a gas cloud is determined by the atomic and/or molecular gas phases, and to compare their predictions to recent observational results. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: 3 pages, 3 figures. To appear in the 64nd Bulletin of the Argentine Astronomical Society

Kinematics of the Local Group gas and galaxies in the Hestia simulations

Authors: Luis Biaus, Sebastián E. Nuza, Cecilia Scannapieco

Abstract: The Local Group (LG) consists of two giant spiral galaxies, the Milky Way (MW) and Andromeda (M31), and several smaller galaxies. The MW and M31 are approaching each other at a radial velocity of about $-109\,$km\,s$^{-1}$. Observational evidence suggests that there is an overall infalling motion of gas and galaxies in the LG, dominated by the dynamics of its two main members. From our perspective… ▽ More The Local Group (LG) consists of two giant spiral galaxies, the Milky Way (MW) and Andromeda (M31), and several smaller galaxies. The MW and M31 are approaching each other at a radial velocity of about $-109\,$km\,s$^{-1}$. Observational evidence suggests that there is an overall infalling motion of gas and galaxies in the LG, dominated by the dynamics of its two main members. From our perspective, this flow imprints a velocity dipole pattern in the sky when Galactic rotation is removed. We investigate the kinematic properties of gas and galaxies in the LG using a suite of high-resolution simulations performed by the {\sc Hestia} (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions. We build sky maps from the local, Galactic and LG standard of rest reference frames. Our findings show that the establishment of a radial velocity dipole near the preferred barycentre direction is a natural outcome of simulation kinematics for material \textit{outside} the MW virial radius after removing galaxy rotation when the relative radial velocity of MW and M31 is similar to the observed value. These results favour a scenario where gas and galaxies stream towards the LG barycentre, producing the observed velocity dipole. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 3 pages, 3 figures. To appear in the 64th Bulletin of the Argentine Astronomical Society. arXiv admin note: substantial text overlap with arXiv:2210.15589

The temporal evolution of gas accretion onto the discs of simulated Milky Way-mass galaxies

Authors: F. G. Iza, S. E. Nuza, C. Scannapieco

Abstract: In the standard model of structure formation, galaxies form in the centre of dark matter haloes that develop as a result of inhomogeneities in the primordial mass distribution of the Universe. Afterwards, galaxies grow by means of continuous accretion of gaseous material stemming from the intergalactic medium, both in diffuse form and through collisions with other systems. After an initial period… ▽ More In the standard model of structure formation, galaxies form in the centre of dark matter haloes that develop as a result of inhomogeneities in the primordial mass distribution of the Universe. Afterwards, galaxies grow by means of continuous accretion of gaseous material stemming from the intergalactic medium, both in diffuse form and through collisions with other systems. After an initial period of violent growth, the gas settles into a rotationally-supported structure where stars are born, giving birth to the stellar disc. The accretion of gaseous material onto the disc plays a fundamental role in its evolution as it can change its dynamical and morphological properties, generating gas flows within the disc. In this work, we use 30 galaxies from the Auriga Project, a set of cosmological magnetohydrodynamical simulations of disc galaxies, to study the temporal dependence of the gas accretion rates, focusing on the inflowing and outflowing fluxes. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 3 pages, 2 figures. To appear in the 64th Bulletin of the Argentine Astronomical Society

Kinematics of the Local Group gas and galaxies in the Hestia simulations

Authors: Luis Biaus, Sebastián E. Nuza, Philipp Richter, Martin Sparre, Cecilia Scannapieco, Mitali Damle, Jenny G. Sorce, Robert J. J. Grand, Elmo Tempel, Noam I. Libeskind, Maan H. Hani

Abstract: We investigate the kinematic properties of gas and galaxies in the Local Group (LG) using high-resolution simulations performed by the {\sc Hestia} (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions consisting of two main spiral galaxies of $\sim 10^{12}$~M$_\odot$, their satellites and minor… ▽ More We investigate the kinematic properties of gas and galaxies in the Local Group (LG) using high-resolution simulations performed by the {\sc Hestia} (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions consisting of two main spiral galaxies of $\sim 10^{12}$~M$_\odot$, their satellites and minor isolated galaxies, all sharing the same large-scale motion within a volume of a few Mpc. We characterise the gas and galaxy kinematics within the simulated LGs, from the perspective of the Sun, to compare with observed trends from recent HST/COS absorption-line observations and LG galaxy data. To analyse the velocity pattern of LG gas and galaxies seen in the observational data, we build sky maps from the local standard of rest, and the galactic and local group barycentre frames. Our findings show that the establishment of a radial velocity dipole at low/high latitudes, near the preferred barycentre direction, is a natural outcome of simulation kinematics for material {\it outside} the Milky Way virial radius after removing galaxy rotation when the two main LG galaxies are approaching. Our results favour a scenario where gas and galaxies stream towards the LG barycentre producing a velocity dipole resembling observations. While our study shows in a qualitative way the global matter kinematics in the LG as part of its on-going assembly, quantitative estimates of gas-flow rates and physical conditions of the LG gas have to await a more detailed modeling of the ionization conditions, which will be presented in a follow-up paper. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

Cosmological gas accretion history onto the stellar discs of Milky Way-like galaxies in the Auriga simulations -- (I) Temporal dependency

Authors: Federico G. Iza, Cecilia Scannapieco, Sebastián E. Nuza, Robert J. J. Grand, Facundo A. Gómez, Volker Springel, Rüdiger Pakmor, Federico Marinacci

Abstract: We use the 30 simulations of the Auriga Project to estimate the temporal dependency of the inflow, outflow and net accretion rates onto the discs of Milky Way-like galaxies. The net accretion rates are found to be similar for all galaxies at early times, increasing rapidly up to $\sim 10~\mathrm{M}_\odot \, \mathrm{yr}^{-1}$. After $\sim 6~\mathrm{Gyr}$ of evolution, however, the net accretion rat… ▽ More We use the 30 simulations of the Auriga Project to estimate the temporal dependency of the inflow, outflow and net accretion rates onto the discs of Milky Way-like galaxies. The net accretion rates are found to be similar for all galaxies at early times, increasing rapidly up to $\sim 10~\mathrm{M}_\odot \, \mathrm{yr}^{-1}$. After $\sim 6~\mathrm{Gyr}$ of evolution, however, the net accretion rates are diverse: in most galaxies, these exhibit an exponential-like decay, but some systems instead present increasing or approximately constant levels up to the present time. An exponential fit to the net accretion rates averaged over the MW analogues yields typical decay time-scale of $7.2~\mathrm{Gyr}$. The analysis of the time-evolution of the inflow and outflow rates, and their relation to the star formation rate (SFR) in the discs, confirms the close connection between these quantities. First, the inflow$/$outflow ratio stays approximately constant, with typical values of $\dot{M}_\mathrm{out}/ \dot{M}_\mathrm{in} \sim 0.75$, indicating that the gas mass involved in outflows is of the order of 25% lower compared to that involved in inflows. A similar behaviour is found for the SFR$/$inflow rate ratio, with typical values between 0.1 and 0.3, and for the outflow rate$/$SFR which varies in the range $3.5$--$5.5$. Our results show that continuous inflow is key to the SFR levels in disc galaxies, and that the star formation activity and the subsequent feedback in the discs is able to produce mass-loaded galaxy winds in the disc-halo interface. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: 21 pages, 15 figures

Journal ref: Mon. Not. Roy. Astron. Soc., 517, 832-852, 2022

The neutron-capture and alpha-elements abundance ratios scatter in old stellar populations. Cosmological simulations of the stellar halo

Authors: Cecilia Scannapieco, Gabriele Cescutti, Cristina Chiappini

Abstract: We investigate the origin of the abundance ratios and scatter of the neutron-capture elements Sr, Ba and Eu in the stellar halo of a Milky Way-mass galaxy formed in a hydrodynamical cosmological simulation, and compare them with those of $α$-elements. For this, we implement a novel treatment for chemical enrichment of Type II supernovae which considers the effects of the rotation of massive stars… ▽ More We investigate the origin of the abundance ratios and scatter of the neutron-capture elements Sr, Ba and Eu in the stellar halo of a Milky Way-mass galaxy formed in a hydrodynamical cosmological simulation, and compare them with those of $α$-elements. For this, we implement a novel treatment for chemical enrichment of Type II supernovae which considers the effects of the rotation of massive stars on the chemical yields and differential enrichment according to the life-times of progenitor stars. We find that differential enrichment has a significant impact on the early enrichment of the interstellar medium which is translated into broader element ratio distributions, particularly in the case of the oldest, most metal-poor stars. We find that the [element/Fe] ratios of the $α-$elements O, Mg and Si have systematically lower scatter compared to the neutron-capture elements ratios Sr, Ba and Eu at [Fe/H]$<-2$, which is $\sim 0.1-0.4$ dex for the former and between $\sim 0.5$ and $1$ dex for the latter. The different scatter levels found for the neutron-capture and $α$-elements is consistent with observations of old stars in the Milky Way. Our model also predicts a high scatter for the [Sr/Ba] ratio, which results from the treatment of the fast-rotating stars and the dependence of the chemical yields on the metallicity, mass and rotational velocities. Such chemical patterns appear naturally if the different ejection times associated to stars of different mass are properly described, without the need to invoke for additional mixing mechanisms or a distinct treatment of the $α-$ and neutron-capture elements. △ Less

Submitted 9 September, 2022; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: MNRAS accepted

Gas accretion onto the disc of a simulated Milky Way-mass galaxy

Authors: F. G. Iza, S. E. Nuza, C. Scannapieco

Abstract: In the standard paradigm of galaxy formation and evolution, the baryonic component of galaxies forms from the collapse and condensation of gas within dark matter haloes, and later grows from continuous accretion of gaseous mass, both in diffuse form and in mergers with other systems. After a first period of rapid and violent halo growth, the gas settles into a rotationally-supported structure, eve… ▽ More In the standard paradigm of galaxy formation and evolution, the baryonic component of galaxies forms from the collapse and condensation of gas within dark matter haloes, and later grows from continuous accretion of gaseous mass, both in diffuse form and in mergers with other systems. After a first period of rapid and violent halo growth, the gas settles into a rotationally-supported structure, eventually giving rise to the formation of a stellar disc. Stars evolve and return chemically-processed gas and energy to the interstellar medium, mainly through Type II supernova explosions. In the disc region, the cosmological accretion of gas combines with the outflows resulting from supernovae, affecting the hydrodynamical and structural properties of the disc and producing gas flows in the vertical and radial directions. In this work, we use a simulation of the Auriga Project, a suite of magneto-hydrodynamical, zoom-in cosmological simulations of Milky Way-like galaxies, to study the temporal and radial dependencies of gas accretion onto the disc. We also investigate the disc evolution, focusing on the inside-out disc formation scenario, which is one of the fundamental hypotheses of chemical evolution models of the Galaxy. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: 3 pages, 3 figures. To appear in the 62nd Bulletin of the Argentine Astronomical Society

Chemical abundance patterns in Local Group galaxies within cosmological simulations

Authors: L. Biaus, C. Scannapieco, S. E. Nuza

Abstract: In the context of the concordance cosmology, structure formation in the Universe is the result of the amplification, by gravitational effects, of small perturbations in the primeval density field. This results in the formation of structures known as dark matter haloes, where gas collapses and forms stars, giving birth to galaxies. Numerical simulations are an important tool in the theoretical stud… ▽ More In the context of the concordance cosmology, structure formation in the Universe is the result of the amplification, by gravitational effects, of small perturbations in the primeval density field. This results in the formation of structures known as dark matter haloes, where gas collapses and forms stars, giving birth to galaxies. Numerical simulations are an important tool in the theoretical study of galaxy formation and evolution. In the present work, we describe the implementation of a chemical enrichment model in a state-of-the-art cosmological simulation of the Local Group. The simulation includes sub-grid models for the most relevant physical processes. We analyze the chemical and morphological evolution of two galaxies with virial masses similar to that of our Milky Way. For each of the stellar components (disc, bulge and halo), we establish links between their formation history and their chemical evolution. We find that $α$-element (O, Mg, Si) enrichment happens at early stages of evolution, as their main producers are short-lived stars which end their lives as type II supernova explosions. There is also a gradual contamination with the rest of the elements as type Ia supernovae and winds of stars in the asymptotic giant branch occur. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 3 pages, 4 figures. To appear in the 62nd Bulletin of the Argentine Astronomical Society

A model for star formation in cosmological simulations of galaxy formation

Authors: Ezequiel Lozano, Cecilia Scannapieco, Sebastian E. Nuza

Abstract: We present a new model to describe the star formation process in galaxies, which includes the description of the different gas phases -- molecular, atomic, and ionized -- together with its metal content. The model, which will be coupled to cosmological simulations of galaxy formation, will be used to investigate the relation between the star formation rate (SFR) and the formation of molecular hydr… ▽ More We present a new model to describe the star formation process in galaxies, which includes the description of the different gas phases -- molecular, atomic, and ionized -- together with its metal content. The model, which will be coupled to cosmological simulations of galaxy formation, will be used to investigate the relation between the star formation rate (SFR) and the formation of molecular hydrogen. The model follows the time evolution of the molecular, atomic and ionized phases in a gas cloud and estimates the amount of stellar mass formed, by solving a set of five coupled differential equations. As expected, we find a positive, strong correlation between the molecular fraction and the initial gas density, which manifests in a positive correlation between the initial gas density and the SFR of the cloud. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 3 pages, 2 figures. To appear in the 62nd Bulletin of the Argentine Astronomical Society

Do galaxies die? Different views from simulations and observations in the Local Universe

Authors: Pablo Corcho-Caballero, Yago Ascasibar, Cecilia Scannapieco

Abstract: For years, the extragalactic community has divided galaxies in two distinct populations. One of them, featuring blue colours, is actively forming stars, while the other is made up of "red-and-dead" objects with negligible star formation. Yet, are these galaxies really dead? Here we would like to highlight that, as previously reported by several independent groups, state-of-the-art cosmological num… ▽ More For years, the extragalactic community has divided galaxies in two distinct populations. One of them, featuring blue colours, is actively forming stars, while the other is made up of "red-and-dead" objects with negligible star formation. Yet, are these galaxies really dead? Here we would like to highlight that, as previously reported by several independent groups, state-of-the-art cosmological numerical simulations predict the existence of a large number of quenched galaxies that have not formed any star over the last few Gyr. In contrast, observational measurements of large galaxy samples in the nearby Universe suggest that even the most passive systems still form stars at some residual level close to $sSFR\sim10^{-12}~\text{yr}^{-1}$. Unfortunately, extremely low star formation poses a challenge for both approaches. We conclude that, at present, the fraction of truly dead galaxies is still an important open question that must be addressed in order to understand galaxy formation and evolution. △ Less

Submitted 29 June, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS. 9 pages, 8 figures

Yule-Simpson's paradox in Galactic Archaeology

Authors: I. Minchev, G. Matijevic, D. W. Hogg, G. Guiglion, M. Steinmetz, F. Anders, C. Chiappini, M. Martig, A. Queiroz, C. Scannapieco

Abstract: Simpson's paradox, or Yule-Simpson effect, arises when a trend appears in different subsets of data but disappears or reverses when these subsets are combined. We describe here seven cases of this phenomenon for chemo-kinematical relations believed to constrain the Milky Way disk formation and evolution. We show that interpreting trends in relations, such as the radial and vertical chemical abunda… ▽ More Simpson's paradox, or Yule-Simpson effect, arises when a trend appears in different subsets of data but disappears or reverses when these subsets are combined. We describe here seven cases of this phenomenon for chemo-kinematical relations believed to constrain the Milky Way disk formation and evolution. We show that interpreting trends in relations, such as the radial and vertical chemical abundance gradients, the age-metallicity relation, and the metallicity-rotational velocity relation (MVR), can lead to conflicting conclusions about the Galaxy past if analyses marginalize over stellar age and/or birth radius. It is demonstrated that the MVR in RAVE giants is consistent with being always strongly negative, when narrow bins of [Mg/Fe] are considered. This is directly related to the negative radial metallicity gradients of stars grouped by common age (mono-age populations) due to the inside out disk formation. The effect of the asymmetric drift can then give rise to a positive MVR trend in high-[alpha/Fe] stars, with a slope dependent on a given survey's selection function and observational uncertainties. We also study the variation of lithium abundance, A(Li), with [Fe/H] of AMBRE:HARPS dwarfs. A strong reversal in the positive A(Li)-[Fe/H] trend of the total sample is found for mono-age populations, flattening for younger groups of stars. Dissecting by birth radius shows strengthening in the positive A(Li)-[Fe/H] trend, shifting to higher [Fe/H] with decreasing birth radius; these observational results suggest new constraints on chemical evolution models. This work highlights the necessity for precise age estimates for large stellar samples covering wide spatial regions. △ Less

Submitted 1 May, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 13 p., 7 Fig., Accepted to MNRAS

Precision Pollution - The effects of enrichment yields and timing on galactic chemical evolution

Authors: Pierre-Antoine Poulhazan, Cecilia Scannapieco, Peter Creasey

Abstract: We present an update to the chemical enrichment component of the smoothed-particle hydrodynamics model for galaxy formation presented in Scannapieco et al. (2005) in order to address the needs of modelling galactic chemical evolution in realistic cosmological environments. Attribution of the galaxy-scale abundance patterns to individual enrichment mechanisms such as the winds from asymptotic giant… ▽ More We present an update to the chemical enrichment component of the smoothed-particle hydrodynamics model for galaxy formation presented in Scannapieco et al. (2005) in order to address the needs of modelling galactic chemical evolution in realistic cosmological environments. Attribution of the galaxy-scale abundance patterns to individual enrichment mechanisms such as the winds from asymptotic giant branch (AGB) stars or the presence of a prompt fraction of Type Ia supernovae is complicated by the interaction between them and gas cooling, subsequent star formation and gas ejection. In this work we address the resulting degeneracies by extending our implementation to a suite of mechanisms that encompasses different IMFs, models for yields from the aforementioned stars, models for the prompt component of the delay-time-distribution (DTDs) for Type Ia SNe and metallicity-dependent gas cooling rates, and then applying these to both isolated initial conditions and cosmological hydrodynamical zoom simulations. We find DTDs with a large prompt fraction (such as the bimodal and power-law models) have, at z=0, similar abundance patterns compared to the low-prompt component time distributions (uniform or wide Gaussian models). However, some differences appear, such as the former having systematically higher [X/Fe] ratios and narrower [O/Fe] distributions compared to the latter, and a distinct evolution of the [Fe/H] abundance. △ Less

Submitted 31 July, 2018; originally announced August 2018.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Gas accretion in Milky Way-like galaxies: temporal and radial dependencies

Authors: Sebastian E. Nuza, Cecilia Scannapieco, Cristina Chiappini, Thiago C. Junqueira, Ivan Minchev, Marie Martig

Abstract: One of the fundamental assumptions of chemical evolution models (CEMs) of the Milky Way (MW) and other spirals is that higher gas accretion rates are expected in the past, and in the inner regions of the Galaxy. This leads to the so-called `inside-out disc formation scenario'. Yet, these are probably the most unconstrained inputs of such models. In the present paper, we aim at investigating these… ▽ More One of the fundamental assumptions of chemical evolution models (CEMs) of the Milky Way (MW) and other spirals is that higher gas accretion rates are expected in the past, and in the inner regions of the Galaxy. This leads to the so-called `inside-out disc formation scenario'. Yet, these are probably the most unconstrained inputs of such models. In the present paper, we aim at investigating these main assumptions by studying how gas is accreted in four simulated MW-like galaxies assembled within the $Λ$CDM scenario. The galaxies were obtained using two different simulation techniques, cosmological setups and initial conditions. Two of them are MW candidates corresponding to the chemodynamical model of Minchev et al. (2013, 2014) (known as MCM) and the Local Group cosmological simulation of Nuza et al. (2014). We investigate vertical and radial gas accretion on to galaxy discs as a function of cosmic time and disc radius. We find that accretion in the MW-like galaxies seem to happen in two distinct phases, namely: an early, more violent period; followed by a subsequent, slowly declining phase. Our simulations seem to give support to the assumption that the amount of gas incorporated into the MW disc exponentially decreases with time, leading to current net accretion rates of $0.6-1\,$M$_\odot\,$yr$^{-1}$. In particular, accretion timescales on to the simulated thin-disc-like structures are within $\sim5-7\,$Gyr, consistent with expectations from CEMs. Moreover, our simulated MW discs are assembled from the inside-out with gas in the inner disc regions accreted in shorter timescales than in external ones, in qualitative agreement with CEMs of the Galaxy. However, this type of growth is not general to all galaxies and it is intimately linked to their particular merger and gas accretion history. △ Less

Submitted 11 December, 2018; v1 submitted 16 May, 2018; originally announced May 2018.

Comments: Replaced to match published version. References and table added; 20 pages, 13 figures, 2 tables. Gas accretion timescale vs. galactocentric radius for simulated MW candidates shown in Table B1

Journal ref: Mon. Not. Roy. Astron. Soc., 482, 3089-3108, 2019

How does the stellar disk of the Milky Way get its gas?

Authors: Sebastian E. Nuza, Cristina Chiappini, Cecilia Scannapieco, Ivan Minchev, Marie Martig, Thiago C. Junqueira

Abstract: In chemodynamical evolution models it is usually assumed that the Milky Way galaxy forms from the inside-out implying that gas inflows onto the disk decrease with galactocentric distance. Similarly, to reproduce differences between chemical abundances of the thick disk and bulge with respect to those of the thin disk, higher accretion fluxes at early times are postulated. By using a suite of Milky… ▽ More In chemodynamical evolution models it is usually assumed that the Milky Way galaxy forms from the inside-out implying that gas inflows onto the disk decrease with galactocentric distance. Similarly, to reproduce differences between chemical abundances of the thick disk and bulge with respect to those of the thin disk, higher accretion fluxes at early times are postulated. By using a suite of Milky Way-like galaxies extracted from cosmological simulations, we investigate the accretion of gas on the simulated stellar disks during their whole evolution. In general, we find that the picture outlined above holds, although the detailed behavior depends on the assembly history of the Galaxy and the complexities inherent to the physics of galaxy formation. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 4 pages, 2 figures. To appear in the Proceedings of IAU Symposium 334,"Rediscovering our Galaxy", Eds. C. Chiappini, I. Minchev, E. Starkenburg & M. Valentini, Cambridge University Press, Cambridge

The SELGIFS data challenge: generating synthetic observations of CALIFA galaxies from hydrodynamical simulations

Authors: Giovanni Guidi, Javier Casado, Yago Ascasibar, Rubén García-Benito, Lluis Galbany, Patricia Sánchez-Blázquez, Sebastian F. Sánchez, F. Fabián Rosales-Ortega, Cecilia Scannapieco

Abstract: In this work we present a set of synthetic observations that mimic the properties of the Integral Field Spectroscopy (IFS) survey CALIFA, generated using radiative transfer techniques applied to hydrodynamical simulations of galaxies in a cosmological context. The simulated spatially-resolved spectra include stellar and nebular emission, kinematic broadening of the lines, and dust extinction and s… ▽ More In this work we present a set of synthetic observations that mimic the properties of the Integral Field Spectroscopy (IFS) survey CALIFA, generated using radiative transfer techniques applied to hydrodynamical simulations of galaxies in a cosmological context. The simulated spatially-resolved spectra include stellar and nebular emission, kinematic broadening of the lines, and dust extinction and scattering. The results of the radiative transfer simulations have been post-processed to reproduce the main properties of the CALIFA V500 and V1200 observational setups. The data has been further formatted to mimic the CALIFA survey in terms of field of view size, spectral range and sampling. We have included the effect of the spatial and spectral Point Spread Functions affecting CALIFA observations, and added detector noise after characterizing it on a sample of 20 galaxies. The simulated datacubes are suited to be analyzed by the same algorithms used on real IFS data. In order to provide a benchmark to compare the results obtained applying IFS observational techniques to our synthetic datacubes, and test the calibration and accuracy of the analysis tools, we have computed the spatially-resolved properties of the simulations. Hence, we provide maps derived directly from the hydrodynamical snapshots or the noiseless spectra, in a way that is consistent with the values recovered by the observational analysis algorithms. Both the synthetic observations and the product datacubes are public and can be found in the collaboration website http://astro.ft.uam.es/selgifs/data_challenge/. △ Less

Submitted 8 June, 2018; v1 submitted 24 October, 2016; originally announced October 2016.

Comments: 16 pages, 9 figures, 8 tables. Accepted for publication in MNRAS

On the nature of diffuse ionized gas in galaxies -- I The contribution of dust scattering to diffuse line emission

Authors: Y. Ascasibar, G. Guidi, J. Casado, C. Scannapieco, A. I. Díaz

Abstract: In this work, we investigate the contribution of dust scattering to the diffuse H-alpha emission observed in nearby galaxies. As initial conditions for the spatial distribution of HII regions, gas, and dust, we take three Milky Way-like galaxies from state-of-the-art cosmological hydrodynamical simulations that implement different prescriptions for star formation, feedback, and chemical enrichment… ▽ More In this work, we investigate the contribution of dust scattering to the diffuse H-alpha emission observed in nearby galaxies. As initial conditions for the spatial distribution of HII regions, gas, and dust, we take three Milky Way-like galaxies from state-of-the-art cosmological hydrodynamical simulations that implement different prescriptions for star formation, feedback, and chemical enrichment. Radiative transfer has been solved a posteriori, using the publicly-available Monte Carlo code Sunrise to take into account dust absorption and scattering of the H-alpha photons, originating exclusively from the HII regions. No contribution from recombinations in the diffuse ionized gas (DIG) component is explicitly or implicitly included in our model. Our main result is that the flux arising from scattered light is of the order of 1-2 per cent of the H-alpha flux coming directly from the HII regions. Building upon previous studies, we conclude that the DIG contributes lass than 50 per cent of the total H-alpha emission. △ Less

Submitted 26 February, 2016; originally announced February 2016.

Comments: 7 pages, 3 figures, submitted to MNRAS

Physical properties of galaxies: toward a consistent comparison between hydrodynamical simulations and SDSS

Authors: Giovanni Guidi, Cecilia Scannapieco, C. Jakob Walcher, Anna Gallazzi

Abstract: We study the effects of applying observational techniques to derive the properties of simulated galaxies, with the aim of making an unbiased comparison between observations and simulations. For our study, we used fifteen galaxies simulated in a cosmological context using three different feedback and chemical enrichment models, and compared their z=0 properties with data from the Sloan Digital Sky… ▽ More We study the effects of applying observational techniques to derive the properties of simulated galaxies, with the aim of making an unbiased comparison between observations and simulations. For our study, we used fifteen galaxies simulated in a cosmological context using three different feedback and chemical enrichment models, and compared their z=0 properties with data from the Sloan Digital Sky Survey (SDSS). We show that the physical properties obtained directly from the simulations without post-processing can be very different to those obtained mimicking observational techniques. In order to provide simulators a way to reliably compare their galaxies with SDSS data, for each physical property that we studied - colours, magnitudes, gas and stellar metallicities, mean stellar ages and star formation rates - we give scaling relations that can be easily applied to the values extracted from the simulations. These scalings have in general a high correlation, except for the galaxy mean stellar ages and gas oxygen metallicities. Our simulated galaxies are photometrically similar to galaxies in the blue cloud/green valley, but in general they appear older, passive and with lower metal content compared to most of the spirals in SDSS. As a careful assessment of the agreement/disagreement with observations is the primary test of the baryonic physics implemented in hydrodynamical codes, our study shows that considering the observational biases in the derivation of the galaxies' properties is of fundamental importance to decide on the failure/success of a galaxy formation model. △ Less

Submitted 20 July, 2016; v1 submitted 19 February, 2016; originally announced February 2016.

Comments: 19 pages, 18 figures, accepted for publication by MNRAS

Biases and systematics in the observational derivation of galaxy properties: comparing different techniques on synthetic observations of simulated galaxies

Authors: Giovanni Guidi, Cecilia Scannapieco, C. Jakob Walcher

Abstract: We study the sources of biases and systematics in the derivation of galaxy properties of observational studies, focusing on stellar masses, star formation rates, gas/stellar metallicities, stellar ages and magnitudes/colors. We use hydrodynamical cosmological simulations of galaxy formation, for which the real quantities are known, and apply observational techniques to derive the observables. We a… ▽ More We study the sources of biases and systematics in the derivation of galaxy properties of observational studies, focusing on stellar masses, star formation rates, gas/stellar metallicities, stellar ages and magnitudes/colors. We use hydrodynamical cosmological simulations of galaxy formation, for which the real quantities are known, and apply observational techniques to derive the observables. We also make an analysis of biases that are relevant for a proper comparison between simulations and observations. For our study, we post-process the simulation outputs to calculate the galaxies' spectral energy distributions (SEDs) using Stellar Population Synthesis models and also generating the fully-consistent far UV-submillimeter wavelength SEDs with the radiative transfer code SUNRISE. We compared the direct results of simulations with the observationally-derived quantities obtained in various ways, and found that systematic differences in all studied galaxy properties appear, which are caused by: (1) purely observational biases, (2) the use of mass-weighted/luminosity-weighted quantities, with preferential sampling of more massive/luminous regions, (3) the different ways to construct the template of models when a fit to the spectra is performed, and (4) variations due to the use of different calibrations, most notably in the cases of the gas metallicities and star formation rates. Our results show that large differences can appear depending on the technique used to derive galaxy properties. Understanding these differences is of primary importance both for simulators, to allow a better judgement on similarities/differences with observations, and for observers, to allow a proper interpretation of the data. △ Less

Submitted 14 September, 2015; v1 submitted 1 July, 2015; originally announced July 2015.

Comments: 22 pages, 14 figures, 2 tables, MNRAS accepted

The Milky Way and Andromeda galaxies in a constrained hydrodynamical simulation: morphological evolution

Authors: Cecilia Scannapieco, Peter Creasey, Sebastian E. Nuza, Gustavo Yepes, Stefan Gottloeber, Matthias Steinmetz

Abstract: We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-o… ▽ More We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs. △ Less

Submitted 2 March, 2015; originally announced March 2015.

Comments: 12 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 577, A3 (2015)

On the Formation of Galactic Thick Disks

Authors: I. Minchev, M. Martig, D. Streich, C. Scannapieco, R. S. de Jong, M. Steinmetz

Abstract: Recent spectroscopic observations in the Milky Way suggest that the chemically defined thick disk (stars with high [alpha/Fe] ratios and thus old) has a significantly smaller scale-length than the thin disk. This is in apparent contradiction with observations of external edge-on galaxies, where the thin and thick components have comparable scale-lengths. Moreover, while observed disks do not flare… ▽ More Recent spectroscopic observations in the Milky Way suggest that the chemically defined thick disk (stars with high [alpha/Fe] ratios and thus old) has a significantly smaller scale-length than the thin disk. This is in apparent contradiction with observations of external edge-on galaxies, where the thin and thick components have comparable scale-lengths. Moreover, while observed disks do not flare (scale-height does not increase with radius), numerical simulations suggest that disk flaring is unavoidable, resulting from both environmental effects and secular evolution. Here we address these problems by studying two different suites of simulated galactic disks formed in the cosmological context. We show that the scale-heights of coeval populations always increase with radius. However, the total population can be decomposed morphologically into thin and thick disks, which do not flare. We relate this to the disk inside-out formation, where younger populations have increasingly larger scale-lengths and flare at progressively larger radii. In this new picture, thick disks are composed of the imbedded flares of mono-age stellar populations. Assuming that disks form inside out, we predict that morphologically defined thick disks must show a decrease in age (or [alpha/Fe] ratios) with radius and that coeval populations should always flare. This also explains the observed inversion in the metallicity and [alpha/Fe] gradients for stars away from the disk midplane in the Milky Way. The results of this work are directly linked to, and can be seen as evidence of, inside-out disk growth. △ Less

Submitted 6 April, 2015; v1 submitted 23 February, 2015; originally announced February 2015.

Comments: 5 pages, 2 figures, Accepted for publication in ApJ Letters

The Effect of Environment on Milky Way-mass galaxies in a Constrained Simulation of the Local Group

Authors: Peter Creasey, Cecilia Scannapieco, Sebastian E. Nuza, Gustavo Yepes, Stefan Gottloeber, Matthias Steinmetz

Abstract: In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a sufficient volume to compare these galaxies environmental densities ranging from the field to that of the Local Gro… ▽ More In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a sufficient volume to compare these galaxies environmental densities ranging from the field to that of the Local Group (LG). This is particularly relevant as it has been shown that, quite generally, galaxy properties depend intimately upon their environment, most prominently when galaxies in clusters are compared to those in the field. For galaxies in loose groups such as the LG, however, environmental effects have been less clear. We consider the galaxy's environmental density in spheres of 1200 kpc (comoving) and find that whilst environment does not appear to directly affect morphology, there is a positive trend with star formation rates. This enhancement in star formation occurs systematically for galaxies in higher density environments, regardless whether they are part of the LG or in filaments. Our simulations suggest that the richer environment at Mpc-scales may help replenish the star-forming gas, allowing higher specific star formation rates in galaxies such as the MW. △ Less

Submitted 13 January, 2015; originally announced January 2015.

Comments: 6 pages, 4 figures, accepted to ApJL

Low-metallicity stellar halo populations as tracers of dark matter haloes

Authors: Patricia B. Tissera, Cecilia Scannapieco

Abstract: We analyse the density profiles of the stellar halo populations in eight Milky-Way mass galaxies, simulated within the $Λ$-Cold Dark Matter scenario. We find that accreted stars can be well-fitted by an Einasto profile, as well as any subsample defined according to metallicity. We detect a clear correlation between the Einasto fitting parameters of the low-metallicity stellar populations and those… ▽ More We analyse the density profiles of the stellar halo populations in eight Milky-Way mass galaxies, simulated within the $Λ$-Cold Dark Matter scenario. We find that accreted stars can be well-fitted by an Einasto profile, as well as any subsample defined according to metallicity. We detect a clear correlation between the Einasto fitting parameters of the low-metallicity stellar populations and those of the dark matter haloes. The correlations for stars with [Fe/H]$<-3$ allow us to predict the shape of the dark matter profiles within residuals of $\sim 10 $ per cent, in case the contribution from in situ stars remains small. Using Einasto parameters estimated for the stellar halo of the Milky Way and assuming the later formed with significant contributions from accreted low-mass satellite, our simulations predict $α\sim 0.15 $ and $r_2 \sim 15$ kpc for its dark matter profile. These values, combined with observed estimations of the local dark matter density, yield an enclosed dark matter mass at $\sim 8$ kpc in the range $3.9 - 6.7 \times 10^{10}$ M$_{\odot}$, in agreement with recent observational results. These findings suggest that low-metallicity stellar haloes could store relevant information on the DM haloes. Forthcoming observations would help us to further constrain our models and predictions. △ Less

Submitted 22 July, 2014; originally announced July 2014.

Comments: 5 pages,3 figures, MNRAS Letters accepted

The distribution of gas in the Local Group from constrained cosmological simulations: the case for Andromeda and the Milky Way galaxies

Authors: Sebastian E. Nuza, Florencia Parisi, Cecilia Scannapieco, Philipp Richter, Stefan Gottloeber, Matthias Steinmetz

Abstract: We study the gas distribution in the Milky Way and Andromeda using a constrained cosmological simulation of the Local Group (LG) within the context of the CLUES (Constrained Local UniversE Simulations) project. We analyse the properties of gas in the simulated galaxies at $z=0$ for three different phases: `cold', `hot' and HI, and compare our results with observations. The amount of material in th… ▽ More We study the gas distribution in the Milky Way and Andromeda using a constrained cosmological simulation of the Local Group (LG) within the context of the CLUES (Constrained Local UniversE Simulations) project. We analyse the properties of gas in the simulated galaxies at $z=0$ for three different phases: `cold', `hot' and HI, and compare our results with observations. The amount of material in the hot halo ($M_{hot}\approx 4-5\times10^{10}\,$M$_{\odot}$), and the cold ($M_{cold}(r\lesssim10\,$kpc$)\approx10^{8}\,$M$_{\odot}$) and HI ($M_{HI}(r\lesssim50\,$kpc$)\approx 3-4\times10^8\,$M$_{\odot}$) components display a reasonable agreement with observations. We also compute the accretion/ejection rates together with the HI (radial and all-sky) covering fractions. The integrated HI accretion rate within $r=50\,$kpc gives $\sim$$0.2-0.3\,$M$_{\odot}\,$yr$^{-1}$, i.e. close to that obtained from high-velocity clouds in the Milky Way. We find that the global accretion rate is dominated by hot material, although ionized gas with $T\lesssim10^5\,$K can contribute significantly too. The $net$ accretion rates of $all$ material at the virial radii are $6-8\,$M$_{\odot}\,$yr$^{-1}$. At $z=0$, we find a significant gas excess between the two galaxies, as compared to any other direction, resulting from the overlap of their gaseous haloes. In our simulation, the gas excess first occurs at $z\sim1$, as a consequence of the kinematical evolution of the LG. △ Less

Submitted 19 June, 2014; v1 submitted 28 March, 2014; originally announced March 2014.

Comments: Typos corrected. Replaced to match published version. 23 pages, 16 figures, 4 tables

Journal ref: Mon. Not. Roy. Astron. Soc., 441, 2593-2612, 2014

A new stellar chemo-kinematic relation reveals the merger history of the Milky Way disc

Authors: I. Minchev, C. Chiappini, M. Martig, M. Steinmetz, R. S. de Jong, C. Boeche, C. Scannapieco, T. Zwitter, R. F. G. Wyse, J. J. Binney, J. Bland-Hawthorn, O. Bienayme, B. Famaey, K. C. Freeman, B. K. Gibson, E. K. Grebel, G. Gilmore, A. Helmi, G. Kordopatis, Y. S. Lee, U. Munari, J. F. Navarro, Q. A. Parker, A. C. Quillen, W. A. Reid , et al. (5 additional authors not shown)

Abstract: The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemica… ▽ More The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemical and kinematical information, in that it decreases strongly for stars with [Mg/Fe] > 0.4 dex (i.e., those that formed in the first Gyr of the Galaxy's life). These findings can be explained by perturbations from massive mergers in the early Universe, which have affected more strongly the outer parts of the disc, and the subsequent radial migration of stars with cooler kinematics from the inner disc. Similar reversed trends in velocity dispersion are also found for different metallicity subpopulations. Our results suggest that the Milky Way disc merger history can be recovered by relating the observed chemo-kinematic relations to the properties of past merger events. △ Less

Submitted 25 December, 2013; v1 submitted 18 October, 2013; originally announced October 2013.

Comments: 9 pages, 5 figures, Accepted to ApJL. Interpretation section expanded, different Fig.5, typos corrected, references added. Conclusions unchanged

The RAVE survey: the Galactic escape speed and the mass of the Milky Way

Authors: T. Piffl, C. Scannapieco, J. Binney, M. Steinmetz, R. -D. Scholz, M. E. K. Williams, R. S. de Jong, G. Kordopatis, G. Matijevic, O. Bienayme, J. Bland-Hawthorn, C. Boeche, K. Freeman, B. Gibson, G. Gilmore, E. K. Grebel, A. Helmi, U. Munari, J. F. Navarro, Q. Parker, W. A. Reid, G. Seabroke, F. Watson, R. F. G. Wyse, T. Zwitter

Abstract: We construct new estimates on the Galactic escape speed at various Galactocentric radii using the latest data release of the Radial Velocity Experiment (RAVE DR4). Compared to previous studies we have a database larger by a factor of 10 as well as reliable distance estimates for almost all stars. Our analysis is based on the statistical analysis of a rigorously selected sample of 90 high-velocity… ▽ More We construct new estimates on the Galactic escape speed at various Galactocentric radii using the latest data release of the Radial Velocity Experiment (RAVE DR4). Compared to previous studies we have a database larger by a factor of 10 as well as reliable distance estimates for almost all stars. Our analysis is based on the statistical analysis of a rigorously selected sample of 90 high-velocity halo stars from RAVE and a previously published data set. We calibrate and extensively test our method using a suite of cosmological simulations of the formation of Milky Way-sized galaxies. Our best estimate of the local Galactic escape speed, which we define as the minimum speed required to reach three virial radii $R_{340}$, is $533^{+54}_{-41}$ km/s (90% confidence) with an additional 5% systematic uncertainty, where $R_{340}$ is the Galactocentric radius encompassing a mean over-density of 340 times the critical density for closure in the Universe. From the escape speed we further derive estimates of the mass of the Galaxy using a simple mass model with two options for the mass profile of the dark matter halo: an unaltered and an adiabatically contracted Navarro, Frenk & White (NFW) sphere. If we fix the local circular velocity the latter profile yields a significantly higher mass than the un-contracted halo, but if we instead use the statistics on halo concentration parameters in large cosmological simulations as a constraint we find very similar masses for both models. Our best estimate for $M_{340}$, the mass interior to $R_{340}$ (dark matter and baryons), is $1.3^{+0.4}_{-0.3} \times 10^{12}$ M$_\odot$ (corresponding to $M_{200} = 1.6^{+0.5}_{-0.4} \times 10^{12}$ M$_\odot$). This estimate is in good agreement with recently published independent mass estimates based on the kinematics of more distant halo stars and the satellite galaxy Leo I. △ Less

Submitted 13 November, 2013; v1 submitted 17 September, 2013; originally announced September 2013.

Comments: 16 pages, 15 figures; accepted for publication in Astronomy & Astrophysics

Stellar haloes in Milky-Way mass galaxies: From the inner to the outer haloes

Authors: P. Tissera, T. Beers, D. Carollo, C. Scannapieco

Abstract: We present a comprehensive study of the chemical properties of the stellar haloes of Milky-Way mass galaxies, analysing the transition between the inner to the outer haloes. We find the transition radius between the relative dominance of the inner-halo and outer-halo stellar populations to be ~15-20 kpc for most of our haloes, similar to that inferred for the Milky Way from recent observations. Wh… ▽ More We present a comprehensive study of the chemical properties of the stellar haloes of Milky-Way mass galaxies, analysing the transition between the inner to the outer haloes. We find the transition radius between the relative dominance of the inner-halo and outer-halo stellar populations to be ~15-20 kpc for most of our haloes, similar to that inferred for the Milky Way from recent observations. While the number density of stars in the simulated inner-halo populations decreases rapidly with distance, the outer-halo populations contribute about 20-40 per cent in the fiducial solar neighborhood, in particular at the lowest metallicities. We have determined [Fe/H] profiles for our simulated haloes; they exhibit flat or mild gradients, in the range [-0.002, -0.01 ] dex/kpc. The metallicity distribution functions exhibit different features, reflecting the different assembly history of the individual stellar haloes. We find that stellar haloes formed with larger contributions from massive subgalactic systems have steeper metallicity gradients. Very metal-poor stars are mainly contributed to the halo systems by lower-mass satellites. There is a clear trend among the predicted metallicity distribution functions that a higher fraction of low-metallicity stars are found with increasing radius. These properties are consistent with the range of behaviours observed for stellar haloes of nearby galaxies. △ Less

Submitted 23 January, 2014; v1 submitted 13 September, 2013; originally announced September 2013.

Comments: 11 pages, 6 figures. Accepted MNRAS. Revised version after referee's comments

Towards a more realistic population of bright spiral galaxies in cosmological simulations

Authors: Michael Aumer, Simon White, Thorsten Naab, Cecilia Scannapieco

Abstract: We present an update to the multiphase SPH galaxy formation code by Scannapieco et al. We include a more elaborate treatment of the production of metals, cooling rates based on individual element abundances, and a scheme for the turbulent diffusion of metals. Our SN feedback model now transfers energy to the ISM in kinetic and thermal form, and we include a prescription for the effects of radiatio… ▽ More We present an update to the multiphase SPH galaxy formation code by Scannapieco et al. We include a more elaborate treatment of the production of metals, cooling rates based on individual element abundances, and a scheme for the turbulent diffusion of metals. Our SN feedback model now transfers energy to the ISM in kinetic and thermal form, and we include a prescription for the effects of radiation pressure from massive young stars on the ISM. We calibrate our new code on the well studied Aquarius haloes and then use it to simulate a sample of 16 galaxies with halo masses between 1x10^11 and 3x10^12 M_sun. In general, the stellar masses of the sample agree well with the stellar mass to halo mass relation inferred from abundance matching techniques for redshifts z=0-4. There is however a tendency to overproduce stars at z>4 and to underproduce them at z<0.5 in the least massive haloes. Overly high SFRs at z<1 for the most massive haloes are likely connected to the lack of AGN feedback in our model. The simulated sample also shows reasonable agreement with observed star formation rates, sizes, gas fractions and gas-phase metallicities at z=0-3. Remaining discrepancies can be connected to deviations from predictions for star formation histories from abundance matching. At z=0, the model galaxies show realistic morphologies, stellar surface density profiles, circular velocity curves and stellar metallicities, but overly flat metallicity gradients. 15 out of 16 of our galaxies contain disk components with kinematic disk fraction ranging between 15 and 65 %. The disk fraction depends on the time of the last destructive merger or misaligned infall event. Considering the remaining shortcomings of our simulations we conclude that even higher kinematic disk fractions may be possible for LambdaCDM haloes with quiet merger histories, such as the Aquarius haloes. △ Less

Submitted 16 July, 2013; v1 submitted 4 April, 2013; originally announced April 2013.

Comments: 26 pages, 20 figures, accepted for publication in MNRAS

Stellar haloes of simulated Milky Way-like galaxies: Chemical and kinematic properties

Authors: P. B. Tissera, C. Scannapieco, T. Beers, D. Carollo

Abstract: We investigate the chemical and kinematic properties of the diffuse stellar haloes of six simulated Milky Way-like galaxies from the Aquarius Project. Binding energy criteria are adopted to defined two dynamically distinct stellar populations: the diffuse inner and outer haloes, which comprise different stellar sub-populations with particular chemical and kinematic characteristics. Our simulated i… ▽ More We investigate the chemical and kinematic properties of the diffuse stellar haloes of six simulated Milky Way-like galaxies from the Aquarius Project. Binding energy criteria are adopted to defined two dynamically distinct stellar populations: the diffuse inner and outer haloes, which comprise different stellar sub-populations with particular chemical and kinematic characteristics. Our simulated inner- and outer-halo stellar populations have received contributions from debris stars (formed in sub-galactic systems while they were outside the virial radius of the main progenitor galaxies) and endo-debris stars (those formed in gas-rich sub-galactic systems inside the dark matter haloes). The inner haloes possess an additional contribution from disc-heated stars in the range $\sim 3 - 30 %$, with a mean of $\sim 20% $. Disc-heated stars might exhibit signatures of kinematical support, in particular among the youngest ones. Endo-debris plus disc-heated stars define the so-called \insitu stellar populations. In both the inner- and outer-halo stellar populations, we detect contributions from stars with moderate to low [$α$/Fe] ratios, mainly associated with the endo-debris or disc-heated sub-populations. The observed abundance gradients in the inner-halo regions are influenced by both the level of chemical enrichment and the relative contributions from each stellar sub-population. Steeper abundance gradients in the inner-halo regions are related to contributions from the disc-heated and endo-debris stars, which tend to be found at lower binding energies than debris stars. (Abridged). △ Less

Submitted 19 April, 2013; v1 submitted 7 January, 2013; originally announced January 2013.

Comments: 12 pages, 6 figures. Accepted version. To appear in MNRAS

Bars in hydrodynamical cosmological simulations

Authors: Cecilia Scannapieco, Evangelia Athanassoula

Abstract: We study the properties of two bars formed in fully cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies. In one case, the bar formed in a system with disc, bulge and halo components and is relatively strong and long, as could be expected for a system where the spheroid strongly influences the evolution. The second bar is less strong, shorter, and formed in a galaxy… ▽ More We study the properties of two bars formed in fully cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies. In one case, the bar formed in a system with disc, bulge and halo components and is relatively strong and long, as could be expected for a system where the spheroid strongly influences the evolution. The second bar is less strong, shorter, and formed in a galaxy with no significant bulge component. We study the strength and length of the bars, the stellar density profiles along and across the bars and the velocity fields in the bar region. We compare them with the results of dynamical (idealised) simulations and with observations, and find, in general, a good agreement, although we detect some important differences as well. Our results show that more or less realistic bars can form naturally in a $Λ$CDM cosmology, and open up the possibility to study the bar formation process in a more consistent way than previously done, since the host galaxies grow, accrete matter and significantly evolve during the formation and evolution of the bar. △ Less

Submitted 22 June, 2012; originally announced June 2012.

Comments: 5 pages, 4 figures, MNRAS accepted

The Aquila comparison Project: The Effects of Feedback and Numerical Methods on Simulations of Galaxy Formation

Authors: C. Scannapieco, M. Wadepuhl, O. H. Parry, J. F. Navarro, A. Jenkins, V. Springel, R. Teyssier, E. Carlson, H. M. P. Couchman, R. A. Crain, C. Dalla Vecchia, C. S. Frenk, C. Kobayashi, P. Monaco, G. Murante, T. Okamoto, T. Quinn, J. Schaye, G. S. Stinson, T. Theuns, J. Wadsley, S. D. M. White, R. Woods

Abstract: We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we… ▽ More We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge. △ Less

Submitted 19 April, 2012; v1 submitted 1 December, 2011; originally announced December 2011.

Comments: 28 pages, 16 figures, MNRAS accepted. High resolution figures and additional images can be found at http://www.aip.de/People/cscannapieco/aquila/

Chemical signatures of formation processes in the stellar populations of simulated galaxies

Authors: Patricia B. Tissera, Simon D. M. White, Cecilia Scannapieco

Abstract: We study the chemical properties of the stellar populations in eight simulations of the formation of Milky-Way mass galaxies in a LCDM Universe. Our simulations include metal-dependent cooling and an explicitly multiphase treatment of the effects on the gas of cooling, enrichment and supernova feedback. We search for correlations between formation history and chemical abundance patterns. Differing… ▽ More We study the chemical properties of the stellar populations in eight simulations of the formation of Milky-Way mass galaxies in a LCDM Universe. Our simulations include metal-dependent cooling and an explicitly multiphase treatment of the effects on the gas of cooling, enrichment and supernova feedback. We search for correlations between formation history and chemical abundance patterns. Differing contributions to spheroids and discs from in situ star formation and from accreted populations are reflected in differing chemical properties. Discs have younger stellar populations, with most stars forming in situ and with low alpha-enhancement from gas which never participated in a galactic outflow. Up to 15 per cent of disc stars can come from accreted satellites. These tend to be alpha-enhanced, older and to have larger velocity dispersions than the in situ population. Inner spheroids have old, metal-rich and alpha-enhanced stars which formed primarily in situ, more than 40 per cent from material recycled through earlier galactic winds. Few accreted stars are found in the inner spheroid unless a major merger occurred recently. Such stars are older, more metal-poor and more alpha-enhanced than the in situ population. Stellar haloes tend to have low metallicity and high alpha-enhancement. The outer haloes are made primarily of accreted stars. Their mean metallicity and alpha-enhancement reflect the masses of the disrupted satellites where they formed: more massive satellites typically have higher [Fe/H] and lower [alpha/Fe]. Surviving satellites have distinctive chemical patterns which reflect their extended, bursty star formation histories. These produce lower alpha-enhancement at given metallicity than in the main galaxy, in agreement with observed trends in the Milky Way. △ Less

Submitted 26 October, 2011; originally announced October 2011.

Comments: 17 pages, 13 figures, accepted MNRAS

The GALEX Arecibo SDSS Survey. IV. Baryonic Mass-Velocity-Size Relations of Massive Galaxies

Authors: Barbara Catinella, Guinevere Kauffmann, David Schiminovich, Jenna Lemonias, Cecilia Scannapieco, Jing Wang, Silvia Fabello, Cameron Hummels, Sean M. Moran, Ronin Wu, Andrew P. Cooper, Riccardo Giovanelli, Martha P. Haynes, Timothy M. Heckman, Amélie Saintonge

Abstract: We present dynamical scaling relations for a homogeneous and representative sample of ~500 massive galaxies, selected only by stellar mass (>10^10 Msun) and redshift (0.025<z<0.05) as part of the ongoing GALEX Arecibo SDSS Survey. We compare baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations for this sample, and investigate how galaxies scatter around the best fits obtained for pruned s… ▽ More We present dynamical scaling relations for a homogeneous and representative sample of ~500 massive galaxies, selected only by stellar mass (>10^10 Msun) and redshift (0.025<z<0.05) as part of the ongoing GALEX Arecibo SDSS Survey. We compare baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations for this sample, and investigate how galaxies scatter around the best fits obtained for pruned subsets of disk-dominated and bulge-dominated systems. The BFJ relation is significantly less scattered than the BTF when the relations are applied to their maximum samples, and is not affected by the inclination problems that plague the BTF. Disk-dominated, gas-rich galaxies systematically deviate from the BFJ relation defined by the spheroids. We demonstrate that by applying a simple correction to the stellar velocity dispersions that depends only on the concentration index of the galaxy, we are able to bring disks and spheroids onto the same dynamical relation -- in other words, we obtain a generalized BFJ relation that holds for all the galaxies in our sample, regardless of morphology, inclination or gas content, and has a scatter smaller than 0.1 dex. We find that disks and spheroids are offset in the stellar dispersion-size relation, and that the offset is removed when corrected dispersions are used instead. The generalized BFJ relation represents a fundamental correlation between the global dark matter and baryonic content of galaxies, which is obeyed by all (massive) systems regardless of morphology. [abridged] △ Less

Submitted 19 October, 2011; originally announced October 2011.

Comments: 20 pages, 15 figures. Accepted for publication in MNRAS. GASS publications and released data can be found at http://www.mpa-garching.mpg.de/GASS/index.php

Formation history, structure and dynamics of discs and spheroids in simulated Milky Way mass galaxies

Authors: Cecilia Scannapieco, Simon D. M White, Volker Springel, Patricia B. Tissera

Abstract: We study the stellar discs and spheroids in eight simulations of galaxy formation within Milky Way-mass haloes in a Lambda Cold Dark Matter cosmology. A first paper in this series concentrated on disc properties. Here, we extend this analysis to study how the formation history, structure and dynamics of discs and spheroids relate to the assembly history and structure of their haloes. We find that… ▽ More We study the stellar discs and spheroids in eight simulations of galaxy formation within Milky Way-mass haloes in a Lambda Cold Dark Matter cosmology. A first paper in this series concentrated on disc properties. Here, we extend this analysis to study how the formation history, structure and dynamics of discs and spheroids relate to the assembly history and structure of their haloes. We find that discs are generally young, with stars spanning a wide range in stellar age: the youngest stars define thin discs and have near-circular orbits, while the oldest stars form thicker discs which rotate ~2 times slower than the thin components, and have 2-3 times larger velocity dispersions. Unlike the discs, spheroids form early and on short time-scales, and are dominated by velocity dispersion. We find great variety in their structure. The inner regions are bar- or bulge-like, while the extended outer haloes are rich in complex non-equilibrium structures such as stellar streams, shells and clumps. Our discs have very high in-situ fractions, i.e. most of their stars formed in the disc itself. Nevertheless, there is a non-negligible contribution (~15 percent) from satellites that are accreted on nearly coplanar orbits. The inner regions of spheroids also have relatively high in-situ fractions, but 65-85 percent of their outer stellar population is accreted. We analyse the circular velocities, rotation velocities and velocity dispersions of our discs and spheroids, both for gas and stars, showing that the dynamical structure is complex as a result of the non-trivial interplay between cooling and SN heating. △ Less

Submitted 3 May, 2011; originally announced May 2011.

Comments: Submitted to MNRAS. 20 pages, 18 figures

Local Group Dwarf Galaxies: Nature And Nurture

Authors: Till Sawala, Cecilia Scannapieco, Simon White

Abstract: We investigate the formation and evolution of dwarf galaxies in a high resolution, hydrodynamical cosmological simulation of a Milky Way sized halo and its environment. Our simulation includes gas cooling, star formation, supernova feedback, metal enrichment and UV heating. In total, 90 satellites and more than 400 isolated dwarf galaxies are formed in the simulation, allowing a systematic study o… ▽ More We investigate the formation and evolution of dwarf galaxies in a high resolution, hydrodynamical cosmological simulation of a Milky Way sized halo and its environment. Our simulation includes gas cooling, star formation, supernova feedback, metal enrichment and UV heating. In total, 90 satellites and more than 400 isolated dwarf galaxies are formed in the simulation, allowing a systematic study of the internal and environmental processes that determine their evolution. We find that 95% of satellite galaxies are gas-free at z=0, and identify three mechanisms for gas loss: supernova feedback, tidal stripping, and photo-evaporation due to re-ionization. Gas-rich satellite galaxies are only found with total masses above ~ 5x10^9 solar masses. In contrast, for isolated dwarf galaxies, a total mass of ~ 10^9 solar masses constitutes a sharp transition; less massive galaxies are predominantly gas-free at z=0, more massive, isolated dwarf galaxies are often able to retain their gas. In general, we find that the total mass of a dwarf galaxy is the main factor which determines its star formation, metal enrichment, and its gas content, but that stripping may explain the observed difference in gas content between field dwarf galaxies and satellites with total masses close to 10^9 solar masses. We also find that a morphological transformation via tidal stripping of infalling, luminous dwarf galaxies whose dark matter is less concentrated than their stars, cannot explain the high total mass-light ratios of the faint dwarf spheroidal galaxies. △ Less

Submitted 23 March, 2011; originally announced March 2011.

Comments: 18 pages, 13 figures, submitted to MNRAS

What is the (Dark) Matter with Dwarf Galaxies?

Authors: Till Sawala, Qi Guo, Cecilia Scannapieco, Adrian Jenkins, Simon D. M. White

Abstract: We present cosmological hydrodynamical simulations of the formation of dwarf galaxies in a representative sample of haloes extracted from the Millennium-II Simulation. Our six haloes have a z = 0 mass of ~10^10 solar masses and show different mass assembly histories which are reflected in different star formation histories. We find final stellar masses in the range 5 x 10^7 - 10^8 solar masses, co… ▽ More We present cosmological hydrodynamical simulations of the formation of dwarf galaxies in a representative sample of haloes extracted from the Millennium-II Simulation. Our six haloes have a z = 0 mass of ~10^10 solar masses and show different mass assembly histories which are reflected in different star formation histories. We find final stellar masses in the range 5 x 10^7 - 10^8 solar masses, consistent with other published simulations of galaxy formation in similar mass haloes. Our final objects have structures and stellar populations consistent with dwarf elliptical and dwarf irregular galaxies. However, in a Lambda CDM universe, 10^10 solar mass haloes must typically contain galaxies with much lower stellar mass than our simulated objects if they are to match observed galaxy abundances. The dwarf galaxies formed in our own and all other current hydrodynamical simulations are more than an order of magnitude more luminous than expected for haloes of this mass. We discuss the significance and possible implications of this result. △ Less

Submitted 7 December, 2010; v1 submitted 3 March, 2010; originally announced March 2010.

Comments: 12 pages, 4 figures, MNRAS accepted

An observer's view of simulated galaxies: disc-to-total ratios, bars, and (pseudo-)bulges

Authors: Cecilia Scannapieco, Dimitri A. Gadotti, Patrik Jonsson, Simon D. M. White

Abstract: We use cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies to study the relative importance of the main stellar components, i.e., discs, bulges, and bars, at redshift zero. The main aim of this work is to understand if estimates of the structural parameters of these components determined from kinematics (as is usually done in simulations) agree well with those obtai… ▽ More We use cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies to study the relative importance of the main stellar components, i.e., discs, bulges, and bars, at redshift zero. The main aim of this work is to understand if estimates of the structural parameters of these components determined from kinematics (as is usually done in simulations) agree well with those obtained using a photometric bulge/disc/bar decomposition (as done in observations). To perform such a comparison, we have produced synthetic observations of the simulation outputs with the Monte-Carlo radiative transfer code SUNRISE and used the BUDDA code to make 2D photometric decompositions of the resulting images (in the i and g bands). We find that the kinematic disc-to-total ratio (D/T) estimates are systematically and significantly lower than the photometric ones. While the maximum D/T ratios obtained with the former method are of the order of 0.2, they are typically >0.4, and can be as high as 0.7, according to the latter. The photometric decomposition shows that many of the simulated galaxies have bars, with Bar/T ratios in the range 0.2-0.4, and that bulges have in all cases low Sersic indices, resembling observed pseudo-bulges instead of classical ones. Simulated discs, bulges and bars generally have similar (g-i) colours, which are in the blue tail of the distribution of observed colours. This is not due to the presence of young stars, but rather to low metallicities and poor gas content in the simulated galaxies, which makes dust extinction low. Photometric decompositions thus match the component ratios usually quoted for spiral galaxies better than kinematic decompositions, but the shift is insufficient to make the simulations consistent with observed late-type systems. △ Less

Submitted 11 June, 2010; v1 submitted 27 January, 2010; originally announced January 2010.

Comments: MNRAS Letters accepted, revised version, conclusions unchanged

Dark matter response to galaxy formation

Authors: Patricia B. Tissera, Simon D. M. White, Susana Pedrosa, Cecilia Scannapieco

Abstract: We have resimulated the six galaxy-sized haloes of the Aquarius Project including metal-dependent cooling, star formation and supernova feedback. This allows us to study not only how dark matter haloes respond to galaxy formation, but also how this response is affected by details of halo assembly history. In agreement with previous work, we find baryon condensation to lead to increased dark matt… ▽ More We have resimulated the six galaxy-sized haloes of the Aquarius Project including metal-dependent cooling, star formation and supernova feedback. This allows us to study not only how dark matter haloes respond to galaxy formation, but also how this response is affected by details of halo assembly history. In agreement with previous work, we find baryon condensation to lead to increased dark matter concentration. Dark matter density profiles differ substantially in shape from halo to halo when baryons are included, but in all cases the velocity dispersion decreases monotonically with radius. Some haloes show an approximately constant dark matter velocity anisotropy with $ β\approx 0.1-02$, while others retain the anisotropy structure of their baryon-free versions. Most of our haloes become approximately oblate in their inner regions, although a few retain the shape of their dissipationless counterparts. Pseudo-phase-space densities are described by a power law in radius of altered slope when baryons are included. The shape and concentration of the dark matter density profiles are not well reproduced by published adiabatic contraction models. The significant spread we find in the density and kinematic structure of our haloes appears related to differences in their formation histories. Such differences already affect the final structure in baryon-free simulations, but they are reinforced by the inclusion of baryons, and new features are produced. The details of galaxy formation need to be better understood before the inner dark matter structure of galaxies can be used to constrain cosmological models or the nature of dark matter. △ Less

Submitted 2 April, 2010; v1 submitted 12 November, 2009; originally announced November 2009.

Comments: 14 pages, 9 figures. Accepted MNRAS. Revised version includes discussion on resolution effects and minor changes.

The joint evolution of baryons and dark matter haloes

Authors: Susana E. Pedrosa, Patricia B. Tissera, Cecilia Scannapieco

Abstract: We have studied the dark matter (DM) distribution in a approx 10^12 h^-1 M_sun mass halo extracted from a simulation consistent with the concordance cosmology, where the physics regulating the transformation of gas into stars was allowed to change producing galaxies with different morphologies. The presence of baryons produces the concentration of the DM halo with respect to its corresponding di… ▽ More We have studied the dark matter (DM) distribution in a approx 10^12 h^-1 M_sun mass halo extracted from a simulation consistent with the concordance cosmology, where the physics regulating the transformation of gas into stars was allowed to change producing galaxies with different morphologies. The presence of baryons produces the concentration of the DM halo with respect to its corresponding dissipationless run, but we found that this response does not only depend on the amount of baryons gathered in the central region but also on the way they have been assembled. DM and baryons affect each other in a complex way so the formation history of a galaxy plays an important role on its final total mass distribution. Supernova (SN) feedback regulates the star formation and triggers galactic outflows not only in the central galaxy but also in its satellites. Our results suggest that, as the effects of SN feedback get stronger, satellites get less massive and can even be more easily disrupted by dynamical friction, transferring less angular momentum. We found indications that this angular momentum could be acquired not only by the outer part of the DM halo but also by the inner ones and by the stellar component in the central galaxy. The latter effect produces stellar migration which contributes to change the inner potential well, probably working against further DM contraction. As a consequence of the action of these processes, when the halo hosts a galaxy with an important disc structure formed by smooth gas accretion, it is more concentrated than when it hosts a spheroidal system which experienced more massive mergers and interactions. (abridged) △ Less

Submitted 30 November, 2009; v1 submitted 22 October, 2009; originally announced October 2009.

Comments: 16 pages, 16 figures, accepted to MNRAS

The impact of baryons on dark matter haloes

Authors: Susana E. Pedrosa, Patricia B. Tissera, Cecilia Scannapieco

Abstract: We analyse the dark matter (DM) distribution in a approx 10^12 M_sun halo extracted from a simulation consistent with the concordance cosmology, where the physics regulating the transformation of gas into stars was allowed to change producing galaxies with different morphologies. Although the DM profiles get more concentrated as baryons are collected at the centre of the haloes compared to a pur… ▽ More We analyse the dark matter (DM) distribution in a approx 10^12 M_sun halo extracted from a simulation consistent with the concordance cosmology, where the physics regulating the transformation of gas into stars was allowed to change producing galaxies with different morphologies. Although the DM profiles get more concentrated as baryons are collected at the centre of the haloes compared to a pure dynamical run, the total baryonic mass alone is not enough to fully predict the reaction of the DM profile. We also note that baryons affect the DM distribution even outside the central regions. Those systems where the transformation of gas into stars is regulated by Supernova (SN) feedback, so that significant disc structures are able to form, are found to have more concentrated dark matter profiles than a galaxy which has efficiently transformed most of its baryons into stars at early times. The accretion of satellites is found to be associated with an expansion of the dark matter profiles, triggered by angular momentum transfer from the incoming satellites. As the impact of SN feedback increases, the satellites get less massive and are even strongly disrupted before getting close to the main structure causing less angular momentum transfer. Our findings suggest that the response of the DM halo is driven by the history of assembly of baryons into a galaxy along their merger tree. △ Less

Submitted 12 February, 2009; originally announced February 2009.

Comments: 5 pages, 5 figures, accepted for publication in MNRAS

Formation of Isolated Dwarf Galaxies with Feedback

Authors: Till Sawala, Cecilia Scannapieco, Umberto Maio, Simon D. M. White

Abstract: We present results of high resolution hydrodynamical simulations of the formation and evolution of dwarf galaxies. Our simulations start from cosmological initial conditions at high redshift. They include metal-dependent cooling, star formation, feedback from type II and type Ia supernovae and UV background radiation, with physical recipes identical to those applied in a previous study of Milky… ▽ More We present results of high resolution hydrodynamical simulations of the formation and evolution of dwarf galaxies. Our simulations start from cosmological initial conditions at high redshift. They include metal-dependent cooling, star formation, feedback from type II and type Ia supernovae and UV background radiation, with physical recipes identical to those applied in a previous study of Milky Way type galaxies. We find that a combination of feedback and the cosmic UV background results in the formation of galaxies with properties similar to the Local Group dwarf spheroidals, and that their effect is strongly moderated by the depth of the gravitational potential. Taking this into account, our models naturally reproduce the observed luminosities and metallicities. The final objects have halo masses between 2.3x10^8 and 1.1x10^9 solar masses, mean velocity dispersions between 6.5 and 9.7 kms-1, stellar masses ranging from 5x10^5 to 1.2x10^7 solar masses, median metallicities between [Fe/H] = -1.8 and -1.1, and half-light radii of the order of 200 to 300 pc, all comparable with Local Group dwarf spheroidals. Our simulations also indicate that the dwarf spheroidal galaxies observed today lie near a halo mass threshold around 10^9 solar masses, in agreement with stellar kinematic data, where supernova feedback not only suffices to completely expel the interstellar medium and leave the residual gas-free, but where the combination of feedback, UV radiation and self-shielding establishes a dichotomy of age distributions similar to that observed in the Milky Way and M31 satellites. △ Less

Submitted 23 November, 2009; v1 submitted 10 February, 2009; originally announced February 2009.

Comments: 17 pages, 14 figures. MNRAS accepted

The Formation and Survival of Discs in a Lambda-CDM Universe

Authors: Cecilia Scannapieco, Simon D. M. White, Volker Springel, Patricia B. Tissera

Abstract: We study the formation of galaxies in a Lambda-CDM Universe using high resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to… ▽ More We study the formation of galaxies in a Lambda-CDM Universe using high resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at z = 0, none of our galaxies contain a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The z = 0 spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at z = 0 nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at z > ~2, regardless of their z = 0 morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic Lambda-CDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models. △ Less

Submitted 3 April, 2009; v1 submitted 4 December, 2008; originally announced December 2008.

Comments: 14 pages, 10 figures, mn2e.cls. MNRAS in press, updated to match published version

Effects of Supernova Feedback on the Formation of Galaxies

Authors: Cecilia Scannapieco, Patricia B. Tissera, Simon D. M. White, Volker Springel

Abstract: We study the effects of Supernova (SN) feedback on the formation of galaxies using hydrodynamical simulations in a Lambda-CDM cosmology. We use an extended version of the code GADGET-2 which includes chemical enrichment and energy feedback by Type II and Type Ia SN, metal-dependent cooling and a multiphase model for the gas component. We focus on the effects of SN feedback on the star formation… ▽ More We study the effects of Supernova (SN) feedback on the formation of galaxies using hydrodynamical simulations in a Lambda-CDM cosmology. We use an extended version of the code GADGET-2 which includes chemical enrichment and energy feedback by Type II and Type Ia SN, metal-dependent cooling and a multiphase model for the gas component. We focus on the effects of SN feedback on the star formation process, galaxy morphology, evolution of the specific angular momentum and chemical properties. We find that SN feedback plays a fundamental role in galaxy evolution, producing a self-regulated cycle for star formation, preventing the early consumption of gas and allowing disks to form at late times. The SN feedback model is able to reproduce the expected dependence on virial mass, with less massive systems being more strongly affected. △ Less

Submitted 21 August, 2008; v1 submitted 20 August, 2008; originally announced August 2008.

Comments: To appear in "The Galaxy Disk in Cosmological Context"; Proceedings of IAU254; 9-13 June 2008; Copenhagen; v2: typo corrected; uses iaus.cls

Effects of SN Feedback on the Dark Matter Distribution

Authors: Susana E. Pedrosa, Patricia B. Tissera, Cecilia Scannapieco

Abstract: We use cosmological simulations to study the effects of supernova (SN) feedback on the dark matter distribution in galaxies. We simulate the formation of a Milky-Way type galaxy using a version of the SPH code GADGET2 which includes chemical enrichment and energy feedback by SN, a multiphase model for the gas component and metal-dependent cooling. We analyse the impact of the main three input SN… ▽ More We use cosmological simulations to study the effects of supernova (SN) feedback on the dark matter distribution in galaxies. We simulate the formation of a Milky-Way type galaxy using a version of the SPH code GADGET2 which includes chemical enrichment and energy feedback by SN, a multiphase model for the gas component and metal-dependent cooling. We analyse the impact of the main three input SN feedback parameters on the amplitude and shape of the dark matter density profiles, focusing on the inner regions of the halo. In order to test the dependence of the results on the halo mass, we simulated a scale-down version of this system. First results of this ongoing work show that the dark matter distribution is affected by the feedback, through the redistribution of the baryons. Our findings suggest that the response of the dark matter halo could be the result of a combination of several physical parameters such as the amount of stellar mass formed at the centre, its shape, and probably the bursty characteristics of the star formation rate. As expected, we find that the dark matter haloes of small galaxies are more sensitive to SN feedback. Higher resolution simulations are being performed to test for numerical effects. △ Less

Submitted 16 July, 2008; originally announced July 2008.

Comments: 6 pages, 5 figures, Proceedings of IAU Symposium 254 "The Galaxy disk in a cosmological context", Copenhagen, June 2008

Effects of Supernova Feedback on the Formation of Galaxy Disks

Authors: Cecilia Scannapieco, Patricia B. Tissera, Simon D. M. White, Volker Springel

Abstract: We use cosmological simulations in order to study the effects of supernova (SN) feedback on the formation of a Milky Way-type galaxy of virial mass ~10^12 M_sun/h. We analyse a set of simulations run with the code described by Scannapieco et al. (2005, 2006), where we have tested our star formation and feedback prescription using isolated galaxy models. Here we extend this work by simulating the… ▽ More We use cosmological simulations in order to study the effects of supernova (SN) feedback on the formation of a Milky Way-type galaxy of virial mass ~10^12 M_sun/h. We analyse a set of simulations run with the code described by Scannapieco et al. (2005, 2006), where we have tested our star formation and feedback prescription using isolated galaxy models. Here we extend this work by simulating the formation of a galaxy in its proper cosmological framework, focusing on the ability of the model to form a disk-like structure in rotational support. We find that SN feedback plays a fundamental role in the evolution of the simulated galaxy, efficiently regulating the star formation activity, pressurizing the gas and generating mass-loaded galactic winds. These processes affect several galactic properties such as final stellar mass, morphology, angular momentum, chemical properties, and final gas and baryon fractions. In particular, we find that our model is able to reproduce extended disk components with high specific angular momentum and a significant fraction of young stars. The galaxies are also found to have significant spheroids composed almost entirely of stars formed at early times. We find that most combinations of the input parameters yield disk-like components, although with different sizes and thicknesses, indicating that the code can form disks without fine-tuning the implemented physics. We also show how our model scales to smaller systems. By analysing simulations of virial masses 10^9 M_sun/h and 10^10 M_sun/h, we find that the smaller the galaxy, the stronger the SN feedback effects. △ Less

Submitted 24 July, 2008; v1 submitted 23 April, 2008; originally announced April 2008.

Comments: 13 pages, 10 figures, mne2.cls, MNRAS, replaced with accepted version

The host galaxies of long-duration GRBs in a cosmological hierarchical scenario

Authors: S. E. Nuza, P. B. Tissera, L. J. Pellizza, D. G. Lambas, C. Scannapieco, M. E. De Rossi

Abstract: We developed a Monte Carlo code to generate long-duration gamma ray burst (LGRB) events within cosmological hydrodynamical simulations consistent with the concordance model. As structure is assembled, LGRBs are generated in the substructure that formed galaxies today. We adopted the collapsar model so that LGRBs are produced by single, massive stars at the final stage of their evolution. We foun… ▽ More We developed a Monte Carlo code to generate long-duration gamma ray burst (LGRB) events within cosmological hydrodynamical simulations consistent with the concordance model. As structure is assembled, LGRBs are generated in the substructure that formed galaxies today. We adopted the collapsar model so that LGRBs are produced by single, massive stars at the final stage of their evolution. We found that the observed properties of the LGRB host galaxies (HGs) are reproduced if LGRBs are also required to be generated by low metallicity stars. The low metallicity condition imposed on the progenitor stars of LGRBs selects a sample of HGs with mean gas abundances of 12 + log O/H \~ 8.6. For z<1 the simulated HGs of low metallicity LGRB progenitors tend to be faint, slow rotators with high star formation efficiency, compared with the general galaxy population, in agreement with observations. At higher redshift, our results suggest that larger systems with high star formation activity could also contribute to the generation of LGRBs from low metallicity progenitors since the fraction of low metallicity gas available for star formation increases for all systems with look-back time. Under the hypothesis of our LGRB model, our results support the claim that LGRBs could be unbiased tracers of star formation at high redshifts. △ Less

Submitted 5 February, 2007; v1 submitted 3 November, 2006; originally announced November 2006.

Comments: Final revised version with minor changes. 9 pages, 9 figures, mn2e.cls. To appear in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.375:665-672,2007

Clues for the origin of the fundamental metallicity relations. I: The hierarchical building up of the structure

Authors: Maria E. De Rossi, Patricia B. Tissera, Cecilia Scannapieco

Abstract: We analyse the evolutionary history of galaxies formed in a hierarchical scenario consistent with the concordance $Λ$-CDM model focusing on the study of the relation between their chemical and dynamical properties. Our simulations consistently describe the formation of the structure and its chemical enrichment within a cosmological context. Our results indicate that the luminosity-metallicity (L… ▽ More We analyse the evolutionary history of galaxies formed in a hierarchical scenario consistent with the concordance $Λ$-CDM model focusing on the study of the relation between their chemical and dynamical properties. Our simulations consistently describe the formation of the structure and its chemical enrichment within a cosmological context. Our results indicate that the luminosity-metallicity (LZR) and the stellar mass-metallicity (MZR) relations are naturally generated in a hierarchical scenario. Both relations are found to evolve with redshift. In the case of the MZR, the estimated evolution is weaker than that deduced from observational works by approximately 0.10 dex. We also determine a characteristic stellar mass, $M_c \approx 3 \times 10^{10} M_{\odot}$, which segregates the simulated galaxy population into two distinctive groups and which remains unchanged since $z\sim 3$, with a very weak evolution of its metallicity content. The value and role played by $M_c$ is consistent with the characteristic mass estimated from the SDSS galaxy survey by Kauffmann et al. (2004). Our findings suggest that systems with stellar masses smaller than $M_c$ are responsible for the evolution of this relation at least from $ z\approx 3$. Larger systems are stellar dominated and have formed more than 50 per cent of their stars at $z \ge 2$, showing very weak evolution since this epoch. We also found bimodal metallicity and age distributions from $z\sim3$, which reflects the existence of two different galaxy populations. Although SN feedback may affect the properties of galaxies and help to shape the MZR, it is unlikely that it will significantly modify $M_c$ since, from $z=3$ this stellar mass is found in systems with circular velocities larger than $100 \kms$. △ Less

Submitted 5 October, 2006; v1 submitted 8 September, 2006; originally announced September 2006.

Comments: 17 pages, 13 figures. Minor changes to match accepted version. Accepted October 3 MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.374:323-336,2007

Galaxy pairs in cosmological simulations: effects of interactions on colours and chemical abundances

Authors: M. Josefa Perez, Patricia B. Tissera, Cecilia Scannapieco, Diego G. Lambas, Maria E. De Rossi, .

Abstract: We perform an statistical analysis of galaxies in pairs in a Lambda-CDM scenario by using the chemical GADGET-2 of Scannapieco et al. (2005) in order to study the effects of galaxy interactions on colours and metallicities. We find that galaxy-galaxy interactions can produce a bimodal colour distribution with galaxies with significant recent star formation activity contributing mainly to blue co… ▽ More We perform an statistical analysis of galaxies in pairs in a Lambda-CDM scenario by using the chemical GADGET-2 of Scannapieco et al. (2005) in order to study the effects of galaxy interactions on colours and metallicities. We find that galaxy-galaxy interactions can produce a bimodal colour distribution with galaxies with significant recent star formation activity contributing mainly to blue colours. In the simulations, the colours and the fractions of recently formed stars of galaxies in pairs depend on environment more strongly than those of galaxies without a close companion, suggesting that interactions play an important role in galaxy evolution. If the metallicity of the stellar populations is used as the chemical indicator, we find that the simulated galaxies determine luminosity-metallicity and stellar mass-metallicity relations which do not depend on the presence of a close companion. However, in the case of the luminosity-metallicity relation, at a given level of enrichment, we detect a systematic displacement of the relation to brighter magnitudes for active star forming systems. Regardless of relative distance and current level of star formation activity, galaxies in pairs have stellar populations with higher level of enrichment than galaxies without a close companion. In the case of the gas component, this is no longer valid for galaxies in pairs with passive star formation which only show an excess of metals for very close pair members, consequence of an important recent past star formation activity. (Abridged). △ Less

Submitted 4 May, 2006; originally announced May 2006.

Comments: Submitted to A&A

Feedback and metal enrichment in cosmological SPH simulations - II. A multiphase model with supernova energy feedback

Authors: Cecilia Scannapieco, Patricia B. Tissera, Simon D. M. White, Volker Springel

Abstract: We have developed a new scheme to treat a multiphase interstellar medium in smoothed particle hydrodynamics simulations of galaxy formation. This scheme can represent a co-spatial mixture of cold and hot ISM components, and is formulated without scale-dependent parameters. It is thus particularly suited to studies of cosmological structure formation where galaxies with a wide range of masses for… ▽ More We have developed a new scheme to treat a multiphase interstellar medium in smoothed particle hydrodynamics simulations of galaxy formation. This scheme can represent a co-spatial mixture of cold and hot ISM components, and is formulated without scale-dependent parameters. It is thus particularly suited to studies of cosmological structure formation where galaxies with a wide range of masses form simultaneously. We also present new algorithms for energy and heavy element injection by supernovae, and show that together these schemes can reproduce several important observed effects in galaxy evolution. Both in collapsing systems and in quiescent galaxies our codes can reproduce the Kennicutt relation between the surface densities of gas and of star formation. Strongly metal-enhanced winds are generated in both cases with ratios of mass-loss to star formation which are similar to those observed. This leads to a self-regulated cycle for star formation activity. The overall impact of feedback depends on galaxy mass. Star formation is suppressed at most by a factor of a few in massive galaxies, but in low-mass systems the effects can be much larger, giving star formation an episodic, bursty character. The larger the energy fraction assumed available in feedback, the more massive the outflows and the lower the final stellar masses. Winds from forming disks are collimated perpendicular to the disk plane, reach velocities up to 1000 km/s, and efficiently transport metals out of the galaxies. The asymptotically unbound baryon fraction drops from >95 per cent to ~30 per cent from the least to the most massive of our idealised galaxies, but the fraction of all metals ejected with this component exceeds 60 per cent regardless of mass. Such winds could plausibly enrich the intergalactic medium to observed levels. △ Less

Submitted 23 August, 2006; v1 submitted 25 April, 2006; originally announced April 2006.

Comments: 18 pages, 13 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.371:1125,2006