Showing 1–39 of 39 results for author: Di Cintio, A

The role of AGN feedback on the evolution of dwarf galaxies from cosmological simulations: SMBHs suppress star formation in low-mass galaxies

Authors: Elena Arjona-Galvez, Arianna Di Cintio, Robert J. J. Grand

Abstract: Aims. Recent observational studies suggest that feedback from active galactic nuclei (AGNs) may play an important role in the formation and evolution of low-mass dwarf galaxies, an issue that has received little attention from a theoretical perspective. Methods. We investigate this using two sets of 12 cosmological magneto-hydrodynamical simulations of the formation of dwarf galaxies: one set usin… ▽ More Aims. Recent observational studies suggest that feedback from active galactic nuclei (AGNs) may play an important role in the formation and evolution of low-mass dwarf galaxies, an issue that has received little attention from a theoretical perspective. Methods. We investigate this using two sets of 12 cosmological magneto-hydrodynamical simulations of the formation of dwarf galaxies: one set using a version of the AURIGA galaxy formation physics model including AGN feedback and a parallel set with AGN feedback turned off. Results. We show that the full-physics AGN runs satisfactorily reproduce several scaling relations, including the M_ BH-M_*, M_ BH-sigma_* and the baryonic Tully-Fisher relation. We find that the global star formation (SF) of galaxies run with AGN is reduced compared to the one in which AGN has been turned off, suggesting that this type of feedback is a viable way of suppressing SF in low-mass dwarfs, even though none of our galaxies is completely quenched by z=0. Furthermore, we found a tight correlation between the median SF rates and the black-hole-to-stellar mass ratio (M_BH/M_star) in our simulated dwarfs. SF is suppressed due to gas heating in the vicinity of the AGN: less HI gas is available in AGN runs, though the total amount of gas is preserved across the two settings within each galaxy. This indicates that the main effect of AGN feedback in our dwarfs is to heat up and push the gas away from the galaxy's centre rather than expelling it completely. Finally, we show that the two galaxies harbouring the largest BHs have suffered a considerable (up to 65%) reduction in their central dark matter density, pinpointing the role of AGNs in determining the final dark matter mass distribution within dwarf galaxies. This pilot paper highlights the importance of modelling AGN feedback at the lowest mass scales and the impact this can have on dwarf galaxy evolution. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 13 pages, 10 figures and 2 tables Submitted to A&A

Constraining gas metal mixing strength in simulations using observations of the Milky Way's disc

Authors: Jorge Sarrato-Alós, Christopher Brook, Arianna Di Cintio

Abstract: This work explores the mixing rate of metals in the interstellar medium (ISM), comparing observational constraints from our solar neighbourhood to high resolution cosmological hydrodynamical simulations of Milky Way (MW)-like galaxies. The mixing rate, described by the coefficient C, is varied in simulations between 0 and 0.05, with resultant simulated galaxies compared to observations of metallic… ▽ More This work explores the mixing rate of metals in the interstellar medium (ISM), comparing observational constraints from our solar neighbourhood to high resolution cosmological hydrodynamical simulations of Milky Way (MW)-like galaxies. The mixing rate, described by the coefficient C, is varied in simulations between 0 and 0.05, with resultant simulated galaxies compared to observations of metallicity dispersion in young star clusters, HII regions and neutral gas in the disc of the MW. A value of C between 0.003125 and 0.0125 is found to self-consistently match a range of observables, with a best estimate of C=0.0064$\pm$0.0004. We demonstrate that the relationship between metal dispersion in young stars, HII regions and neutral gas, versus the coefficient C, can be described by a power law. These constrained mixing rates infer a comparatively well mixed ISM in the solar neighbourhood, at odds with some recent observations that have reported a highly inhomogeneous ISM. The degree of mixing suggested by this work is lower than what often employed in many hydrodynamical simulations. Our results have implications for studying the metallicity distribution of stars as well as of gas in the interstellar and circumgalactic media. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 7 pages, 3 figures, accepted for publication in MNRAS

On the anti-correlation between pericentric distance and inner dark matter density of Milky Way's dwarf spheroidal galaxies

Authors: Salvador Cardona-Barrero, Giuseppina Battaglia, Carlo Nipoti, Arianna Di Cintio

Abstract: An anti-correlation between the central density of the dark matter halo ($ρ_{150,\ {\rm DM}}$) and the pericentric distances ($r_{p}$) of the Milky Way's (MW's) dwarf spheroidal galaxies (dSphs) has been reported in the literature. The existence and origin of such anti-correlation is however controversial, one possibility being that only the densest dSphs can survive the tidal field towards the ce… ▽ More An anti-correlation between the central density of the dark matter halo ($ρ_{150,\ {\rm DM}}$) and the pericentric distances ($r_{p}$) of the Milky Way's (MW's) dwarf spheroidal galaxies (dSphs) has been reported in the literature. The existence and origin of such anti-correlation is however controversial, one possibility being that only the densest dSphs can survive the tidal field towards the centre of our Galaxy. In this work, we place particular emphasis on quantifying the statistical significance of such anti-correlation, by using available literature data in order to explore its robustness under different assumptions on the MW gravitational potential, and for various derivations of $ρ_{150}$ and $r_{p}$. We consider models in which the MW is isolated and has a low ($8.8\times10^{11}\,M_{\odot}$) and high ($1.6\times10^{12}\, M_{\odot}$) halo mass, respectively, as well as configurations in which the MW's potential is perturbed by a Large Magellanic Cloud (LMC) infall. We find that, while data generally support models in which the dSphs' central DM density decreases as a function of their pericentric radius, this anti-correlation is statistically significant at $3σ$ level only in $\sim$12$\%$ of the combinations of $ρ_{150}$ and $r_{p}$ explored. Moreover, including the impact of the LMC's infall onto the MW weakens or even washes away this anti-correlation, with respect to models in which the MW is isolated. Our results suggest that the strength and existence of such anti-correlation is still debatable: exploring it with high-resolution simulations including baryonic physics and different DM flavours will help us to understand its emergence. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: 10 pages, 3 figures. Accepted for publication in MNRAS

The undiscovered ultra-diffuse galaxies of the Local Group

Authors: Oliver Newton, Arianna Di Cintio, Salvador Cardona-Barrero, Noam I. Libeskind, Yehuda Hoffman, Alexander Knebe, Jenny Sorce, Matthias Steinmetz, Elmo Tempel

Abstract: Ultra-diffuse galaxies (UDGs) are attractive candidates to probe cosmological models and test theories of galaxy formation at low masses; however, they are difficult to detect because of their low surface brightness. In the Local Group a handful of UDGs have been found to date, most of which are satellites of the Milky Way and M31, and only two are isolated galaxies. It is unclear whether so few U… ▽ More Ultra-diffuse galaxies (UDGs) are attractive candidates to probe cosmological models and test theories of galaxy formation at low masses; however, they are difficult to detect because of their low surface brightness. In the Local Group a handful of UDGs have been found to date, most of which are satellites of the Milky Way and M31, and only two are isolated galaxies. It is unclear whether so few UDGs are expected. We address this by studying the population of UDGs formed in hydrodynamic constrained simulations of the Local Group from the HESTIA suite. For a Local Group with a total enclosed mass $M_{\rm LG}\!\left(< 2.5\, {\rm Mpc}\right)=8\times10^{12}\, {\rm M_\odot}$, we predict that there are $12\pm3$ isolated UDGs (68% confidence) with stellar masses $10^6 \leq M_\ast\, /\, {\rm M_\odot} < 10^9,$ and effective radii $R_{\rm e} \geq 1.5\, {\rm kpc}$, within 2.5 Mpc of the Local Group, of which $2^{+2}_{-1}$ (68% confidence) are detectable in the footprint of the Sloan Digital Sky Survey (SDSS). Accounting for survey incompleteness, we find that almost the entire population of UDGs in the Local Group field would be observable in a future all-sky survey with a depth similar to the SDSS, the Dark Energy Survey, or the Legacy Survey of Space and Time. Our results suggest that there is a population of UDGs in the Local Group awaiting discovery. △ Less

Submitted 28 March, 2023; v1 submitted 9 December, 2022; originally announced December 2022.

Comments: 9 pages including acknowledgements and references. Contains 5 figures and two tables. A repository of reduced data and scripts to reproduce the results in this manuscript is available at https://github.com/Musical-Neutron/lg_field_udgs and archived at https://doi.org/10.5281/zenodo.7757685 . Accepted by ApJL

A probabilistic deep learning model to distinguish cusps and cores in dwarf galaxies

Authors: J. Expósito-Márquez, C. B. Brook, M. Huertas-Company, A. Di Cintio, A. V. Macciò, R. J. J. Grand, G. Battaglia

Abstract: Numerical simulations within a cold dark matter (DM) cosmology form halos whose density profiles have a steep inner slope (`cusp'), yet observations of galaxies often point towards a flat central `core'. We develop a convolutional mixture density neural network model to derive a probability density function (PDF) of the inner density slopes of DM halos. We train the network on simulated dwarf gala… ▽ More Numerical simulations within a cold dark matter (DM) cosmology form halos whose density profiles have a steep inner slope (`cusp'), yet observations of galaxies often point towards a flat central `core'. We develop a convolutional mixture density neural network model to derive a probability density function (PDF) of the inner density slopes of DM halos. We train the network on simulated dwarf galaxies from the NIHAO and AURIGA projects, which include both DM cusps and cores: line-of-sight velocities and 2D spatial distributions of their stars are used as inputs to obtain a PDF representing the probability of predicting a specific inner slope. The model recovers accurately the expected DM profiles: $\sim$82$\%$ of the galaxies have a derived inner slope within $\pm$0.1 of their true value, while $\sim$98$\%$ within $\pm$0.3. We apply our model to four Local Group dwarf spheroidal galaxies and find results consistent with those obtained with the Jeans modelling based code GravSphere: the Fornax dSph has a strong indication of possessing a central DM core, Carina and Sextans have cusps (although the latter with large uncertainties), while Sculptor shows a double peaked PDF indicating that a cusp is preferred, but a core can not be ruled out. Our results show that simulation-based inference with neural networks provide a innovative and complementary method for the determination of the inner matter density profiles in galaxies, which in turn can help constrain the properties of the elusive DM. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 14 pages, 11 figures, submitted to MNRAS. comments welcome

Metallicity profiles of Ultra Diffuse Galaxies in NIHAO simulations

Authors: S. Cardona-Barrero, A. Di Cintio, G. Battaglia, A. V. Macciò, S. Taibi

Abstract: Supernovae feedback driven expansion has proven to be a viable mechanism to explain the average properties of Ultra Diffuse Galaxies (UDGs) such as the sizes, colors, mass and internal kinematics. Here, we explore the origin of stellar metallicity gradients in feedback driven simulated UDGs from the NIHAO project and compare them with the observed distribution of metallicity gradients of both Loca… ▽ More Supernovae feedback driven expansion has proven to be a viable mechanism to explain the average properties of Ultra Diffuse Galaxies (UDGs) such as the sizes, colors, mass and internal kinematics. Here, we explore the origin of stellar metallicity gradients in feedback driven simulated UDGs from the NIHAO project and compare them with the observed distribution of metallicity gradients of both Local Group dwarfs as well as of the recently observed UDG DF44. Simulated UDGs display a large variety of metallicity profiles, showing flat to negative gradients, similarly to what is observed in LG dwarfs, while DF44 data suggest a flat to positive gradient. The variety of metallicity gradients in simulations is set by the interplay between the radius at which star formation occurs and the subsequent supernovae feedback driven stellar redistribution: rotation supported systems tend to have flat metallicity profiles while dispersion supported galaxies show negative and steep profiles. Our results suggest that UDGs are not peculiar in what regards their metallicity gradients, when compared to regular dwarfs. Desirably, a larger observational sample of UDGs' gradients shall be available in the future, in order to test our predictions. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 13 pages, 6+3 figures

A shallow dark matter halo in Ultra Diffuse Galaxy AGC 242019: are UDGs structurally similar to low surface brightness galaxies?

Authors: Chris B. Brook, Arianna Di Cintio, Andrea V. Maccio, Marvin Blank

Abstract: A central question regarding Ultra Diffuse Galaxies (UDGs) is whether they are a separate category to Low Surface Brightness (LSB) galaxies, or just their natural continuation towards low stellar masses. In this letter, we show that the rotation curve of the gas rich UDG AGC 242019 is well fit by a dark matter halo with inner slope that asymptotes to -0.54, and that such fit provides a concentrati… ▽ More A central question regarding Ultra Diffuse Galaxies (UDGs) is whether they are a separate category to Low Surface Brightness (LSB) galaxies, or just their natural continuation towards low stellar masses. In this letter, we show that the rotation curve of the gas rich UDG AGC 242019 is well fit by a dark matter halo with inner slope that asymptotes to -0.54, and that such fit provides a concentration parameter that matches theoretical expectations. This finding, together with previously works in which shallow inner profiles are derived for UDGs, shows that the structural properties of these galaxies are like other observed LSBs. UDGs show slowly rising rotation curves and this favours formation scenarios in which internal processes, such as SNae driven gas outflows, are acting to modify UDGs profiles. △ Less

Submitted 3 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJL

Hydrographic variability and biomass fluctuations of European anchovy (Engraulis encrasicolus) in the Central Mediterranean Sea: Monetary estimations and impacts on fishery from Lagrangian analysis

Authors: Antonio Di Cintio, Marco Torri, Federico Falcini, Raffaele Corrado, Guglielmo Lacorata, Angela Cuttitta, Bernardo Patti, Rosalia Santoleri

Abstract: During the last decades, scientific community has been investigating both biological and hydrographic processes that affect fisheries. Such an interdisciplinary and synergic approach is nowadays giving a fundamental contribution, in particular, in connecting the dots between hydrographic phenomena and biomass variability and distribution of small pelagic fish. Here we estimate impacts of hydrograp… ▽ More During the last decades, scientific community has been investigating both biological and hydrographic processes that affect fisheries. Such an interdisciplinary and synergic approach is nowadays giving a fundamental contribution, in particular, in connecting the dots between hydrographic phenomena and biomass variability and distribution of small pelagic fish. Here we estimate impacts of hydrographic fluctuations on small pelagic fishery, focusing on the inter-annual variability that characterizes connectivity between spawning and recruiting areas for the European anchovy (Engraulis encrasicolus, Linnaues 1758), in the Northern side of the Sicily Channel (Mediterranean Sea). Results show that coastal transport dynamics of a specific year largely affect the biomass recorded the following year. Our work, moreover, quantifies the specific monetary impacts on landings of European anchovy fishery due to hydrodynamics variability, connecting biomass fluctuations with fishery economics in a highly dynamic and exploited marine environment as the Sicily Channel. In particular, we build a model that attributes a monetary value to the hydrographic phenomena (i.e., cross-shore vs. alongshore eggs and larvae transport), registered in the FAO Geographical Sub-Area (GSA) 16 (Southern Sicily). This allows us to provide a monetary estimation of catches, derived from different transport dynamics. Our results highlight the paramount importance that hydrographic phenomena can have over the socio-economic performance of a fishery. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Origin of stellar prolate rotation in a cosmologically simulated faint dwarf galaxy

Authors: Salvador Cardona-Barrero, Giuseppina Battaglia, Arianna Di Cintio, Yves Revaz, Pascale Jablonka

Abstract: Stellar prolate rotation in dwarf galaxies is rather uncommon, with only two known galaxies in the Local Group showing such feature (Phoenix and And II). Cosmological simulations show that in massive early-type galaxies prolate rotation likely arises from major mergers. However, the origin of such kinematics in the dwarf galaxies regime has only been explored using idealized simulations. Here we m… ▽ More Stellar prolate rotation in dwarf galaxies is rather uncommon, with only two known galaxies in the Local Group showing such feature (Phoenix and And II). Cosmological simulations show that in massive early-type galaxies prolate rotation likely arises from major mergers. However, the origin of such kinematics in the dwarf galaxies regime has only been explored using idealized simulations. Here we made use of hydrodynamical cosmological simulations of dwarfs galaxies with stellar mass between $3\times10^5$ and $5\times10^8$ M$_{\odot}$ to explore the formation of prolate rotators. Out of $27$ dwarfs, only one system showed clear rotation around the major axis, whose culprit is a major merger at $z=1.64$, which caused the transition from an oblate to a prolate configuration. Interestingly, this galaxy displays a steep metallicity gradient, reminiscent of the one measured in Phoenix and And II: this is the outcome of the merger event that dynamically heats old, metal-poor stars, and of the centrally concentrated residual star formation. Major mergers in dwarf galaxies offer a viable explanation for the formation of such peculiar systems, characterized by steep metallicity gradients and prolate rotation. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 5 pages, 3 figures

Pericentric passage-driven star formation in satellite galaxies and their hosts: CLUES from Local Group simulations

Authors: Arianna Di Cintio, Robert Mostoghiu, Alexander Knebe, Julio Navarro

Abstract: Local Group satellite galaxies show a wide diversity of star formation histories (SFHs) whose origin is yet to be fully understood. Using hydrodynamical simulations from the Constrained Local UniversE project, we study the SFHs of satellites of Milky Way-like galaxies in a cosmological context: while in the majority of the cases the accretion onto their host galaxy causes the satellites to lose th… ▽ More Local Group satellite galaxies show a wide diversity of star formation histories (SFHs) whose origin is yet to be fully understood. Using hydrodynamical simulations from the Constrained Local UniversE project, we study the SFHs of satellites of Milky Way-like galaxies in a cosmological context: while in the majority of the cases the accretion onto their host galaxy causes the satellites to lose their gas, with a subsequent suppression in SF, in about 25$\%$ of our sample we observe a clear enhancement of SF after infall. Peaks in SF clearly correlate with the satellite pericentric passage around its host and, in one case, with a satellite-satellite interaction. We identify two key ingredients that result in enhanced SF after infall: galaxies must enter the host's virial radius with a reservoir of cold gas $M_{\rm gas,inf}/M_{\rm vir,inf}\gtrsim 10^{-2}$ and with a minimum pericentric distance $\gtrsim$10 kpc (mean distance $\sim$50 kpc for the full sample), in order to form new stars due to compression of cold gas at pericentric passage. On the other hand, satellites that infall with little gas or whose pericentric distance is too small, have their gas ram-pressure stripped and subsequent SF quenched. The pericentric passage of satellites likewise correlates with SF peaks in their hosts, suggesting that this mechanism induces bursts of SF in satellites and central galaxies alike, in agreement with recent studies of our Galaxy's SFH. Our findings can explain the recently reported multiple stellar populations observed in dwarf galaxies such as Carina and Fornax, and should be taken into account in semi-analytic models of galaxy formation and satellite quenching. △ Less

Submitted 9 June, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: 15 pages, 12 figures, matches accepted MNRAS version. Expanded sec 3.2. New figures 2 and 3 show comparison between simulated and observed SFHs of satellites and isolated dwarfs

NIHAO XXIV: Rotation or pressure supported systems? Simulated Ultra Diffuse Galaxies show a broad distribution in their stellar kinematics

Authors: Salvador Cardona-Barrero, Arianna Di Cintio, Christopher B. A. Brook, Tomas Ruiz-Lara, Michael A. Beasley, Jesus Falcón-Barroso, Andrea V. Macciò

Abstract: In recent years a new window on galaxy evolution opened, thanks to the increasing discovery of galaxies with a low surface brightness, such as Ultra Diffuse Galaxies (UDGs). The formation mechanism of these systems is still a much debated question, and so are their kinematical properties. In this work, we address this topic by analyzing the stellar kinematics of isolated UDGs formed in the hydrody… ▽ More In recent years a new window on galaxy evolution opened, thanks to the increasing discovery of galaxies with a low surface brightness, such as Ultra Diffuse Galaxies (UDGs). The formation mechanism of these systems is still a much debated question, and so are their kinematical properties. In this work, we address this topic by analyzing the stellar kinematics of isolated UDGs formed in the hydrodynamical simulation suite NIHAO. We construct projected line-of-sight velocity and velocity dispersion maps to compute the projected specific angular momentum, $λ_{\rm R}$, to characterize the kinematical support of the stars in these galaxies. We found that UDGs cover a broad distribution, ranging from dispersion to rotation supported galaxies, with similar abundances in both regimes. The degree of rotation support of simulated UDGs correlates with several properties such as galaxy morphology, higher HI fractions and larger effective radii with respect to the dispersion supported group, while the dark matter halo spin and mass accretion history are similar amongst the two populations. We demonstrate that the alignment of the infalling baryons into the protogalaxy at early $z$ is the principal driver of the $z$=0 stellar kinematic state: pressure supported isolated UDGs form via mis-aligned gas accretion while rotation supported ones build-up their baryons in an ordered manner. Accounting for random inclination effects, we predict that a comprehensive survey will find nearly half of field UDGs to have rotationally supported stellar disks. △ Less

Submitted 16 July, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Self-Interacting Dark Matter and the Delay of Super-Massive Black Hole Growth

Authors: Akaxia Cruz, Andrew Pontzen, Marta Volonteri, Thomas R. Quinn, Michael Tremmel, Alyson M. Brooks, Nicole N. Sanchez, Ferah Munshi, Arianna Di Cintio

Abstract: Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find t… ▽ More Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find that SMBH formation is suppressed in the early universe due to SIDM coring. SMBH-SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyrs of their evolution. Nonetheless, by $z = 0.8$ the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the $M_{BH}-M_{*}$ relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model. △ Less

Submitted 17 April, 2020; originally announced April 2020.

The mass of our Galaxy from satellite proper motions in the Gaia era

Authors: T. K. Fritz, A. Di Cintio, G. Battaglia, C. Brook, S. Taibi

Abstract: We use Gaia DR2 systemic proper motions of 45 satellite galaxies to constrain the mass of the Milky Way using the scale free mass estimator of Watkins et al. (2010). We first determine the anisotropy parameter $β$, and the tracer satellites' radial density index $γ$ to be $β$=$-0.67^{+0.45}_{-0.62}$ and $γ=2.11\pm0.23$. When we exclude possible former satellites of the Large Magellanic Cloud, the… ▽ More We use Gaia DR2 systemic proper motions of 45 satellite galaxies to constrain the mass of the Milky Way using the scale free mass estimator of Watkins et al. (2010). We first determine the anisotropy parameter $β$, and the tracer satellites' radial density index $γ$ to be $β$=$-0.67^{+0.45}_{-0.62}$ and $γ=2.11\pm0.23$. When we exclude possible former satellites of the Large Magellanic Cloud, the anisotropy changes to $β$=$-0.21^{+0.37}_{-0.51}$. We find that the index of the Milky Way's gravitational potential $α$, which is dependent on the mass itself, is the parameter with the largest impact on the mass determination. Via comparison with cosmological simulations of Milky Way-like galaxies, we carried out a detailed analysis of the estimation of the observational uncertainties and their impact on the mass estimator. We found that the mass estimator is biased when applied naively to the satellites of simulated Milky Way halos. Correcting for this bias, we obtain for our Galaxy a mass of $0.58^{+0.15}_{-0.14}\times10^{12}$M$_\odot$ within 64 kpc, as computed from the inner half of our observational sample, and $1.43^{+0.35}_{-0.32}\times10^{12}$M$_\odot$ within 273 kpc, from the full sample; this latter value extrapolates to a virial mass of $M_\mathrm{vir\,Δ=97}$=$1.51^{+0.45}_{-0.40} \times 10^{12}M_{\odot}$ corresponding to a virial radius of R$_\mathrm{vir}$=$308\pm29$ kpc. This value of the Milky Way mass lies in-between other mass estimates reported in the literature, from various different methods. △ Less

Submitted 10 April, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 17 pages, 10 figures, 2 tables. Accepted for publication in MNRAS

NIHAO XXI: The emergence of Low Surface Brightness galaxies

Authors: Arianna Di Cintio, Chris B. Brook, Andrea V. Macciò, Aaron A. Dutton, Salvador Cardona-Barrero

Abstract: The existence of galaxies with a surface brightness $μ$ lower than the night sky has been known since three decades. Yet, their formation mechanism and emergence within a $\rmΛCDM$ universe has remained largely undetermined. For the first time, we investigated the origin of Low Surface Brightness (LSB) galaxies with M$_{\star}$$\sim$10$^{9.5-10}$M$_{\odot}$, which we are able to reproduce within h… ▽ More The existence of galaxies with a surface brightness $μ$ lower than the night sky has been known since three decades. Yet, their formation mechanism and emergence within a $\rmΛCDM$ universe has remained largely undetermined. For the first time, we investigated the origin of Low Surface Brightness (LSB) galaxies with M$_{\star}$$\sim$10$^{9.5-10}$M$_{\odot}$, which we are able to reproduce within hydrodynamical cosmological simulations from the NIHAO suite. The simulated and observed LSBs share similar properties, having large HI reservoir, extended star formation histories and effective radii, low Sérsic index and slowly rising rotation curves. The formation mechanism of these objects is explored: simulated LSBs form as a result of co-planar co-rotating mergers and aligned accretion of gas at early times, while perpendicular mergers and mis-aligned gas accretion result in higher $μ$ galaxies by $z$=0. The larger the merger, the stronger the correlation between merger orbital configuration and final $μ$. While the halo spin parameter is consistently high in simulated LSB galaxies, the impact of halo concentration, feedback-driven gas outflows and merger time only plays a minor-to-no role in determining $μ$. Interestingly, the formation scenario of such `classical' LSBs differs from the one of less massive, M$_{\star}$$\sim$10$^{7-9}$M$_{\odot}$, Ultra-Diffuse Galaxies, the latter resulting from the effects of SNae driven gas outflows: a M$_{\star}$ of $\sim$10$^9$M$_{\odot}$ thus represents the transition regime between a feedback-dominated to an angular momentum-dominated formation scenario in the LSB realm. Observational predictions are offered regarding spatially resolved star formation rates through LSB discs: these, together with upcoming surveys, can be used to verify the proposed emergence scenario of LSB galaxies. △ Less

Submitted 9 April, 2019; v1 submitted 24 January, 2019; originally announced January 2019.

Comments: MNRAS Accepted 2019 April 02; in original form 2019 January 24. 15 pages,17 figures. main results in Fig 12 and 14. Added appendix A

Dark Matter Response to Galaxy Assembly History

Authors: M. Celeste Artale, Susana E. Pedrosa, Patricia B. Tissera, Pedro Cataldi, Arianna Di Cintio

Abstract: Aims: It is well known that the presence of baryons affects the dark matter host haloes. Exploring the galaxy assembly history together with the dark matter haloes properties through time can provide a way to measure these effects. Methods: We study the properties of four Milky Way mass dark matter haloes from the Aquarius project during their assembly history, between $z = 0 - 4$. In this work, w… ▽ More Aims: It is well known that the presence of baryons affects the dark matter host haloes. Exploring the galaxy assembly history together with the dark matter haloes properties through time can provide a way to measure these effects. Methods: We study the properties of four Milky Way mass dark matter haloes from the Aquarius project during their assembly history, between $z = 0 - 4$. In this work, we use the SPH run from Scannapieco et al. (2009) and the dark matter only counterpart as case studies. To asses the robustness of our findings, we compare them with one of the haloes run using a moving-mesh technique and different sub-grid scheme. Results: Our results show that the cosmic evolution of the dark matter halo profiles depends on the assembly history of the baryons. We find that the dark matter profiles do not significantly change with time, hence they become stable, when the fraction of baryons accumulated in the central regions reaches 80 percent of its present mass within the virial radius. Furthermore, the mass accretion history shows that the haloes that assembled earlier are those that contain a larger amount of baryonic mass aforetime, which in turn allows the dark matter halo profiles to reach a stable configuration earlier. For the SPH haloes, we find that the specific angular momentum of the dark matter particles within the five percent of the virial radius at z = 0, remains approximately constant from the time at which 60 percent of the stellar mass is gathered. We explore different theoretical and empirical models for the contraction of the haloes through redshift. A model to better describe the contraction of the haloes through redshift evolution must depend on the stellar mass content in the inner regions. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: 17 pages, 14 figures. Accepted for publication in A&A. Comments are welcome

Journal ref: A&A 622, A197 (2019)

Formation of ultra-diffuse galaxies in the field and in galaxy groups

Authors: Fangzhou Jiang, Avishai Dekel, Jonathan Freundlich, Aaron J. Romanowsky, Aaron Dutton, Andrea Maccio, Arianna Di Cintio

Abstract: We study ultra-diffuse galaxies (UDGs) in zoom in cosmological simulations, seeking the origin of UDGs in the field versus galaxy groups. We find that while field UDGs arise from dwarfs in a characteristic mass range by multiple episodes of supernova feedback (Di Cintio et al. 2017), group UDGs may also form by tidal puffing up and they become quiescent by ram-pressure stripping. The field and gro… ▽ More We study ultra-diffuse galaxies (UDGs) in zoom in cosmological simulations, seeking the origin of UDGs in the field versus galaxy groups. We find that while field UDGs arise from dwarfs in a characteristic mass range by multiple episodes of supernova feedback (Di Cintio et al. 2017), group UDGs may also form by tidal puffing up and they become quiescent by ram-pressure stripping. The field and group UDGs share similar properties, independent of distance from the group centre. Their dark-matter haloes have ordinary spin parameters and centrally dominant dark-matter cores. Their stellar components tend to have a prolate shape with a Sersic index n~1 but no significant rotation. Ram pressure removes the gas from the group UDGs when they are at pericentre, quenching star formation in them and making them redder. This generates a colour/star-formation-rate gradient with distance from the centre, as observed in clusters. We find that ~20 per cent of the field UDGs that fall into a massive halo survive as satellite UDGs. In addition, normal field dwarfs on highly eccentric orbits can become UDGs near pericentre due to tidal puffing up, contributing about half of the group-UDG population. We interpret our findings using simple toy models, showing that gas stripping is mostly due to ram pressure rather than tides. We estimate that the energy deposited by tides in the bound component of a satellite over one orbit can cause significant puffing up provided that the orbit is sufficiently eccentric. △ Less

Submitted 26 November, 2018; originally announced November 2018.

Comments: 20 pages, 13 figures, submitted to MNRAS

CLUES about M33: the reversed radial stellar age gradient in the outskirts of Triangulum galaxy

Authors: Robert Mostoghiu, Arianna Di Cintio, Alexander Knebe, Noam I. Libeskind, Ivan Minchev, Chris Brook

Abstract: HST/ACS observations along the major axis of M33 show that the mean age of its stars decreases with increasing distance from the galaxy center. Such a behavior is consistent with an inside-out growth of the disc. However, in the outermost observed field, at r$\simeq$11.6 kpc, a reversal of this gradient is detected, with old stars found in high percentages beyond this radius. In this work we inves… ▽ More HST/ACS observations along the major axis of M33 show that the mean age of its stars decreases with increasing distance from the galaxy center. Such a behavior is consistent with an inside-out growth of the disc. However, in the outermost observed field, at r$\simeq$11.6 kpc, a reversal of this gradient is detected, with old stars found in high percentages beyond this radius. In this work we investigate the origin of such a reversal in stellar age gradient, by using a simulated M33 analogue from the Constrained Local UniversE Simulations (CLUES). The simulated M33 is similar to the observed one in terms of mass, rotation velocity, surface brightness and, similar to what has been reported in observations, shows a stellar age turnaround at large radii. We demonstrate that this reversal is mostly a result of stellar accretion from old satellite galaxies and, to a lesser extent, of stellar migration of in-situ stars. The old accreted stars, with formation times tf < 4 Gyrs, are kinematically hot and can be differentiated from the in-situ stars by their high velocity dispersion and the fact that they do not have rotationally-supported orbits. In the future, obtaining kinematic information of the stars in the outskirt of M33 will help to verify this scenario. △ Less

Submitted 9 August, 2018; originally announced August 2018.

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

Journal ref: MNRAS 480, Issue 4 (2018) 4455-4467

Stellar Feedback and the Energy Budget of Late-Type Galaxies: Missing Baryons and Core Creation

Authors: Harley Katz, Harry Desmond, Federico Lelli, Stacy McGaugh, Arianna Di Cintio, Chris Brook, James Schombert

Abstract: In a $Λ$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems.… ▽ More In a $Λ$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems. Assuming that haloes initially accrete the cosmological baryon fraction, we show that the maximum energy available from supernovae is typically not enough to completely eject these "missing baryons" from a halo, but it is often sufficient to heat them to the virial temperature. The energy available from supernovae has the same scaling with galaxy mass as the energy needed to heat or eject the "missing baryons", indicating that the coupling efficiency of the feedback to the ISM may be constant with galaxy virial mass. We further find that the energy available from supernova feedback is always enough to convert a primordial cusp into a core and has magnitude consistent with what is required to heat the "missing baryons" to the virial temperature. Taking a census of the baryon content of galaxies with ${\rm 10^9<M_{vir}/M_{\odot}<10^{12}}$ reveals that $\sim86\%$ of baryons are likely to be in a hot phase surrounding the galaxies and possibly observable in the X-ray, $\sim7\%$ are in the form of cold gas, and $\sim7\%$ are in stars. △ Less

Submitted 2 August, 2018; originally announced August 2018.

Comments: 16 pages, 8 figures, MNRAS accepted

A distance of 13 Mpc resolves the claimed anomalies of the galaxy lacking dark matter

Authors: Ignacio Trujillo, Michael A. Beasley, Alejandro Borlaff, Eleazar R. Carrasco, Arianna Di Cintio, Mercedes Filho, Matteo Monelli, Mireia Montes, Javier Roman, Tomas Ruiz-Lara, Jorge Sanchez Almeida, David Valls-Gabaud, Alexandre Vazdekis

Abstract: The claimed detection of a diffuse galaxy lacking dark matter represents a possible challenge to our understanding of the properties of these galaxies and galaxy formation in general. The galaxy, already identified in photographic plates taken in the summer of 1976 at the UK 48-in Schmidt telescope, presents normal distance-independent properties (e.g. colour, velocity dispersion of its globular c… ▽ More The claimed detection of a diffuse galaxy lacking dark matter represents a possible challenge to our understanding of the properties of these galaxies and galaxy formation in general. The galaxy, already identified in photographic plates taken in the summer of 1976 at the UK 48-in Schmidt telescope, presents normal distance-independent properties (e.g. colour, velocity dispersion of its globular clusters). However, distance-dependent quantities are at odds with those of other similar galaxies, namely the luminosity function and sizes of its globular clusters, mass-to-light ratio and dark matter content. Here we carry out a careful analysis of all extant data and show that they consistently indicate a much shorter distance (13 Mpc) than previously indicated (20 Mpc). With this revised distance, the galaxy appears to be a rather ordinary low surface brightness galaxy (R_e=1.4+-0.1 kpc; M*=6.0+-3.6x10^7 Msun) with plenty of room for dark matter (the fraction of dark matter inside the half mass radius is >75% and M_halo/M*>20) corresponding to a minimum halo mass >10^9 Msun. At 13 Mpc, the luminosity and structural properties of the globular clusters around the object are the same as those found in other galaxies. △ Less

Submitted 15 March, 2019; v1 submitted 26 June, 2018; originally announced June 2018.

Comments: MNRAS, in press. This version includes an extra new independent distance indicator confirming the 13 Mpc distance. New tests also reject the possibility that blends are causing a misidentification of the location of the TRGB

On the early evolution of Local Group dwarf galaxy types: star formation and supernova feedback

Authors: Jose R. Bermejo-Climent, Giuseppina Battaglia, Carme Gallart, Arianna Di Cintio, Chris B. Brook, Luis Cicuendez, Matteo Monelli, Ryan Leaman, Lucio Mayer, Jorge Penarrubia, Justin I. Read

Abstract: According to star formation histories (SFHs), Local Group dwarf galaxies can be broadly classified in two types: those forming most of their stars before $z=2$ (${\it fast}$) and those with more extended SFHs (${\it slow}$). The most precise SFHs are usually derived from deep but not very spatially extended photometric data; this might alter the ratio of old to young stars when age gradients are p… ▽ More According to star formation histories (SFHs), Local Group dwarf galaxies can be broadly classified in two types: those forming most of their stars before $z=2$ (${\it fast}$) and those with more extended SFHs (${\it slow}$). The most precise SFHs are usually derived from deep but not very spatially extended photometric data; this might alter the ratio of old to young stars when age gradients are present. Here we correct for this effect and derive the mass formed in stars by $z=2$ for a sample of 16 Local Group dwarf galaxies. We explore early differences between ${\it fast}$ and ${\it slow}$ dwarfs, and evaluate the impact of internal feedback by supernovae (SN) on the baryonic and dark matter (DM) component of the dwarfs. ${\it Fast}$ dwarfs assembled more stellar mass at early times and have larger amounts of DM within the half-light radius than ${\it slow}$ dwarfs. By imposing that ${\it slow}$ dwarfs cannot have lost their gas by $z=2$, we constrain the maximum coupling efficiency of SN feedback to the gas and to the DM to be $\sim$10%. We find that internal feedback alone appears insufficient to quench the SFH of ${\it fast}$ dwarfs by gas deprivation, in particular for the fainter systems. Nonetheless, SN feedback can core the DM halo density profiles relatively easily, producing cores of the sizes of the half-light radius in ${\it fast}$ dwarfs by $z=2$ with very low efficiencies. Amongst the "classical" Milky Way satellites, we predict that the smallest cores should be found in Draco and Ursa Minor, while Sculptor and Fornax should host the largest ones. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: 15 pages, 3 figures, 2 tables. Accepted for publication in MNRAS

Spectroscopic characterisation of the stellar content of ultra diffuse galaxies

Authors: T. Ruiz-Lara, M. A. Beasley, J. Falcón-Barroso, J. Román, F. Pinna, C. Brook, A. Di Cintio, I. Martín-Navarro, I. Trujillo, A. Vazdekis

Abstract: Understanding the peculiar properties of Ultra Diffuse Galaxies (UDGs) via spectroscopic analysis is a challenging task requiring very deep observations and exquisite data reduction. In this work we perform one of the most complete characterisations of the stellar component of UDGs to date using deep optical spectroscopic data from OSIRIS at GTC. We measure radial and rotation velocities, star for… ▽ More Understanding the peculiar properties of Ultra Diffuse Galaxies (UDGs) via spectroscopic analysis is a challenging task requiring very deep observations and exquisite data reduction. In this work we perform one of the most complete characterisations of the stellar component of UDGs to date using deep optical spectroscopic data from OSIRIS at GTC. We measure radial and rotation velocities, star formation histories (SFH) and mean population parameters, such as ages and metallicities, for a sample of five UDG candidates in the Coma cluster. From the radial velocities, we confirm the Coma membership of these galaxies. We find that their rotation properties, if detected at all, are compatible with dwarf-like galaxies. The SFHs of the UDG are dominated by old (~ 7 Gyr), metal-poor ([M/H] ~ -1.1) and alpha-enhanced ([Mg/Fe] ~ 0.4) populations followed by a smooth or episodic decline which halted ~ 2 Gyr ago, possibly a sign of cluster-induced quenching. We find no obvious correlation between individual SFH shapes and any UDG morphological properties. The recovered stellar properties for UDGs are similar to those found for DDO44, a local UDG analogue resolved into stars. We conclude that the UDGs in our sample are extended dwarfs whose properties are likely the outcome of both internal processes, such as bursty SFHs and/or high-spin haloes, as well as environmental effects within the Coma cluster. △ Less

Submitted 6 July, 2018; v1 submitted 16 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS

NIHAO XIV: Reproducing the observed diversity of dwarf galaxy rotation curve shapes in LCDM

Authors: Isabel M. Santos-Santos, Arianna Di Cintio, Chris B. Brook, Andrea Macciò, Aaron Dutton, Rosa Domínguez-Tenreiro

Abstract: The significant diversity of rotation curve (RC) shapes in dwarf galaxies has recently emerged as a challenge to LCDM: in dark matter (DM) only simulations, DM halos have a universal cuspy density profile that results in self-similar RC shapes. We compare RC shapes of simulated galaxies from the NIHAO project with observed galaxies from the homogeneous SPARC dataset. The DM halos of the NIHAO gala… ▽ More The significant diversity of rotation curve (RC) shapes in dwarf galaxies has recently emerged as a challenge to LCDM: in dark matter (DM) only simulations, DM halos have a universal cuspy density profile that results in self-similar RC shapes. We compare RC shapes of simulated galaxies from the NIHAO project with observed galaxies from the homogeneous SPARC dataset. The DM halos of the NIHAO galaxies can expand to form cores, with the degree of expansion depending on their stellar-to-halo mass ratio. By means of the $V_{\rm 2kpc}-V_{\rm Rlast}$ relation (where $V_{\rm Rlast}$ is the outermost measured rotation velocity), we show that both the average trend and the scatter in RC shapes of NIHAO galaxies are in reasonable agreement with SPARC: this represents a significant improvement compared to simulations that do not result in DM core formation, suggesting that halo expansion is a key process in matching the diversity of dwarf galaxy RCs. Note that NIHAO galaxies can reproduce even the extremely slowly rising RCs of IC 2574 and UGC 5750. Revealingly, the range where observed galaxies show the highest diversity corresponds to the range where core formation is most efficient in NIHAO simulations, 50$<V_{\rm Rlast}$\km s$^{-1}<$100. A few observed galaxies in this range cannot be matched by any NIHAO RC nor by simulations that predict a universal halo profile. Interestingly, the majority of these are starbursts or emission-line galaxies, with steep RCs and small effective radii. Such galaxies represent an interesting observational target providing new clues to the process/viability of cusp-core transformation, the relationship between starburst and inner potential well, and the nature of DM. △ Less

Submitted 16 October, 2017; v1 submitted 13 June, 2017; originally announced June 2017.

Comments: 13 pages, 6 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society 9 October 2017

A rumble in the dark: signatures of self-interacting dark matter in Super-Massive Black Hole dynamics and galaxy density profiles

Authors: Arianna Di Cintio, Michael Tremmel, Fabio Governato, Andrew Pontzen, Jesús Zavala, Alexander Bastidas Fry, Alyson Brooks, Mark Vogelsberger

Abstract: We explore for the first time the effect of self-interacting dark matter (SIDM) on the dark matter (DM) and baryonic distribution in massive galaxies formed in hydrodynamical cosmological simulations, including explicit baryonic physics treatment. A novel implementation of Super-Massive Black Hole (SMBH) formation and evolution is used, as in Tremmel et al.(2015, 2016), allowing to explicitly foll… ▽ More We explore for the first time the effect of self-interacting dark matter (SIDM) on the dark matter (DM) and baryonic distribution in massive galaxies formed in hydrodynamical cosmological simulations, including explicit baryonic physics treatment. A novel implementation of Super-Massive Black Hole (SMBH) formation and evolution is used, as in Tremmel et al.(2015, 2016), allowing to explicitly follow SMBH dynamics at the center of galaxies. A high SIDM constant cross-section is chosen, $σ$=10 $\rm cm^2/gr$, to amplify differences from CDM models. Milky Way-like galaxies form a shallower DM density profile in SIDM than they do in CDM, with differences already at 20 kpc scales. This demonstrates that even for the most massive spirals the effect of SIDM dominates over the adiabatic contraction due to baryons. Strikingly, the dynamics of SMBHs differs in the SIDM and reference CDM case. SMBHs in massive spirals have sunk to the centre of their host galaxy in both the SIDM and CDM run, while in less massive galaxies about 80$\%$ of the SMBH population is off-centered in the SIDM case, as opposed to the CDM case in which $\sim$90$\%$ of SMBHs have reached their host's centre. SMBHs are found as far as $\sim$9 kpc away from the centre of their host SIDM galaxy. This difference is due to the increased dynamical friction timescale caused by the lower DM density in SIDM galaxies compared to CDM, resulting in 'core stalling'. This pilot work highlights the importance of simulating in a full hydrodynamical context different DM models combined to SMBH physics to study their influence on galaxy formation. △ Less

Submitted 16 May, 2017; v1 submitted 16 January, 2017; originally announced January 2017.

Comments: 10 pages, 6 figures, matches accepted version (MNRAS, April 27 2017)

NIHAO XI: Formation of Ultra-Diffuse Galaxies by outflows

Authors: Arianna Di Cintio, Chris B. Brook, Aaron A. Dutton, Andrea V. Macciò, Aura C. Obreja, Avishai Dekel

Abstract: We address the origin of Ultra-Diffuse Galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and whether they are failed $\rm L_{\star}$ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the NIHAO project, we show that UDG analogues form naturally in medium-mass h… ▽ More We address the origin of Ultra-Diffuse Galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and whether they are failed $\rm L_{\star}$ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the NIHAO project, we show that UDG analogues form naturally in medium-mass haloes due to episodes of gas outflows associated with star formation. The simulated UDGs live in isolated haloes of masses $10^{10-11}\rm M_{\odot}$, have stellar masses of $10^{7-8.5}\rm M_{\odot}$, effective radii larger than 1 kpc and dark matter cores. They show a broad range of colors, an average Sérsic index of 0.83, a typical distribution of halo spin and concentration, and a non-negligible HI gas mass of $10^{7-9}\rm M_{\odot}$, which correlates with the extent of the galaxy. Gas availability is crucial to the internal processes that form UDGs: feedback driven gas outflows, and subsequent dark matter and stellar expansion, are the key to reproduce faint, yet unusually extended, galaxies. This scenario implies that UDGs represent a dwarf population of low surface brightness galaxies and should exist in the field. The largest isolated UDGs should contain more HI gas than less extended dwarfs of similar $\rm M_{\star}$. △ Less

Submitted 22 November, 2016; v1 submitted 3 August, 2016; originally announced August 2016.

Comments: matches accepted version, MNRAS Letter 2016-10-12

Testing Feedback-Modified Dark Matter Haloes with Galaxy Rotation Curves: Estimation of Halo Parameters and Consistency with $Λ$CDM

Authors: Harley Katz, Federico Lelli, Stacy S. McGaugh, Arianna Di Cintio, Chris B. Brook, James M. Schombert

Abstract: Cosmological $N$-body simulations predict dark matter (DM) haloes with steep central cusps (e.g. NFW, Navarro et al. 1996). This contradicts observations of gas kinematics in low-mass galaxies that imply the existence of shallow DM cores. Baryonic processes such as adiabatic contraction and gas outflows can, in principle, alter the initial DM density profile, yet their relative contributions to th… ▽ More Cosmological $N$-body simulations predict dark matter (DM) haloes with steep central cusps (e.g. NFW, Navarro et al. 1996). This contradicts observations of gas kinematics in low-mass galaxies that imply the existence of shallow DM cores. Baryonic processes such as adiabatic contraction and gas outflows can, in principle, alter the initial DM density profile, yet their relative contributions to the halo transformation remain uncertain. Recent high resolution, cosmological hydrodynamic simulations (Di Cintio et al. 2014, DC14) predict that inner density profiles depend systematically on the ratio of stellar to DM mass (M$_*$/M$_{\text{halo}}$). Using a Markov Chain Monte Carlo approach, we test the NFW and the M$_*$/M$_{\text{halo}}$-dependent DC14 halo models against a sample of 147 galaxy rotation curves from the new {\it Spitzer} Photometry and Accurate Rotation Curves (SPARC) data set. These galaxies all have extended H{\small I} rotation curves from radio interferometry as well as accurate stellar mass density profiles from near-infrared photometry. The DC14 halo profile provides markedly better fits to the data compared to the NFW profile. Unlike NFW, the DC14 halo parameters found in our rotation curve fits naturally fall within two standard deviations of the mass-concentration relation predicted by $Λ$CDM and the stellar mass-halo mass relation inferred from abundance matching with few outliers. Halo profiles modified by baryonic processes are therefore more consistent with expectations from $Λ$ cold dark matter ($Λ$CDM) cosmology and provide better fits to galaxy rotation curves across a wide range of galaxy properties than do halo models that neglect baryonic physics. Our results offer a solution to the decade long cusp-core discrepancy. △ Less

Submitted 11 January, 2017; v1 submitted 19 May, 2016; originally announced May 2016.

Comments: 23 Pages, 18 Figures, MNRAS Accepted

NIHAO IX: the role of gas inflows and outflows in driving the contraction and expansion of cold dark matter haloes

Authors: Aaron A. Dutton, Andrea V. Macciò, Avishai Dekel, Liang Wang, Gregory S. Stinson, Aura Obreja, Arianna Di Cintio, Chris B. Brook, Tobias Buck, Xi Kang

Abstract: We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with… ▽ More We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(Ω_b/Ω_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes. △ Less

Submitted 5 July, 2016; v1 submitted 17 May, 2016; originally announced May 2016.

Comments: 20 pages, 15 figures, accepted to MNRAS

The mass discrepancy acceleration relation in a $Λ$CDM context

Authors: Arianna Di Cintio, Federico Lelli

Abstract: The mass discrepancy acceleration relation (MDAR) describes the coupling between baryons and dark matter (DM) in galaxies: the ratio of total-to-baryonic mass at a given radius anti-correlates with the acceleration due to baryons. The MDAR has been seen as a challenge to the $Λ$CDM galaxy formation model, while it can be explained by Modified Newtonian Dynamics. In this Letter we show that the MDA… ▽ More The mass discrepancy acceleration relation (MDAR) describes the coupling between baryons and dark matter (DM) in galaxies: the ratio of total-to-baryonic mass at a given radius anti-correlates with the acceleration due to baryons. The MDAR has been seen as a challenge to the $Λ$CDM galaxy formation model, while it can be explained by Modified Newtonian Dynamics. In this Letter we show that the MDAR arises in a $Λ$CDM cosmology once observed galaxy scaling relations are taken into account. We build semi-empirical models based on $Λ$CDM haloes, with and without the inclusion of baryonic effects, coupled to empirically motivated structural relations. Our models can reproduce the MDAR: specifically, a mass-dependent density profile for DM haloes can fully account for the observed MDAR shape, while a universal profile shows a discrepancy with the MDAR of dwarf galaxies with $\rm M^{\star}$$<$$\rm10^{9.5}M_{\odot}$, a further indication suggesting the existence of DM cores. Additionally, we reproduce slope and normalization of the baryonic Tully-Fisher relation (BTFR) with 0.17 dex scatter. These results imply that in $Λ$CDM (i) the MDAR is driven by structural scaling relations of galaxies and DM density profile shapes, and (ii) the baryonic fractions determined by the BTFR are consistent with those inferred from abundance-matching studies. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: 5 pages,4 figures. Accepted for publication in MNRAS Letter

The distribution of mass components in simulated disc galaxies

Authors: Isabel M. Santos-Santos, Chris B. Brook, Greg Stinson, Arianna Di Cintio, James Wadsley, Rosa Domínguez-Tenreiro, Stefan Gottlöber, Gustavo Yepes

Abstract: Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V… ▽ More Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V$_{\rm flat}$ is the circular velocity at the flat part of the rotation curve, and match galaxy scaling relations, are able to recover the observed relations between the distributions of stars, gas and dark matter over the radial range for which we have observational rotation curve data. Furthermore, we explicitly match the observed baryonic to halo mass relation for the first time with simulated galaxies. We discuss our results in the context of the baryon cycle that is inherent in these simulations, and with regards to the effect of baryonic processes on the distribution of dark matter. △ Less

Submitted 13 October, 2015; v1 submitted 8 October, 2015; originally announced October 2015.

Comments: 8 pages, 7 pdf figures. Accepted for publication in MNRAS on 2015 October 06

NIHAO IV: Core creation and destruction in dark matter density profiles across cosmic time

Authors: Edouard Tollet, Andrea V. Macciò, Aaron A. Dutton, Greg S. Stinson, Liang Wang, Camilla Penzo, Thales A. Gutcke, Tobias Buck, Xi Kang, Chris Brook, Arianna Di Cintio, Ben W. Keller, James Wadsley

Abstract: We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter densi… ▽ More We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $α$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $α$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $α$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation. △ Less

Submitted 6 May, 2019; v1 submitted 13 July, 2015; originally announced July 2015.

Comments: 11 pages, 13 figures. Corrected typo in table 2 (middle row) with respect to the version published in MNRAS

Signatures of Dark Matter Halo Expansion in Galaxy Populations

Authors: Chris B. Brook, Arianna Di Cintio

Abstract: Dark matter cores within galaxy haloes can be formed by energy feedback from star forming regions: an energy balance suggests that the maximum core formation efficiency arises in galaxies with M$_{\star}\sim10^{8.5}$M$_{\odot}$. We show that a model population of galaxies, in which the density profile has been modified by such baryonic feedback, is able to explain the observed galaxy velocity func… ▽ More Dark matter cores within galaxy haloes can be formed by energy feedback from star forming regions: an energy balance suggests that the maximum core formation efficiency arises in galaxies with M$_{\star}\sim10^{8.5}$M$_{\odot}$. We show that a model population of galaxies, in which the density profile has been modified by such baryonic feedback, is able to explain the observed galaxy velocity function and Tully-Fisher relations significantly better than a model in which a universal cuspy density profile is assumed. Alternative models, namely warm or self-interacting dark matter, also provide a better match to these observed relations than a universal profile model does, but make different predictions for how halo density profiles vary with mass compared to the baryonic feedback case. We propose that different core formation mechanisms may be distinguished based on the imprint they leave on galaxy populations over a wide range of mass. Within the current observational data we find evidence of the expected signatures of the mass dependence of core formation generated by baryonic feedback. △ Less

Submitted 30 July, 2015; v1 submitted 2 June, 2015; originally announced June 2015.

Comments: 11 pages, 9 figures, accepted in MNRAS

Expanded haloes, abundance matching and too-big-to-fail in the Local Group

Authors: Chris B. Brook, Arianna Di Cintio

Abstract: Observed kinematical data of 40 Local Group (LG) members are used to derive the dark matter halo mass of such galaxies. Haloes are selected from the theoretically expected LG mass function and two different density profiles are assumed, a standard universal cuspy model and a mass dependent profile which accounts for the effects of baryons in modifying the dark matter distribution within galaxies.… ▽ More Observed kinematical data of 40 Local Group (LG) members are used to derive the dark matter halo mass of such galaxies. Haloes are selected from the theoretically expected LG mass function and two different density profiles are assumed, a standard universal cuspy model and a mass dependent profile which accounts for the effects of baryons in modifying the dark matter distribution within galaxies. The resulting relations between stellar and halo mass are compared with expectations from abundance matching. Using a universal cuspy profile, the ensemble of LG galaxies is fit in relatively low mass haloes, leaving "dark" many massive haloes of \mhalo$\gtrsim$10$^{10}$\msun: this reflects the "too big to fail" problem and results in a \mstar-\mhalo\ relation that differs from abundance matching predictions. Moreover, the star formation efficiency of isolated LG galaxies increases with decreasing halo mass when adopting a cuspy model. By contrast, using the mass dependent density profile, dwarf galaxies with \mstar$\gtrsim$10$^{6}$\msun are assigned to more massive haloes, which have a central cored distribution of dark matter: the "too big to fail" problem is alleviated, the resultant \mstar-\mhalo\ relation follows abundance matching predictions down to the completeness limit of current surveys, and the star formation efficiency of isolated members decreases with decreasing halo mass, in agreement with theoretical expectations. Finally, the cusp/core space of LG galaxies is presented, providing a framework to understand the non-universality of their density profiles. △ Less

Submitted 15 April, 2015; v1 submitted 14 October, 2014; originally announced October 2014.

Comments: Accepted in MNRAS, 15 pages 7 figures. section 3.3 has been added after reports

A mass-dependent density profile for dark matter haloes including the influence of galaxy formation

Authors: Arianna Di Cintio, Chris B. Brook, Aaron A. Dutton, Andrea V. Macciò, Greg S. Stinson, Alexander Knebe

Abstract: We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range… ▽ More We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity Mstar/Mhalo constrains the inner ($γ$) and outer ($β$) slopes of dark matter density, and the sharpness of transition between the slopes($α$), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs. This mass dependent density profile can be directly applied to rotation curve data of observed galaxies and to semi analytic galaxy formation models as a significant improvement over the commonly used NFW profile. △ Less

Submitted 23 April, 2014; originally announced April 2014.

Comments: 10 pages, 4 figures. Accepted for publication in MNRAS

Galaxy formation

Authors: Joseph Silk, Arianna Di Cintio, Irina Dvorkin

Abstract: Galaxy formation is at the forefront of observation and theory in cosmology. An improved understanding is essential for improving our knowledge both of the cosmological parameters, of the contents of the universe, and of our origins. In these lectures intended for graduate students, galaxy formation theory is reviewed and confronted with recent observational issues. In Lecture 1, the following top… ▽ More Galaxy formation is at the forefront of observation and theory in cosmology. An improved understanding is essential for improving our knowledge both of the cosmological parameters, of the contents of the universe, and of our origins. In these lectures intended for graduate students, galaxy formation theory is reviewed and confronted with recent observational issues. In Lecture 1, the following topics are presented: star formation considerations, including IMF, star formation efficiency and star formation rate, the origin of the galaxy luminosity function, and feedback in dwarf galaxies. In Lecture 2, we describe formation of disks and massive spheroids, including the growth of supermassive black holes, negative feedback in spheroids, the AGN-star formation connection, star formation rates at high redshift and the baryon fraction in galaxies. △ Less

Submitted 30 November, 2013; originally announced December 2013.

Comments: Lectures given at Post-Planck Cosmology, Ecole de Physique des Houches, Les Houches, July 8-Aug 2, 2013, eds. B. Wandelt, C. Deffayet, P. Peter, to be published by Oxford University Press, and New Horizons for Observational Cosmology, International School of Physics Enrico Fermi, Varenna, July 1-6, 2013, eds. A. Melchiorri, A. Cooray, E. Komatsu, to be published by the Italian Society of Physics

The stellar-to-halo mass relation for Local Group galaxies

Authors: C. B. Brook, A. Di Cintio, A. Knebe, S. Gottlöber, Y. Hoffman, G. Yepes, S. Garrison-Kimmel

Abstract: We contend that a single power law halo mass distribution is appropriate for direct matching to the stellar masses of observed Local Group dwarf galaxies, allowing the determination of the slope of the stellar mass-halo mass relation for low mass galaxies. Errors in halo masses are well defined as the Poisson noise of simulated local group realisations, which we determine using constrained local u… ▽ More We contend that a single power law halo mass distribution is appropriate for direct matching to the stellar masses of observed Local Group dwarf galaxies, allowing the determination of the slope of the stellar mass-halo mass relation for low mass galaxies. Errors in halo masses are well defined as the Poisson noise of simulated local group realisations, which we determine using constrained local universe simulations (CLUES). For the stellar mass range 10$^7$<M*<10$^8$M$_\odot$, for which we likely have a complete census of observed galaxies, we find that the stellar mass-halo mass relation follows a power law with slope of 3.1, significantly steeper than most values in the literature. The steep relation between stellar and halo masses indicates that Local Group dwarf galaxies are hosted by dark matter halos with a small range of mass. Our methodology is robust down to the stellar mass to which the census of observed Local Group galaxies is complete, but the significant uncertainty in the currently measured slope of the stellar-to halo mass relation will decrease dramatically if the Local Group completeness limit was $10^{6.5}$M$\odot$ or below, highlighting the importance of pushing such limit to lower masses and larger volumes. △ Less

Submitted 1 February, 2014; v1 submitted 21 November, 2013; originally announced November 2013.

Comments: ApJL Accepted version. Comments to cbabrook@gmail.com

The dependence of dark matter profiles on the stellar to halo mass ratio: a prediction for cusps vs cores

Authors: Arianna Di Cintio, Chris B. Brook, Andrea V. Maccio, Greg S. Stinson, Alexander Knebe, Aaron A. Dutton, James Wadsley

Abstract: We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks… ▽ More We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks. The main result is a clear dependence of the inner slope of the DM density profile, α in ρr^α, on the ratio between stellar-to-halo mass (M*/Mhalo). This relation is independent of the stellar feedback scheme, allowing a prediction for cusp vs core formation. When M*/Mhalo is low, ~0.01%, energy from stellar feedback is insufficient to significantly alter the inner DM density and the galaxy retains a cuspy profile. At higher M*/Mhalo, feedback drives the expansion of the DM and generates cored profiles. The flattest profiles form where M*/Mhalo~0.5%. Above this ratio, stars formed in the central regions deepen the gravitational potential enough to oppose this supernova-driven expansion process, resulting in smaller cores and cuspier profiles. Combining the dependence of α on M*/Mhalo with the abundance matching relation between M* and Mhalo provides a prediction for how α varies with M*. Further, using the Tully-Fisher relation allows a prediction for the dependence of the DM inner slope on the observed rotation velocity of galaxies. The most cored galaxies are expected to have Vrot~50km/s, with α decreasing for more massive disc galaxies: spirals with Vrot~150km/s have central slopes α<-0.8, approaching the NFW profile. This novel prediction for the dependence of α on disc galaxy mass can be tested using current observational data sets, and can be applied to theoretical modeling of mass profiles and populations of disc galaxies. △ Less

Submitted 7 October, 2013; v1 submitted 4 June, 2013; originally announced June 2013.

Comments: accepted version, MNRAS in press. results unchanged,expanded discussion in sec 3.1,3.2,conclusions

Cold versus Warm Dark Matter simulations of a galaxy group

Authors: Noam I. Libeskind, Arianna Di Cintio, Alexander Knebe, Gustavo Yepes, Stefan Gottloeber, Matthias Steinmetz, Yehuda Hoffman, Luis A. Martinez-Vaquero

Abstract: The differences between cold (CDM) and warm (WDM) dark matter in the formation of a group of galaxies is examined by running two identical simulations where in the WDM case the initial power spectrum has been altered to mimic a 1keV dark matter particle. The CDM initial conditions were constrained to reproduce at z = 0 the correct local environment within which a "Local Group" (LG) of galaxies may… ▽ More The differences between cold (CDM) and warm (WDM) dark matter in the formation of a group of galaxies is examined by running two identical simulations where in the WDM case the initial power spectrum has been altered to mimic a 1keV dark matter particle. The CDM initial conditions were constrained to reproduce at z = 0 the correct local environment within which a "Local Group" (LG) of galaxies may form. Two significant differences between the two simulations are found. While in the CDM case a group of galaxies that resembles the real LG forms, the WDM run fails to reproduce a viable LG, instead forming a diffuse group which is still expanding at z = 0. This is surprising since, due to the suppression of small scale power in its power spectrum, WDM is naively expected to only affect the collapse of small haloes and not necessarily the dynamics on a scale of a group of galaxies. Furthermore the concentration of baryons in halo center's is greater in CDM than in WDM and the properties of the disks differ. △ Less

Submitted 23 May, 2013; originally announced May 2013.

Comments: 13 Pages, 5 figures, Accepted for publication in Publications of the Astronomical Society of Australia special issue on Warm Dark Matter

Size matters: the non-universal density profile of subhaloes in SPH simulations and implications for the Milky Way's dSphs

Authors: Arianna Di Cintio, Alexander Knebe, Noam I. Libeskind, Chris Brook, Gustavo Yepes, Stefan Gottloeber, Yehuda Hoffman

Abstract: We use dark matter only and full hydrodynamical Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group to study the density profile of subhaloes of the simulated Milky Way and Andromeda galaxies. We show that the Einasto model provides the best description of the subhaloes' density profile, as opposed to the more commonly used NFW profile or any generalisation of it. We… ▽ More We use dark matter only and full hydrodynamical Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group to study the density profile of subhaloes of the simulated Milky Way and Andromeda galaxies. We show that the Einasto model provides the best description of the subhaloes' density profile, as opposed to the more commonly used NFW profile or any generalisation of it. We further find that the Einasto shape parameter \nEin\ is strongly correlated with the total subhalo mass, pointing towards the notion of a non-universality of the subhaloes' density profile. We observe that the effect of mass loss due to tidal stripping, in both the dark matter only and the hydrodynamical run, is the reduction of the shape parameter \nEin\ between the infall and the present time. Assuming now that the dSphs of our Galaxy follow the Einasto profile and using the maximum and minimum values of \nEin\ from our hydrodynamical simulation as a gauge, we can improve the observational constraints on the \Rmax-\Vmax\ pairs obtained for the brightest satellite galaxies of the Milky Way. When considering only the subhaloes with $-13.2\lesssim M_V\lesssim-8.8$, i.e. the range of luminosity of the classical dwarfs, we find that all our simulated objects are consistent with the observed dSphs if their haloes follow the Einasto model with $1.6\lesssim n_{\rm E} \lesssim5.3$. The numerically motivated Einasto profile for the observed dSphs will alleviate the recently presented "massive failures" problem. △ Less

Submitted 25 February, 2013; v1 submitted 2 April, 2012; originally announced April 2012.

Comments: 11 pages, 5 figs, MNRAS accepted version, figures added, section 4 expanded

Applying scale-free mass estimators to the Local Group in Constrained Local Universe Simulations

Authors: Arianna Di Cintio, Alexander Knebe, Noam I. Libeskind, Yehuda Hoffman, Gustavo Yepes, Stefan Gottloeber

Abstract: We use the recently proposed scale-free mass estimators to determine the masses of the Milky Way (MW) and Andromeda (M31) galaxy in a dark matter only Constrained Local UniversE Simulation (CLUES). While these mass estimators work rather well for isolated spherical host systems, we examine here their applicability to a simulated binary system with a unique satellite population similar to the obser… ▽ More We use the recently proposed scale-free mass estimators to determine the masses of the Milky Way (MW) and Andromeda (M31) galaxy in a dark matter only Constrained Local UniversE Simulation (CLUES). While these mass estimators work rather well for isolated spherical host systems, we examine here their applicability to a simulated binary system with a unique satellite population similar to the observed satellites of MW and M31. We confirm that the scale-free estimators work also very well in our simulated Local Group galaxies with the right number of satellites which follow the observed radial distribution. In the isotropic case and under the assumption that the satellites are tracking the total gravitating mass, the power-law index of the radial satellite distribution $N(<r)\propto r^{3-γ}$ is directly related to the host's mass profile $M(<r)\propto r^{1-α}$ as $α=γ-2$. The use of this relation for any given $γ$ leads to highly accurate mass estimations which is a crucial point for observer, since they do not know a priori the mass profile of the MW and M31 haloes. We discuss possible bias in the mass estimators and conclude that the scale-free mass estimators can be satisfactorily applied to the real MW and M31 system. △ Less

Submitted 30 March, 2012; originally announced April 2012.

Comments: 14 pages, 6 figures, 6 tables. Accepted in MNRAS 2012 March 29. Received 2012 March 29; in original form 2011 September 29

Too small to succeed? Lighting up massive dark matter subhaloes of the Milky Way

Authors: Arianna Di Cintio, Alexander Knebe, Noam I. Libeskind, Gustavo Yepes, Stefan Gottloeber, Yehuda Hoffman

Abstract: Using Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group in a cosmological context we investigate the recently highlighted problem that the majority of the most massive dark subhaloes of the Milky Way are too dense to host any of its bright satellites. In particular, we examine the influence of baryonic processes and find that they leave a twofold effect on the rela… ▽ More Using Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group in a cosmological context we investigate the recently highlighted problem that the majority of the most massive dark subhaloes of the Milky Way are too dense to host any of its bright satellites. In particular, we examine the influence of baryonic processes and find that they leave a twofold effect on the relation between the peak of the rotation curve and its position (Vmax and Rmax). Satellites with a large baryon fraction experience adiabatic contraction thus decreasing Rmax while leaving Vmax more or less unchanged. Subhaloes with smaller baryon fractions undergo a decrease in Vmax possibly due to outflows of material. Furthermore, the situation of finding subhaloes in simulations that lie outside the confidence interval for possible hosts of the bright MW dwarf spheroidals, appears to be far more prominent in cosmologies with a high sigma8 normalisation and depends on the mass of the host. We conclude that the problem cannot be simply solved by including baryonic processes and hence demands further investigations. △ Less

Submitted 25 July, 2011; originally announced July 2011.

Comments: 6 pages, 3 figures. Accepted as MNRAS Letter 2011 July 21. Received 2011 July 20; in original form 2011 June 6