Showing 1–47 of 47 results for author: Abadi, M G

One size fits all: Insights on Extrinsic Thermal Absorption from Similarity of Supernova Remnant Radio Continuum Spectra

Authors: M. G. Abadi, G. Castelletti, L. Supan, N. E. Kassim, J. W. Lazio

Abstract: Typically, integrated radio frequency continuum spectra of supernova remnants (SNRs) exhibit a power-law form due to their synchrotron emission. In numerous cases, these spectra show an exponential turnover, long assumed due to thermal free-free absorption in the interstellar medium. We use a compilation of Galactic radio continuum SNR spectra, with and without turnovers, to constrain the distribu… ▽ More Typically, integrated radio frequency continuum spectra of supernova remnants (SNRs) exhibit a power-law form due to their synchrotron emission. In numerous cases, these spectra show an exponential turnover, long assumed due to thermal free-free absorption in the interstellar medium. We use a compilation of Galactic radio continuum SNR spectra, with and without turnovers, to constrain the distribution of the absorbing ionized gas. We introduce a novel parameterization of SNR spectra in terms of a characteristic frequency v_*, which depends both on the absorption turnover frequency and the power-law slope. Normalizing to v_* and to the corresponding flux density, S_*, we demonstrate that the stacked spectra of our sample reveal a similarity in behavior with low scatter (r.m.s. ~15%), and a unique exponential drop-off fully consistent with the predictions of a free-free absorption process. Observed SNRs, whether exhibiting spectral turnovers or not, appear to be spatially well mixed in the Galaxy without any evident segregation between them. Moreover, their Galactic distribution does not show a correlation with general properties such as heliocentric distance or Galactic longitude, as might have been expected if the absorption were due to a continuous distribution of ionized gas. However, it naturally arises if the absorbers are discretely distributed, as suggested by early low-frequency observations. Modeling based on HII regions tracking Galactic spiral arms successfully reproduces the patchy absorption observed to date. While more extensive statistical datasets should yield more precise spatial models of the absorbing gas distribution, our present conclusion regarding its inhomogeneity will remain robust. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: Accepted for A&A. Comments welcome. 11 pages, 11 figures, 1 table

Environmental effects on associations of dwarf galaxies

Authors: C. Yamila Yaryura, Mario G. Abadi, Stefan Gottlöber, Noam I. Libeskind, Sofía A. Cora, Andrés N. Ruiz, Cristian A. Vega-Martínez, Gustavo Yepes

Abstract: We study the properties of associations of dwarf galaxies and their dependence on the environment. Associations of dwarf galaxies are extended systems composed exclusively of dwarf galaxies, considering as dwarf galaxies those galaxies less massive than $M_{\star, \rm max} = 10^{9.0}$ ${\rm M}_{\odot}\,h^{-1}$. We identify these particular systems using a semi-analytical model of galaxy formation… ▽ More We study the properties of associations of dwarf galaxies and their dependence on the environment. Associations of dwarf galaxies are extended systems composed exclusively of dwarf galaxies, considering as dwarf galaxies those galaxies less massive than $M_{\star, \rm max} = 10^{9.0}$ ${\rm M}_{\odot}\,h^{-1}$. We identify these particular systems using a semi-analytical model of galaxy formation coupled to a dark matter only simulation in the $Λ$ Cold Dark Matter cosmological model. To classify the environment, we estimate eigenvalues from the tidal field of the dark matter particle distribution of the simulation. We find that the majority, two thirds, of associations are located in filaments ($ \sim 67$ per cent), followed by walls ($ \sim 26 $ per cent), while only a small fraction of them are in knots ($ \sim 6 $ per cent) and voids ($ \sim 1 $ per cent). Associations located in more dense environments present significantly higher velocity dispersion than those located in less dense environments, evidencing that the environment plays a fundamental role in their dynamical properties. However, this connection between velocity dispersion and the environment depends exclusively on whether the systems are gravitational bound or unbound, given that it disappears when we consider associations of dwarf galaxies that are gravitationally bound. Although less than a dozen observationally detected associations of dwarf galaxies are currently known, our results are predictions on the eve of forthcoming large surveys of galaxies, which will enable these very particular systems to be identified and studied. △ Less

Submitted 26 July, 2023; originally announced July 2023.

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

DS+: a method for the identification of cluster substructures

Authors: Jose A. Benavides, Andrea Biviano, Mario G. Abadi

Abstract: The study of cluster substructures is important for the determination of the cluster dynamical status, assembly history, and the evolution of cluster galaxies, and it allows to set of constraints on the nature of dark matter and cosmological parameters. We present and test DS+, a new method for the identification and characterization of group-sized substructures in clusters. Our new method is base… ▽ More The study of cluster substructures is important for the determination of the cluster dynamical status, assembly history, and the evolution of cluster galaxies, and it allows to set of constraints on the nature of dark matter and cosmological parameters. We present and test DS+, a new method for the identification and characterization of group-sized substructures in clusters. Our new method is based on the projected positions and line-of-sight velocities of cluster galaxies, and it is an improvement and extension of the traditional method of Dressler & Shectman (1988). We test it on cluster-size cosmological halos extracted from the IllustrisTNG simulations, with virial masses $\rm{14 \lesssim \log (M_{200}/M_{\odot}) \lesssim 14.6}$, that contain on average $\sim 190$ galaxies. We also present an application of our method on a real data set, the Bullet cluster. DS+ is able to identify $\sim 80\%$ of real group galaxies as members of substructures, and at least 60\% of the galaxies assigned to substructures belong to real groups. The physical properties of the real groups are significantly correlated with those of the corresponding detected substructures, albeit with significant scatter, and overestimated on average. Application of the DS+ method to the Bullet cluster confirms the presence and main properties of the high-speed collision and identifies other substructures along the main cluster axis. DS+ proves to be a reliable method for the identification of substructures in clusters. The method is made freely available to the community as a Python code. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: 13 pages, 8 figures, accepted for publication in A&A

Journal ref: A&A 669, A147 (2023)

Modeling globular clusters in the TNG50 simulation: predictions from dwarfs to giants

Authors: Jessica E. Doppel, Laura V. Sales, Dylan Nelson, Annalisa Pillepich, Mario G. Abadi, Eric W. Peng, Federico Marinacci, Jill Naiman, Paul Torrey, Mark Vogelsberger, Rainer Weinberger, Lars Hernquist

Abstract: We present a post-processing catalog of globular clusters (GCs) for the $39$ most massive groups and clusters in the TNG50 simulation of the IlllustrisTNG project (virial masses $M_{200} =[5\times 10^{12} \rm - 2 \times 10^{14}$] M$_{\odot}$). We tag GC particles to all galaxies with stellar mass $M_* \geq 5\times10^6$ M$_{\odot}$, and we calibrate their masses to reproduce the observed power-law… ▽ More We present a post-processing catalog of globular clusters (GCs) for the $39$ most massive groups and clusters in the TNG50 simulation of the IlllustrisTNG project (virial masses $M_{200} =[5\times 10^{12} \rm - 2 \times 10^{14}$] M$_{\odot}$). We tag GC particles to all galaxies with stellar mass $M_* \geq 5\times10^6$ M$_{\odot}$, and we calibrate their masses to reproduce the observed power-law relation between GC mass and halo mass for galaxies with $M_{200} \geq 10^{11}$ M$_{\odot}$ (corresponding to $M_* \sim 10^9$ $M_{\odot}$). Here we explore whether an extrapolation of this $M_{\rm GC}$-$M_{200}$ relation to lower-mass dwarfs is consistent with current observations. We find a good agreement between our predicted number and specific frequency of GCs in dwarfs with $\rm M_*=[5 \times 10^6 \rm - 10^9]$ M$_{\odot}$ and observations. Moreover, we predict a steep decline in the GC occupation fraction for dwarfs with $M_*<10^9$ M$_{\odot}$ which agrees well with current observational constraints. This declining occupation fraction is due to a combination of tidal stripping in all dwarfs plus a stochastic sampling of the GC mass function for dwarfs with $M_* < 10^{7.5}$ M$_{\odot}$. Our simulations also reproduce available constraints on the abundance of intra-cluster GCs in Virgo and Centaurus A. These successes provide support to the hypothesis that the $M_{\rm GC}$-$M_{200}$ relation holds, albeit with more scatter, all the way down to the regime of classical dwarf spheroidals in these environments. Our GC catalogs are publicly available as part of the IllustrisTNG data release. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 16 pages, 8 figures (main body of paper). Resubmitted to MNRAS after positive referee report. The globular cluster catalogs for the 39 most massive groups and clusters in TNG50 are available for download from the IllustrisTNG public data release at www.tng-project.org/doppel22

Origin and evolution of ultra-diffuse galaxies in different environments

Authors: Jose A. Benavides, Laura V. Sales, Mario. G. Abadi, Federico Marinacci, Mark Vogelsberger, Lars Hernquist

Abstract: We study the formation of ultra-diffuse galaxies (UDGs) using the cosmological hydrodynamical simulation TNG50 of the Illustris-TNG suite. We define UDGs as dwarf galaxies in the stellar mass range $\rm{7.5 \leq log (M_{\star} / M_{\odot}) \leq 9 }$ that are in the $5\%$ most extended tail of the simulated mass-size relation. This results in a sample of UDGs with half-mass radii… ▽ More We study the formation of ultra-diffuse galaxies (UDGs) using the cosmological hydrodynamical simulation TNG50 of the Illustris-TNG suite. We define UDGs as dwarf galaxies in the stellar mass range $\rm{7.5 \leq log (M_{\star} / M_{\odot}) \leq 9 }$ that are in the $5\%$ most extended tail of the simulated mass-size relation. This results in a sample of UDGs with half-mass radii $\rm{r_{h \star } \gtrsim 2 \ kpc}$ and surface brightness between $\rm{24.5}$ and $\rm{28 \ mag \ arcsec^{-2}}$, similar to definitions of UDGs in observations. The large cosmological volume in TNG50 allows for a comparison of UDGs properties in different environments, from the field to galaxy clusters with virial mass $\rm{M_{200} \sim 2 \times 10^{14} ~ M_{\odot}}$. All UDGs in our sample have dwarf-mass haloes ($\rm{M_{200}\sim 10^{11} ~ M_{\odot} }$) and show the same environmental trends as normal dwarfs: field UDGs are star-forming and blue while satellite UDGs are typically quiescent and red. The TNG50 simulation predicts UDGs that populate preferentially higher spin haloes and more massive haloes at fixed $\rm{M_{\star}}$ compared to non-UDG dwarfs. This applies also to most satellite UDGs, which are actually ``born" UDGs in the field and infall into groups and clusters without significant changes to their size. We find, however, a small subset of satellite UDGs ($\lesssim 10 \%$) with present-day stellar size a factor $\geq 1.5$ larger than at infall, confirming that tidal effects, particularly in the lower mass dwarfs, are also a viable formation mechanism for some of these dwarfs, although subdominant in this simulation. △ Less

Submitted 12 April, 2023; v1 submitted 15 September, 2022; originally announced September 2022.

Comments: 17 pages, 14 figures, accepted to MNRAS

The Next Generation Virgo Cluster Survey. XXXIII. Stellar Population Gradients in the Virgo Cluster Core Globular Cluster System

Authors: Youkyung Ko, Eric W. Peng, Patrick Côté, Laura Ferrarese, Chengze Liu, Alessia Longobardi, Ariane Lançon, Roberto P. Muñoz, Thomas H. Puzia, Karla A. Alamo-Martínez, Laura V. Sales, Felipe Ramos-Almendares, Mario G. Abadi, Myung Gyoon Lee, Ho Seong Hwang, Nelson Caldwell, John P. Blakeslee, Alessandro Boselli, Jean-Charles Cuillandre, Pierre-Alain Duc, Susana Eyheramendy, Puragra Guhathakurta, Stephen Gwyn, Andrés Jordán, Sungsoon Lim , et al. (2 additional authors not shown)

Abstract: We present a study of the stellar populations of globular clusters (GCs) in the Virgo Cluster core with a homogeneous spectroscopic catalog of 692 GCs within a major axis distance $R_{\rm maj} = $ 840 kpc from M87. We investigate radial and azimuthal variations in the mean age, total metallicity, [Fe/H], and $α$-element abundance, of blue (metal-poor) and red (metal-rich) GCs using their co-added… ▽ More We present a study of the stellar populations of globular clusters (GCs) in the Virgo Cluster core with a homogeneous spectroscopic catalog of 692 GCs within a major axis distance $R_{\rm maj} = $ 840 kpc from M87. We investigate radial and azimuthal variations in the mean age, total metallicity, [Fe/H], and $α$-element abundance, of blue (metal-poor) and red (metal-rich) GCs using their co-added spectra. We find that the blue GCs have a steep radial gradient in [Z/H] within $R_{\rm maj} =$ 165 kpc, with roughly equal contributions from [Fe/H] and [$α$/Fe], and flat gradients beyond. By contrast, the red GCs show a much shallower gradient in [Z/H], which is entirely driven by [Fe/H]. We use GC-tagged Illustris simulations to demonstrate an accretion scenario where more massive satellites (with more metal- and $α$-rich GCs) sink further into the central galaxy than less massive ones, and where the gradient flattening occurs because of the low GC occupation fraction of low-mass dwarfs disrupted at larger distances. The dense environment around M87 may also cause the steep [$α$/Fe] gradient of the blue GCs, mirroring what is seen in the dwarf galaxy population. The progenitors of red GCs have a narrower mass range than those of blue GCs, which makes their gradients shallower. We also explore spatial inhomogeneity in GC abundances, finding that the red GCs to the northwest of M87 are slightly more metal-rich. Future observations of GC stellar population gradients will be useful diagnostics of halo merger histories. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 22 pages, 8 figures, 4 tables, accepted for publication in ApJ

Numerical simulations of bar formation in the Local Group

Authors: Ornela F. Marioni, Mario G. Abadi, Stefan Gottlöber, Gustavo Yepes

Abstract: More than 50 per cent of present-day massive disc galaxies show a rotating stellar bar. Their formation and dynamics have been widely studied both numerically and observationally. Although numerical simulations in the $Λ$CDM cosmological framework predict the formation of such stellar components, there seems to be a tension between theoretical and observational results. Simulated bars are typicall… ▽ More More than 50 per cent of present-day massive disc galaxies show a rotating stellar bar. Their formation and dynamics have been widely studied both numerically and observationally. Although numerical simulations in the $Λ$CDM cosmological framework predict the formation of such stellar components, there seems to be a tension between theoretical and observational results. Simulated bars are typically larger in size and have slower pattern speed than observed ones. We study the formation and evolution of barred galaxies, using two $Λ$CDM zoom-in hydrodynamical simulations of the CLUES project that follow the evolution of a cosmological Local Group-like volume. We found that our simulated bars, at $z = 0$, are both shorter and faster rotators than previous ones found in other studies on cosmological simulations alleviating the tension mentioned above. These bars match the short tail-end of the observed bar length distribution. In agreement with previous numerical works, we find that bars form in those systems where the disc self-gravity is dominant over the dark matter halo, making them unstable against bar formation. Our bars developed in the last 3-4 Gyr until they achieve their current length and strength; as bars grow, their lengths increase while their rotation speeds decrease. Despite this slowdown, at redshift $z = 0$ their rotation speeds and size match well the observational data. △ Less

Submitted 13 January, 2022; originally announced January 2022.

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

Quiescent Ultra-diffuse galaxies in the field originating from backsplash orbits

Authors: José A. Benavides, Laura V. Sales, Mario. G. Abadi, Annalisa Pillepich, Dylan Nelson, Federico Marinacci, Michael Cooper, Ruediger Pakmor, Paul Torrey, Mark Vogelsberger, Lars Hernquist

Abstract: Ultra-diffuse galaxies (UDGs) are the lowest-surface brightness galaxies known, with typical stellar masses of dwarf galaxies but sizes similar to larger galaxies like the Milky Way. The reason for their extended sizes is debated, with suggested internal processes like angular momentum, feedback or mergers versus external mechanisms or a combination of both. Observationally, we know that UDGs are… ▽ More Ultra-diffuse galaxies (UDGs) are the lowest-surface brightness galaxies known, with typical stellar masses of dwarf galaxies but sizes similar to larger galaxies like the Milky Way. The reason for their extended sizes is debated, with suggested internal processes like angular momentum, feedback or mergers versus external mechanisms or a combination of both. Observationally, we know that UDGs are red and quiescent in groups and clusters while their counterparts in the field are blue and star-forming. This dichotomy suggests environmental effects as main culprit. However, this scenario is challenged by recent observations of isolated quiescent UDGs in the field. Here we use $Λ$CDM cosmological hydrodynamical simulation to show that isolated quenched UDGs are formed as backsplash galaxies that were once satellites of another galactic, group or cluster halo but are today a few Mpc away from them. These interactions, albeit brief, remove the gas and tidally strip the outskirts of the dark matter haloes of the now quenched seemingly-isolated UDGs, which are born as star-forming field UDGs occupying dwarf-mass dark matter haloes. Quiescent UDGs may therefore be found in non-negligible numbers in filaments and voids, bearing the mark of past interactions as stripped outer haloes devoid of dark matter and gas compared to dwarfs with similar stellar content. △ Less

Submitted 3 September, 2021; originally announced September 2021.

Comments: Author's version for the main article and supplementary information (13 pages, 10 figures). Accepted for publication in Nature Astronomy

Journal ref: Nature Astronomy 2021

Globular clusters as tracers of the dark matter content of dwarfs in galaxy clusters

Authors: Jessica E. Doppel, Laura V. Sales, Julio F. Navarro, Mario G. Abadi, Eric W. Peng, Elisa Toloba

Abstract: Globular clusters (GCs) are often used to estimate the dark matter content of galaxies, especially dwarfs, where other kinematic tracers are lacking. These estimates typically assume spherical symmetry and dynamical equilibrium, assumptions that may not hold for the sparse GC population of dwarfs in galaxy clusters. We use a catalog of GCs tagged onto the Illustris simulation to study the accuracy… ▽ More Globular clusters (GCs) are often used to estimate the dark matter content of galaxies, especially dwarfs, where other kinematic tracers are lacking. These estimates typically assume spherical symmetry and dynamical equilibrium, assumptions that may not hold for the sparse GC population of dwarfs in galaxy clusters. We use a catalog of GCs tagged onto the Illustris simulation to study the accuracy of GC-based mass estimates. We focus on galaxies in the stellar mass range 10$^{8} - 10^{11.8}$ M$_{\odot}$ identified in $9$ simulated Virgo-like clusters. Our results indicate that mass estimates are, on average, quite accurate in systems with GC numbers $N_{\rm GC} \geq 10$ and where the uncertainty of individual GC line-of-sight velocities is smaller than the inferred velocity dispersion, $σ_{\rm GC}$. In cases where $N_{\rm GC} \leq 10$, however, biases may result depending on how $σ_{\rm GC}$ is computed. We provide calibrations that may help alleviate these biases in methods widely used in the literature. As an application, we find a number of dwarfs with $M_{*} \sim 10^{8.5}\, M_{\odot}$ (comparable to the ultradiffuse galaxy DF2, notable for the low $σ_{GC}$ of its $10$ GCs) with $σ_{\rm GC} \sim 7$ - $15\; \rm km \rm s^{-1}$. These DF2 analogs correspond to relatively massive systems at their infall time ($M_{200} \sim 1$ - $3 \times 10^{11}$ $M_{\odot}$) which have retained only $3$-$17$ GCs and have been stripped of more than 95$\%$ of their dark matter. Our results suggest that extreme tidal mass loss in otherwise normal dwarf galaxies may be a possible formation channel for ultradiffuse objects like DF2. △ Less

Submitted 10 February, 2021; v1 submitted 13 October, 2020; originally announced October 2020.

Comments: 19 pages, 15 figures. Accepted to MNRAS Dec. 11 2020

Associations of dwarf galaxies in a $Λ$CDM Universe

Authors: C. Y. Yaryura, M. G. Abadi, S. Gottlober, N. I. Libeskind, S. A. Cora, A. N. Ruiz, C. A. Vega-Martínez, Gustavo Yepes, Peter Behroozi

Abstract: Associations of dwarf galaxies are loose systems composed exclusively of dwarf galaxies. These systems were identified in the Local Volume for the first time more than thirty years ago. We study these systems in the cosmological framework of the $Λ$ Cold Dark Matter ($Λ$CDM) model. We consider the Small MultiDark Planck simulation and populate its dark matter haloes by applying the semi-analytic m… ▽ More Associations of dwarf galaxies are loose systems composed exclusively of dwarf galaxies. These systems were identified in the Local Volume for the first time more than thirty years ago. We study these systems in the cosmological framework of the $Λ$ Cold Dark Matter ($Λ$CDM) model. We consider the Small MultiDark Planck simulation and populate its dark matter haloes by applying the semi-analytic model of galaxy formation SAG. We identify galaxy systems using a friends of friends algorithm with a linking length equal to $b=0.4 \,{\rm Mpc}\,h^{-1}$, to reproduce the size of dwarf galaxy associations detected in the Local Volume. Our samples of dwarf systems are built up removing those systems that have one (or more) galaxies with stellar mass larger than a maximum threshold $M_{\rm max}$. We analyse three different samples defined by ${\rm log}_{10}(M_{\rm max}[{\rm M}_{\odot}\,h^{-1}]) = 8.5, 9.0$ and $9.5$. On average, our systems have typical sizes of $\sim 0.2\,{\rm Mpc}\,h^{-1}$, velocity dispersion of $\sim 30 {\rm km\,s^{-1}} $ and estimated total mass of $\sim 10^{11} {\rm M}_{\odot}\,h^{-1}$. Such large typical sizes suggest that individual members of a given dwarf association reside in different dark matter haloes and are generally not substructures of any other halo. Indeed, in more than 90 per cent of our dwarf systems their individual members inhabit different dark matter haloes, while only in the remaining 10 per cent members do reside in the same halo. Our results indicate that the $Λ$CDM model can naturally reproduce the existence and properties of dwarf galaxies associations without much difficulty. △ Less

Submitted 12 October, 2020; originally announced October 2020.

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

A unified scenario for the origin of spiral and elliptical galaxy structural scaling laws

Authors: I. Ferrero, J. F. Navarro, M. G. Abadi, J. A. Benavides, D. Mast

Abstract: Elliptical (E) and spiral (S) galaxies follow tight, but different, scaling laws that link their stellar masses, radii, and characteristic velocities. Mass and velocity, for example, scale tightly in spirals with little dependence on galaxy radius (the "Tully-Fisher relation"; TFR). On the other hand, ellipticals appear to trace a 2D surface in size-mass-velocity space (the "Fundamental Plane"; FP… ▽ More Elliptical (E) and spiral (S) galaxies follow tight, but different, scaling laws that link their stellar masses, radii, and characteristic velocities. Mass and velocity, for example, scale tightly in spirals with little dependence on galaxy radius (the "Tully-Fisher relation"; TFR). On the other hand, ellipticals appear to trace a 2D surface in size-mass-velocity space (the "Fundamental Plane"; FP). Over the years, a number of studies have attempted to understand these empirical relations, usually in terms of variations of the virial theorem for E galaxies and in terms of the scaling relations of dark matter halos for spirals. We use Lambda cold dark matter (LCDM) cosmological hydrodynamical simulations to show that the scaling relations of both ellipticals and spirals arise as the result of (i) a tight galaxy mass-dark halo mass relation and (ii) the self-similar mass profile of cold dark matter halos. In this interpretation, E and S galaxies of a given stellar mass inhabit halos of similar masses, and their different scaling laws result from the varying amounts of dark matter enclosed within their luminous radii. This scenario suggests a new galaxy distance indicator applicable to galaxies of all morphologies and provides simple and intuitive explanations for long-standing puzzles, such as why the TFR is independent of surface brightness, or what causes the "tilt" in the FP. Our results provide strong support for the predictions of LCDM in the strongly non-linear regime, as well as guidance for further improvements to cosmological simulations of galaxy formation. △ Less

Submitted 3 March, 2021; v1 submitted 8 September, 2020; originally announced September 2020.

Comments: 17 pages, 13 figures. Accepted for publication in A&A

Journal ref: A&A 648, A124 (2021)

Accretion of Galaxy Groups into Galaxy Clusters

Authors: José. A. Benavides, Laura. V. Sales, Mario. G. Abadi

Abstract: We study the role of group infall in the assembly and dynamics of galaxy clusters in $Λ$CDM. We select $10$ clusters with virial mass $M_{\rm 200} \sim 10^{14} \, M_{\odot}$ from the cosmological hydrodynamical simulation Illustris and follow their galaxies with stellar mass $M_{\star} \geq 1.5 \times 10^8 \, M_{\odot}$. A median of $\sim 38\%$ of surviving galaxies at $z=0$ are accreted as part o… ▽ More We study the role of group infall in the assembly and dynamics of galaxy clusters in $Λ$CDM. We select $10$ clusters with virial mass $M_{\rm 200} \sim 10^{14} \, M_{\odot}$ from the cosmological hydrodynamical simulation Illustris and follow their galaxies with stellar mass $M_{\star} \geq 1.5 \times 10^8 \, M_{\odot}$. A median of $\sim 38\%$ of surviving galaxies at $z=0$ are accreted as part of groups and did not infall directly from the field, albeit with significant cluster-to-cluster scatter. The evolution of these galaxy associations is quick, with observational signatures of their common origin eroding rapidly in $1$-$3$ Gyr after infall. Substructure plays a dominant role in fostering the conditions for galaxy mergers to happen, even within the cluster environment. Integrated over time, we identify (per cluster) an average of $17 \pm 6$ mergers that occur in infalling galaxy associations, of which $7 \pm 3$ occur well within the virial radius of their cluster hosts. The number of mergers show large dispersion from cluster to cluster, with our most massive system having $42$ mergers above our mass cut-off. These mergers, which are typically gas rich for dwarfs and a combination of gas rich and gas poor for $M_{\star} \sim 10^{11} \, M_{\odot}$, may contribute significantly within $Λ$CDM to the formation of specific morphologies, such as lenticulars (S0) and blue compact dwarfs in groups and clusters. △ Less

Submitted 2 September, 2020; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: 12 pages, 11 figures. Accepted for publication in MNRAS

Properties of subhalos in the interacting dark matter scenario

Authors: Ángeles Moliné, Jascha A. Schewtschenko, Miguel A. Sánchez-Conde, Alejandra Aguirre-Santaella, Sofía A. Cora, Mario G. Abadi

Abstract: One possible and natural derivation from the collisionless cold dark matter (CDM) standard cosmological framework is the assumption of the existence of interactions between dark matter (DM) and photons or neutrinos. Such possible interacting dark matter (IDM) model would imply a suppression of small-scale structures due to a large collisional damping effect, even though the weakly interacting mass… ▽ More One possible and natural derivation from the collisionless cold dark matter (CDM) standard cosmological framework is the assumption of the existence of interactions between dark matter (DM) and photons or neutrinos. Such possible interacting dark matter (IDM) model would imply a suppression of small-scale structures due to a large collisional damping effect, even though the weakly interacting massive particle (WIMP) can still be the DM candidate. Because of this, IDM models can help alleviate alleged tensions between standard CDM predictions and observations at small mass scales. In this work, we investigate the properties of DM halo substructure or subhalos formed in a high-resolution cosmological N-body simulation specifically run within these alternative models. We also run its CDM counterpart, which allowed us to compare subhalo properties in both cosmologies. We show that, in the lower mass range covered by our simulation runs, both subhalo concentrations and abundances are systematically lower in IDM compared to the CDM scenario. Yet, as in CDM, we find that median IDM subhalo concentration values increase towards the innermost regions of their hosts for same mass subhalos. Also similarly to CDM, we find IDM subhalos to be more concentrated than field halos of the same mass. Our work has a direct application on studies aimed at the indirect detection of DM where subhalos are expected to boost the DM signal of their host halos significantly. From our results, we conclude that the role of halo substructure in DM searches will be less important in interacting scenarios than in CDM, but is nevertheless far from being negligible. △ Less

Submitted 9 October, 2019; v1 submitted 29 July, 2019; originally announced July 2019.

Comments: 16 pages, 6 figures and 2 tables. Invited contribution submitted to Galaxies, Special Issue "The Role of Halo Substructure in Gamma-Ray Dark Matter Searches" v2: minor changes, matches the published version

Report number: HALO-SPECIAL_ISSUE/2020/09

Journal ref: Galaxies 2019, 7(4), 80

Simulating the Spatial Distribution and Kinematics of Globular Clusters within Galaxy Clusters in Illustris

Authors: Felipe Ramos-Almendares, Laura V. Sales, Mario G. Abadi, Jessica E. Doppel, Hernan Muriel, Eric W. Peng

Abstract: We study the assembly of globular clusters (GCs) in 9 galaxy clusters using the cosmological simulation Illustris. GCs are tagged to individual galaxies at their infall time. The tidal removal of GCs from their galaxies and the distribution of the GCs within the cluster is later followed self-consistently by the simulation. The method relies on the simple assumption of a single power-law relation… ▽ More We study the assembly of globular clusters (GCs) in 9 galaxy clusters using the cosmological simulation Illustris. GCs are tagged to individual galaxies at their infall time. The tidal removal of GCs from their galaxies and the distribution of the GCs within the cluster is later followed self-consistently by the simulation. The method relies on the simple assumption of a single power-law relation between halo mass (M_vir) and mass in GCs (M_GC) as found in observations. We find that the GCs specific frequency $S_N$ as a function of V-band magnitude naturally reproduces the observed "U"-shape, due to the combination of a power law M_GC-M_vir relation and the non-linear M_*-M_vir relation from the simulation. Additional scatter in the $S_N$ values are traced back to galaxies with early infall times due to the evolution in the M_*-M_vir relation with redshift. GCs that have been tidally removed from their galaxies form today the intra-cluster component from which about ~60% were brought in by galaxies that orbit today within the cluster potential. The remaining "orphan" GCs are contributed by satellite galaxies with a wide range of stellar masses that are fully tidally disrupted at z=0. This intra-cluster component is a good dynamical tracer of the dark matter potential. As a consequence of the accreted nature of most intra-cluster GCs, their orbits are fairly radial with a predicted orbital anisotropy β>= 0.5. However, local tangential motions may appear as a consequence of localized substructure, providing a possible interpretation to the β<0 values suggested in observations of M87. △ Less

Submitted 21 February, 2020; v1 submitted 27 June, 2019; originally announced June 2019.

Comments: 13 pages, 10 figures, accepted in MNRAS

Ring Galaxies in the EAGLE Hydrodynamical Simulations

Authors: Ahmed Elagali, Claudia D. P. Lagos, O. Ivy Wong, Lister Staveley-Smith, James W. Trayford, Matthieu Schaller, Tiantian Yuan, Mario G. Abadi

Abstract: We study the formation and evolution of ring galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. We use the largest reference model Ref-L100N1504, a cubic cosmological volume of 100 comoving megaparsecs on a side, to identify and characterise these systems through cosmic time. The number density of ring galaxies in EAGLE is in broad agreement with the obs… ▽ More We study the formation and evolution of ring galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. We use the largest reference model Ref-L100N1504, a cubic cosmological volume of 100 comoving megaparsecs on a side, to identify and characterise these systems through cosmic time. The number density of ring galaxies in EAGLE is in broad agreement with the observations. The vast majority of ring galaxies identified in EAGLE (83 per cent) have an interaction origin, i.e., form when one or more companion galaxies drop-through a disk galaxy. The remainder (17 per cent) have very long-lived ring morphologies (> 2 Gyr) and host strong bars. Ring galaxies are HI rich galaxies, yet display inefficient star formation activity and tend to reside in the green valley particularly at z > 0.5. This inefficiency is mainly due to the low pressure and metallicity of their interstellar medium (ISM) compared with the ISM of similar star-forming galaxies. We find that the interaction(s) is responsible for decreasing the ISM pressure by causing the ISM gas to flow from the inner regions to the outer disk, where the ring feature forms. At a fixed radius, the star formation efficiency of ring galaxies is indistinguishable from their star-forming counterparts, and thus the main reason for their integrated lower efficiency is the different gas surface density profiles. Since galaxy morphologies are not used to tune the parameters in hydrodynamical simulations, the experiment performed here demonstrates the success of the current numerical models in EAGLE. △ Less

Submitted 7 September, 2018; v1 submitted 22 July, 2018; originally announced July 2018.

Comments: Accepted for publication in MNRAS!

Intra-cluster Globular Clusters in a Simulated Galaxy Cluster

Authors: Felipe Ramos-Almendares, Mario G. Abadi, Hernán Muriel, Valeria Coenda

Abstract: Using a cosmological dark matter simulation of a galaxy-cluster halo, we follow the temporal evolution of its globular cluster population. To mimic the red and blue globular cluster populations, we select at high redshift $(z\sim 1)$ two sets of particles from individual galactic halos constrained by the fact that, at redshift $z=0$, they have density profiles similar to observed ones. At redshift… ▽ More Using a cosmological dark matter simulation of a galaxy-cluster halo, we follow the temporal evolution of its globular cluster population. To mimic the red and blue globular cluster populations, we select at high redshift $(z\sim 1)$ two sets of particles from individual galactic halos constrained by the fact that, at redshift $z=0$, they have density profiles similar to observed ones. At redshift $z=0$, approximately 60\% of our selected globular clusters were removed from their original halos building up the intra-cluster globular cluster population, while the remaining 40\% are still gravitationally bound to their original galactic halos. Since the blue population is more extended than the red one, the intra-cluster globular cluster population is dominated by blue globular clusters, with a relative fraction that grows from 60\% at redshift $z=0$ up to 83\% for redshift $z\sim 2$. In agreement with observational results for the Virgo galaxy cluster, the blue intra-cluster globular cluster population is more spatially extended than the red one, pointing to a tidally disrupted origin. △ Less

Submitted 14 December, 2017; originally announced December 2017.

Comments: accepted for publication in ApJ

Barred galaxies in the EAGLE cosmological hydrodynamical simulation

Authors: David G. Algorry, Julio F. Navarro, Mario G. Abadi, Laura V. Sales, Richard G. Bower, Robert A. Crain, Claudio Dalla Vecchia, Carlos S. Frenk, Matthieu Schaller, Joop Schaye, Tom Theuns

Abstract: We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift… ▽ More We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift = 1.3, on timescales that depend sen- sitively on the strength of the pattern. Strong bars develop relatively quickly (in a few Gyr, = 10 disc rotation periods) in systems that are disc dominated, gas poor, and have declining rotation curves. Weak bars develop more slowly in systems where the disc is less gravitation- ally important, and are still growing at z = 0. Unbarred galaxies are comparatively gas-rich discs whose rotation speeds do not exceed the maximum circular velocity of the halos they inhabit. Bar lengths compare favourably with observations, ranging from 0.2 to 0.8 times the radius containing 90% of the stars. Bars slow down remarkably quickly as they grow, causing the inner regions of the surrounding dark halo to expand. At z = 0 strong bars have corota- tion radii roughly ten times the bar length. Such slow bars are inconsistent with the few cases where pattern speeds have been measured or inferred observationally, a discrepancy that, if confirmed, might prove a challenge for disc galaxy formation in LCDM. △ Less

Submitted 19 September, 2016; originally announced September 2016.

Comments: 11 pages, 10 figures

Size matters: abundance matching, galaxy sizes, and the Tully-Fisher relation in EAGLE

Authors: Ismael Ferrero, Julio F. Navarro, Mario G. Abadi, Laura V. Sales, Richard G. Bower, Robert A. Crain, Carlos S. Frenk, Matthieu Schaller, Joop Schaye, Tom Theuns

Abstract: The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM)… ▽ More The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM) models, which predict the galaxy-halo mass relation to be non-monotonic and rapidy evolving. We study the TFR of luminous spirals and its relation to AM using the EAGLE set of $Λ$CDM cosmological simulations. Matching both relations requires disk sizes to satisfy constraints given by the concentration of halos and their response to galaxy assembly. EAGLE galaxies approximately match these constraints and show a tight mass-velocity scaling that compares favourably with the observed TFR. The TFR is degenerate to changes in galaxy formation efficiency and the mass-size relation; simulations that fail to match the galaxy stellar mass function may fit the observed TFR if galaxies follow a different mass-size relation. The small scatter in the simulated TFR results because, at fixed halo mass, galaxy mass and rotation speed correlate strongly, scattering galaxies along the main relation. EAGLE galaxies evolve with lookback time following approximately the prescriptions of AM models and the observed mass-size relation of bright spirals, leading to a weak TFR evolution consistent with observation out to $z=1$. $Λ$CDM models that match both the abundance and size of galaxies as a function of stellar mass have no difficulty reproducing the observed TFR and its evolution. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Comments: 12 pages, 7 figures, 1 table, submitted to MNRAS

The low-mass end of the baryonic Tully-Fisher relation

Authors: Laura V. Sales, Julio F. Navarro, Kyle Oman, Azadeh Fattahi, Ismael Ferrero, Mario G. Abadi, Richard Bower, Robert A. Crain, Carlos S. Frenk, Till Sawala, Matthieu Schaller, Joop Schaye, Tom Theuns, Simon D. M. White

Abstract: The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be exp… ▽ More The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation. △ Less

Submitted 5 February, 2016; originally announced February 2016.

Comments: 10 pages, 7 figures, submitted to MNRAS

The low mass end of the neutral gas mass and velocity width functions of galaxies in $Λ$CDM

Authors: C. Y. Yaryura, A. Helmi, M. G. Abadi, E. Starkenburg

Abstract: We use the high-resolution Aquarius cosmological dark matter simulations coupled to the semi-analytic model by Starkenburg et al. (2013) to study the HI content and velocity width properties of field galaxies at the low mass end in the context of $Λ$CDM. We compare our predictions to the observed ALFALFA survey HI mass and velocity width functions, and find very good agreement without fine-tuning,… ▽ More We use the high-resolution Aquarius cosmological dark matter simulations coupled to the semi-analytic model by Starkenburg et al. (2013) to study the HI content and velocity width properties of field galaxies at the low mass end in the context of $Λ$CDM. We compare our predictions to the observed ALFALFA survey HI mass and velocity width functions, and find very good agreement without fine-tuning, when considering central galaxies. Furthermore, the properties of the dark matter halos hosting galaxies, characterised by their peak velocity and circular velocity at 2 radial disk scalelengths overlap perfectly with the inferred values from observations. This suggests that our galaxies are placed in the right dark matter halos, and consequently at face value, we do not find any discrepancy with the predictions from the $Λ$CDM model. Our analysis indicates that previous tensions, apparent when using abundance matching models, arise because this technique cannot be straightforwardly applied for objects with masses $M_{vir} < 10^{10} M_{\odot}$. △ Less

Submitted 27 January, 2016; originally announced January 2016.

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

Mergers and the outside-in formation of dwarf spheroidals

Authors: Alejandro Benítez-Llambay, Julio F. Navarro, Mario G. Abadi, Stefan Gottloeber, Gustavo Yepes, Yehuda Hoffman, Matthias Steinmetz

Abstract: We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few example… ▽ More We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few examples where stars form in two populations distinct in age in order to elucidate the origin of these gradients. The spatial distributions of the two components reflect their diverse origin; the old stellar component is assembled through mergers, but the young population forms largely in situ. The older component results from a first episode of star formation that begins early but is quickly shut off by the combined effects of stellar feedback and reionization. The younger component forms when a late accretion event adds gas and reignites star formation. The effect of mergers is to disperse the old stellar population, increasing their radius and decreasing their central density relative to the young population. We argue that dwarf-dwarf mergers offer a plausible scenario for the formation of systems with multiple distinct populations and, more generally, for the origin of age and metallicity gradients in dwarf spheroidals. △ Less

Submitted 19 November, 2015; originally announced November 2015.

Comments: 10 pages, 8 figures, Accepted for publication in MNRAS

The Imprint of Reionization on the Star Formation Histories of Dwarf Galaxies

Authors: Alejandro Benítez-Llambay, Julio F. Navarro, Mario G. Abadi, Stefan Gottloeber, Gustavo Yepes, Yehuda Hoffman, Matthias Steinmetz

Abstract: We explore the impact of cosmic reionization on nearby isolated dwarf galaxies using a compilation of SFHs estimated from deep HST data and a cosmological hydrodynamical simulation of the Local Group. The nearby dwarfs show a wide diversity of star formation histories; from ancient systems that have largely completed their star formation $\sim 10$ Gyr ago to young dwarfs that have formed the major… ▽ More We explore the impact of cosmic reionization on nearby isolated dwarf galaxies using a compilation of SFHs estimated from deep HST data and a cosmological hydrodynamical simulation of the Local Group. The nearby dwarfs show a wide diversity of star formation histories; from ancient systems that have largely completed their star formation $\sim 10$ Gyr ago to young dwarfs that have formed the majority of their stars in the past $\sim 5$ Gyr to two-component systems characterized by the overlap of comparable numbers of old and young stars. Taken as an ensemble, star formation in nearby dwarfs dips to lower-than-average rates at intermediate times ($4<t$/Gyr $<8$), a feature that we trace in the simulation to the effects of cosmic reionization. Reionization heats the gas and drives it out of the shallow potential wells of low mass halos, affecting especially those below a sharp mass threshold that corresponds to a virial temperature of $\sim 2 \times 10^4 $ $\mathrm{K}$ at $z_{\rm reion}$. The loss of baryons leads to a sharp decline in the star forming activity of early-collapsing systems, which, compounded by feedback from early star formation, empties halos of gas and leaves behind systems where a single old stellar component prevails. In halos below the threshold at $z_{\rm reion}$, reionization heating leads to a delay in the onset of star formation that lasts until the halo grows massive enough to allow some of the remaining gas to cool and form stars. Young stellar components therefore dominate in dwarfs whose halos assemble late and thus form few stars before reionization. Two-component systems may be traced to late mergers of individual examples of the two aforementioned cases. The relative dearth of intermediate-age stars in nearby dwarfs might thus be the clearest signature yet identified of the imprint of cosmic reionization on the star formation history of dwarf galaxies. △ Less

Submitted 28 May, 2014; v1 submitted 21 May, 2014; originally announced May 2014.

Comments: 13 pages, 1 table, 10 figures. Submitted to MNRAS. Find supplementary movies at http://www.iate.oac.uncor.edu/~alejandro/sfrreion-web/ or at http://www.astro.uvic.ca/~abll/sfrreion-web/ ; Dwarf galaxy sample updated to remove possible satellites of M81, M82, and NGC 2403. Main results and conclusions unchanged

Counterrotating Stars in Simulated Galaxy Disks

Authors: David G. Algorry, Julio F. Navarro, Mario G. Abadi, Laura V. Sales, Matthias Steinmetz, Franziska Piontek

Abstract: Counterrotating stars in disk galaxies are a puzzling dynamical feature whose origin has been ascribed to either satellite accretion events or to disk instabilities triggered by deviations from axisymmetry. We use a cosmological simulation of the formation of a disk galaxy to show that counterrotating stellar disk components may arise naturally in hierarchically-clustering scenarios even in the ab… ▽ More Counterrotating stars in disk galaxies are a puzzling dynamical feature whose origin has been ascribed to either satellite accretion events or to disk instabilities triggered by deviations from axisymmetry. We use a cosmological simulation of the formation of a disk galaxy to show that counterrotating stellar disk components may arise naturally in hierarchically-clustering scenarios even in the absence of merging. The simulated disk galaxy consists of two coplanar, overlapping stellar components with opposite spins: an inner counterrotating bar-like structure made up mostly of old stars surrounded by an extended, rotationally-supported disk of younger stars. The opposite-spin components originate from material accreted from two distinct filamentary structures which at turn around, when their net spin is acquired, intersect delineating a "V"-like structure. Each filament torques the other in opposite directions; the filament that first drains into the galaxy forms the inner counterrotating bar, while material accreted from the other filament forms the outer disk. Mergers do not play a substantial role and most stars in the galaxy are formed in situ; only 9% of all stars are contributed by accretion events. The formation scenario we describe here implies a significant age difference between the co- and counterrotating components, which may be used to discriminate between competing scenarios for the origin of counterrotating stars in disk galaxies. △ Less

Submitted 5 November, 2013; originally announced November 2013.

Comments: 7 pages, 7 figures. Accepted for publication in MNRAS

Dwarf Galaxies and the Cosmic Web

Authors: Alejandro Benitez-Llambay, Julio F. Navarro, Mario G. Abadi, Stefan Gottloeber, Gustavo Yepes, Yehuda Hoffman, Matthias Steinmetz

Abstract: We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the… ▽ More We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the process. The dark matter content of these halos is unaffected but their gas can be efficiently removed by ram-pressure. The loss of gas is especially pronounced in low-mass halos due to their lower binding energy and has a dramatic effect on the star formation history of affected systems. This "cosmic web stripping" may help to explain the scarcity of dwarf galaxies compared with the numerous low-mass halos expected in ΛCDM and the large diversity of star formation histories and morphologies characteristic of faint galaxies. Although our results are based on a single high-resolution simulation, it is likely that the hydrodynamical interaction of dwarf galaxies with the cosmic web is a crucial ingredient so far missing from galaxy formation models. △ Less

Submitted 2 November, 2012; originally announced November 2012.

Comments: Submitted to ApJL. 6 pages, 4 figures. A set of movies showing the interaction between dwarf galaxies and the Cosmic Web can be found at mirror 1 http://www.astro.uvic.ca/~mario/dwarf-web/ or at mirror 2 http://www.iate.oac.uncor.edu/~alejandro/dwarf-web/ . Comments are welcomed

Effects of Baryon Mass Loss on Profiles of Large Galactic Dark Matter Haloes

Authors: Cinthia Ragone-Figueroa, Gian Luigi Granato, Mario G. Abadi

Abstract: We perform controlled numerical experiments to asses the effect of baryon mass loss on the inner structure of large galactic dark matter haloes. This mass expulsion is intended to mimic both the supernovae and AGN feedbacks, as well as the evolution of stellar populations. This study is meant in particular for precursors of massive Early Type Galaxies, wherein strong AGN feedback (often dubbed "QS… ▽ More We perform controlled numerical experiments to asses the effect of baryon mass loss on the inner structure of large galactic dark matter haloes. This mass expulsion is intended to mimic both the supernovae and AGN feedbacks, as well as the evolution of stellar populations. This study is meant in particular for precursors of massive Early Type Galaxies, wherein strong AGN feedback (often dubbed "QSO mode" in galaxy formation models) has been proposed to remove on a short timescale, of the order of a few dynamical times, a substantial fraction of their baryons. In a previous paper we evaluated the observational consequences (size increase) of this process on the galactic structure (Ragone-Figueroa & Granato 2011). Here we focus on the distribution of dark matter in the galactic region. It is shown that the inner region of the DM halo expands and its density profile flattens by a sizeable amount, with little dependence on the expulsion timescale. We also evaluate the effect of the commonly made approximation of treating the baryonic component as a potential that changes in intensity without any variation in shape. This approximation leads to some underestimates of the halo expansion and its slope flattening. We conclude that cuspy density profiles in ETGs could be difficult to reconcile with an effective AGN (or stellar) feedback during the evolution of these systems. △ Less

Submitted 17 April, 2012; v1 submitted 7 February, 2012; originally announced February 2012.

Comments: 9 Pages, 5 figures, MNRAS accepted. Final version with minor revision

The dark matter halos of dwarf galaxies: a challenge for the LCDM paradigm?

Authors: Ismael Ferrero, Mario G. Abadi, Julio F. Navarro, Laura V. Sales, Sebastian Gurovich

Abstract: The cold dark matter halo mass function is much steeper than the galaxy stellar mass function on galactic and subgalactic scales. This difference is usually reconciled by assuming that the galaxy formation efficiency drops sharply with decreasing halo mass, so that virtually no dwarf galaxies form in halos less massive than ~ 10^10 M_sun. In turn, this implies that, at any given radius, the dark m… ▽ More The cold dark matter halo mass function is much steeper than the galaxy stellar mass function on galactic and subgalactic scales. This difference is usually reconciled by assuming that the galaxy formation efficiency drops sharply with decreasing halo mass, so that virtually no dwarf galaxies form in halos less massive than ~ 10^10 M_sun. In turn, this implies that, at any given radius, the dark mass enclosed by a galaxy must exceed a certain minimum. We use rotation curves of dwarf galaxies compiled from the literature to explore whether their enclosed mass is consistent with these constraints. We find that almost one half of the dwarfs in our sample with stellar mass between 10^6-10^7 M_sun are at odds with this restriction: either they live in halos with masses substantially below 10^10 M_sun or there is a mechanism capable of reducing the dark mass enclosed by some of the faintest dwarfs. Neither possibility is easily accommodated within the standard LCDM scenario. Extending galaxy formation to halos well below 10^10 M_sun would lead to large numbers of dwarf galaxies in excess of current estimates; at the same time, the extremely low stellar mass of the systems involved makes it unlikely that baryonic effects can reduce their dark matter content. Resolving this challenge seems to require new insights into dwarf galaxy formation, or perhaps a radical revision of the prevailing paradigm. △ Less

Submitted 11 September, 2012; v1 submitted 28 November, 2011; originally announced November 2011.

Comments: 8 pages, 4 Figures, MNRAS accepted

Journal ref: Monthly Notices of the Royal Astronomical Society, 425: 2817-2823 (2012)

A Sagittarius-Induced Origin for the Monoceros Ring

Authors: Leo Michel-Dansac, Mario G. Abadi, Julio F. Navarro, Matthias Steinmetz

Abstract: The Monoceros ring is a collection of stars in nearly-circular orbits at roughly 18 kpc from the Galactic center. It may have originated (i) as the response of the disc to perturbations excited by satellite companions or (ii) from the tidal debris of a disrupted dwarf galaxy. The metallicity of Monoceros stars differs from that of disc stars at comparable Galactocentric distances, an observation t… ▽ More The Monoceros ring is a collection of stars in nearly-circular orbits at roughly 18 kpc from the Galactic center. It may have originated (i) as the response of the disc to perturbations excited by satellite companions or (ii) from the tidal debris of a disrupted dwarf galaxy. The metallicity of Monoceros stars differs from that of disc stars at comparable Galactocentric distances, an observation that disfavours the first scenario. On the other hand, circular orbits are difficult to accommodate in the tidal-disruption scenario, since it requires a satellite which at the time of disruption was itself in a nearly circular orbit. Such satellite could not have formed at the location of the ring and, given its low mass, dynamical friction is unlikely to have played a major role in its orbital evolution. We search cosmological simulations for low-mass satellites in nearly-circular orbits and find that they result, almost invariably, from orbital changes induced by collisions with more massive satellites: the radius of the circular orbit thus traces the galactocentric distance of the collision. Interestingly, the Sagittarius dwarf, one of the most luminous satellites of the Milky Way, is in a polar orbit that crosses the Galactic plane at roughly the same Galactocentric distance as Monoceros. We use idealized simulations to demonstrate that an encounter with Sagittarius might well have led to the circularization and subsequent tidal demise of the progenitor of the Monoceros ring. △ Less

Submitted 9 May, 2011; v1 submitted 22 February, 2011; originally announced February 2011.

Comments: 6 pages, 4 figures, to match version published in MNRAS Letters (http://onlinelibrary.wiley.com/doi/10.1111/j.1745-3933.2011.01035.x/abstract)

Through Thick and Thin: Kinematic and Chemical Components in the Solar Neighbourhood

Authors: Julio F. Navarro, Mario G. Abadi, Kim A. Venn, Kenneth C. Freeman, Borja Anguiano

Abstract: We search for the existence of chemically-distinct stellar components in the solar neighbourhood using published data. Extending earlier work, we show that when the abundances of Fe, alpha elements, and the r-process element Eu are considered, stars separate neatly into two groups that delineate the traditional thin and thick disk of the Milky Way. The group akin to the thin disk is traced by star… ▽ More We search for the existence of chemically-distinct stellar components in the solar neighbourhood using published data. Extending earlier work, we show that when the abundances of Fe, alpha elements, and the r-process element Eu are considered, stars separate neatly into two groups that delineate the traditional thin and thick disk of the Milky Way. The group akin to the thin disk is traced by stars with [Fe/H]>-0.7 and alpha/Fe<0.2. The thick disk-like group overlaps the thin disk in [Fe/H] but has higher abundances of alpha elements and Eu. Stars in the range -1.5<[Fe/H]<-0.7 with low [alpha/Fe] ratios, however, seem to belong to a separate, dynamically-cold, non-rotating component that we associate with tidal debris. The kinematically-hot stellar halo dominates the sample for [Fe/H]<-1.5. These results suggest that it may be possible to define the main dynamical components of the solar neighbourhood using only their chemistry, an approach with a number of interesting consequences. The kinematics of thin disk stars is then independent of metallicity: their average rotation speed remains roughly constant in the range -0.7<[Fe/H]<+0.4, a result that argues against radial migration having played a substantial role in the evolution of the thin disk. The velocity dispersion of stars assigned to the thin disk is also independent of [Fe/H], implying that the familiar increase in velocity dispersion with decreasing metallicity is the result of the increased prevalence of the thick disk at lower metallicities, rather than of the sustained operation of a dynamical heating mechanism. The substantial overlap in [Fe/H] and, probably, stellar age, of the various components might affect other reported trends in the properties of stars in the solar neighbourhood. △ Less

Submitted 29 October, 2010; v1 submitted 31 August, 2010; originally announced September 2010.

Comments: 8 pages, 5 figures, some clarifications after referee report. Conclusions unchanged

Ram pressure stripping in a galaxy formation model. I. A novel numerical approach

Authors: Tomás E. Tecce, Sofía A. Cora, Patricia B. Tissera, Mario G. Abadi, Claudia del P. Lagos

Abstract: We develop a new numerical approach to describe the action of ram pressure stripping (RPS) within a semi-analytic model of galaxy formation and evolution which works in combination with non-radiative hydrodynamical simulations of galaxy clusters. The new feature in our method is the use of the gas particles to obtain the kinematical and thermodynamical properties of the intragroup and intracluster… ▽ More We develop a new numerical approach to describe the action of ram pressure stripping (RPS) within a semi-analytic model of galaxy formation and evolution which works in combination with non-radiative hydrodynamical simulations of galaxy clusters. The new feature in our method is the use of the gas particles to obtain the kinematical and thermodynamical properties of the intragroup and intracluster medium (ICM). This allows a self-consistent estimation of the RPS experienced by satellite galaxies. We find that the ram pressure (RP) in the central regions of clusters increases approximately one order of magnitude between z = 1 and 0, consistent with the increase in the density of the ICM. The mean RP experienced by galaxies within the virial radius increases with decreasing redshift. In clusters with virial masses M_vir ~10^15 h^-1 M_Sun, over 50 per cent of satellite galaxies have experienced RP ~10^(-11) h^2 dyn cm^-2 or higher for z <= 0.5. In smaller clusters (M_vir ~10^14 h^-1 M_Sun) the mean RP are approximately one order of magnitude lower at all redshifts. RPS has a strong effect on the cold gas content of galaxies for all cluster masses. At z = 0, over 70 per cent of satellite galaxies within the virial radius are completely depleted of cold gas. For the more massive clusters the fraction of depleted galaxies is already established at z ~ 1, whereas for the smaller clusters this fraction increases appreciably between z = 1 and 0. This indicates that the rate at which the cold gas is stripped depends on the virial mass of the host cluster. Compared to our new approach, the use of an analytic profile to describe the ICM results in an overestimation of the RP larger than 50 per cent for z >= 0.5. △ Less

Submitted 13 October, 2010; v1 submitted 28 June, 2010; originally announced June 2010.

Comments: 15 pages, 8 figures. Accepted for publication in MNRAS. Minor changes to match published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 408, Issue 4, pp. 2008-2021

The Sphericalization of Dark Matter Halos by Galaxy Disks

Authors: Stelios Kazantzidis, Mario G. Abadi, Julio F. Navarro

Abstract: Cosmological simulations indicate that cold dark matter (CDM) halos should be triaxial. Verifying observationally this theoretical prediction is, however, less than straightforward because the assembly of galaxies is expected to modify the halo shapes and to render them more axisymmetric. We use a suite of N-body simulations to investigate quantitatively the effect of the growth of a central disk… ▽ More Cosmological simulations indicate that cold dark matter (CDM) halos should be triaxial. Verifying observationally this theoretical prediction is, however, less than straightforward because the assembly of galaxies is expected to modify the halo shapes and to render them more axisymmetric. We use a suite of N-body simulations to investigate quantitatively the effect of the growth of a central disk galaxy on the shape of triaxial dark matter halos. As expected, the halo responds to the presence of the disk by becoming more spherical. The net effect depends only weakly on the orientation of the disk relative to the halo principal axes or the timescale of disk assembly, but strongly on the overall gravitational importance of the disk. Our results show that exponential disks whose contribution peaks at less than ~50% of their circular velocity are unable to modify noticeably the shape of the gravitational potential of their surrounding halos. Many dwarf and low surface brightness galaxies are expected to be in this regime, and therefore their detailed kinematics could be used to probe halo triaxiality, one of the basic predictions of the CDM paradigm. We argue that the complex disk kinematics of the dwarf galaxy NGC 2976 might be the reflection of a triaxial halo. Such signatures of halo triaxiality should be common in galaxies where the luminous component is subdominant. △ Less

Submitted 17 August, 2010; v1 submitted 2 June, 2010; originally announced June 2010.

Comments: 7 pages, 4 figures, LaTeX (uses emulateapj.cls). Replaced to match the version accepted for publication in ApJL. This version includes an additional figure in response to referee. Conclusions remain unchanged

Journal ref: Astrophys.J. 720 (2010) L62-L66

Orbital Eccentricity as a probe of Thick Disk Formation Scenarios

Authors: Laura V. Sales, Amina Helmi, Mario G. Abadi, Chris B. Brook, Facundo A. Gomez, Rok Roskar, Victor P. Debattista, Elisa House, Matthias Steinmetz, Alvaro Villalobos

Abstract: We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected f… ▽ More We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected for the heating, migration and gas-rich merging scenarios, while the eccentricity distribution is broader and shifted towards higher values for the accretion model. These differences can be traced back to whether the bulk of the stars in each case is formed 'in-situ' or is 'accreted', and are robust to the peculiarities of each model. A simple test based on the eccentricity distribution of nearby thick disk stars may thus help elucidate the dominant formation mechanism of the Galactic thick disk. △ Less

Submitted 23 September, 2009; v1 submitted 21 September, 2009; originally announced September 2009.

Comments: 5 pages, 3 figures, accepted to MNRAS Letters. Typo corrected

Galaxy-Induced Transformation of Dark Matter Halos

Authors: Mario G. Abadi, Julio F. Navarro, Mark Fardal, Arif Babul, Matthias Steinmetz

Abstract: We use N-body/gasdynamical LambdaCDM cosmological simulations to examine the effect of the assembly of a central galaxy on the shape and mass profile of its dark halo. Two series of simulations are compared; one that follows only the evolution of the dark matter component and a second one where a baryonic component is added. These simulations include radiative cooling but neglect star formation… ▽ More We use N-body/gasdynamical LambdaCDM cosmological simulations to examine the effect of the assembly of a central galaxy on the shape and mass profile of its dark halo. Two series of simulations are compared; one that follows only the evolution of the dark matter component and a second one where a baryonic component is added. These simulations include radiative cooling but neglect star formation and feedback, leading most baryons to collect at the halo center in a disk which is too small and too massive when compared with typical spiral. This unrealistic model allows us, nevertheless, to gauge the maximum effect that galaxies may have in transforming their dark halos. We find that the shape of the halo becomes more axisymmetric: halos are transformed from triaxial into essentially oblate systems, with well-aligned isopotential contours of roughly constant flattening (c/a ~ 0.85). Halos always contract as a result of galaxy assembly, but the effect is substantially less pronounced than predicted by the "adiabatic contraction" hypothesis. The reduced contraction helps to reconcile LambdaCDM halos with constraints on the dark matter content inside the solar circle and should alleviate the long-standing difficulty of matching simultaneously the scaling properties of galaxy disks and the luminosity function. The halo contraction is also less pronounced than found in earlier simulations, a disagreement that suggests that halo contraction is not solely a function of the initial and final distribution of baryons. Not only how much baryonic mass has been deposited at the center of a halo matters, but also the mode of its deposition. It might prove impossible to predict the halo response without a detailed understanding of a galaxy's assembly history. (Abriged) △ Less

Submitted 14 February, 2009; originally announced February 2009.

Comments: 11 pages and 9 figures

An Alternative Origin for Hypervelocity Stars

Authors: Mario G. Abadi, Julio F. Navarro, Matthias Steinmetz

Abstract: Halo stars with unusually high radial velocity ("hypervelocity" stars, or HVS) are thought to be stars unbound to the Milky Way that originate from the gravitational interaction of stellar systems with the supermassive black hole at the Galactic center. We examine the latest HVS compilation and find peculiarities that are unexpected in this black hole-ejection scenario. For example, a large frac… ▽ More Halo stars with unusually high radial velocity ("hypervelocity" stars, or HVS) are thought to be stars unbound to the Milky Way that originate from the gravitational interaction of stellar systems with the supermassive black hole at the Galactic center. We examine the latest HVS compilation and find peculiarities that are unexpected in this black hole-ejection scenario. For example, a large fraction of HVS cluster around the constellation of Leo and share a common travel time of $\sim 100$-200 Myr. Furthermore, their velocities are not really extreme if, as suggested by recent galaxy formation models, the Milky Way is embedded within a $2.5\times 10^{12} h^{-1} M_{\odot}$ dark halo with virial velocity of $\sim 220$ km/s. In this case, the escape velocity at $\sim 50$ kpc would be $\sim 600$ km/s and very few HVS would be truly unbound. We use numerical simulations to show that disrupting dwarf galaxies may contribute halo stars with velocities up to and sometimes exceeding the nominal escape speed of the system. These stars are arranged in a thinly-collimated outgoing ``tidal tail'' stripped from the dwarf during its latest pericentric passage. We speculate that some HVS may therefore be tidal debris from a dwarf recently disrupted near the center of the Galaxy. In this interpretation, the angular clustering of HVS results because from our perspective the tail is seen nearly ``end on'', whereas the common travel time simply reflects the fact that these stars were stripped simultaneously from the dwarf during a single pericentric passage. This proposal is eminently falsifiable, since it makes a number of predictions that are distinct from the black-hole ejection mechanism and that should be testable with improved HVS datasets. △ Less

Submitted 28 November, 2008; v1 submitted 8 October, 2008; originally announced October 2008.

Comments: 4 pages, 4 figures. Replacement to match version accepted to ApJL

Cosmic Ménage à Trois: The Origin of Satellite Galaxies On Extreme Orbits

Authors: L. V. Sales, J. F. Navarro, M. G. Abadi, M. Steinmetz

Abstract: We examine the orbits of satellite galaxies identified in a suite of N-body/gasdynamical simulations of the formation of $L_*$ galaxies in a LCDM universe. Most satellites follow conventional orbits; after turning around, they accrete into their host halo and settle on orbits whose apocentric radii are steadily eroded by dynamical friction. However, a number of outliers are also present, we find… ▽ More We examine the orbits of satellite galaxies identified in a suite of N-body/gasdynamical simulations of the formation of $L_*$ galaxies in a LCDM universe. Most satellites follow conventional orbits; after turning around, they accrete into their host halo and settle on orbits whose apocentric radii are steadily eroded by dynamical friction. However, a number of outliers are also present, we find that ~1/3 of satellites identified at $z=0$ are on unorthodox orbits, with apocenters that exceed their turnaround radii. This population of satellites on extreme orbits consists typically of the faint member of a satellite pair that has been ejected onto a highly-energetic orbit during its first approach to the primary. Since the concurrent accretion of multiple satellite systems is a defining feature of hierarchical models of galaxy formation, we speculate that this three-body ejection mechanism may be the origin of (i) some of the newly discovered high-speed satellites around M31 (such as Andromeda XIV); (ii) some of the distant fast-receding Local Group members, such as Leo I; and (iii) the oddly isolated dwarf spheroidals Cetus and Tucana in the outskirts of the Local Group. Our results suggest that care must be exercised when using the orbits of the most weakly bound satellites to place constraints on the total mass of the Local Group. △ Less

Submitted 24 May, 2007; v1 submitted 13 April, 2007; originally announced April 2007.

Comments: 10 pages, 6 figures, MNRAS in press. Accepted version with minor changes. Version with high resolution figures available at: http://www.astro.uvic.ca/~lsales/SatPapers/SatPapers.html

Journal ref: Mon.Not.Roy.Astron.Soc.379:1475-1483,2007

Satellites of Simulated Galaxies: survival, merging, and their relation to the dark and stellar halos

Authors: L. V. Sales, J. F. Navarro, M. G. Abadi, M. Steinmetz

Abstract: We study the population of satellite galaxies formed in a suite of N-body/gasdynamical simulations of galaxy formation in a LCDM universe. We find little spatial or kinematic bias between the dark matter and the satellite population. The velocity dispersion of the satellites is a good indicator of the virial velocity of the halo: σ_{sat}/V_{vir}=0.9 +/- 0.2. Applied to the Milky Way and M31 this… ▽ More We study the population of satellite galaxies formed in a suite of N-body/gasdynamical simulations of galaxy formation in a LCDM universe. We find little spatial or kinematic bias between the dark matter and the satellite population. The velocity dispersion of the satellites is a good indicator of the virial velocity of the halo: σ_{sat}/V_{vir}=0.9 +/- 0.2. Applied to the Milky Way and M31 this gives V_{vir}^{MW}=109 +/- 22$ km/s and V_{vir}^{M31} = 138 +/- 35 km/s, respectively, substantially lower than the rotation speed of their disk components. The detailed kinematics of simulated satellites and dark matter are also in good agreement. By contrast, the stellar halo of the simulated galaxies is kinematically and spatially distinct from the population of surviving satellites. This is because the survival of a satellite depends on mass and on time of accretion; surviving satellites are biased toward low-mass systems that have been recently accreted by the galaxy. Our results support recent proposals for the origin of the systematic differences between stars in the Galactic halo and in Galactic satellites: the elusive ``building blocks'' of the Milky Way stellar halo were on average more massive, and were accreted (and disrupted) earlier than the population of dwarfs that has survived self-bound until the present. △ Less

Submitted 24 May, 2007; v1 submitted 13 April, 2007; originally announced April 2007.

Comments: 13 pages, 11 figures, MNRAS in press. Accepted version with minor changes. Version with high resolution figures available at: http://www.astro.uvic.ca/~lsales/SatPapers/SatPapers.html

Journal ref: Mon.Not.Roy.Astron.Soc.379:1464-1474,2007

Stars beyond Galaxies: The Origin of Extended Luminous Halos around Galaxies

Authors: Mario G. Abadi, Julio F. Navarro, Matthias Steinmetz

Abstract: (Abridged) We use numerical simulations to investigate the origin and structure of the luminous halos that surround isolated galaxies. These stellar structures extend out to several hundred kpc away from a galaxy, and consist of stars shed by merging subunits during the many accretion events that characterize the hierarchical assembly of galaxies. Such origin suggests that outer luminous halos a… ▽ More (Abridged) We use numerical simulations to investigate the origin and structure of the luminous halos that surround isolated galaxies. These stellar structures extend out to several hundred kpc away from a galaxy, and consist of stars shed by merging subunits during the many accretion events that characterize the hierarchical assembly of galaxies. Such origin suggests that outer luminous halos are ubiquitous and that they should appear as an excess of light over extrapolations of the galaxy's inner profile beyond its traditional luminous radius. The mass profile of the accreted stellar component is well approximated by a model where the logarithmic slope steepens monotonically with radius; from -3 at the luminous edge of the galaxy to -4 or steeper near the virial radius of the system. Such spatial distribution is consistent with that of Galactic and M31 globular clusters, suggesting that many of the globulars were brought in by accretion events, in a manner akin to the classic Searle-Zinn scenario. The outer stellar spheroid is supported by a velocity dispersion tensor with a substantial and radially increasing radial anisotropy. These properties distinguish the stellar halo from the dark matter component, which is more isotropic in velocity space, as well as from some tracers of the outer spheroid such as satellite galaxies. Most stars in the outer halo formed in progenitors that have since merged with the central galaxy; very few stars in the halo are contributed by satellites that survive as self-bound entities at the present. These features are in reasonable agreement with recent observations of the outer halo of the MW, of M31, and of other isolated spirals, and suggest that all of these systems underwent an early period of active merging, as envisioned in hierarchical models of galaxy formation. △ Less

Submitted 27 June, 2005; originally announced June 2005.

Comments: Submitted to MNRAS, 13 pages, 12 figures

Journal ref: Mon.Not.Roy.Astron.Soc.365:747-758,2006

Pieces of the puzzle: Ancient substructure in the Galactic disk

Authors: Amina Helmi, J. F. Navarro, B. Nordstrom, J. Holmberg, M. G. Abadi, M. Steinmetz

Abstract: We search for signatures of past accretion events in the Milky Way in the recently published catalogue by Nordstrom et al.(2004), containing accurate spatial and kinematic information as well as metallicities for 13240 nearby stars. To optimize our strategy, we use numerical simulations to characterize the properties of the debris from disrupted satellites. We find that stars with a common proge… ▽ More We search for signatures of past accretion events in the Milky Way in the recently published catalogue by Nordstrom et al.(2004), containing accurate spatial and kinematic information as well as metallicities for 13240 nearby stars. To optimize our strategy, we use numerical simulations to characterize the properties of the debris from disrupted satellites. We find that stars with a common progenitor should show distinct correlations between their orbital parameters; in particular, between the apocentre A and pericentre P, as well as their z-angular momentum (L_z). In the APL-space, such stars are expected to cluster around regions of roughly constant eccentricity. The APL space for the Nordstrom catalogue exhibits a wealth of substructure, much of which can be linked to dynamical perturbations induced by spiral arms and the Galactic bar. However, our analysis also reveals a statistically significant excess of stars on orbits of common (moderate) eccentricity, analogous to the pattern expected for merger debris. Besides being dynamically peculiar, the 274 stars in these substructures have very distinct metallicity and age distributions, providing further evidence of their extra-Galactic provenance. It is possible to identify among these stars, three coherent Groups with characteristic ages and metallicities, that, in all likelihood, correspond to the remains of disrupted satellites. [abridged] △ Less

Submitted 3 November, 2005; v1 submitted 19 May, 2005; originally announced May 2005.

Comments: 18 pages, 21 figures, MNRAS in press. Accepted version with high-resolution figures available from http://www.astro.rug.nl/~ahelmi/research/puzzle.html

Journal ref: Mon.Not.Roy.Astron.Soc. 365 (2006) 1309-1323

Accretion relicts in the solar neighbourhood: debris from omegaCen's parent galaxy

Authors: Andres Meza, Julio F. Navarro, Mario G. Abadi, Matthias Steinmetz

Abstract: We use numerical simulations to investigate the orbital characteristics of tidal debris from satellites whose orbits are dragged into the plane of galactic disks by dynamical friction before disruption. We find that these satellites may deposit a significant fraction of their stars into the disk components of a galaxy, and use our results to motivate the search for accretion relicts in samples o… ▽ More We use numerical simulations to investigate the orbital characteristics of tidal debris from satellites whose orbits are dragged into the plane of galactic disks by dynamical friction before disruption. We find that these satellites may deposit a significant fraction of their stars into the disk components of a galaxy, and use our results to motivate the search for accretion relicts in samples of metal-poor disk stars in the vicinity of the Sun. Satellites disrupted on very eccentric orbits coplanar with the disk are expected to shed stars in "trails" of distinct orbital energy and angular momentum during each pericentric passage. To an observer located between the pericenter and apocenter of such orbits, these trails would show as distinct groupings of stars with low vertical velocity and a broad, symmetric, often double-peaked distribution of Galactocentric radial velocities. One group of stars with these characteristics stands out in available compilations of nearby metal-poor stars. These stars have specific angular momenta similar to that of the globular cluster omegaCen, long hypothesized to be the nucleus of a dwarf galaxy disrupted by the Milky Way tidal field. In addition to their kindred kinematics, stars in the omegaCen group share distinct chemical abundance characteristics, and trace a well-defined track in the [alpha/Fe] versus [Fe/H] plane, consistent with simple closed-box enrichment models and a protracted star formation history. The dynamical and chemical coherence of this group suggests that it consists of stars that once belonged to the dwarf that brought omegaCen into the Galaxy. The presence of this and other "tidal relicts" in the solar neighbourhood suggest an extra-Galactic origin for the presence of nearby stars with odd kinematics and chemistry. △ Less

Submitted 25 January, 2005; v1 submitted 30 August, 2004; originally announced August 2004.

Comments: 11 pages, 11 color figures, new title and minor changes to match MNRAS accepted version

Journal ref: Mon.Not.Roy.Astron.Soc. 359 (2005) 93-103

Tidal Torques and the Orientation of Nearby Disk Galaxies

Authors: Julio F. Navarro, Mario G. Abadi, Matthias Steinmetz

Abstract: We use numerical simulations to investigate the orientation of the angular momentum axis of disk galaxies relative to their surrounding large scale structure. We find that this is closely related to the spatial configuration at turnaround of the material destined to form the galaxy, which is often part of a coherent two-dimensional slab criss-crossed by filaments. The rotation axis is found to a… ▽ More We use numerical simulations to investigate the orientation of the angular momentum axis of disk galaxies relative to their surrounding large scale structure. We find that this is closely related to the spatial configuration at turnaround of the material destined to form the galaxy, which is often part of a coherent two-dimensional slab criss-crossed by filaments. The rotation axis is found to align very well with the intermediate principal axis of the inertia momentum tensor at this time. This orientation is approximately preserved during the ensuing collapse, so that the rotation axis of the resulting disk ends up lying on the plane traced by the protogalactic material at turnaround. This suggests a tendency for disks to align themselves so that their rotation axis is perpendicular to the minor axis of the structure defined by surrounding matter. One example of this trend is provided by our own Galaxy, where the Galactic plane is almost at right angles with the supergalactic plane (SGP) drawn by nearby galaxies; indeed, the SGP latitude of the North Galactic Pole is just 6 degrees. We have searched for a similar signature in catalogs of nearby disk galaxies, and find a significant excess of edge-on spirals (for which the orientation of the disk rotation axis may be determined unambiguously) highly inclined relative to the SGP. This result supports the view that disk galaxies acquire their angular momentum as a consequence of early tidal torques acting during the expansion phase of the protogalactic material. △ Less

Submitted 9 August, 2004; v1 submitted 21 May, 2004; originally announced May 2004.

Comments: 5 pages, 2 figures, accepted for publication in ApJL

Journal ref: Astrophys.J. 613 (2004) L41-L44

Simulations of Galaxy Formation in a Lambda CDM Universe II: The Fine Structure of Simulated Galactic Disks

Authors: Mario G. Abadi, Julio F. Navarro, Matthias Steinmetz, Vincent R. Eke

Abstract: We present a detailed analysis of the dynamical properties of a simulated disk galaxy assembled hierarchically in the Lambda CDM cosmogony. At z=0, two distinct dynamical components are identified on the basis of the orbital parameters of stars in the galaxy: a slowly rotating, centrally concentrated spheroid and a disk-like component largely supported by rotation. These components are also reco… ▽ More We present a detailed analysis of the dynamical properties of a simulated disk galaxy assembled hierarchically in the Lambda CDM cosmogony. At z=0, two distinct dynamical components are identified on the basis of the orbital parameters of stars in the galaxy: a slowly rotating, centrally concentrated spheroid and a disk-like component largely supported by rotation. These components are also recognized in the surface brightness profile of the galaxy, which can be very well approximated by the superposition of an R^{1/4} spheroid and an exponential disk. However, neither does the dynamically-identified spheroid follow de Vaucouleurs' law nor is the disk purely exponential, a resultwhich calls for caution when estimating the importance of the disk from traditional photometric decomposition techniques. The disk may be further decomposed into a thin, dynamically cold component with stars on nearly circularorbits and a hotter, thicker component with orbital parameters transitional between the thin disk and the spheroid. The spheroid is old, and has essentiallyno stars younger than the time elapsed since the last major accretion event ~8 Gyr ago. The majority of thin disk stars, form after the merging activity is over, although a significant fraction ~15% of thin-disk stars are old enough to predate the last major merger event. This unexpected population of old disk stars consists mainly of the tidal debris of satellites whose orbital plane was coincident with the disk and whose orbits were circularized by dynamical friction prior to full disruption. Our results highlight the role of satellite accretion events in shaping the disk and the spheroidal component and reveal some of the clues to the assembly process of a galaxy preserved in the detailed dynamics of old stellar populations. △ Less

Submitted 11 December, 2002; originally announced December 2002.

Comments: 15 pages, 14 figures, submitted to ApJ

Journal ref: Astrophys.J.597:21-34,2003

Simulations of Galaxy Formation in a Lambda CDM Universe I: Dynamical and Photometric Properties of a Simulated Disk Galaxy

Authors: Mario G. Abadi, Julio F. Navarro, Matthias Steinmetz, Vincent R. Eke

Abstract: We analyze the properties of a disk galaxy simulated with unprecedented numerical resolution in the Lambda CDM cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z~1 which lead to the formation of two distinct dynamical components: a spheroid and a disk. The surface brightness profile is very well approximated by the sup… ▽ More We analyze the properties of a disk galaxy simulated with unprecedented numerical resolution in the Lambda CDM cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z~1 which lead to the formation of two distinct dynamical components: a spheroid and a disk. The surface brightness profile is very well approximated by the superposition of an R^{1/4} spheroid and an exponential disk. The surface brightness profile is remarkably similar to that of Sab galaxy UGC615, but the simulated galaxy rotates significantly faster and has a declining rotation curve dominated by the spheroid near the center. The decline in circular velocity is at odds with observation and results from the high concentration of the dark matter and baryonic components, as well as from the relatively high mass-to-light ratio of the stars in the simulation. The simulated galaxy lies ~ 1 mag off the I-band Tully-Fisher relation of late-type spirals, but seems to be in reasonable agreement with Tully-Fisher data on S0 galaxies. The angular momentum of the luminous component is an order of magnitude lower than that of late-type spirals of similar rotation speed. This reflects the dominance of the slowly-rotating, dense spheroidal component. The disk component, on the other hand, has properties rather similar to those of late-type spirals. This suggests that a different form of feedback than adopted in this simulation is required to inhibit the efficient collapse and cooling of gas at high redshift that leads to the formation of the spheroid. Reconciling the properties of disk galaxies with the early collapse and high merging rates characteristic of hierarchical scenarios such as Lambda CDM remains a challenging, yet so far elusive, proposition. △ Less

Submitted 14 November, 2002; originally announced November 2002.

Comments: 15 pages, 13 figures, submitted to ApJ

Journal ref: Astrophys.J.591:499-514,2003

The Power Spectrum of Flux Limited X-Ray Galaxy Cluster Surveys

Authors: Ariel Zandivarez, Mario G. Abadi, Diego G. Lambas

Abstract: We compute the redshift space power spectrum of two X-ray cluster samples: the X-ray Brightest Abell Cluster Sample (XBACS) and the Brightest Cluster Sample (BCS) using the method developed by Feldman, Kaiser & Peacock. The power spectrums derived for these samples are in agreement with determinations of other optical and X-ray cluster samples. For XBACS we find the largest power spectrum amplit… ▽ More We compute the redshift space power spectrum of two X-ray cluster samples: the X-ray Brightest Abell Cluster Sample (XBACS) and the Brightest Cluster Sample (BCS) using the method developed by Feldman, Kaiser & Peacock. The power spectrums derived for these samples are in agreement with determinations of other optical and X-ray cluster samples. For XBACS we find the largest power spectrum amplitude expected given the high richness of this sample ($R \ge 2$). In the range $0.05 \uk < k < 0.4 \uk$ the power spectrum shows a power law behavior $P(k)\propto k^{n}$ with an index $n\simeq-1.2$. In a similar range $0.04 \uk < k < 0.3 \uk $ BCS power spectrum has a smaller amplitude with index $n\simeq-1.0$. We do not find significant evidence for a peak at $k \simeq 0.05 \uk$ suggesting that claims such of feature detections in some cluster samples could relay on artificial inhomogeneities of the data. We compare our results with power spectrum predictions derived by Moscardini et al. within current cosmological models (LCDM and OCDM). For XBACS we find that both models underestimate the amplitude of the power spectrum but for BCS there is reasonably good agreement at $k\gsim 0.03 \uk$ for both models. △ Less

Submitted 22 March, 2001; originally announced March 2001.

Comments: 9 pages (LateX, mn.sty), 9 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.326:147,2001

Self-Similar Shocked Accretion of Collisional Gas with Radiative Cooling

Authors: M. G. Abadi, R. G. Bower, J. F. Navarro

Abstract: We describe similarity solutions that characterize the collapse of collisional gas onto scalefree perturbations in an Einstein-de Sitter universe. We consider the effects of radiative cooling and derive self-similar solutions under the assumption that the cooling function is a power law of density and temperature, $Λ(T,ρ) \propto ρ^{3/2} T$. We use these results to test the ability of Smooth Par… ▽ More We describe similarity solutions that characterize the collapse of collisional gas onto scalefree perturbations in an Einstein-de Sitter universe. We consider the effects of radiative cooling and derive self-similar solutions under the assumption that the cooling function is a power law of density and temperature, $Λ(T,ρ) \propto ρ^{3/2} T$. We use these results to test the ability of Smooth Particle Hydrodynamics (SPH) techniques to follow the collapse and accretion of shocked, rapidly cooling gas in a cosmological context. Our SPH code reproduces the analytical results very well in cases that include or exclude radiative cooling. No substantial deviations from the predicted central mass accretion rates or from the temperature, density, and velocity profiles are observed in well resolved regions inside the shock radius. This test problem lends support to the reliability of SPH techniques to model the complex process of galaxy formation. △ Less

Submitted 26 April, 1999; originally announced April 1999.

Comments: 9 pages, 5 postscript figures, MNRAS submitted

Ram Pressure Stripping of Spiral Galaxies in Clusters

Authors: Mario G. Abadi, Ben Moore, Richard G. Bower

Abstract: We use 3-dimensional SPH/N-BODY simulations to study ram pressure stripping of gas from spiral galaxies orbiting in clusters. We find that the analytic expectation of Gunn & Gott (1972) relating the gravitational restoring force provided by the disk to the ram pressure force, provides a good approximation to the radius that gas will be stripped from a galaxy. However, at small radii it is also i… ▽ More We use 3-dimensional SPH/N-BODY simulations to study ram pressure stripping of gas from spiral galaxies orbiting in clusters. We find that the analytic expectation of Gunn & Gott (1972) relating the gravitational restoring force provided by the disk to the ram pressure force, provides a good approximation to the radius that gas will be stripped from a galaxy. However, at small radii it is also important to consider the potential provided by the bulge component. A spiral galaxy passing through the core of a rich cluster such as Coma, will have its gaseous disk truncated to $\sim 4$ kpc, thus losing $\sim 80%$ of its diffuse gas mass. The timescale for this to occur is a fraction of a crossing time $\sim 10^7$ years. Galaxies orbiting within poorer clusters, or inclined to the direction of motion through the intra-cluster medium will lose significantly less gas. We conclude that ram-pressure alone is insufficient to account for the rapid and widespread truncation of star-formation observed in cluster galaxies, or the morphological transformation of Sab's to S0's that is necessary to explain the Butcher-Oemler effect. △ Less

Submitted 31 March, 1999; v1 submitted 29 March, 1999; originally announced March 1999.

Comments: 8 pages, 7 figures, to be published in MNRAS. Levels added/corrected on figures 3, 4 and 5

Journal ref: Mon.Not.Roy.Astron.Soc.308:947-954,1999

Correlation length of X-ray brightest Abell clusters

Authors: Mario G. Abadi, Diego G. Lambas, Hernán Muriel

Abstract: We compute the cluster auto-correlation function $ξ_{cc}(r)$ of an X-ray flux limited sample of Abell clusters (XBACs, \cite{ebe}). For the total XBACs sample we find a power-law fit $ξ_{cc}=(r/r_0)^γ$ with $r_0=21.1$ Mpc h$^{-1}$and $γ=-1.9$ consistent with the results of $R \ge 1 $ Abell clusters. We also analyze $ξ_{cc}(r)$ for subsamples defined by different X-ray luminosity thresholds where… ▽ More We compute the cluster auto-correlation function $ξ_{cc}(r)$ of an X-ray flux limited sample of Abell clusters (XBACs, \cite{ebe}). For the total XBACs sample we find a power-law fit $ξ_{cc}=(r/r_0)^γ$ with $r_0=21.1$ Mpc h$^{-1}$and $γ=-1.9$ consistent with the results of $R \ge 1 $ Abell clusters. We also analyze $ξ_{cc}(r)$ for subsamples defined by different X-ray luminosity thresholds where we find a weak tendency of larger values of $r_0$ with increasing X-ray luminosity although with a low statistical significance. In the different subsamples analyzed we find $21 < r_0 < 35 $ Mpc h$^{-1}$ and $-1.9< γ< -1.6$. Our analysis suggests that cluster X-ray luminosities may be used for a reliable confrontation of cluster spatial distribution properties in models and observations. △ Less

Submitted 19 June, 1998; originally announced June 1998.

Comments: Accepted for publication in Astrophysical Journal

Cluster-Galaxy Correlations in CDM Models

Authors: M. E. Merchan, M. G. Abadi, D. G. Lambas, C. A. Valotto

Abstract: We study the ability of COBE-normalized CDM models to reproduce observed properties of the distribution of galaxies and clusters using N-body numerical simulations. We analyze the galaxy-galaxy and cluster-galaxy two-point correlation functions, $ξ_{gg}$ and $ξ_{cg}$, in open ($Ω_{0}=0.4, Ω_Λ=0, σ_8=0.75$), and flat ($Ω_{0}=0.3, Ω_Λ=0.7, σ_8=1.05$) CDM models which both reproduce the observed ab… ▽ More We study the ability of COBE-normalized CDM models to reproduce observed properties of the distribution of galaxies and clusters using N-body numerical simulations. We analyze the galaxy-galaxy and cluster-galaxy two-point correlation functions, $ξ_{gg}$ and $ξ_{cg}$, in open ($Ω_{0}=0.4, Ω_Λ=0, σ_8=0.75$), and flat ($Ω_{0}=0.3, Ω_Λ=0.7, σ_8=1.05$) CDM models which both reproduce the observed abundances of rich clusters of galaxies. To compare models with observations we compute projected cross-correlation functions $ω_{gg}$ and $ω_{cg}$ to derive the corresponding $ξ_{gg}$ and $ξ_{cg}$. We use target galaxies selected from Las Campanas Redshift Survey, target clusters selected from the APM Cluster Survey and tracer galaxies from the Edinburgh Durham Sky Survey catalog. △ Less

Submitted 20 November, 1997; originally announced November 1997.

Comments: 13 pages, 6 figures, to appear in: The Astrophysical Journal

Analysis of Galaxy Formation with Hydrodynamics

Authors: P. B. Tissera, D. G. Lambas, M. G. Abadi

Abstract: We present a hydrodynamical code based on the Smooth Particle Hydrodynamics technique implemented in an AP3M code aimed at solving the hydrodynamical and gravitational equations in a cosmological frame. We analyze the ability of the code to reproduce standard tests and perform numerical simulations to study the formation of galaxies in a typical region of a CDM model. These numerical simulations… ▽ More We present a hydrodynamical code based on the Smooth Particle Hydrodynamics technique implemented in an AP3M code aimed at solving the hydrodynamical and gravitational equations in a cosmological frame. We analyze the ability of the code to reproduce standard tests and perform numerical simulations to study the formation of galaxies in a typical region of a CDM model. These numerical simulations include gas and dark matter particles and take into account physical processes such as shock waves, radiative cooling, and a simplified model of star formation. Several observed properties of normal galaxies such as $M_{gas}/M_{total}$ ratios, the luminosity function and the Tully-Fisher relation are analyzed within the limits imposed by numerical resolution. △ Less

Submitted 28 March, 1996; originally announced March 1996.

Comments: 21 pages, 2 postscript tables. Submitted MNRAS 04.03.96