• ### The Sausage Globular Clusters(1805.00453)

May 1, 2018 astro-ph.CO, astro-ph.GA
We analyse the structure of the Milky Way's globular cluster (GC) system in action space using the Gaia Data Release 2 catalogue. There is a characteristic energy $E_{\rm crit} = -1.6 \times 10^5$ km$^2$s$^{-2}$. The in situ objects, such as the bulge/disc clusters, all lie below this energy, the accreted objects (such as the young halo clusters) above. There are 12 old halo GCs that have $E > E_{\rm crit}$ and they all betray evidence of recent accretion. Eight of the high energy, old halo GCs are strongly clustered in azimuthal and vertical action, yet are strung out like beads on a chain at extreme radial action. They are the Sausage Globular Clusters, expelled from a massive dwarf galaxy ($\sim 5 \times 10^{10} M_\odot$) on a strongly radial orbit that merged with the Milky Way at a redshift $z\lesssim 3$. The Sausage GCs are strongly radially anisotropic ($\beta \sim 0.95$), move on orbits that are all highly eccentric ($e \gtrsim 0.85$) and form a track in the age-metallicity plane distinct from the bulk of the Milky Way GCs.
• We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.
• ### The Shards of $\omega$ Centauri(1804.07050)

April 23, 2018 astro-ph.GA
We use the SDSS-Gaia catalogue to search for material stripped off the atypical globular cluster $\omega$ Centauri. As actions are conserved under slow changes of the potential, this allows identification of groups of stars with a common accretion history. We devise a method to assess the significance of halo substructures based on their clustering in action space, using metallicity as a secondary check. This is validated against smooth models and numerically constructed stellar halos. We identify a number of new substructures in the SDSS-Gaia catalogue, including 7 high significance, high energy and retrograde ones. Using a simple model of the accretion of the progenitor of the $\omega$ Centauri, we argue for the possible association of up to 4 of our new substructures (labelled Rg1, Rg3, Rg4 and Rg6) with this event. This sets a minimum mass of $5 \times 10^8 M_\odot$ for the progenitor, so as to bring $\omega$ Centauri to its current location in action -- energy space. Our proposal can be tested by high resolution spectroscopy of the candidates to look for the unusual abundance patterns possessed by $\omega$ Centauri stars.
• ### Co-formation of the Galactic disc and the stellar halo(1802.03414)

Feb. 9, 2018 astro-ph.GA
Using a large sample of Main Sequence stars with 7-D measurements supplied by Gaia and SDSS, we study the kinematic properties of the local (within ~10 kpc from the Sun) stellar halo. We demonstrate that the halo's velocity ellipsoid evolves strongly with metallicity. At the low [Fe/H] end, the orbital anisotropy (the amount of motion in the radial direction compared to the tangential one) is mildly radial with 0.2<beta<0.4. However, for stars with [Fe/H]>-1.7 we measure extreme values of beta~0.9. Across the metallicity range considered, i.e. -3<[Fe/H]<-1, the stellar halo's spin is minimal, at the level of 20<v_theta (km/s) <30. Using a suite of cosmological zoom-in simulations of halo formation, we deduce that the observed acute anisotropy is inconsistent with the continuous accretion of dwarf satellites. Instead, we argue, the stellar debris in the inner halo were deposited in a major accretion event by a satellite with Mvir>10^10 Msun around the epoch of the Galactic disc formation, i.e. between 8 and 11 Gyr ago. The radical halo anisotropy is the result of the dramatic radialisation of the massive progenitor's orbit, amplified by the action of the growing disc.
• ### Unmixing the Galactic Halo with RR Lyrae tagging(1711.07485)

Nov. 20, 2017 astro-ph.GA
We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously.
• ### The first all-sky view of the Milky Way stellar halo with Gaia+2MASS RR Lyrae(1707.03833)

Nov. 2, 2017 astro-ph.GA
We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo ($\text{R}<28$ kpc). We study the morphology of the stellar halo by analysing the RR Lyrae distribution with parametric and non-parametric techniques. Taking advantage of the uniform all-sky coverage, we test halo models more sophisticated than usually considered in the literature, such as those with varying flattening, tilt and/or offset of the halo with respect to the Galactic disc. A consistent picture emerges: the inner halo is well reproduced by a smooth distribution of stars settled on triaxial ellipsoids. The minor axis is perpendicular to the Milky Ways disc, while the major axis is misaligned by $\sim20^{\circ}$ from the Galactic Y axis. The elongation along the major axis is mild ($\text{p}=1.27$), and the vertical flattening is shown to evolve from a squashed state with $\text{q}\approx0.57$ in the centre to a more spherical $\text{q}\approx0.75$ at the outer edge of our dataset. The density slope is well approximated by a single power-law with exponent $\alpha=-2.96$. Within the range probed, we see no significant evidence for a change of the radial density slope, out of the plane tilt or an offset of the halo with respect to the Galaxy's centre.
• ### The Gaia-ESO Survey: Matching Chemo-Dynamical Simulations to Observations of the Milky Way(1709.01523)

Sept. 5, 2017 astro-ph.GA
The typical methodology for comparing simulated galaxies with observational surveys is usually to apply a spatial selection to the simulation to mimic the region of interest covered by a comparable observational survey sample. In this work we compare this approach with a more sophisticated post-processing in which the observational uncertainties and selection effects (photometric, surface gravity and effective temperature) are taken into account. We compare a `solar neighbourhood analogue' region in a model Milky Way-like galaxy simulated with RAMSES-CH with fourth release Gaia-ESO survey data. We find that a simple spatial cut alone is insufficient and that observational uncertainties must be accounted for in the comparison. This is particularly true when the scale of uncertainty is large compared to the dynamic range of the data, e.g. in our comparison, the [Mg/Fe] distribution is affected much more than the more accurately determined [Fe/H] distribution. Despite clear differences in the underlying distributions of elemental abundances between simulation and observation, incorporating scatter to our simulation results to mimic observational uncertainty produces reasonable agreement. The quite complete nature of the Gaia-ESO survey means that the selection function has minimal impact on the distribution of observed age and metal abundances but this would become increasingly more important for surveys with narrower selection functions.
• Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.
• Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.
• ### Gaia data release 1, the photometric data(1612.02952)

Dec. 9, 2016 astro-ph.IM
Context. This paper presents an overview of the photometric data that are part of the first Gaia data release. Aims. The principles of the processing and the main characteristics of the Gaia photometric data are presented. Methods. The calibration strategy is outlined briefly and the main properties of the resulting photometry are presented. Results. Relations with other broadband photometric systems are provided. The overall precision for the Gaia photometry is shown to be at the milli-magnitude level and has a clear potential to improve further in future releases.
• ### The Gaia-ESO Survey: Lithium depletion in the Gamma Velorum cluster and inflated radii in low-mass pre-main-sequence stars(1609.07150)

Sept. 22, 2016 astro-ph.SR
We show that non-magnetic models for the evolution of pre-main-sequence (PMS) stars *cannot* simultaneously describe the colour-magnitude diagram (CMD) and the pattern of lithium depletion seen in the cluster of young, low-mass stars surrounding $\gamma^2$ Velorum. The age of 7.5+/-1 Myr inferred from the CMD is much younger than that implied by the strong Li depletion seen in the cluster M-dwarfs and the Li depletion occurs at much redder colours than predicted. The epoch at which a star of a given mass depletes its Li and the surface temperature of that star are both dependent on its radius. We demonstrate that if the low-mass stars have radii ~10 per cent larger at a given mass and age, then both the CMD and Li depletion pattern of the Gamma Vel cluster are explained at a common age of 18-21 Myr. This radius inflation could be produced by some combination of magnetic suppression of convection and extensive cool starspots. Models that incorporate radius inflation suggest that PMS stars similar to those in the Gamma Vel cluster, in the range 0.2<M/Msun<0.7, are at least a factor of two older and ~7 per cent cooler than previously thought and that their masses are much larger (by >30 per cent) than inferred from conventional, non-magnetic models in the Hertzsprung-Russell diagram. Systematic changes of this size may be of great importance in understanding the evolution of young stars, disc lifetimes and the formation of planetary systems.
• ### At the survey limits: discovery of the Aquarius 2 dwarf galaxy in the VST ATLAS and the SDSS data(1605.05338)

Aug. 15, 2016 astro-ph.GA
We announce the discovery of the Aquarius~2 dwarf galaxy, a new distant satellite of the Milky Way, detected on the fringes of the VST ATLAS and the SDSS surveys. The object was originally identified as an overdensity of Red Giant Branch stars, but chosen for subsequent follow-up based on the presence of a strong Blue Horizontal Branch, which was also used to measure its distance of $\sim 110$ kpc. Using deeper imaging from the IMACS camera on the 6.5m Baade and spectroscopy with DEIMOS on Keck, we measured the satellite's half-light radius $5.1\pm 0.8$ arcmin, or $\sim 160$ pc at this distance, and its stellar velocity dispersion of $5.4^{+3.4}_{-0.9}$ km s$^{-1}$. With $\mu=30.2$ mag arcsec$^{-2}$ and $M_V=-4.36$, the new satellite lies close to two important detection limits: one in surface brightness; and one in luminosity at a given distance, thereby making Aquarius~2 one of the hardest dwarfs to find.
• ### Chemodynamic subpopulations of the Carina dwarf galaxy(1601.01589)

Jan. 7, 2016 astro-ph.GA
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.
• ### The Gaia-ESO Survey: chemical signatures of rocky accretion in a young solar-type star(1509.00933)

Sept. 11, 2015 astro-ph.SR
It is well known that newly formed planetary systems undergo processes of orbital reconfiguration and planetary migration. As a result, planets or protoplanetary objects may accrete onto the central star, being fused and mixed into its external layers. If the accreted mass is sufficiently high and the star has a sufficiently thin convective envelope, such events may result in a modification of the chemical composition of the stellar photosphere in an observable way, enhancing it with elements that were abundant in the accreted mass. The recent Gaia-ESO Survey observations of the 10-20 Myr old Gamma Velorum cluster have enabled identifying a star that is significantly enriched in iron with respect to other cluster members. In this Letter we further investigate the abundance pattern of this star, showing that its abundance anomaly is not limited to iron, but is also present in the refractory elements, whose overabundances are correlated with the condensation temperature. This finding strongly supports the hypothesis of a recent accretion of rocky material.
• ### An extended view of the Pisces Overdensity from the SCUSS survey(1508.01272)

Sept. 7, 2015 astro-ph.GA
SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.
• ### The Gaia-ESO Survey: Characterisation of the [alpha/Fe] sequences in the Milky Way discs(1507.08066)

July 29, 2015 astro-ph.GA
We investigate, using the Gaia-ESO Survey internal Data-Release 2, the properties of the double sequence of the Milky Way discs (defined chemically as the high-alpha and low-alpha populations), and discuss their compatibility with discs defined by other means such as metallicity, kinematics or positions. This investigation uses two different approaches: in velocity space for stars located in the extended Solar neighbourhood, and in chemical space for stars at different ranges of Galactocentric radii and heights from the plane. The separation we find in velocity space allows us to investigate, in a novel manner, the extent in metallicity of each of the two sequences, identifying them with the two discs, without making any assumption about the shape of their metallicity distribution functions. Then, using the separation in chemical space, we characterise the spatial variation of the slopes of the [alpha/Fe] - [Fe/H] sequences for the thick and thin discs and the way in which the relative proportions of the two discs change across the Galaxy. We find that the thick disc (high-alpha sequence), extends up to [Fe/H]~ +0.2 and the thin disc (low-alpha sequence), at least down to [Fe/H]~ -0.8. Radial and vertical gradients in alpha-abundances are found for the thin disc, with mild spatial variations in its [alpha/Fe] - [Fe/H] paths, whereas for the thick disc we do not detect any such spatial variations. The small variations in the spatial [alpha/Fe] - [Fe/H] paths of the thin disc do not allow us to distinguish between formation models of this structure. On the other hand, the lack of radial gradients and [alpha/Fe] - [Fe/H] variations for the thick disc indicate that the mechanism responsible for the mixing of the metals in the young Galaxy (e.g. radial stellar migration or turbulent gaseous disc) was more efficient before the (present) thin disc started forming.
• ### The Catalina Surveys Periodic Variable Star Catalog(1405.4290)

May 16, 2014 astro-ph.SR
We present ~47,000 periodic variables found during the analysis of 5.4 million variable star candidates within a 20,000 square degree region covered by the Catalina Surveys Data Release-1 (CSDR1). Combining these variables with type-ab RR Lyrae from our previous work, we produce an on-line catalog containing periods, amplitudes, and classifications for ~61,000 periodic variables. By cross-matching these variables with those from prior surveys, we find that > 90% of the ~8,000 known periodic variables in the survey region are recovered. For these sources we find excellent agreement between our catalog and prior values of luminosity, period and amplitude, as well as classification. We investigate the rate of confusion between objects classified as contact binaries and type-c RR Lyrae (RRc's) based on periods, colours, amplitudes, metalicities, radial velocities and surface gravities. We find that no more than few percent of these variables in these classes are misidentified. By deriving distances for this clean sample of ~5,500 RRc's, we trace the path of the Sagittarius tidal streams within the Galactic halo. Selecting 146 outer-halo RRc's with SDSS radial velocities, we confirm the presence of a coherent halo structure that is inconsistent with current N-body simulations of the Sagittarius tidal stream. We also find numerous long-period variables that are very likely associated within the Sagittarius tidal streams system. Based on the examination of 31,000 contact binary light curves we find evidence for two subgroups exhibiting irregular lightcurves. One subgroup presents significant variations in mean brightness that are likely due to chromospheric activity. The other subgroup shows stable modulations over more than a thousand days and thereby provides evidence that the O'Connell effect is not due to stellar spots.
• ### Precession of the Sagittarius stream(1301.7069)

Dec. 10, 2013 astro-ph.GA
Using a variety of stellar tracers -- blue horizontal branch stars, main-sequence turn-off stars and red giants -- we follow the path of the Sagittarius (Sgr) stream across the sky in Sloan Digital Sky Survey data. Our study presents new Sgr debris detections, accurate distances and line-of-sight velocities that together help to shed new light on the puzzle of the Sgr tails. For both the leading and the trailing tail, we trace the points of their maximal extent, or apo-centric distances, and find that they lie at $R^L$ = 47.8 $\pm$ 0.5 kpc and $R^T$ = 102.5 $\pm$ 2.5 kpc respectively. The angular difference between the apo-centres is 93.2 $\pm$ 3.5 deg, which is smaller than predicted for logarithmic haloes. Such differential orbital precession can be made consistent with models of the Milky Way in which the dark matter density falls more quickly with radius. However, currently, no existing Sgr disruption simulation can explain the entirety of the observational data. Based on its position and radial velocity, we show that the unusually large globular cluster NGC 2419 can be associated with the Sgr trailing stream. We measure the precession of the orbital plane of the Sgr debris in the Milky Way potential and show that, surprisingly, Sgr debris in the primary (brighter) tails evolves differently to the secondary (fainter) tails, both in the North and the South.
• ### Automated search for galactic star clusters in large multiband surveys: I. Discovery of 15 new open clusters in the Galactic anticenter region(0709.1275)

May 9, 2008 astro-ph
Aims: According to some estimations, there are as many as 100000 open clusters in the Galaxy, but less than 2000 of them have been discovered, measured, and cataloged. We plan to undertake data mining of multiwavelength surveys to find new star clusters. Methods: We have developed a new method to search automatically for star clusters in very large stellar catalogs, which is based on convolution with density functions. We have applied this method to a subset of the Two Micron All Sky Survey catalog toward the Galactic anticenter. We also developed a method to verify whether detected stellar groups are real star clusters, which tests whether the stars that form the spatial density peak also fall onto a single isochrone in the color-magnitude diagram. By fitting an isochrone to the data, we estimate at the same time the main physical parameters of a cluster: age, distance, color excess. Results: For the present paper, we carried out a detailed analysis of 88 overdensity peaks detected in a field of $16\times16$ degrees near the Galactic anticenter. From this analysis, 15 overdensities were confirmed to be new open clusters and the physical and structural parameters were determined for 12 of them; 10 of them were previously suspected to be open clusters by Kronberger (2006) and Froebrich (2007). The properties were also determined for 13 yet-unstudied known open clusters, thus almost tripling the sample of open clusters with studied parameters in the anticenter. The parameters determined with this method showed a good agreement with published data for a set of well-known clusters.