-
Testing strong lensing subhalo detection with a cosmological simulation
Authors:
Qiuhan He,
James Nightingale,
Andrew Robertson,
Aristeidis Amvrosiadis,
Shaun Cole,
Carlos S. Frenk,
Richard Massey,
Ran Li,
Nicola C. Amorisco,
R. Benton Metcalf,
Xiaoyue Cao,
Amy Etherington
Abstract:
Strong gravitational lensing offers a compelling test of the cold dark matter paradigm, as it allows for subhaloes with masses of $\sim10^{9}$ M$_\odot$ and below to be detected. We test commonly-used techniques for detecting subhaloes superposed in images of strongly lensed galaxies. For the lens we take a simulated galaxy in a $\sim10^{13}$ M$_\odot$ halo grown in a high-resolution cosmological…
▽ More
Strong gravitational lensing offers a compelling test of the cold dark matter paradigm, as it allows for subhaloes with masses of $\sim10^{9}$ M$_\odot$ and below to be detected. We test commonly-used techniques for detecting subhaloes superposed in images of strongly lensed galaxies. For the lens we take a simulated galaxy in a $\sim10^{13}$ M$_\odot$ halo grown in a high-resolution cosmological hydrodynamical simulation, which we view from two different directions. Though the resolution is high, we note the simulated galaxy still has an artificial core which adds additional complexity to the baryon dominated region. To remove particle noise, we represent the projected galaxy mass distribution by a series of Gaussian profiles which precisely capture the features of the projected galaxy. We first model the lens mass as a (broken) power-law density profile and then search for small haloes. Of the two projections, one has a regular elliptical shape, while the other has distinct deviations from an elliptical shape. For the former, the broken power-law model gives no false positives and correctly recovers the mass of the superposed small halo, but for the latter we find false positives and the inferred halo mass is overestimated by $\sim4-5$ times. We then use a more complex model in which the lens mass is decomposed into stellar and dark matter components. In this case, we show that we can capture the simulated galaxy's complex projected structures and correctly infer the input small halo.
△ Less
Submitted 26 September, 2022; v1 submitted 21 February, 2022;
originally announced February 2022.
-
Automated galaxy-galaxy strong lens modelling: no lens left behind
Authors:
Amy Etherington,
James W. Nightingale,
Richard Massey,
XiaoYue Cao,
Andrew Robertson,
Nicola C. Amorisco,
Aristeidis Amvrosiadis,
Shaun Cole,
Carlos S. Frenk,
Qiuhan He,
Ran Li,
Sut-Ieng Tam
Abstract:
The distribution of dark and luminous matter can be mapped around galaxies that gravitationally lens background objects into arcs or Einstein rings. New surveys will soon observe hundreds of thousands of galaxy lenses, and current, labour-intensive analysis methods will not scale up to this challenge. We instead develop a fully automatic, Bayesian method which we use to fit a sample of 59 lenses i…
▽ More
The distribution of dark and luminous matter can be mapped around galaxies that gravitationally lens background objects into arcs or Einstein rings. New surveys will soon observe hundreds of thousands of galaxy lenses, and current, labour-intensive analysis methods will not scale up to this challenge. We instead develop a fully automatic, Bayesian method which we use to fit a sample of 59 lenses imaged by the Hubble Space Telescope in uniform conditions. We set out to \textit{leave no lens behind} and focus on ways in which automated fits fail in a small handful of lenses, describing adjustments to the pipeline that allows us to infer accurate lens models. Our pipeline ultimately fits {\em all} 59 lenses in our sample, with a high success rate key because catastrophic outliers would bias large samples with small statistical errors. Machine Learning techniques might further improve the two most difficult steps: subtracting foreground lens light and initialising a first, approximate lens model. After that, increasing model complexity is straightforward. We find a mean $\sim1\%$ measurement precision on the measurement of the Einstein radius across the lens sample which {\em does not degrade with redshift} up to at least $z=0.7$ -- in stark contrast to other techniques used to study galaxy evolution, like stellar dynamics. Our \texttt{PyAutoLens} software is open source, and is also installed in the Science Data Centres of the ESA Euclid mission.
△ Less
Submitted 12 January, 2023; v1 submitted 18 February, 2022;
originally announced February 2022.
-
Implications for Galaxy Formation Models from Observations of Globular Clusters around Ultra-Diffuse Galaxies
Authors:
Teymoor Saifollahi,
Dennis Zaritsky,
Ignacio Trujillo,
Reynier F. Peletier,
Johan H. Knapen,
Nicola Amorisco,
Michael A. Beasley,
Richard Donnerstein
Abstract:
We present an analysis of Hubble Space Telescope observations of globular clusters (GCs) in six ultra-diffuse galaxies (UDGs) in the Coma cluster, a sample that represents UDGs with large effective radii ($R_{\rm e}$), and use the results to evaluate competing formation models. We eliminate two significant sources of systematic uncertainty in the determination of the number of GCs, $N_{\rm GC}$ by…
▽ More
We present an analysis of Hubble Space Telescope observations of globular clusters (GCs) in six ultra-diffuse galaxies (UDGs) in the Coma cluster, a sample that represents UDGs with large effective radii ($R_{\rm e}$), and use the results to evaluate competing formation models. We eliminate two significant sources of systematic uncertainty in the determination of the number of GCs, $N_{\rm GC}$ by using sufficiently deep observations that (i) reach the turnover of the GC luminosity function and (ii) provide a sufficient number of GCs with which to measure the GC number radial distribution. We find that $N_{\rm GC}$ for these galaxies is on average $\sim$20, which implies an average total mass, $M_{\rm total}$, $\sim$ $10^{11}$ $M_{\odot}$ when applying the relation between $N_{\rm GC}$ and $M_{\rm total}$. This value of $N_{\rm GC}$ lies at the upper end of the range observed for dwarf galaxies of the same stellar mass and is roughly a factor of two larger than the mean. The GC luminosity function, radial profile and average colour are more consistent with those observed for dwarf galaxies than with those observed for the more massive ($L^*$) galaxies, while both the radial and azimuthal GC distributions closely follow those of the stars in the host galaxy. Finally, we discuss why our observations, specifically the GC number and GC distribution around these six UDGs, pose challenges for several of the currently favoured UDG formation models.
△ Less
Submitted 3 February, 2022; v1 submitted 27 January, 2022;
originally announced January 2022.
-
Cold dark matter subhaloes at arbitrarily low masses
Authors:
Nicola C. Amorisco
Abstract:
A defining prediction of the cold dark matter (CDM) cosmological model is the existence of a very large population of low mass haloes, down to planet-size masses. However, their fate as they are accreted onto haloes many orders of magnitude more massive remains fundamentally uncertain. A number of numerical explorations have found subhaloes to be very resilient to tides, but resolution limits make…
▽ More
A defining prediction of the cold dark matter (CDM) cosmological model is the existence of a very large population of low mass haloes, down to planet-size masses. However, their fate as they are accreted onto haloes many orders of magnitude more massive remains fundamentally uncertain. A number of numerical explorations have found subhaloes to be very resilient to tides, but resolution limits make it difficult to explore tidal evolution at arbitrarily low masses. What are the structural properties of heavily stripped subhaloes? Do tidal effects destroy low-mass CDM subhaloes? Here we focus on cosmologically motivated subhaloes orbiting CDM hosts, and show that subhaloes of any initial mass can be stripped to arbitrarily small mass fractions. We show that previous numerical results can be reproduced by a simple model that describes tidal evolution as a progressive `peeling' in energy space and subsequent re-virialization. Tidal heating can effectively be neglected at all masses, because its importance i) does not increase for subhaloes of decreasing mass at accretion, and ii) it decreases as stripping proceeds. This allows us to predict analytically the structural properties of subhaloes with any initial mass and arbitrary degrees of mass loss. Under the hypotheses that CDM haloes have centrally divergent density profiles and approximately isotropic phase space distributions, our results prove that subhaloes of very low masses are a robust prediction of the CDM model: as for haloes, CDM subhalo populations extend and are abundant down to very low masses.
△ Less
Submitted 1 November, 2021;
originally announced November 2021.
-
Systematic errors induced by the elliptical power-law model in galaxy-galaxy strong lens modeling
Authors:
Xiaoyue Cao,
Ran Li,
J. W. Nightingale,
Richard Massey,
Andrew Robertson,
Carlos S. Frenk,
Aristeidis Amvrosiadis,
Nicola C. Amorisco,
Qiuhan He,
Amy Etherington,
Shaun Cole,
Kai Zhu
Abstract:
The elliptical power-law (EPL) model of the mass in a galaxy is widely used in strong gravitational lensing analyses. However, the distribution of mass in real galaxies is more complex. We quantify the biases due to this model mismatch by simulating and then analysing mock {\it Hubble Space Telescope} imaging of lenses with mass distributions inferred from SDSS-MaNGA stellar dynamics data. We find…
▽ More
The elliptical power-law (EPL) model of the mass in a galaxy is widely used in strong gravitational lensing analyses. However, the distribution of mass in real galaxies is more complex. We quantify the biases due to this model mismatch by simulating and then analysing mock {\it Hubble Space Telescope} imaging of lenses with mass distributions inferred from SDSS-MaNGA stellar dynamics data. We find accurate recovery of source galaxy morphology, except for a slight tendency to infer sources to be more compact than their true size. The Einstein radius of the lens is also robustly recovered with 0.1% accuracy, as is the global density slope, with 2.5% relative systematic error, compared to the 3.4% intrinsic dispersion. However, asymmetry in real lenses also leads to a spurious fitted `external shear' with typical strength, $γ_{\rm ext}=0.015$. Furthermore, time delays inferred from lens modelling without measurements of stellar dynamics are typically underestimated by $\sim$5%. Using such measurements from a sub-sample of 37 lenses would bias measurements of the Hubble constant $H_0$ by $\sim$9%. Although this work is based on a particular set of MaNGA galaxies, and the specific value of the detected biases may change for another set of strong lenses, our results strongly suggest the next generation cosmography needs to use more complex lens mass models.
△ Less
Submitted 21 November, 2021; v1 submitted 27 October, 2021;
originally announced October 2021.
-
Dwarf stellar haloes: a powerful probe of small-scale galaxy formation and the nature of dark matter
Authors:
Alis J. Deason,
Sownak Bose,
Azadeh Fattahi,
Nicola C. Amorisco,
Wojciech Hellwing,
Carlos S. Frenk
Abstract:
We use N-body cosmological simulations and empirical galaxy models to study the merger history of dwarf-mass galaxies (with M_halo~10^10 M_Sun). Our input galaxy models describe the stellar mass-halo mass relation, and the galaxy occupation fraction. The number of major and minor mergers depends on the type of dark matter; in particular, minor mergers are greatly suppressed in warm dark matter mod…
▽ More
We use N-body cosmological simulations and empirical galaxy models to study the merger history of dwarf-mass galaxies (with M_halo~10^10 M_Sun). Our input galaxy models describe the stellar mass-halo mass relation, and the galaxy occupation fraction. The number of major and minor mergers depends on the type of dark matter; in particular, minor mergers are greatly suppressed in warm dark matter models. In addition, the number of mergers that bring in stars is strongly dependent on the galaxy occupation model. For example, minor mergers are negligible for stellar halo growth in models with a high mass threshold for galaxy formation (i.e. 10^9.3 M_Sun at z=0). Moreover, this threshold for galaxy formation can also determine the relative difference (if any) between the stellar haloes of satellite and field dwarfs. Using isolated simulations of dwarf-dwarf mergers, we show that the relative frequency of major and minor mergers predict very different stellar haloes: Typically, "intermediate" dark matter merger ratios (~1:5) maximise the growth of distant stellar haloes. We discuss the observability of dwarf stellar haloes and find that the surface brightness of these features are incredibly faint. However, when several dwarfs are stacked together models that form particularly rich stellar haloes could be detectable. Finally, we show that stellar streams in the Galactic halo overlapping in phase-space with known dwarf satellites are likely remnants of their stripped stellar haloes. The mere existence of dwarf stellar haloes can already put constraints on some small-scale models, and thus observational probes should be a high priority.
△ Less
Submitted 1 January, 2022; v1 submitted 11 October, 2021;
originally announced October 2021.
-
Galaxy-galaxy strong lens perturbations: line-of-sight haloes versus lens subhaloes
Authors:
Qiuhan He,
Ran Li,
Carlos S. Frenk,
James Nightingale,
Shaun Cole,
Nicola C. Amorisco,
Richard Massey,
Andrew Robertson,
Amy Etherington,
Aristeidis Amvrosiadis,
Xiaoyue Cao
Abstract:
We rederive the number density of intervening line-of-sight haloes relative to lens subhaloes in galaxy-galaxy strong lensing observations, where these perturbers can generate detectable image fluctuations. Previous studies have calculated the detection limit of a line-of-sight small-mass dark halo by comparing the lensing deflection angles it would cause, to those caused by a subhalo within the l…
▽ More
We rederive the number density of intervening line-of-sight haloes relative to lens subhaloes in galaxy-galaxy strong lensing observations, where these perturbers can generate detectable image fluctuations. Previous studies have calculated the detection limit of a line-of-sight small-mass dark halo by comparing the lensing deflection angles it would cause, to those caused by a subhalo within the lens. However, this overly simplifies the difference in observational consequences between a subhalo and a line-of-sight halo. Furthermore, it does not take into account degeneracies between an extra subhalo and the uncertain properties of the main lens. More in keeping with analyses of real-world observations, we regard a line-of-sight halo as detectable only if adding it to a smooth model generates a statistically significant improvement in the reconstructed image. We find that the number density of detectable line-of-sight perturbers has been overestimated by as much as a factor of two in the previous literature. For typical lensing geometries and configurations, very deep imaging is sensitive to twice as many line-of-sight perturbers as subhaloes, but moderate depth imaging is sensitive to only slightly more line-of-sight perturbers than subhaloes.
△ Less
Submitted 26 September, 2022; v1 submitted 9 October, 2021;
originally announced October 2021.
-
Halo concentration strengthens dark matter constraints in galaxy-galaxy strong lensing analyses
Authors:
Nicola C. Amorisco,
James Nightingale,
Qiuhan He,
Aristeidis Amvrosiadis,
Xiaoyue Cao,
Shaun Cole,
Amy Etherington,
Carlos S. Frenk,
Ran Li,
Richard Massey,
Andrew Robertson
Abstract:
A defining prediction of the cold dark matter (CDM) cosmological model is the existence of a very large population of low-mass haloes. This population is absent in models in which the dark matter particle is warm (WDM). These alternatives can, in principle, be distinguished observationally because halos along the line-of-sight can perturb galaxy-galaxy strong gravitational lenses. Furthermore, the…
▽ More
A defining prediction of the cold dark matter (CDM) cosmological model is the existence of a very large population of low-mass haloes. This population is absent in models in which the dark matter particle is warm (WDM). These alternatives can, in principle, be distinguished observationally because halos along the line-of-sight can perturb galaxy-galaxy strong gravitational lenses. Furthermore, the WDM particle mass could be deduced because the cut-off in their halo mass function depends on the mass of the particle. We systematically explore the detectability of low-mass haloes in WDM models by simulating and fitting mock lensed images. Contrary to previous studies, we find that halos are harder to detect when they are either behind or in front of the lens. Furthermore, we find that the perturbing effect of haloes increases with their concentration: detectable haloes are systematically high-concentration haloes, and accounting for the scatter in the mass-concentration relation boosts the expected number of detections by as much as an order of magnitude. Haloes have lower concentration for lower particle masses and this further suppresses the number of detectable haloes beyond the reduction arising from the lower halo abundances alone. Taking these effects into account can make lensing constraints on the value of the mass function cut-off at least an order of magnitude more stringent than previously appreciated.
△ Less
Submitted 31 August, 2021;
originally announced September 2021.
-
Giant cold satellites from low-concentration haloes
Authors:
Nicola C. Amorisco
Abstract:
The dwarf satellite galaxies of the Milky Way Crater II and Antlia II have uncommonly low dynamical mass densities, due to their large size and low velocity dispersion. Previous work have failed to identify formation scenarios within the $Λ$CDM framework and have invoked cored dark matter haloes, processed by tides. I show that the tidal evolution of $Λ$CDM NFW haloes is richer than previously rec…
▽ More
The dwarf satellite galaxies of the Milky Way Crater II and Antlia II have uncommonly low dynamical mass densities, due to their large size and low velocity dispersion. Previous work have failed to identify formation scenarios within the $Λ$CDM framework and have invoked cored dark matter haloes, processed by tides. I show that the tidal evolution of $Λ$CDM NFW haloes is richer than previously recognised: tidal heating causes the innermost regions of haloes that fall short of the mass-concentration relation to expand significantly, resulting in the formation of giant, kinematically cold satellites like Crater II and Antlia II. Furthermore, while the satellite is reaching apocenter, extra-tidal material can cause an even more inflated appearance. When present, as likely for the larger Antlia II, nominally unbound material can be recognised thanks to its somewhat hotter kinematics and line-of-sight velocity gradient. Contrary to other formation scenarios, Crater II and Antlia II may well have experienced very little mass loss, as in fact hinted by their observed metallicity. If indeed a satellite of NGC1052, tidal evolution of a low-concentration halo may similarly have led to the formation of NGC1052-DF2.
△ Less
Submitted 16 January, 2019;
originally announced January 2019.
-
First constraints on Fuzzy Dark Matter from the dynamics of stellar streams in the Milky Way
Authors:
Nicola C. Amorisco,
A. Loeb
Abstract:
We present a novel method to constrain the mass of ultra-light bosons as the dark matter using stellar streams formed by disrupting Globular Clusters in the Milky Way. The turbulent density field of Fuzzy Dark Matter (FDM) haloes results in perturbations and dynamical heating of thin streams. Using numerical simulations based on an effective model, we explore the magnitude of this phenomenon and s…
▽ More
We present a novel method to constrain the mass of ultra-light bosons as the dark matter using stellar streams formed by disrupting Globular Clusters in the Milky Way. The turbulent density field of Fuzzy Dark Matter (FDM) haloes results in perturbations and dynamical heating of thin streams. Using numerical simulations based on an effective model, we explore the magnitude of this phenomenon and show that this is observable for the range of axion masses $m_a$ that is interesting for alleviating the `small-scale problems' of $Λ$CDM. We derive an analytical model for the thickening of thin stellar streams and obtain an early conservative lower limit for the boson mass of $m_a > 1.5 \times 10^{-22}~$eV, using pre-Gaia literature data for six Milky Way streams and after marginalizing over physical parameters. This demonstrates the great promise for using this novel dynamical method as a fully independent probe of FDM, to complement results based on Lyman-$α$ forest data.
△ Less
Submitted 1 August, 2018;
originally announced August 2018.
-
Globular cluster populations and the kinematical fingerprints of minor mergers
Authors:
Nicola C. Amorisco
Abstract:
We use Monte Carlo $Λ$CDM assembly histories, minor-merger N-body simulations and empirical relations between halo mass and the globular cluster (GC) abundance to study the kinematical properties of halo GCs in massive galaxies, with $M_{vir}(z=0)=10^{13.5}M_\odot$. While the accreted stellar halo is dominated by the contributions of massive satellites, we show that satellites with low virial mass…
▽ More
We use Monte Carlo $Λ$CDM assembly histories, minor-merger N-body simulations and empirical relations between halo mass and the globular cluster (GC) abundance to study the kinematical properties of halo GCs in massive galaxies, with $M_{vir}(z=0)=10^{13.5}M_\odot$. While the accreted stellar halo is dominated by the contributions of massive satellites, we show that satellites with low virial mass (i.e. low satellite-to-host virial mass ratio, VMRs) are important contributors to the population of accreted GCs. The relative contribution of accretion events with low VMRs is highest for the halo population of blue GCs and gradually decreases for red GCs and accreted stars. As a consequence of the reduced efficiency of dynamical friction on minor mergers, our populations of cosmologically accreted blue GCs are systematically more spatially extended and have higher velocity dispersions than for red GCs, in agreement with observations. For the same reason, assembly histories including a higher fraction of minor mergers result in more spatially extended GC populations. GC line-of-sight velocity distributions featuring negative values of the kurtosis $κ$, as recently observed, are ubiquitous in our models. Therefore, $κ<0$ is not at odds with an accretion scenario, and in fact a fingerprint of the important contribution of minor mergers. However, our populations of accreted GCs remain mostly radially biased, with profiles of the anisotropy parameter $β$ that are mildly radial in the center ($β(r<10~{\rm kpc})\sim0.2$) and strongly radially anisotropic at large galactocentric distances ($β(r>30~{\rm kpc})\gtrsim0.6$), for both red and blue populations.
△ Less
Submitted 26 October, 2018; v1 submitted 2 February, 2018;
originally announced February 2018.
-
Halo Substructure in the SDSS-Gaia Catalogue : Streams and Clumps
Authors:
G. C. Myeong,
N. W. Evans,
V. Belokurov,
N. C. Amorisco,
S. Koposov
Abstract:
We use the SDSS-Gaia Catalogue to identify six new pieces of halo substructure. SDSS-Gaia is an astrometric catalogue that exploits SDSS data release 9 to provide first epoch photometry for objects in the Gaia source catalogue. We use a version of the catalogue containing $245\,316$ stars with all phase space coordinates within a heliocentric distance of $\sim 10$ kpc. We devise a method to assess…
▽ More
We use the SDSS-Gaia Catalogue to identify six new pieces of halo substructure. SDSS-Gaia is an astrometric catalogue that exploits SDSS data release 9 to provide first epoch photometry for objects in the Gaia source catalogue. We use a version of the catalogue containing $245\,316$ stars with all phase space coordinates within a heliocentric distance of $\sim 10$ kpc. We devise a method to assess the significance of halo substructures based on their clustering in velocity space. The two most substantial structures are multiple wraps of a stream which has undergone considerable phase mixing (S1, with 94 members) and a kinematically cold stream (S2, with 61 members). The member stars of S1 have a median position of ($X,Y,Z$) = ($8.12, -0.22, 2.75$) kpc and a median metallicity of [Fe/H] $= -1.78$. The stars of S2 have median coordinates ($X,Y,Z$) = ($8.66, 0.30, 0.77$) kpc and a median metallicity of [Fe/H] $= -1.91$. They lie in velocity space close to some of the stars in the stream reported by Helmi et al. (1999). By modelling, we estimate that both structures had progenitors with virial masses $\approx 10^{10} M_\odot$ and infall times $\gtrsim 9$ Gyr ago. Using abundance matching, these correspond to stellar masses between $10^6$ and $10^7 M_\odot$. These are somewhat larger than the masses inferred through the mass-metallicity relation by factors of 5 to 15. Additionally, we identify two further substructures (S3 and S4 with 55 and 40 members) and two clusters or moving groups (C1 and C2 with 24 and 12) members. In all 6 cases, clustering in kinematics is found to correspond to clustering in both configuration space and metallicity, adding credence to the reliability of our detections.
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
The virial mass distribution of ultra-diffuse galaxies in clusters and groups
Authors:
Nicola C. Amorisco
Abstract:
We use the observed abundances of ultra-diffuse galaxies (UDGs) in clusters and groups and ΛCDM subhalo mass functions to put constraints on the distribution {of present-day halo masses of satellite} UDGs. If all of the most massive subhaloes in the cluster host a UDG, UDGs occupy all subhaloes with \log M_{sub}/M_\odot\gtrsim11. For a model in which the efficiency of UDG formation is higher aroun…
▽ More
We use the observed abundances of ultra-diffuse galaxies (UDGs) in clusters and groups and ΛCDM subhalo mass functions to put constraints on the distribution {of present-day halo masses of satellite} UDGs. If all of the most massive subhaloes in the cluster host a UDG, UDGs occupy all subhaloes with \log M_{sub}/M_\odot\gtrsim11. For a model in which the efficiency of UDG formation is higher around some characteristic halo mass, higher fractions of massive UDGs require larger spreads in the UDG mass distribution. In a cluster with a virial mass of 10^{15}M_\odot, the 90% upper limit for the fraction of UDGs with \log M_{sub}/M_\odot>12 is 7%, occupying 70% of all cluster subhaloes above the same mass. To reproduce the observed abundances, however, the mass distribution of satellite UDGs has to be broad, with >30% having \log M_{sub}/M_\odot<10.9. This strongly supports that UDGs are part of a continuous distribution in which a majority are hosted by low mass haloes. The abundance of satellite UDGs may fall short of the linear relation with the cluster/group mass M_{host} in low-mass hosts, \log M_{host}/M_\odot\sim 12. Characterising these deviations -- or the lack thereof -- will allow for stringent constraints on the UDG virial mass distribution.
△ Less
Submitted 18 January, 2018; v1 submitted 5 December, 2017;
originally announced December 2017.
-
Galaxies with Shells in the Illustris Simulation: Metallicity Signatures
Authors:
Ana-Roxana Pop,
Annalisa Pillepich,
Nicola C. Amorisco,
Lars Hernquist
Abstract:
Stellar shells are low surface brightness arcs of overdense stellar regions, extending to large galactocentric distances. In a companion study, we identified 39 shell galaxies in a sample of 220 massive ellipticals ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_\odot$) from the Illustris cosmological simulation. We used stellar history catalogs to trace the history of each individual s…
▽ More
Stellar shells are low surface brightness arcs of overdense stellar regions, extending to large galactocentric distances. In a companion study, we identified 39 shell galaxies in a sample of 220 massive ellipticals ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_\odot$) from the Illustris cosmological simulation. We used stellar history catalogs to trace the history of each individual star particle inside the shell substructures, and we found that shells in high-mass galaxies form through mergers with massive satellites (stellar mass ratios $μ_{\mathrm{stars}}\gtrsim1:10$). Using the same sample of shell galaxies, the current study extends the stellar history catalogs in order to investigate the metallicity of stellar shells around massive galaxies. Our results indicate that outer shells are often times more metal-rich than the surrounding stellar material in a galaxy's halo. For a galaxy with two different satellites forming $z=0$ shells, we find a significant difference in the metallicity of the shells produced by each progenitor. We also find that shell galaxies have higher mass-weighted logarithmic metallicities ([Z/H]) at $2$-$4\,\mathrm{R}_{\mathrm{eff}}$ compared to galaxies without shells. Our results indicate that observations comparing the metallicities of stars in tidal features, such as shells, to the average metallicities in the stellar halo can provide information about the assembly histories of galaxies.
△ Less
Submitted 4 August, 2017;
originally announced August 2017.
-
Formation and Incidence of Shell Galaxies in the Illustris Simulation
Authors:
Ana-Roxana Pop,
Annalisa Pillepich,
Nicola C. Amorisco,
Lars Hernquist
Abstract:
Shells are low surface brightness tidal debris that appear as interleaved caustics with large opening angles, often situated on both sides of the galaxy center. In this paper, we study the incidence and formation processes of shell galaxies in the cosmological gravity+hydrodynamics Illustris simulation. We identify shells at redshift z=0 using stellar surface density maps, and we use stellar histo…
▽ More
Shells are low surface brightness tidal debris that appear as interleaved caustics with large opening angles, often situated on both sides of the galaxy center. In this paper, we study the incidence and formation processes of shell galaxies in the cosmological gravity+hydrodynamics Illustris simulation. We identify shells at redshift z=0 using stellar surface density maps, and we use stellar history catalogs to trace the birth, trajectory and progenitors of each individual star particle contributing to the tidal feature. Out of a sample of the 220 most massive galaxies in Illustris ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_{\odot}$), $18\%\pm3\%$ of the galaxies exhibit shells. This fraction increases with increasing mass cut: higher mass galaxies are more likely to have stellar shells. Furthermore, the fraction of massive galaxies that exhibit shells decreases with increasing redshift. We find that shell galaxies observed at redshift $z=0$ form preferentially through relatively major mergers ($\gtrsim$1:10 in stellar mass ratio). Progenitors are accreted on low angular momentum orbits, in a preferred time-window between $\sim$4 and 8 Gyrs ago. Our study indicates that, due to dynamical friction, more massive satellites are allowed to probe a wider range of impact parameters at accretion time, while small companions need almost purely radial infall trajectories in order to produce shells. We also find a number of special cases, as a consequence of the additional complexity introduced by the cosmological setting. These include galaxies with multiple shell-forming progenitors, satellite-of-satellites also forming shells, or satellites that fail to produce shells due to multiple major mergers happening in quick succession.
△ Less
Submitted 21 August, 2018; v1 submitted 19 June, 2017;
originally announced June 2017.
-
Deadly dark matter cusps vs faint and extended star clusters: Eridanus II and Andromeda XXV
Authors:
Nicola C. Amorisco
Abstract:
The recent detection of two faint and extended star clusters in the central regions of two Local Group dwarf galaxies, Eridanus II and Andromeda XXV, raises the question of whether clusters with such low densities can survive the tidal field of cold dark matter haloes with central density cusps. Using both analytic arguments and a suite of collisionless N-body simulations, I show that these cluste…
▽ More
The recent detection of two faint and extended star clusters in the central regions of two Local Group dwarf galaxies, Eridanus II and Andromeda XXV, raises the question of whether clusters with such low densities can survive the tidal field of cold dark matter haloes with central density cusps. Using both analytic arguments and a suite of collisionless N-body simulations, I show that these clusters are extremely fragile and quickly disrupted in the presence of central cusps $ρ\sim r^{-α}$ with $α\gtrsim 0.2$. Furthermore, the scenario in which the clusters where originally more massive and sank to the center of the halo requires extreme fine tuning and does not naturally reproduce the observed systems. In turn, these clusters are long lived in cored haloes, whose central regions are safe shelters for $α\lesssim 0.2$. The only viable scenario for hosts that have preserved their primoridal cusp to the present time is that the clusters formed at rest at the bottom of the potential, which is easily tested by measurement of the clusters proper velocity within the host. This offers means to readily probe the central density profile of two dwarf galaxies as faint as $L_V\sim5\times 10^5 L_\odot$ and $L_V\sim6\times10^4 L_\odot$, in which stellar feedback is unlikely to be effective.
△ Less
Submitted 17 May, 2017; v1 submitted 20 April, 2017;
originally announced April 2017.
-
The accreted stellar halo as a window on halo assembly in L* galaxies
Authors:
Nicola C. Amorisco
Abstract:
Theory and observations agree that the accreted stellar halos (ASHs) of Milky Way-like galaxies display significant scatter. I take advantage of this stochasticity to invert the link between halo assembly history (HAH) and ASH, using mock ASHs corresponding to 750 $Λ$CDM HAHs, sharing a final virial mass of $M_{h}(z=0)=10^{12.25}M_\odot$. Hosts with poor/rich ASHs assemble following orthogonal gro…
▽ More
Theory and observations agree that the accreted stellar halos (ASHs) of Milky Way-like galaxies display significant scatter. I take advantage of this stochasticity to invert the link between halo assembly history (HAH) and ASH, using mock ASHs corresponding to 750 $Λ$CDM HAHs, sharing a final virial mass of $M_{h}(z=0)=10^{12.25}M_\odot$. Hosts with poor/rich ASHs assemble following orthogonal growth-patterns. Hosts with rich ASHs experience accretion events (AEs) with high virial mass ratios (HVMRs, $M_s/M_h\gtrsim 0.1$) at $0.5\lesssim z_{infall}\lesssim1.5$, in a phase of fast growth. This maximizes the accreted stellar mass under the condition these satellites are disrupted by $z=0$. At similar times, hosts with poor ASHs grow slowly through minor mergers, with only very recent HVMR AEs: this results in a globally more abundant satellite population and in distinctive surviving massive satellites (stellar mass $\log M_{s,*}/M_\odot\gtrsim 9$). Several properties of the Milky Way are in agreement with the predictions of this framework for hosts with poor, concentrated ASHs, including: i) the recent infall of Sagittarius and Magellanic Clouds, ii) the likely higher-than-average concentration of its dark halo, iii) the signatures of fast chemical enrichment of a sizable fraction of its halo stellar populations.
△ Less
Submitted 27 March, 2017; v1 submitted 10 January, 2017;
originally announced January 2017.
-
The globular cluster systems of 54 Coma ultra-diffuse galaxies: statistical constraints from HST data
Authors:
N. C. Amorisco,
A. Monachesi,
A. Agnello,
S. D. M. White
Abstract:
We use data from the HST Coma Cluster Treasury program to assess the richness of the Globular Cluster Systems (GCSs) of 54 Coma ultra-diffuse galaxies (UDGs), 18 of which have a half-light radius exceeding 1.5 kpc. We use a hierarchical Bayesian method tested on a large number of mock datasets to account consistently for the high and spatially varying background counts in Coma. These include both…
▽ More
We use data from the HST Coma Cluster Treasury program to assess the richness of the Globular Cluster Systems (GCSs) of 54 Coma ultra-diffuse galaxies (UDGs), 18 of which have a half-light radius exceeding 1.5 kpc. We use a hierarchical Bayesian method tested on a large number of mock datasets to account consistently for the high and spatially varying background counts in Coma. These include both background galaxies and intra-cluster GCs (ICGCs), which are disentangled from the population of member GCs in a probabilistic fashion. We find no candidate for a GCS as rich as that of the Milky Way, our sample has GCSs typical of dwarf galaxies. For the standard relation between GCS richness and halo mass, 33 galaxies have a virial mass $M_{vir}\leq10^{11}M_\odot$ at 90% probability. Only three have $M_{vir}>10^{11}M_\odot$ with the same confidence. The mean colour and spread in colour of the UDG GCs are indistinguishable from those of the abundant population of ICGCs. The majority of UDGs in our sample are consistent with the relation between stellar mass and GC richness of 'normal' dwarf galaxies. Nine systems, however, display GCSs that are richer by a factor of 3 or more (at 90% probability). Six of these have sizes $\lesssim1.4$ kpc. Our results imply that the physical mechanisms responsible for the extended size of the UDGs and for the enhanced GC richness of some cluster dwarfs are at most weakly correlated.
△ Less
Submitted 11 January, 2018; v1 submitted 5 October, 2016;
originally announced October 2016.
-
Exploring the connection between stellar halo profiles and accretion histories in $L_*$ galaxies
Authors:
Nicola C. Amorisco
Abstract:
I use a library of controlled minor merger N-body simulations, a particle tagging technique and Monte Carlo generated $Λ$CDM accretion histories to study the highly stochastic process of stellar deposition onto the accreted stellar halos (ASHs) of $L_*$ galaxies. I explore the main physical mechanisms that drive the connection between the accretion history and the density profile of the ASH. I fin…
▽ More
I use a library of controlled minor merger N-body simulations, a particle tagging technique and Monte Carlo generated $Λ$CDM accretion histories to study the highly stochastic process of stellar deposition onto the accreted stellar halos (ASHs) of $L_*$ galaxies. I explore the main physical mechanisms that drive the connection between the accretion history and the density profile of the ASH. I find that: i) through dynamical friction, more massive satellites are more effective at delivering their stars deeper into the host; ii) as a consequence, ASHs feature a negative gradient between radius and the local mass-weighed virial satellite-to-host mass ratio; iii) in $L_*$ galaxies, most ASHs feature a density profile that steepens towards sharper logarithmic slopes at increasing radii, though with significant halo-to-halo scatter; iv) the ASHs with the largest total ex-situ mass are such because of the chance accretion of a small number of massive satellites (rather than of a large number of low-mass ones).
△ Less
Submitted 30 July, 2016;
originally announced August 2016.
-
Gaps in globular cluster streams: giant molecular clouds can cause them too
Authors:
Nicola C. Amorisco,
Facundo A. Gòmez,
Simona Vegetti,
Simon D. M. White
Abstract:
As a result of their internal dynamical coherence, thin stellar streams formed by disrupting globular clusters (GCs) can act as detectors of dark matter (DM) substructure in the Galactic halo. Perturbations induced by close flybys amplify into detectable density gaps, providing a probe both of the abundance and of the masses of DM subhaloes. Here, we use N-body simulations to show that the Galacti…
▽ More
As a result of their internal dynamical coherence, thin stellar streams formed by disrupting globular clusters (GCs) can act as detectors of dark matter (DM) substructure in the Galactic halo. Perturbations induced by close flybys amplify into detectable density gaps, providing a probe both of the abundance and of the masses of DM subhaloes. Here, we use N-body simulations to show that the Galactic population of giant molecular clouds (GMCs) can also produce gaps (and clumps) in GC streams, and so may confuse the detection of DM subhaloes. We explore the cases of streams analogous to the observed Palomar 5 and GD1 systems, quantifying the expected incidence of structure caused by GMC perturbations. Deep observations should detect such disturbances regardless of the substructure content of the Milky Way's halo. Detailed modelling will be needed to demonstrate that any detected gaps or clumps were produced by DM subhaloes rather than by molecular clouds.
△ Less
Submitted 26 July, 2016; v1 submitted 8 June, 2016;
originally announced June 2016.
-
Ultra-diffuse galaxies: the high-spin tail of the abundant dwarf galaxy population
Authors:
N. C. Amorisco,
A. Loeb
Abstract:
Recent observations have revealed the existence of an abundant population of faint, low surface brightness (SB) galaxies, which appear to be numerous and ubiquitous in nearby galaxy clusters, including the Virgo, Coma and Fornax clusters. With median stellar masses of dwarf galaxies, these ultra-diffuse galaxies (UDGs) have unexpectedly large sizes, corresponding to a mean SB of…
▽ More
Recent observations have revealed the existence of an abundant population of faint, low surface brightness (SB) galaxies, which appear to be numerous and ubiquitous in nearby galaxy clusters, including the Virgo, Coma and Fornax clusters. With median stellar masses of dwarf galaxies, these ultra-diffuse galaxies (UDGs) have unexpectedly large sizes, corresponding to a mean SB of $24\lesssim\langleμ_e\rangle_r\ {\rm mag}^{-1} {\rm arcsec}^2\lesssim27$ within the effective radius. We show that the UDG population represents the tail of galaxies formed in dwarf-sized haloes with higher-than-average angular momentum. By adopting the standard model of disk formation -- in which the size of galaxies is set by the spin of the halo -- we recover both the abundance of UDGs as a function of the host cluster mass and the distribution of sizes within the UDG population. According to this model, UDGs are not failed $L_*$ galaxies, but genuine dwarfs, and their low SB is not uniquely connected to the harsh cluster environment. We therefore expect a correspondingly abundant population of UDGs in the field, with likely different morphologies and colours.
△ Less
Submitted 30 March, 2016; v1 submitted 1 March, 2016;
originally announced March 2016.
-
Chemodynamic subpopulations of the Carina dwarf galaxy
Authors:
G. Kordopatis,
N. C. Amorisco,
N. W. Evans,
G. Gilmore,
S. E. Koposov
Abstract:
We study the chemodynamical properties of the Carina dwarf spheroidal by combining an intermediate spectroscopic resolution dataset of more than 900 red giant and red clump stars, with high-precision photometry to derive the atmospheric parameters, metallicities and age estimates for our targets. Within the red giant branch population, we find evidence for the presence of three distinct stellar su…
▽ More
We study the chemodynamical properties of the Carina dwarf spheroidal by combining an intermediate spectroscopic resolution dataset of more than 900 red giant and red clump stars, with high-precision photometry to derive the atmospheric parameters, metallicities and age estimates for our targets. Within the red giant branch population, we find evidence for the presence of three distinct stellar sub-populations with different metallicities, spatial distributions, kinematics and ages. As in the Fornax and Sculptor dwarf spheroidals, the subpopulation with the lowest average metallicity is more extended and kinematically hotter than all other populations. However, we identify an inversion in the parallel ordering of metallicity, kinematics and characteristic length scale in the two most metal rich subpopulations, which therefore do not contribute to a global negative chemical gradient. Contrary to common trends in the chemical properties with radius, the metal richest population is more extended and mildly kinematically hotter than the main component of intermediate metallicity. More investigations are required to ascertain the nature of this inversion, but we comment on the mechanisms that might have caused it.
△ Less
Submitted 7 January, 2016;
originally announced January 2016.
-
Contributions to the accreted stellar halo: an atlas of stellar deposition
Authors:
N. C. Amorisco
Abstract:
The accreted component of stellar halos is composed of the contributions of several satellites, falling onto their host with their different masses, at different times, on different orbits. This work uses a suite of idealised, collisionless N-body simulations of minor mergers and a particle tagging technique to understand how these different ingredients shape each contribution to the accreted halo…
▽ More
The accreted component of stellar halos is composed of the contributions of several satellites, falling onto their host with their different masses, at different times, on different orbits. This work uses a suite of idealised, collisionless N-body simulations of minor mergers and a particle tagging technique to understand how these different ingredients shape each contribution to the accreted halo, in both density and kinematics. I find that more massive satellites deposit their stars deeper into the gravitational potential of the host, with a clear segregation enforced by dynamical friction. Earlier accretion events contribute more to the inner regions of the halo; more concentrated subhaloes sink deeper through increased dynamical friction. The orbital circularity of the progenitor at infall is only important for low-mass satellites: dynamical friction efficiently radialises the most massive minor mergers erasing the imprint of the infall orbit for satellite-to-host virial mass ratios $\gtrsim1/20$. The kinematics of the stars contributed by each satellite is also ordered with satellite mass: low-mass satellites contribute fast-moving populations, in both ordered rotation and radial velocity dispersion. In turn, contributions by massive satellites have lower velocity dispersion and lose their angular momentum to dynamical friction, resulting in a strong radial anisotropy.
△ Less
Submitted 2 September, 2016; v1 submitted 27 November, 2015;
originally announced November 2015.
-
A dwarf galaxy's transformation and a massive galaxy's edge: autopsy of kill and killer in NGC 1097
Authors:
N. C. Amorisco,
D. Martinez-Delgado,
J. Schedler
Abstract:
(abridged) We present a dynamical analysis of the extended stellar stream encircling NGC 1097. Within a statistical framework, we model its surface brightness using mock streams as in Amorisco (2015) and deep imaging data from the CHART32 telescope (Stellar Tidal Stream Survey). We reconstruct the post-infall evolution of the progenitor, which has experienced 3 pericentric passages and lost more t…
▽ More
(abridged) We present a dynamical analysis of the extended stellar stream encircling NGC 1097. Within a statistical framework, we model its surface brightness using mock streams as in Amorisco (2015) and deep imaging data from the CHART32 telescope (Stellar Tidal Stream Survey). We reconstruct the post-infall evolution of the progenitor, which has experienced 3 pericentric passages and lost more than 2 orders of magnitude in mass. At infall, $5.4\pm0.6$ Gyr ago, the progenitor was a disky dwarf with mass of $\log_{10}[m(<3.4\pm1 {\rm kpc})/ M_\odot]=10.35\pm0.25$. We illustrate how the 90$^\circ$ turn in the stream, identifying the `dog leg', is the signature of the progenitor's prograde rotation. Today, the remnant is a nucleated dwarf, with a LOS velocity of $v_{\rm p, los}^{\rm obs}=-30\pm 30$ kms$^{-1}$, and a luminosity of $3.3\times 10^7 L_{V,\odot}$ (Galianni et al. 2010). Our independent analysis predicts $v_{\rm p, los}=-51^{-17}_{+14}$ kms$^{-1}$, and measures $\log_{10}(m/ M_\odot)=7.4^{+0.6}_{-0.8}$, so that the compact nucleus is soon becoming a low-luminosity UCD. We find that NGC 1097 has a mass of $M_{200}=1.8^{+0.5}_{-0.4} \times 10^{12}\; M_{\odot}$, and its concentration $c_{200}=6.7^{+2.4}_{-1.3}$ is in agreement with LCDM. The stream is described almost down to the noise in a spherical host potential, we find this would not be possible if the halo was substantially triaxial at large radii. Its morphology shows that the slope of the total density profile bends from an inner $γ(r_{\rm peri})=1.5\pm0.15$. The progenitor's orbit reaches $r_{\rm apo}=150\pm 15$ kpc, more than a half of the virial radius of the host, so that, for the first time on an individual extragalactic halo, we measure the outer density slope, $γ(0.6r_{200,c})=3.9\pm0.5$. This demonstrates the promise of the newborn field of detailed, statistical modelling of extragalactic tidal streams.
△ Less
Submitted 14 April, 2015;
originally announced April 2015.
-
On feathers, bifurcations and shells: the dynamics of tidal streams across the mass scale
Authors:
N. C. Amorisco
Abstract:
I present an organic description of the regimes of collisionless tidal streams and define the orderings between the physical quantities that shape their morphology. Three fundamental dichotomies are identified in the form of dimensionless inequalities. These govern i) the speed of the stream's growth, ii) its internal coherence, iii) its thickness or opening angles. The mechanisms that regulate su…
▽ More
I present an organic description of the regimes of collisionless tidal streams and define the orderings between the physical quantities that shape their morphology. Three fundamental dichotomies are identified in the form of dimensionless inequalities. These govern i) the speed of the stream's growth, ii) its internal coherence, iii) its thickness or opening angles. The mechanisms that regulate such main properties are analysed. The slope of the host's density profile influences the speed of the stream's growth, in both length and width, as steeper profiles enhance differential streaming. Internal coherence is the requirement for the appearance of substructure in tidal debris, and I concentrate on the `feathering' typical of GC streams. Overdensities are associated with minima in the relative streaming velocity of the stream members. For streams with high circularity, these are caused by the epicyclic oscillations of stars; however, for highly non-circular progenitor's orbits, substructure is caused by the oscillating differences in energy and actions with which material is shed at different orbital phases of the progenitor. This modulation results in different streaming speeds: the streakline of material shed between two successive apocentric passages is folded along its length, pulled at its centre by the faster streaming of particles released near pericenter, which are therefore more widely scattered. When the stream is coherent enough, this mechanism is potentially capable of generating a bimodal profile in the density distributions of the longer wraps of more massive progenitors, which I dub `bifurcations'. The conditions for internal coherence are explored and I comment on the cases of Palomar 5, Willman 1, the Anticenter and Sagittarius' streams. Analytical methods are accompanied by numerical experiments, performed using a purposely built generative model, also presented here.
△ Less
Submitted 16 March, 2015; v1 submitted 1 October, 2014;
originally announced October 2014.
-
The remnant of a merger between two dwarf galaxies in Andromeda II
Authors:
N. C. Amorisco,
N. W. Evans,
G. van de Ven
Abstract:
Driven by gravity, massive structures like galaxies and clusters of galaxies are believed to grow continuously through hierarchical merging and accretion of smaller systems. Observational evidence of accretion events is provided by the coherent stellar streams crossing the outer haloes of massive galaxies, such as the Milky Way or Andromeda. At similar mass-scales, around $10^{11}$ solar masses in…
▽ More
Driven by gravity, massive structures like galaxies and clusters of galaxies are believed to grow continuously through hierarchical merging and accretion of smaller systems. Observational evidence of accretion events is provided by the coherent stellar streams crossing the outer haloes of massive galaxies, such as the Milky Way or Andromeda. At similar mass-scales, around $10^{11}$ solar masses in stars, further observational evidence of merging activity is also ample. Mergers of lower-mass galaxies are expected within the hierarchical process of galaxy formation, but have hitherto not been seen for galaxies with less than about $10^9$ solar masses in stars. Here, we report the kinematic detection of a stellar stream in one of the satellite galaxies of Andromeda, the dwarf spheroidal galaxy Andromeda II, which has a mass of only $10^7$ solar masses in stars. The properties of the stream show that we are observing the remnant of a merger between two dwarf galaxies. This had a dramatic influence on the dynamics of the remnant, which is now rotating around its projected major axis. The stellar stream in Andromeda II illustrates the scale-free character of the formation of galaxies, down to the lowest galactic mass scales.
△ Less
Submitted 20 February, 2014;
originally announced February 2014.
-
Dark matter cores in the Fornax and Sculptor dwarf galaxies: joining halo assembly and detailed star formation histories
Authors:
Nicola C. Amorisco,
Jesus Zavala,
Thomas J. L. de Boer
Abstract:
We combine the detailed Star Formation Histories of the Fornax and Sculptor dwarf Spheroidals with the Mass Assembly History of their dark matter (DM) halo progenitors to estimate if the energy deposited by Supernova type II (SNeII) is sufficient to create a substantial DM core. Assuming the efficiency of energy injection of the SNeII into DM particles is $ε_{\rm gc}=0.05$, we find that a single e…
▽ More
We combine the detailed Star Formation Histories of the Fornax and Sculptor dwarf Spheroidals with the Mass Assembly History of their dark matter (DM) halo progenitors to estimate if the energy deposited by Supernova type II (SNeII) is sufficient to create a substantial DM core. Assuming the efficiency of energy injection of the SNeII into DM particles is $ε_{\rm gc}=0.05$, we find that a single early episode, $z \gtrsim z_{\rm infall}$, that combines the energy of all SNeII due to explode over 0.5 Gyr, is sufficient to create a core of several hundred parsecs in both Sculptor and Fornax. Therefore, our results suggest that it is energetically plausible to form cores in Cold Dark Matter (CDM) halos via early episodic gas outflows triggered by SNeII. Furthermore, based on CDM merger rates and phase-space density considerations, we argue that the probability of a subsequent complete regeneration of the cusp is small for a substantial fraction of dwarf-size haloes.
△ Less
Submitted 28 January, 2014; v1 submitted 23 September, 2013;
originally announced September 2013.
-
The core size of the Fornax dwarf Spheroidal
Authors:
N. C. Amorisco,
A. Agnello,
N. W. Evans
Abstract:
We exploit the detection of three distinct stellar subpopulations in the red giant branch of the Fornax dwarf Spheroidal to probe its density distribution. This allows us to resolve directly the evolution with radius of the dark matter mass profile. We find that a cored dark matter halo provides a perfect fit to the data, being consistent with all three stellar populations well within 1-sigma, and…
▽ More
We exploit the detection of three distinct stellar subpopulations in the red giant branch of the Fornax dwarf Spheroidal to probe its density distribution. This allows us to resolve directly the evolution with radius of the dark matter mass profile. We find that a cored dark matter halo provides a perfect fit to the data, being consistent with all three stellar populations well within 1-sigma, and for the first time we are able to put constraints on the core size of such a halo. With respect to previous work, we do not strengthen the statistical exclusion of a dark matter cusp in Fornax, but we find that Navarro-Frenk-White haloes would be required to have unrealistically large scale radii in order to be compatible with the data, hence low values of the concentration parameter. We are then forced to conclude that the Fornax dwarf Spheroidal sits within a dark matter halo having a constant density core, with a core size of between 0.6 and 1.8 kpc.
△ Less
Submitted 11 October, 2012;
originally announced October 2012.
-
A Troublesome Past: Chemodynamics of the Fornax dwarf spheroidal
Authors:
N. C. Amorisco,
N. W. Evans
Abstract:
We present compelling evidence for the complexity of the Fornax dwarf spheroidal. By disentangling three different stellar subpopulations in its red giant branch, we are able to study in detail the dependence between kinematics and metallicity. A well-defined ordering in velocity dispersion, spatial concentration, and metallicity is evident in the subpopulations. We also present evidence for a sig…
▽ More
We present compelling evidence for the complexity of the Fornax dwarf spheroidal. By disentangling three different stellar subpopulations in its red giant branch, we are able to study in detail the dependence between kinematics and metallicity. A well-defined ordering in velocity dispersion, spatial concentration, and metallicity is evident in the subpopulations. We also present evidence for a significant misalignment between the angular momentum vectors of the old and intermediate-age populations. According to the HST measurement of Fornax's proper motion, this corresponds to counter-rotation. These ingredients are used to construct a novel evolutionary history of the Fornax dwarf spheroidal, characterized as a late merger of a bound pair.
△ Less
Submitted 3 August, 2012; v1 submitted 28 June, 2012;
originally announced June 2012.
-
Line Profiles from Discrete Kinematic Data
Authors:
N. C. Amorisco,
N. W. Evans
Abstract:
We develop a method to extract the shape information of line profiles from discrete kinematic data. The Gauss-Hermite expansion, which is widely used to describe the line of sight velocity distributions extracted from absorption spectra of elliptical galaxies, is not readily applicable to samples of discrete stellar velocity measurements, accompanied by individual measurement errors and probabilit…
▽ More
We develop a method to extract the shape information of line profiles from discrete kinematic data. The Gauss-Hermite expansion, which is widely used to describe the line of sight velocity distributions extracted from absorption spectra of elliptical galaxies, is not readily applicable to samples of discrete stellar velocity measurements, accompanied by individual measurement errors and probabilities of membership. We introduce two parameter families of probability distributions describing symmetric and asymmetric distortions of the line profiles from Gaussianity. These are used as the basis of a maximum likelihood estimator to quantify the shape of the line profiles. Tests show that the method outperforms a Gauss-Hermite expansion for discrete data, with a lower limit for the relative gain of approx 2 for sample sizes N approx 800. To ensure that our methods can give reliable descriptions of the shape, we develop an efficient test to assess the statistical quality of the obtained fit. As an application, we turn our attention to the discrete velocity datasets of the dwarf spheroidals of the Milky Way. In Sculptor, Carina and Sextans the symmetric deviations are consistent with velocity distributions more peaked than Gaussian. In Fornax, instead, there is an evolution in the symmetric deviations of the line profile from a peakier to more flat-topped distribution on moving outwards. These results suggest a radially biased orbital structure for the outer parts of Sculptor, Carina and Sextans. On the other hand, tangential anisotropy is favoured in Fornax. This is all consistent with a picture in which Fornax may have had a different evolutionary history to Sculptor, Carina and Sextans.
△ Less
Submitted 28 June, 2012; v1 submitted 23 April, 2012;
originally announced April 2012.
-
Dark Matter Cores and Cusps: The Case of Multiple Stellar Populations in Dwarf Spheroidals
Authors:
N. C. Amorisco,
N. W. Evans
Abstract:
A number of dwarf spheroidal (dSph) galaxies are known to contain a more extended, metal-poor population with a flattish velocity dispersion profile, and a more concentrated, metal-rich population with a velocity dispersion declining with radius. The two populations can be modelled with Michie-King distribution functions (DFs) in the isothermal and in the sharply truncated limits, respectively. We…
▽ More
A number of dwarf spheroidal (dSph) galaxies are known to contain a more extended, metal-poor population with a flattish velocity dispersion profile, and a more concentrated, metal-rich population with a velocity dispersion declining with radius. The two populations can be modelled with Michie-King distribution functions (DFs) in the isothermal and in the sharply truncated limits, respectively. We argue that the truncation of the metal-rich population can be traced back to the spatial distribution of the star forming gas. Suppose D is the exponent of the first non-constant term in the Taylor expansion of the total potential at the center (D=1 for NFW halos, D=2 for cored halos). Then, we show that the ratio of the half-light radii of the populations raised to the D/2 power must be smaller than the ratio of the line-of-sight velocity dispersions. Specialising to the case of the Sculptor dSph, we develop a technique to fit simultaneously both populations with Michie-King DFs. This enables us to determine the mass profile of the Sculptor dSph with unprecedented accuracy in the radial range between 0.2 and 1.2 kpc. We show that cored halo models are preferred over cusped halo models, with a likelihood ratio test rejecting NFW models at any significance level higher than 0.05%. Even more worryingly, the best-fit NFW models require concentrations with c<20, which is not in the cosmologically preferred range for dwarf galaxies. We conclude that the kinematics of multiple populations in dSphs provides a substantial new challenge for theories of galaxy formation, with the weight of available evidence strongly against dark matter cusps at the centre.
△ Less
Submitted 21 September, 2011; v1 submitted 6 June, 2011;
originally announced June 2011.
-
Phase Space Models of the Dwarf Spheroidal Galaxies
Authors:
N. C. Amorisco,
N. W. Evans
Abstract:
This paper introduces new phase-space models of dwarf spheroidal galaxies (dSphs). The stellar component has an isotropic, lowered isothermal (or King) distribution function. A physical basis for the isotropization of stellar velocities is given by tidal stirring, whilst the isothermality of the distribution function guarantees the observed flatness of the velocity dispersion profile in the inner…
▽ More
This paper introduces new phase-space models of dwarf spheroidal galaxies (dSphs). The stellar component has an isotropic, lowered isothermal (or King) distribution function. A physical basis for the isotropization of stellar velocities is given by tidal stirring, whilst the isothermality of the distribution function guarantees the observed flatness of the velocity dispersion profile in the inner parts. Our models reproduce the data on the half-light radius and line of sight central velocity dispersion of the dSphs. We show that different dark halo profiles -- whether cored or cusped -- lead to very similar mass estimates within one particular radius, namely 1.7 half-light radii. Deviations between mass measures due to different density profiles are substantially smaller than the uncertainties propagated by the observational errors. We produce a mass measure for each of the Milky Way dSphs and find that the two most massive are the most luminous, namely Sagittarius (~ 2.8 x 10^8 solar masses) and Fornax (~ 1.3 x 10^8 solar masses). The least massive of the Milky Way satellites are Willman 1 (~ 4 x 10^5 solar masses) and Segue 1 (~ 6 x 10^5 solar masses).
△ Less
Submitted 26 January, 2011; v1 submitted 9 September, 2010;
originally announced September 2010.
-
Self-consistent nonspherical isothermal halos embedding zero-thickness disks
Authors:
N. C. Amorisco,
G. Bertin
Abstract:
Disk-halo decompositions of galaxy rotation curves are generally performed in a parametric way. We construct self-consistent models of nonspherical isothermal halos embedding a zero-thickness disk, by assuming that the halo distribution function is a Maxwellian. The method developed here can be used to study other physically-based choices for the halo distribution function and the case of a disk a…
▽ More
Disk-halo decompositions of galaxy rotation curves are generally performed in a parametric way. We construct self-consistent models of nonspherical isothermal halos embedding a zero-thickness disk, by assuming that the halo distribution function is a Maxwellian. The method developed here can be used to study other physically-based choices for the halo distribution function and the case of a disk accompanied by a bulge. In a preliminary investigation we note the existence of a fine tuning between the scalelengths R_Ω and h, respectively characterizing the rise of the rotation curve and the luminosity profile of the disk, which surprisingly applies to both high surface brightness and low surface brightness galaxies. This empirical correlation identifies a much stronger conspiracy than the one required by the smoothness and flatness of the rotation curve (disk-halo conspiracy). The self-consistent models are characterized by smooth and flat rotation curves for very different disk-to-halo mass ratios, hence suggesting that conspiracy is not as dramatic as often imagined. For a typical rotation curve, with asymptotically flat rotation curve at V_{\infty} (the precise value of which can also be treated as a free parameter), and a typical density profile of the disk, self-consistent models are characterized by two dimensionless parameters, which correspond to the dimensional scales (the disk mass-to-light ratio M/L and the halo central density) of standard disk-halo decompositions. We show that if the rotation curve is decomposed by means of our self-consistent models, the disk-halo degeneracy is removed and typical rotation curves are fitted by models that are below the maximum-disk prescription. Similar results are obtained from a study of NGC 3198. Finally, we quantify the flattening of the spheroidal halo, which is significant, especially on the scale of the visible disk.
△ Less
Submitted 18 May, 2010;
originally announced May 2010.
-
Prominent spiral arms in the gaseous outer galaxy disks
Authors:
G. Bertin,
N. C. Amorisco
Abstract:
Context: Several spiral galaxies, as beautifully exhibited by the case of NGC 6946, display a prominent large-scale spiral structure in their gaseous outer disk. Such structure is often thought to pose a dynamical puzzle, because grand-design spiral structure is traditionally interpreted as the result of density waves carried mostly in the stellar disk. Aims. Here we argue that the outer spiral…
▽ More
Context: Several spiral galaxies, as beautifully exhibited by the case of NGC 6946, display a prominent large-scale spiral structure in their gaseous outer disk. Such structure is often thought to pose a dynamical puzzle, because grand-design spiral structure is traditionally interpreted as the result of density waves carried mostly in the stellar disk. Aims. Here we argue that the outer spiral arms in the cold gas outside the bright optical disk actually have a natural interpretation as the manifestation of the mechanism that excites grand-design spiral structure in the main, star-dominated body of the disk: the excitation is driven by angular momentum transport to the outer regions, through trailing density waves outside the corotation circle that can penetrate beyond the Outer Lindblad Resonance in the gaseous component of the disk.
Methods: Because of conservation of the density wave action, these outgoing waves are likely to become more prominent in the outer disk and eventually reach non-linear amplitudes. To calculate the desired amplitude profiles, we make use of the theory of dispersive waves.
Results: If the conditions beyond the optical radius allow for an approximate treatment in terms of a linear theory, we show that fitting the observed amplitude profiles leads to a quantitative test on the density of the disk material and thus on the dark matter distribution in the outer parts of the galaxy.
Conclusions: This study is thus of interest to the general problem of the disk-halo decomposition of rotation curves.
△ Less
Submitted 16 December, 2009;
originally announced December 2009.