Satellite Stars Research Articles

Metallicity distributions of halo stars: do they trace the Galactic accretion history?

The standard cosmological scenario predicts a hierarchical formation for galaxies. Many substructures have been found in the Galactic halo, usually identified as clumps in kinematic spaces, like the energy-angular momentum space ($E-L_z$), under the hypothesis that these quantities should be conserved during the interaction. If these clumps also feature different chemical abundances, such as the metallicity distribution function (MDF), these two arguments together (different kinematic and chemical properties) are often used to motivate their association with distinct and independent merger debris. The aim of this study is to explore to what extent we can couple kinematic characteristics and metallicities of stars in the Galactic halo to reconstruct the accretion history of the Milky Way (MW). In particular, we want to understand whether different clumps in the $E-L_z$ space with different MDFs should be associated with distinct merger debris. We analysed dissipationless, self-consistent, high-resolution N-body simulations of a MW-type galaxy accreting a satellite with a mass ratio of 1:10, with different orbital parameters and different metallicity gradients, which were assigned a posteriori. We confirm that accreted stars from a sim 1:10 mass ratio merger event redistribute in a wide range of $E$ and $L_z$, due to the dynamical friction process, and are thus not associated with a single region. Because satellite stars with different metallicities can be deposited in different regions of the $E-L_z$ space (on average the more metal-rich ones end up more gravitationally bound to the MW), this implies that a single accretion of sim 1:10 can manifest with different MDFs, in different regions of the $E-L_z$ space. Groups of stars with different $E$, $L_z$, and metallicities may be interpreted as originating from different satellite galaxies, but our analysis shows that these interpretations are not physically motivated. In fact, as we show, the coupling of kinematic information with MDFs to reconstruct the accretion history of the MW can bias the reconstructed merger tree towards increasing the number of past accretions and decreasing the masses of the progenitor galaxies.

VPNEP: Detailed characterization of TESS targets around the Northern Ecliptic Pole

Context.We embarked on a high-resolution optical spectroscopic survey of bright Transiting Exoplanet Survey Satellite (TESS) stars around the Northern Ecliptic Pole (NEP), dubbed the Vatican-Potsdam-NEP (VPNEP) survey.Aims.Our NEP coverage comprises ≈770 square degrees with 1067 stars, of which 352 are bona fide dwarf stars and 715 are giant stars, all cooler than spectral type F0 and brighter thanV= 8m.5. Our aim is to characterize these stars for the benefit of future studies in the community.Methods.We analyzed the spectra via comparisons with synthetic spectra. Particular line profiles were analyzed by means of eigenprofiles, equivalent widths, and relative emission-line fluxes (when applicable).Results.TwoR= 200 000 spectra were obtained for each of the dwarf stars with the Vatican Advanced Technology Telescope (VATT) and the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI), with typically threeR= 55 000 spectra obtained for the giant stars with STELLA and the STELLA Echelle Spectrograph (SES). Combined withV-band magnitudes,GaiaEDR3 parallaxes, and isochrones from the Padova and Trieste Stellar Evolutionary Code, the spectra can be used to obtain radial velocities, effective temperatures, gravities, rotational and turbulence broadenings, stellar masses and ages, and abundances for 27 chemical elements, as well as isotope ratios for lithium and carbon, line bisector spans, convective blue-shifts (when feasible), and levels of magnetic activity from Hα, Hβ, and the Ca IIinfrared triplet. In this initial paper, we discuss our analysis tools and biases, presenting our first results from a pilot sub-sample of 54 stars (27 bona-fide dwarf stars observed with VATT+PEPSI and 27 bona-fide giant stars observed with STELLA+SES) and making all reduced spectra available to the community. We carried out a follow-up error analysis, including systematic biases and standard deviations based on a joint target sample for both facilities, as well as a comparison with external data sources.

Machine learning for galactic archaeology: a chemistry-based neural network method for identification of accreted disc stars

ABSTRACT We develop a method (‘Galactic Archaeology Neural Network’, gann) based on neural network models (NNMs) to identify accreted stars in galactic discs by only their chemical fingerprint and age, using a suite of simulated galaxies from the Auriga Project. We train the network on the target galaxy’s own local environment defined by the stellar halo and the surviving satellites. We demonstrate that this approach allows the detection of accreted stars that are spatially mixed into the disc. Two performance measures are defined – recovery fraction of accreted stars, frecov and the probability that a star with a positive (accreted) classification is a true-positive result, P(TP). As the NNM output is akin to an assigned probability (Pa), we are able to determine positivity based on flexible threshold values that can be adjusted easily to refine the selection of presumed-accreted stars. We find that gann identifies accreted disc stars within simulated galaxies, with high frecov and/or high P(TP). We also find that stars in Gaia–Enceladus–Sausage (GES) mass systems are over 50 per cent recovered by our NNMs in the majority (18/24) of cases. Additionally, nearly every individual source of accreted stars is detected at 10 per cent or more of its peak stellar mass in the disc. We also demonstrate that a conglomerated NNM, trained on the halo and satellite stars from all of the Auriga galaxies provides the most consistent results, and could prove to be an intriguing future approach as our observational capabilities expand.

The effect of tides on the Fornax dwarf spheroidal galaxy

Estimates of the mass distribution and dark-matter (DM) content of dwarf spheroidal galaxies (dSphs) are usually derived under the assumption that the effect of the tidal field of the host galaxy is negligible over the radial extent probed by kinematic data-sets. We assess the implications of this assumption in the specific case of the Fornax dSph by means of N-body simulations of a satellite orbiting around the Milky Way. We consider observationally-motivated orbits and we tailor the initial distributions of the satellite's stars and DM to match, at the end of the simulations, the observed structure and kinematics of Fornax. In all our simulations the present-day observable properties of Fornax are not significantly influenced by tidal effects. The DM component is altered by the interaction with the Galactic field (up to 20% of the DM mass within 1.6 kpc is lost), but the structure and kinematics of the stellar component are only mildly affected even in the more eccentric orbit (more than 99% of the stellar particles remain bound to the dwarf). In the simulations that successfully reproduce Fornax's observables, the dark-to-luminous mass ratio within 1.6 kpc is in the range 5-6, and up to 16-18 if measured within 3 kpc.

Planetary Nebulae and their parent stellar populations. Tracing the mass assembly of M87 and Intracluster light in the Virgo cluster core

AbstractThe diffuse extended outer regions of galaxies are hard to study because they are faint, with typical surface brightness of 1% of the dark night sky. We can tackle this problem by using resolved star tracers which remain visible at large distances from the galaxy centers. This article describes the use of Planetary Nebulae as tracers and the calibration of their properties as indicators of the star formation history, mean age and metallicity of the parent stars in the Milky Way and Local Group galaxies. We then report on the results from a deep, extended, planetary nebulae survey in a 0.5 deg2region centered on the brightest cluster galaxy NGC 4486 (M87) in the Virgo cluster core, carried out with SuprimeCam@Subaru and FLAMES-GIRAFFE@VLT. Two planetary nebulae populations are identified out to 150 kpc distance from the center of M87. One population is associated with the M87 halo and the second one with the intracluster light in the Virgo cluster core. They have different line-of-sight velocity and spatial distributions, as well as different planetary nebulae specific frequencies and luminosity functions. The intracluster planetary nebulae in the surveyed region correspond to a luminosity of four times the luminosity of the Large Magellanic Cloud. The M87 halo planetary nebulae trace an older, more metal-rich, parent stellar population. A substructure detected in the projected phase-space of the line-of-sight velocity vs. major axis distance for the M87 halo planetary nebulae provides evidence for the recent accretion event of a satellite galaxy with luminosity twice that of M33. The satellite stars were tidally stripped about 1 Gyr ago, and reached apocenter at a major axis distance of 60–90 kpc from the center of M87. The M87 halo is still growing significantly at the distances where the substructure is detected.

The dual nature of the Milky Way stellar halo

AbstractThe theory of the Milky Way formation, in the framework of the ΛCDM model, predicts galactic stellar halos to be built from multiple accretion events starting from the first structure to collapse in the Universe.Evidences in the past few decades have indicated that the Galactic halo consists of two overlapping structural components, an inner and an outer halo. We provide a set of numerical N-body simulations aimed to study the formation of the outer Milky Way (MW) stellar halo through accretion events between a (bulgeless) MW-like system and a satellite galaxy. After these minor mergers take place, in several orbital configurations, we analyze the signal left by satellite stars in the rotation velocity distribution. The aim is to explore the orbital conditions of the mergers where a signal of retrograde rotation in the outer part of the halo can be obtained, in order to give a possible explanation of the observed rotational properties of the MW stellar halo.Our results show that the dynamical friction has a fundamental role in assembling the final velocity distributions originated by different orbits and that retrograde satellites moving on low inclination orbits deposit more stars in the outer halo regions and therefore can produce the counter-rotating behavior observed in the outer MW halo.

ASSEMBLY OF THE OUTER GALACTIC STELLAR HALO IN THE HIERARCHICAL MODEL

We provide a set of numerical N-body simulations for studying the formation of the outer Milky Ways's stellar halo through accretion events. After simulating minor mergers of prograde and retrograde orbiting satellite halo with a Dark Matter main halo, we analyze the signal left by satellite stars in the rotation velocity distribution. The aim is to explore the orbital conditions where a retrograde signal in the outer part of the halo can be obtained, in order to give a possible explanation of the observed rotational properties of the Milky Way stellar halo. Our results show that, for satellites more massive than $\sim 1/40$ of the main halo, the dynamical friction has a fundamental role in assembling the final velocity distributions resulting from different orbits and that retrograde satellites moving on low inclination orbits deposit more stars in the outer halo regions end therefore can produce the counter-rotating behavior observed in the outer Milky Way halo.

The dynamics of satellite disruption in cold dark matter haloes

We investigate the physical mechanisms of tidal heating and satellite disruption in cold dark matter host haloes using N-body simulations based on cosmological initial conditions. We show the importance of resonant shocks and resonant torques with the host halo to satellite heating. A resonant shock (torque) couples the radial (tangential) motion of a satellite in its orbit to its phase space. For a satellite on a circular orbit, an ILR-like resonance dominates the heating and this heating results in continuous satellite mass loss. We estimate the requirements for simulations to achieve these dynamics using perturbation theory. Both resonant shocks and resonant torques affect satellites on eccentric orbits. We demonstrate that satellite mass loss is an outside-in process in energy space; a satellite's stars and gas are thus protected by their own halo against tidal stripping. We simulate the evolution of a halo similar to the Large Magellanic Cloud (LMC) in our Galactic dark matter halo and conclude that the LMC stars have not yet been stripped. Finally, we present a simple algorithm for estimating the evolution of satellite mass that includes both shock heating and resonant torques.

Disassembling the Galaxy with angle-action coordinates

Angle-action coordinates are used to study the relic of an N-body simulation of a self-gravitating satellite galaxy that was released on a short-period orbit within the disc of the Galaxy. Satellite stars that lie within 1.5 kpc of the Sun are confined to a grid of patches in action space. As the relic phase mixes for longer, the patches become smaller and more numerous. These patches can be seen even when the angle-action coordinates of an erroneous Galactic potential are used, but using the wrong potential displaces them. Diagnostic quantities constructed from the angle coordinates both allow the true potential to be identified, and the relic to be dated. Hence, when the full phase space coordinates of large numbers of solar-neighbourhood stars are known, it should be possible to identify members of particular relics from the distribution of stars in an approximate action space. This would then open up the possibility of determining the time since the relic was disrupted and gaining better knowledge of the Galactic potential. The availability of angle-action coordinates for arbitrary potentials is the key to these developments. The paper includes a brief introduction to the torus technique used to generate them.

The Hierarchical Build‐Up of Massive Galaxies and the Intracluster Light sincez= 1

We use a set of simulation-based models for the dissipationless evolution of galaxies since z = 1 to constrain the fate of accreted satellites embedded in dark matter subhalos. These models assign stellar mass to dark matter halos at z = 1 by relating the observed galaxy stellar mass function (GSMF) to the halo + subhalo mass function monotonically. The evolution of the stellar mass content is then followed using halo merger trees extracted from N-body simulations. Our models are differentiated only in the fate assigned to satellite galaxies once subhalos, within which satellites are embedded, disrupt. These models are confronted with the observed evolution in the massive end of the GSMF, the z ~ 0 brightest cluster galaxy (BCG) to cluster mass relation, and the combined BCG and intracluster light (ICL) luminosity distribution—all observables expected to evolve approximately dissipationlessly since z = 1. The combined observational constraints favor a model in which the vast majority (≳80%) of satellite stars from disrupted subhalos go into the ICL (operationally defined here as light below a surface brightness cut of μi ≈ 23 mag arcsec-2). Conversely, models that leave behind a significant population of satellite galaxies once the subhalo has disrupted are strongly disfavored, as are models that put a significant fraction of satellite stars into the BCG. Our results show that observations of the ICL provide useful and unique constraints on models of galaxy merging and the dissipationless evolution of galaxies in groups and clusters.

Israeli satellite stars

Israeli satellite stars