Search

The ESO public survey VISTA Variables in the V\'ia L\'actea (VVV) surveyed the inner Galactic bulge and the adjacent southern Galactic disk from $2009-2015$. Upon its conclusion, the complementary VVV eXtended (VVVX) survey has expanded both the temporal as well as spatial coverage of the original VVV area, widening it from $562$ to $1700$ sq. deg., as well as providing additional epochs in $JHK_{\rm s}$ filters from $2016-2023$. With the completion of VVVX observations during the first semester of 2023, we present here the observing strategy, a description of data quality and access, and the legacy of VVVX. VVVX took $\sim 2000$ hours, covering about 4% of the sky in the bulge and southern disk. VVVX covered most of the gaps left between the VVV and the VISTA Hemisphere Survey (VHS) areas and extended the VVV time baseline in the obscured regions affected by high extinction and hence hidden from optical observations. VVVX provides a deep $JHK_{\rm s}$ catalogue of $\gtrsim 1.5\times10^9$ point sources, as well as a $K_{\rm s}$ band catalogue of $\sim 10^7$ variable sources. Within the existing VVV area, we produced a $5D$ map of the surveyed region by combining positions, distances, and proper motions of well-understood distance indicators such as red clump stars, RR Lyrae, and Cepheid variables. In March 2023 we successfully finished the VVVX survey observations that started in 2016, an accomplishment for ESO Paranal Observatory upon 4200 hours of observations for VVV+VVVX. The VVV+VVVX catalogues complement those from the Gaia mission at low Galactic latitudes and provide spectroscopic targets for the forthcoming ESO high-multiplex spectrographs MOONS and 4MOST., Comment: 17 pages, 11 figures (+ appendix). Accepted for publication in Astronomy and Astrophysics in section 14: Catalogs and data

Planetary nebulae (PNe) are excellent tracers of the metal-poor haloes of nearby early-type galaxies. They are commonly used to trace spatial distribution and kinematics of the halo and intracluster light at distances of up to 100 Mpcs. The results on the early-type galaxy M105 in the Leo I group represent a benchmark for the quantitative analysis of halo and intragroup light. Since the Leo I group lies at just a 10 Mpc distance, it is at the ideal location to compare results from resolved stellar populations with the homogeneous constraints over a much larger field of view from the PN populations. In M105, we have -- for the first time -- established a direct link between the presence of a metal-poor halo as traced by resolved red-giant branch stars and a PN population with a high specific frequency ($\alpha$-parameter). This confirms our inferences that the high $\alpha$-parameter PN population in the outer halo of M49 in the Virgo Cluster traces the metal-poor halo and intra-group light., Comment: 12 pages, 5 figures, accepted for publication in the proceedings of IAU Symposium 384: Planetary Nebulae: a Universal Toolbox in the Era of Precision Astrophysics

Mass and angular momentum are key parameters of galaxies. Their co-evolution establishes an empirical relation between the specific stellar angular momentum j* and the stellar mass M* that depends on morphology. In this work, we measure j* in a sample of 32 early type galaxies (ETGs) from the ePN.S survey, using full 2D kinematic information out to a mean 6Re. We present lambda and j* profiles. We derive the distribution of these galaxies on the total j*-M* plane and determine the ratio between the stellar j* and the specific angular momentum of the host dark matter halo. The radially extended, 2D kinematic data show that the stellar halos of ETGs do not contain large stellar mass fractions of high j*. The j*-profiles of fast-rotator ETGs are largely converged within the range of the data. For slow rotators, j* is still rising and is estimated to increase beyond 6Re by up to 40%. More than 60% of their halo angular momentum is in misaligned rotation. We find that the ePN.S ETG sample displays the well-known correlation between j*, M*, and morphology: elliptical galaxies have systematically lower j* than similar mass S0 galaxies. However, fast and slow rotators lie on the same relation within errors with the slow rotators falling at the high M* end. A power-law fit to the j*-M* relation gives a slope of 0.55+-0.17 for the S0s and 0.76+-0.23 for the ellipticals, with normalisation about 4 and 9 times lower than spirals, respectively. The estimated retained fraction of angular momentum at 10^10-10^10.5 Msun is 25% for S0s and >10% for ellipticals, and decreases by 1.5 orders of magnitude at M*~10^12 Msun. Our results show that ETGs have substantially lower j* than spiral galaxies with similar M*. Their j* must be lost during their evolution, and/or retained in the hot gas component and the satellite galaxies that have not yet merged with the central galaxy. [abridged], Comment: 26 pages, 14 figures, accepted for publication in A&A

The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022., Comment: 38 pages. A&A in press

In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come., Comment: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

M105 is an early-type galaxy in the nearby Leo I group, the closest galaxy group to contain all galaxy types and therefore an excellent environment to explore the low-mass end of intra-group light (IGL) assembly. We present a new extended kinematic survey of planetary nebulae (PNe) in M105 and the surrounding 30'x30' in the Leo I group with the Planetary Nebula Spectrograph. We use PNe as kinematic tracers of the diffuse stellar light in the halo and IGL and employ Gaussian mixture models to separate contributions from the companion galaxy NGC 3384 and associate PNe with halo and IGL components around M105. We present a catalogue of 314 PNe and firmly associate 93 with NGC 3384 and 169 with M105. The PNe in M105 are further associated with its halo and the surrounding exponential envelope. We construct smooth velocity and velocity dispersion fields and calculate projected rotation, velocity dispersion, and $\lambda_R$ profiles for each component. Halo PNe exhibit declining velocity dispersion and rotation profiles, while the velocity dispersion and rotation of the exponential envelope increase notably at large radii. We identify three regimes with distinct kinematics that are linked to distinct stellar population properties: (i) the rotating core (within $1~R_\mathrm{eff}$) formed in situ and dominated by metal-rich ([M/H]~0) stars likely formed in situ, (ii) the halo from 1 to $7.5~R_\mathrm{eff}$ consisting of intermediate-metallicity stars ([M/H]>-1), either formed in situ or brought in through major mergers, and (iii) the exponential envelope reaching beyond our farthest data point at 16 $R_\mathrm{eff}$, predominately composed of metal-poor ([M/H]<-1) stars. The high velocity dispersion and moderate rotation of the latter are consistent with that measured for dwarf satellite galaxies in the Leo I group, indicating that the exponential envelope traces the transition to the IGL., Comment: Submitted to A&A, comments welcome. 20 pages, 14 figures, abstract abridged

We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a self-consistent dynamical model of the MW defined by the best available parameters for the dark matter halo, the stellar disc and the bulge. For the known (stellar and gas) disc mass (4.5 x 10^10 Msun) and disc mass fraction at R ~ 2.2 R_d (f_d ~ 0.3 - 0.6), the low mass limit does not yield a bar in a Hubble time. The high mass limit adopted here produces a box/peanut bar within about 2 Gyr with the correct mass (~10^10 Msun), size (~5 kpc) and peak pattern speed (~ 40-45 km/s/kpc). In agreement with earlier work, the bar formation timescale scales inversely with f_d (i.e. log [T/Gyr] ~ 0.60/f_d - 0.83 for 1 < f_d < 0.3). The disc radial heating is strong, but, in contrast to earlier claims, we find that disc vertical heating outside of the box/peanut bulge structure is negligible. The synthetic barred MW exhibits long-term stability, except for the slow decline (roughly -2 km/s/kpc/Gyr) of the bar pattern speed, consistent with recent estimates. If our model is indicative of the Milky Way, we estimate that the bar first emerged 3-4 Gyr ago., Comment: We welcome comments and missing references

We use the data-driven method, The Cannon, to bring 21,000 stars from the ARGOS bulge survey, including 10,000 red clump stars, onto the parameter and abundance scales of the cross-Galactic survey, APOGEE, obtaining rms precisions of 0.10 dex, 0.07 dex, 74 K, and 0.18 dex for [Fe/H], [Mg/Fe], Teff, and log(g), respectively. The re-calibrated ARGOS survey - which we refer to as the A2A survey - is combined with the APOGEE survey to investigate the abundance structure of the Galactic bulge. We find X-shaped [Fe/H] and [Mg/Fe] distributions in the bulge that are more pinched than the bulge density, a signature of its disk origin. The mean abundance along the major axis of the bar varies such that the stars are more [Fe/H]-poor and [Mg/Fe]-rich near the Galactic center than in the long bar/outer bulge region. The vertical [Fe/H] and [Mg/Fe] gradients vary between the inner bulge and long bar with the inner bulge showing a flattening near the plane that is absent in the long bar. The [Fe/H]-[Mg/Fe] distribution shows two main maxima, an ``[Fe/H]-poor [Mg/Fe]- rich'' maximum and an ``[Fe/H]-rich [Mg/Fe]-poor'' maximum, that vary in strength with position in the bulge. In particular, the outer long bar close to the Galactic plane is dominated by super-solar [Fe/H], [Mg/Fe]-normal stars. Stars composing the [Fe/H]-rich maximum show little kinematic dependence on [Fe/H], but for lower [Fe/H] the rotation and dispersion of the bulge increase slowly. Stars with [Fe/H]<-1 dex have a very different kinematic structure than stars with higher [Fe/H]. Comparing with recent models for the Galactic boxy-peanut bulge, the abundance gradients and distribution, and the relation between [Fe/H] and kinematics suggest that the stars comprising each maximum have separate disk origins with the ``[Fe/H]-poor [Mg/Fe]-rich'' stars originating from a thicker disk than the ``[Fe/H]-rich [Mg/Fe]-poor'' stars., Comment: 26 pages, 29 figures, revised version as submitted to A&A following referee report

The nuclear stellar disc (NSD) is, together with the nuclear star cluster (NSC) and the central massive black hole, one of the main components in the central parts of our Milky Way. However, until recently, only few studies of the stellar content of the NSD have been obtained due to extreme extinction and stellar crowding. With a dedicated KMOS (VLT, ESO) spectroscopic survey, we study the kinematics and global metallicities of the NSD based on the observations of K/M giant stars. We trace radial velocities and metallicities which were derived based on spectral indices (Na I and CO) along the NSD and compare those with a Galactic Bulge sample of APOGEE (DR16) and data from the NSC. We find that the metallicity distribution function and the fraction of metal-rich and metal-poor stars in the NSD are different from the corresponding distributions and ratios of the NSC and the Galactic Bulge. By tracing the velocity dispersion as a function of metallicity, we clearly see that the NSD is kinematically cool and that the velocity dispersion decreases with increasing metallicity contrary to the inner Bulge sample of APOGEE ($\rm |b| < 4^{o}$). Using molecular gas tracers ($\rm H_{2}CO$, CO(4-3)) of the Central Molecular Zone (CMZ) we find an astonishing agreement between the gas rotation and the rotation of the metal-rich population indicating that the metal-rich stars could have formed from gas in the CMZ. On the other hand, the metal-poor stars show a much slower rotation profile with signs of counter-rotation indicating a different origin of these stars. Coupling kinematics with global metallicities, our results demonstrate that the NSD is chemically and kinematically distinct with respect to the inner Bulge indicating a different formation scenario., Comment: accepted for publication in A\&A

The mass-to-light ratio (M/L) is a key parameter in decomposing galactic rotation curves into contributions from the baryonic components and the dark halo of a galaxy. One direct observational method to determine the disc M/L is by calculating the surface mass density of the disc from the stellar vertical velocity dispersion and the scale height of the disc. Usually, the scale height is obtained from near-IR studies of edge-on galaxies and pertains to the older, kinematically hotter stars in the disc, while the vertical velocity dispersion of stars is measured in the optical band and refers to stars of all ages (up to ~10 Gyr) and velocity dispersions. This mismatch between the scale height and the velocity dispersion can lead to underestimates of the disc surface density and a misleading conclusion of the sub-maximality of galaxy discs. In this paper we present the study of the stellar velocity dispersion of the disc galaxy NGC 6946 using integrated star light and individual planetary nebulae as dynamical tracers. We demonstrate the presence of two kinematically distinct populations of tracers which contribute to the total stellar velocity dispersion. Thus, we are able to use the dispersion and the scale height of the same dynamical population to derive the surface mass density of the disc over a radial extent. We find the disc of NGC 6946 to be closer to maximal with the baryonic component contributing most of the radial gravitational field in the inner parts of the galaxy (Vmax(bar) = 0.76($\pm$0.14)Vmax)., Comment: 17 pages, 14 figures, Accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:1802.00465

Stellar halos in early-type galaxies (ETGs) are shaped by their accretion and merger histories. We use a sample of 1114 ETGs in the TNG100 simulation with stellar masses $10^{10.3}\leq M_{*}/M_\odot\leq 10^{12}$, selected at z=0 within the range of g-r colour and lambda-ellipticity diagram populated by observed ETGs. We study how the rotational support and intrinsic shapes of the stellar halos depend on the fraction of accreted stars, overall and separately by major, minor, and mini mergers. Accretion histories in TNG100 ETGs as well as the radial distributions of ex-situ stars $f_{ex}(R)$ strongly correlate with stellar mass. Low-mass ETGs have characteristic peaked rotation profiles and near-oblate shapes with rounder halos that are completely driven by the in-situ stars. At high $f_{ex}$ major mergers decrease the in-situ peak in rotation velocity, flatten the $V_{*}/\sigma_{*}(R)$ profiles, and increase the triaxiality of the stellar halos. Kinematic transition radii do not trace the transition between in-situ and ex-situ dominated regions, but for systems with $M_{*}>10^{10.6}M_\odot$ the local rotational support decreases with the local ex-situ fraction $f_{ex}(R)$ and triaxiality increases with $f_{ex}$. These correlations are followed by fast and slow rotators alike with a continuous and overlapping sequence of properties. Merger events dynamically couple stars and dark matter: in high mass ETGs and at large radii where $f_{ex}\gtrsim0.5$, both components tend to have similar intrinsic shapes and rotational support, and nearly aligned principal axes and spin directions. Based on these results we suggest that extended photometry and kinematics of massive ETGs ($M_{*}>10^{10.6}M_\odot$) can be used to estimate the local fraction of ex-situ stars and to approximate the intrinsic shapes and rotational support of the co-spatial dark matter component. [abridged], Comment: 22 pages, 17 figures, accepted for publication on A&A update version including a new section on high-redshift compact progenitor galaxies and a new Appendix with with tables reporting relevant physical relations derived throughout this paper

M105 (NGC 3379) is an early-type galaxy in the Leo I group. This group is the nearest group that contains all main galaxy types and can thus be used as a benchmark to study the properties of the intra-group light (IGL) in low-mass groups. We use PNe as discrete stellar tracers of the diffuse light around M105. PNe were identified on the basis of their bright [OIII]5007 AA emission and the absence of a broad-band continuum. We compare the PN number density profile with the galaxy surface-brightness profile decomposed into metallicity components using published HST photometry in two halo fields. We identify 226 PNe candidates within a limiting magnitude of mlim = 28.1 from our Subaru-SuprimeCam imaging, covering 67.6 kpc along the major axis of M105 and the halos of NGC 3384 and NGC 3398. We find an excess of PNe at large radii compared to the stellar surface brightness profile from broad-band surveys. This excess is related to a variation in the luminosity-specific PN number $\alpha$ with radius. The $\alpha$-parameter value of the extended halo is more than 7 times higher than that of the inner halo. We also measure an increase in the slope of the PN luminosity function at fainter magnitudes with radius. We infer that the radial variation of the PN population properties is due to a diffuse population of metal-poor stars ([M/H] < -1.0) following an exponential profile, in addition to the M105 halo. The spatial coincidence between the number density profile of these metal-poor stars and the increase in the $\alpha$-parameter value with radius establishes the missing link between metallicity and the post-AGB phases of stellar evolution. We estimate that the total bolometric luminosity associated with the exponential IGL population is 2.04x10^9 Lsun as a lower limit, corresponding to an IGL fraction of 3.8%. This work sets the stage for kinematic studies of the IGL in low-mass groups., Comment: 22 pages, 18 figures, accepted for publication in Astronomy & Astrophysics. Abridged abstract

We characterize the photometric and kinematic properties of simulated early-type galaxy (ETG) stellar halos, and compare them to observations. We select a sample of ~1200 ETGs in the TNG100 and TNG50 simulations, spanning a stellar mass range of $10^{10.3}-10^{12}M_{\odot}$ and within the range of (g-r) colour and lambda-ellipticity diagram populated by observed ETGs. We determine photometric parameters, intrinsic shapes, and kinematic observables in their extended stellar halos. We study the variation in kinematics from center to halo and connect it to a change in the intrinsic shape of the galaxies. We find that the simulated galaxy sample reproduces the diversity of kinematic properties observed in ETG halos. Simulated fast rotators (FRs) divide almost evenly in one third having flat lambda profiles and high halo rotational support, a third with gently decreasing profiles, and another third with low halo rotation. Slow rotators (SRs) tend to have increased rotation in the outskirts, with half of them exceeding lambda=0.2. For $M_{*}>10^{11.5}M_{\odot}$ halo rotation is unimportant. A similar variety of properties is found for the stellar halo intrinsic shapes. Rotational support and shape are deeply related: the kinematic transition to lower rotational support is accompanied by a change towards rounder intrinsic shape. Triaxiality in the halos of FRs increases outwards and with stellar mass. Simulated SRs have relatively constant triaxiality profiles. Simulated stellar halos show a large variety of structural properties, with quantitative but no clear qualitative differences between FRs and SRs. At the same stellar mass, stellar halo properties show a gradual transition and significant overlap between the two families, despite the clear bimodality in the central regions. This is in agreement with observations of extended photometry and kinematics. [abridged], Comment: accepted for publication in A&A, 25 pages, 22 figures

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare., Comment: Accepted in A&A

In this paper we introduce a new method for analysing Milky Way phase-space which allows us to reveal the imprint left by the Milky Way bar and spiral arms on the stars with full phase-space data in Gaia Data Release 2. The unprecedented quality and extended spatial coverage of these data enable us to discover six prominent stellar density structures in the disc to a distance of 5 kpc from the Sun. Four of these structures correspond to the spiral arms detected previously in the gas and young stars (Scutum-Centaurus, Sagittarius, Local and Perseus). The remaining two are associated with the main resonances of the Milky Way bar where corotation is placed at around 6.2 kpc and the outer Lindblad resonance beyond the Solar radius, at around 9 kpc. For the first time we provide evidence of the imprint left by spiral arms and resonances in the stellar densities not relying on a specific tracer, through enhancing the signatures left by these asymmetries. Our method offers new avenues for studying how the stellar populations in our Galaxy are shaped., Comment: 10 pages, 6 figures, accepted for publication in A&A Letters

We use a sample of 938 red clump giant stars located in the direction of the galactic long bar to study the chemistry of Milky Way bar stars. Kinematically separating stars on bar orbits from stars with inner disc orbits, we find that stars on bar-like orbits are more metal rich with a mean iron abundance of <[Fe/H]>=+0.30 compared to <[Fe/H]>=+0.03 for the inner disc. Spatially selecting bar stars is complicated by a strong vertical metallicity gradient of -1.1dex/kpc, but we find the metallicity distribution varies in a manner consistent with our orbital selection. Our results have two possible interpretations. The first is that the most metal rich stars in the inner Galaxy pre-existed the bar, but were kinematically cold at the time of bar formation and therefore more easily captured onto bar orbits when the bar formed. The second is that the most metal rich stars formed after the bar, either directly onto the bar following orbits or were captured by the bar after their formation., Comment: 7 pages, 6 figures. Accepted by A&A

The second data release of the Gaia mission has revealed, in stellar velocity and action space, multiple ridges, the exact origin of which is still debated. Recently, we demonstrated that a large Galactic bar with pattern speed 39 km/s/kpc does create most of the observed ridges. Among those ridges, the Hercules moving group would then be associated to orbits trapped at the co-rotation resonance of the bar. Here we show that a distinctive prediction of such a model is that the angular momentum of Hercules at the Sun's radius must significantly decrease with increasing Galactocentric azimuth, i.e. when getting closer to the major axis of the bar. We show that such a dependence of the angular momentum of trapped orbits on the azimuth would on the other hand not happen close to the outer Lindblad resonance of a faster bar, unless the orbital distribution is still far from phase-mixed, namely for a bar perturbation younger than ~ 2 Gyr. Using Gaia DR2 and Bayesian distances from the StarHorse code, and tracing the average Galactocentric radial velocity as a function of angular momentum and azimuth, we show that the Hercules angular momentum changes significantly with azimuth as expected for the co-rotation resonance of a dynamically old large bar., Comment: 4 pages, 2 figures, A&A in press

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05., Comment: 9 pages, 7 figures, submitted to A&A

The signatures of the formation and evolution of a galaxy are imprinted in its stars. Their velocities, ages, and chemical compositions present major constraints on models of galaxy formation, and on various processes such as the gas inflows and outflows, the accretion of cold gas, radial migration, and the variability of star formation activity. Understanding the evolution of the Milky Way requires large observational datasets of stars via which these quantities can be determined accurately. This is the science driver of the 4MOST MIlky way Disc And BuLgE High-Resolution (4MIDABLE-HR) survey: to obtain high-resolution spectra at $R \sim 20\,000$ and to provide detailed elemental abundances for large samples of stars in the Galactic disc and bulge. High data quality will allow us to provide accurate spectroscopic diagnostics of two million stellar spectra: precise radial velocities; rotation; abundances of many elements, including those that are currently only accessible in the optical, such as Li, s-, and r-process; and multi-epoch spectra for a sub-sample of stars. Synergies with complementary missions like Gaia and TESS will provide masses, stellar ages and multiplicity, forming a multi-dimensional dataset that will allow us to explore and constrain the origin and structure of the Milky Way., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

We present a kinematic study of a sample of 298 planetary nebulas (PNs) in the outer halo of the central Virgo galaxy M87 (NGC~4486). The line-of-sight velocities of these PNs are used to identify sub-components, to measure the angular momentum content of the smooth M87 halo, and to constrain the orbital distribution of the stars at these large radii. We find that the surface density profile of the ICPNS at $100\, \mathrm{kpc}$ radii has a shallow logarithmic slope, $-\alpha_{\rm ICL}\simeq -0.8$, dominating the light at the largest radii. The angular momentum-related $\lambda_{\mathrm{R}}$ profile for the smooth halo remains in the slow rotator regime, out to 135 kpc average ellipse radius. The velocity dispersione profile is measured to decrease sharply to $\sigma_{\mathrm{halo}} \sim 100\, \mathrm{kms}^{-1}$ at $R_{\rm avg}=135$ kpc. Such a steeply decline of the velocity dispersion and the surface density profile of the smooth halo can be reconciled with the circular velocity curve inferred from assuming hydrostatic equilibrium for the hot X-ray gas if the orbit distribution of the smooth M87 halo changes strongly from approximately isotropic within $R_{\rm avg}\sim 60$ kpc to very radially anisotropic at the largest distances probed. The extended LOSVD of the PNs in the M87 halo allows the identification of several subcomponents: the ICPNs, the "crown" accretion event, and the smooth M87 halo. In galaxies like M87, the presence of these sub-components needs to be taken into account to avoid systematic biases in estimating the total enclosed mass. The dynamical structure inferred from the velocity dispersion profile indicates that the smooth halo of M87 steepens beyond $R_{\rm avg}=60$ kpc and becomes strongly radially anisotropic, and that the velocity dispersion profile is consistent with the X-ray circular velocity curve at these radii without non-thermal pressure effects., Comment: 24 pages. Accepted on A\&A

We continue the analysis of the dataset of our spectroscopic observation campaign of M31, by deriving simple stellar population properties (age metallicity and alpha-elements overabundance) from the measurement of Lick/IDS absorption line indices. We describe their two-dimensional maps taking into account the dust distribution in M31. 80\% of the values of our age measurements are larger than 10 Gyr. The central 100 arcsec of M31 are dominated by the stars of the classical bulge of M31. They are old (11-13 Gyr), metal-rich (as high as [Z/H]~0.35 dex) at the center with a negative gradient outwards and enhanced in alpha-elements ([alpha/Fe]~ 0.28+- 0.01 dex). The bar stands out in the metallicity map, where an almost solar value of [Z/H] (~0.02+-0.01 dex) with no gradient is observed along the bar position angle (55.7 deg) out to 600 arcsec from the center. In contrast, no signature of the bar is seen in the age and [alpha/Fe] maps, that are approximately axisymmetric, delivering a mean age and overabundance for the bar and the boxy-peanut bulge of 10-13 Gyr and 0.25-0.27 dex, respectively. The boxy/peanut-bulge has almost solar metallicity (-0.04+- 0.01 dex). The mass-to-light ratio of the three components is approximately constant at M/LV ~ 4.4-4.7 Msol/Lsol. The disk component at larger distances is made of a mixture of stars, as young as 3-4 Gyr, with solar metallicity and smaller M/LV (~3+-0.1 Msol/Lsol). We propose a two-phase formation scenario for the inner region of M31, where most of the stars of the classical bulge come into place together with a proto-disk, where a bar develops and quickly transforms it into a boxy-peanut bulge. Star formation continues in the bulge region, producing stars younger than 10 Gyr, in particular along the bar, enhancing its metallicity. The disk component appears to build up on longer time-scales., Comment: Language-edited version, Accepted for publication in A&A

M49 (NGC 4472) is the dominant galaxy in subcluster B of the Virgo Cluster, and a benchmark for studying the build-up of the extended halos of brightest group galaxies in the outskirts of galaxy clusters. We investigate the kinematics in the outer halo of M49, look for substructures, and describe the transition to the surrounding intra-group light. As kinematic tracers we use planetary nebulae (PNe), combining kinematics from the extended Planetary Nebula Spectrograph (PN.S) early-type galaxy survey with our recent deep photometric sample. We study the position-velocity-plane for bright and faint PN populations out to 95 kpc radius, and employ a multi-Gaussian model for the velocity distribution to identify stellar populations with distinct kinematics and histories. We report the detection of stellar-kinematic substructure associated with the interaction of M49 with the dwarf irregular galaxy VCC 1249. We find two kinematically distinct PN populations associated with the main M49 halo and the extended intra-group light (IGL). These have velocity dispersions $\sigma_{halo}\simeq 170$ km/s and $\sigma_{IGL} \simeq400$ km/s at 10-80 kpc radii. The overall luminosity profile and velocity dispersion at $\sim80$ kpc are consistent with a flat circular velocity curve extrapolated from X-ray observations. The dispersion of the PNe associated with the IGL joins onto that of the satellite galaxies in subcluster B at $\sim100$ kpc radius. This is the first time that the transition from halo to IGL is observed based on the velocities of individual stars. Therefore the halo of M49, consisting of at least three distinct components, has undergone an extended accretion history within its parent group potential. The blue colours of the IGL component are consistent with a population of stars formed in low-mass galaxies at redshift $\sim0.5$ that has since evolved passively, as suggested by other data., Comment: 10 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

The Galactic Bulge can uniquely be studied from large samples of individual stars, and is therefore of prime importance for understanding the stellar population structure of bulges in general. Here the observational evidence on the kinematics, chemical composition, and ages of Bulge stellar populations based on photometric and spectroscopic data is reviewed. The bulk of Bulge stars are old and span a metallicity range -1.5<~[Fe/H]<~+0.5. Stellar populations and chemical properties suggest a star formation timescale below ~2 Gyr. The overall Bulge is barred and follows cylindrical rotation, and the more metal-rich stars trace a Box/Peanut (B/P) structure. Dynamical models demonstrate the different spatial and orbital distributions of metal-rich and metal-poor stars. We discuss current Bulge formation scenarios based on dynamical, chemical, chemodynamical and cosmological models. Despite impressive progress we do not yet have a successful fully self-consistent chemodynamical Bulge model in the cosmological framework, and we will also need more extensive chrono-chemical-kinematic 3D map of stars to better constrain such models., Comment: 9 figures, 55 pages final version to appear in the Annual Reviews of Astronomy & Astrophysics, volume 56

We present the largest and most spatially extended planetary nebulae (PNe) catalog ever obtained for the Fornax cluster. We measured velocities of 1452 PNe out to 200 kpc in the cluster core using a counter-dispersed slitless spectroscopic technique with data from FORS2 on the VLT. With such extended spatial coverage, we can study separately the stellar halos of some of the cluster main galaxies and the intracluster light. In this second paper of the Fornax Cluster VLT Spectroscopic Survey (FVSS), we identify and classify the emission-line sources, describe the method to select PNe and calculate their coordinates and velocities from the dispersed slitless images. From the PN 2D velocity map we identify stellar streams that are possibly tracing the gravitational interaction of NGC1399 with NGC1404 and NGC1387. We also present the velocity dispersion profile out to $\sim 200$ kpc radii, which shows signatures of a superposition of the bright central galaxy and the cluster potential, with the latter clearly dominating the regions outside R $\sim 1000$" ($\sim 100$ kpc)., Comment: 13 pages, 8 figures, Accepted for publication in MNRAS, companion paper FVSS-I (Pota et al, 2018, submitted)

The decomposition of the rotation curve of galaxies into contribution from the disc and dark halo remains uncertain and depends on the adopted mass to light ratio (M/L) of the disc. Given the vertical velocity dispersion of stars and disc scale height, the disc surface mass density and hence the M/L can be estimated. We address a conceptual problem with previous measurements of the scale height and dispersion. When using this method, the dispersion and scale height must refer to the same population of stars. The scale height is obtained from near-IR studies of edge-on galaxies and is weighted towards older kinematically hotter stars, whereas the dispersion obtained from integrated light in the optical bands includes stars of all ages. We aim to extract the dispersion for the hotter stars, so that it can then be used with the correct scale height to obtain the disc surface mass density. We use a sample of planetary nebulae (PNe) as dynamical tracers in the face-on galaxy NGC 628. We extract two different dispersions from its velocity histogram -- representing the older and younger PNe. We also present complementary stellar absorption spectra in the inner regions of this galaxy and use a direct pixel fitting technique to extract the two components. Our analysis concludes that previous studies, which do not take account of the young disc, underestimate the disc surface mass density by a factor of ~ 2. This is sufficient to make a maximal disc for NGC 628 appear like a submaximal disc., Comment: 24 pages, 18 figures, 7 tables; accepted for publication in MNRAS

In the hierarchical two-phase formation scenario, the extended halos of early type galaxies (ETGs) are expected to have different physical properties from those of the galaxies' central regions. This work aims at characterizing the kinematic properties of ETG halos using planetary nebulae (PNe) as tracers, which allow us to overcome the limitations of absorption line spectroscopy of continuum at low surface brightness. We present two-dimensional velocity and velocity dispersion fields for 33 ETGs, including both fast (FRs) and slow rotators (SRs), making this the largest kinematic survey to-date of extragalactic PNe. The velocity fields extend out to a median 5.6 effective radii (Re), with a range [3Re-13Re] for the PN.S ETGs. We complemented the PN kinematics with absorption line data from the literature. We find that ETGs typically show a kinematic transition between inner regions and halo. Estimated transition radii in units of Re anticorrelate with stellar mass. SRs have increased but still modest rotational support at large radii, while most of the FRs show a decrease in rotation, due to the fading of the stellar disk in the outer, more slowly rotating spheroid. 30% of the FRs are dominated by rotation also at large radii. Most ETGs have flat or slightly falling halo velocity dispersion profiles, but 15% of the sample have steeply falling profiles. 40% of the FRs show kinematic twists, misalignments, or rotation along two axes, indicating that they turn from oblate at the center to triaxial in the halo, consistently with the distribution of isophote twists from extended photometry. ETGs have more diverse kinematic properties in their halos than in the central regions. The observed kinematic transition to the halo and its dependence on stellar mass is consistent with LambdaCDM simulations and supports a two-phase formation scenario. [abridged], Comment: Accepted for publication in A&A, 52 pages, 30 figures

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science., Comment: 17 pages, 9 figures, accepted for publication in Astronomical Notes

The galaxy M49 (NGC 4472) is the brightest early-type galaxy in the Virgo Cluster. It is located in Subcluster B and has an unusually blue, metal-poor outer halo. Planetary nebulae (PNe) are excellent tracers of diffuse galaxy and intragroup light. We present a photometric survey of PNe in the galaxy's extended halo to characterise its PN population, as well as the surrounding intragroup light (IGL) of the Subcluster B. PNe were identified based on their bright [OIII]5007 \AA\ emission and absence of a broad-band continuum. We identify 738 PNe out to a radius of 155 kpc from M49's centre from which we define a complete sample of 624 PNe within a limiting magnitude of m_5007=28.8. Comparing the PN number density to the broad-band stellar surface brightness profile, we find a variation of the PN-specific frequency (alpha-parameter) with radius. The outer halo beyond 60 kpc has a 3.2 times higher alpha-parameter compared to the main galaxy halo, which is likely due to contribution from the surrounding blue IGL. We use the Planetary Nebulae Luminosity Function (PNLF) as an indicator of distance and stellar population. Its slope, which correlates empirically with galaxy type, varies within the inner halo. In the eastern quadrant of M49, the PNLF slope is shallower, indicating an additional localised, bright PN population following an accretion event, likely that of the dwarf irregular galaxy VCC1249. We also determined a distance modulus of mu = 31.29+/-0.08 for M49, corresponding to a physical distance of 18.1+/-0.6 Mpc, which agrees with a recent surface-brightness fluctuations distance. The PN populations in the outer halo of M49 are consistent with the presence of a main Sersic galaxy halo with a slight (B-V) colour gradient of 10${}^{-4}$ mag/arcsec surrounded by intragroup light with a very blue colour of (B-V)=0.25 and a constant surface brightness mu_V=28.0 mag/arcsec${}^2$., Comment: 17 pages, 14 figures, accepted for publication in A&A

In cold dark matter cosmology, the baryonic components of galaxies are thought to be mixed with and embedded in non-baryonic and non-relativistic dark matter, which dominates the total mass of the galaxy and its dark matter halo. In the local Universe, the mass of dark matter within a galactic disk increases with disk radius, becoming appreciable and then dominant in the outer, baryonic regions of the disks of star-forming galaxies. This results in rotation velocities of the visible matter within the disk that are constant or increasing with disk radius. Comparison between the dynamical mass and the sum of stellar and cold gas mass at the peak epoch of galaxy formation, inferred from ancillary data, suggest high baryon factions in the inner, star-forming regions of the disks. Although this implied baryon fraction may be larger than in the local Universe, the systematic uncertainties (stellar initial mass function, calibration of gas masses) render such comparisons inconclusive in terms of the mass of dark matter. Here we report rotation curves for the outer disks of six massive star-forming galaxies, and find that the rotation velocities are not constant, but decrease with radius. We propose that this trend arises because of two main factors: first, a large fraction of the massive, high-redshift galaxy population was strongly baryon dominated, with dark matter playing a smaller part than in the local Universe; and second, the large velocity dispersion in high-redshift disks introduces a substantial pressure term that leads to a decrease in rotation velocity with increasing radius. The effect of both factors appears to increase with redshift. Qualitatively, the observations suggest that baryons in the early Universe efficiently condensed at the centres of dark matter halos when gas fractions were high, and dark matter was less concentrated. [Abridged], Comment: Accepted for publication in Nature. Note press embargo until 1800 London time / 1400 US Eastern Time on 15 March 2017

In decomposing the HI rotation curves of disc galaxies, it is necessary to break a degeneracy between the gravitational fields of the disc and the dark halo by estimating the disc surface density. This is done by combining measurements of the vertical velocity dispersion of the disc with the disc scale height. The vertical velocity dispersion of the discs is measured from absorption lines (near the V-band) of near-face-on spiral galaxies, with the light coming from a mixed population of giants of all ages. However, the scale heights for these galaxies are estimated statistically from near-IR surface photometry of edge-on galaxies. The scale height estimate is therefore dominated by a population of older (> 2 Gyr) red giants. In this paper, we demonstrate the importance of measuring the velocity dispersion for the same older population of stars that is used to estimate the vertical scale height. We present an analysis of the vertical kinematics of K-giants in the solar vicinity. We find the vertical velocity distribution best fit by two components with dispersions of 9.6 +/- 0.5 km/s and 18.6 +/- 1.0 km/s, which we interpret as the dispersions of the young and old disc populations respectively. Combining the (single) measured velocity dispersion of the total young + old disc population (13.0 +/- 0.1 km/s) with the scale height estimated for the older population would underestimate the disc surface density by a factor of ~ 2. Such a disc would have a peak rotational velocity that is only 70% of that for the maximal disc, thus making it appear submaximal.

We present kinematic and photometric evidence for an accretion event in the halo of the cD galaxy M87 in the last Gyr. Using velocities for ~300 planetary nebulas (PNs) in the M87 halo, we identify a chevron-like substructure in the PN phase-space. We implement a probabilistic Gaussian mixture model to identify the PNs that belong to the chevron. From analysis of deep V-band images of M87, we find that the region with the highest density of PNs associated to the chevron, is a crown-shaped substructure in the optical light. We assign a total of N_(PN,sub)=54 to the substructure, which extends over ~50 kpc along the major axis where we also observe radial variations of the ellipticity profile and a colour gradient. The substructure has highest surface brightness in a 20kpc x 60kpc region around 70 kpc in radius. In this region, it causes an increase in surface brightness by >60%. The accretion event is consistent with a progenitor galaxy with a V-band luminosity of L=2.8\pm1.0 x 10^9 L_(sun,V), a colour of (B-V)=0.76\pm0.05, and a stellar mass of M=6.4\pm2.3 x 10^9 M_sun. The accretion of this progenitor galaxy has caused an important modification of the outer halo of M87 in the last Gyr. By itself it is strong evidence that the galaxy's cD halo is growing through the accretion of smaller galaxies as predicted by hierarchical galaxy evolution models., Comment: 4 pages, 2 figures, in press, A&A

Recent observations have discovered the presence of a Box/Peanut or X-shape structure in the Galactic bulge. Such Box/Peanut structures are common in external disc galaxies, and are well-known in N-body simulations where they form following the buckling instability of a bar. From studies of analytical potentials and N-body models it has been claimed in the past that Box/Peanut bulges are supported by "bananas", or x1v1 orbits. We present here a set of N-body models where instead the peanut bulge is mainly supported by brezel-like orbits, allowing strong peanuts to form with short extent relative to the bar length. This shows that stars in the X-shape do not necessarily stream along banana orbits which follow the arms of the X-shape. The brezel orbits are also found to be the main orbital component supporting the peanut shape in our recent Made-to-Measure dynamical models of the Galactic bulge. We also show that in these models the fraction of stellar orbits that contribute to the X-structure account for 40-45% of the stellar mass., Comment: Accepted for publication in MNRAS Letters

We construct dynamical models of the Milky Way's Box/Peanut (B/P) bulge, using the recently measured 3D density of Red Clump Giants (RCGs) as well as kinematic data from the BRAVA survey. We match these data using the NMAGIC Made-to-Measure method, starting with N-body models for barred discs in different dark matter haloes. We determine the total mass in the bulge volume of the RCGs measurement (+-2.2 x +- 1.4 x +- 1.2 kpc) with unprecedented accuracy and robustness to be 1.84 +- 0.07 x10^10 Msun. The stellar mass in this volume varies between 1.25-1.6 x10^10 Msun, depending on the amount of dark matter in the bulge. We evaluate the mass-to-light and mass-to-clump ratios in the bulge and compare them to theoretical predictions from population synthesis models. We find a mass-to-light ratio in the K-band in the range 0.8-1.1. The models are consistent with a Kroupa or Chabrier IMF, but a Salpeter IMF is ruled out for stellar ages of 10 Gyr. To match predictions from the Zoccali IMF derived from the bulge stellar luminosity function requires about 40% or 0.7 x10^10 Msun dark matter in the bulge region. The BRAVA data together with the RCGs 3D density imply a low pattern speed for the Galactic B/P bulge of 25-30 km.s-1.kpc-1. This would place the Galaxy among the slow rotators (R >= 1.5). Finally, we show that the Milky Way's B/P bulge has an off-centred X structure, and that the stellar mass involved in the peanut shape accounts for at least 20% of the stellar mass of the bulge, significantly larger than previously thought., Comment: Accepted for publication in MNRAS

Tidal debris which are rich in HI gas, formed in interacting and merging systems, are suitable laboratories to study star formation outside galaxies. Recently, several such systems were observed, which contained many young star forming regions outside the galaxies. In previous works, we have studied young star forming regions outside galaxies in different systems with optical and/or gaseous tidal debris, all of them with available archive GALEX/UV images, in order to understand how often they occur and in which type of environments. In this paper we searched for star forming regions around the galaxy NGC2865, a shell galaxy which is circled by a ring of HI, with a total mass of 1.2 x 10$^9$ M$_\odot$. Using the Multi-Slit Imaging Spectroscopy Technique with the Gemini telescope, we detected all H$\alpha$ emitting sources in the surroundings of the galaxy NGC2865, down to a flux limit of 10$^{-18}$ erg cm$^{-2}$ s$^{-1}$ \AA$^{-1}$. Together with Near and Far-Ultraviolet flux information we characterize the star formation rates, masses, ages, and metallicities for these HII regions. In total, we found 26 emission-line sources in a 60 $\times$ 60 Kpc field centered over the southeastern tail of the HI gas present around the galaxy NGC2865. Out of the 26 H$\alpha$ emitters, 19 are in the satellite galaxy FGCE 0745 and seven are intergalactic HII regions scattered over the south tail of the HI gas around NGC2865. We found that the intergalactic HII regions are young ($<$200 Myr) with stellar masses in the range 4 X 10$^3$M$_\odot$ to 17x10$^6$ M$_\odot$. These are found in a region of low HI gas density, where the probability of forming stars is expected to be low. For one of the intergalactic HII regions we estimated a solar oxygen abundance of 12 + log(O/H) $\sim$ 8.7. We also were able to estimate the metallicity for the satellite galaxy FGCE0745 to be 12 + log(O/H) ~ 8.0., Comment: 13 pages, 13 figures, Accepted in A&A

We present the stellar atmospheric parameters (effective temperature, surface gravity, overall metallicity), radial velocities, individual abundances and distances determined for 425 561 stars, which constitute the fourth public data release of the RAdial Velocity Experiment (RAVE). The stellar atmospheric parameters are computed using a new pipeline, based on the algorithms of MATISSE and DEGAS. The spectral degeneracies and the 2MASS photometric information are now better taken into consideration, improving the parameter determination compared to the previous RAVE data releases. The individual abundances for six elements (magnesium, aluminum, silicon, titanium, iron and nickel) are also given, based on a special-purpose pipeline which is also improved compared to that available for the RAVE DR3 and Chemical DR1 data releases. Together with photometric information and proper motions, these data can be retrieved from the RAVE collaboration website and the Vizier database., Comment: 40 pages, 36 figures, accepted in AJ

We investigate the diffuse light in the outer regions of the nearby elliptical galaxy M87 in the Virgo cluster, using planetary nebulas (PNs) as tracers. The surveyed areas (0.43 squared degrees) cover M87 up to a radial distance of 150 kpc, in the ransition region between galaxy halo and intracluster light (ICL). All PNs are identified through the on-off band technique using automatic selection criteria based on the distribution of the detected sources in the colour-magnitude diagram and the properties of their point-spread function. We extract a catalogue of 688 objects down to m_5007=28.4, with an estimated residual contamination from foreground stars and background Lyalpha galaxies, which amounts to ~35% of the sample. This is one of the largest extragalactic PN samples in number of candidates, magnitude depth, and radial extent, which allows us to carry out an unprecedented photometric study of the PN population in the outer regions of M87. We find that the logarithmic density profile of the PN distribution is shallower than the surface brightness profile at large radii. This behaviour is consistent with the superposition of two components associated with the halo of M87 and with the ICL, which have different luminosity specific PN numbers, the ICL contributing three times more PNs per unit light. Because of the depth of this survey we are also able to study the shape of the PN luminosity function (PNLF) in the outer regions of M87. We find a slope for the PNLF that is steeper at fainter magnitudes than the standard analytical PNLF formula and adopt a generalised model that treats the slope as a free parameter. Comparing the PNLF of M87 and the M31 bulge, both normalised by the sampled luminosity, the M87 PNLF contains fewer bright PNs and has a steeper slope towards fainter magnitudes., Comment: 16 pages, 13 figures, 5 tables, accepted for publication in A&A

The stellar kinematics of the spheroids and discs of S0 galaxies contain clues to their formation histories. Unfortunately, it is difficult to disentangle the two components and to recover their stellar kinematics in the faint outer parts of the galaxies using conventional absorption line spectroscopy. This paper therefore presents the stellar kinematics of six S0 galaxies derived from observations of planetary nebulae (PNe), obtained using the Planetary Nebula Spectrograph. To separate the kinematics of the two components, we use a maximum-likelihood method that combines the discrete kinematic data with a photometric component decomposition. The results of this analysis reveal that: the discs of S0 galaxies are rotationally supported; however, the amount of random motion in these discs is systematically higher than in comparable spiral galaxies; and the S0s lie around one magnitude below the Tully--Fisher relation for spiral galaxies, while their spheroids lie nearly one magnitude above the Faber--Jackson relation for ellipticals. All of these findings are consistent with a scenario in which spirals are converted into S0s through a process of mild harassment or "pestering," with their discs somewhat heated and their spheroid somewhat enhanced by the conversion process. In such a scenario, one might expect the properties of S0s to depend on environment. We do not see such an effect in this fairly small sample, although any differences would be diluted by the fact that the current location does not necessarily reflect the environment in which the transformation occurred. Similar observations of larger samples probing a broader range of environments, coupled with more detailed modelling of the transformation process to match the wide range of parameters that we have shown can now be measured, should take us from these first steps to the definitive answer as to how S0 galaxies form., Comment: 12 pages, 8 figures, accepted for publication in MNRAS

The origins of S0 galaxies remain obscure, with various mechanisms proposed for their formation, likely depending on environment. These mechanisms would imprint different signatures in the galaxies' stellar kinematics out to large radii, offering a method for distinguishing between them. We aim to study a sample of six S0 galaxies from a range of environments, and use planetary nebulae (PNe) as tracers of their stellar populations out to very large radii, to determine their kinematics in order to understand their origins. Using a special-purpose instrument, the Planetary Nebula Spectrograph, we observe and extract PNe catalogues for these six systems*. We show that the PNe have the same spatial distribution as the starlight, that the numbers of them are consistent with what would be expected in a comparable old stellar population in elliptical galaxies, and that their kinematics join smoothly onto those derived at smaller radii from conventional spectroscopy. The high-quality kinematic observations presented here form an excellent set for studying the detailed kinematics of S0 galaxies, in order to unravel their formation histories. We find that PNe are good tracers of stellar kinematics in these systems. We show that the recovered kinematics are largely dominated by rotational motion, although with significant random velocities in most cases., Comment: 13, 8 figures, accepted in A&A

(Abridged) The Galactic Center (GC) hosts a population of young stars some of which seem to form mutually inclined discs of clockwise and counter clockwise rotating stars. We present a warped disc origin scenario for these stars assuming that an initially flat accretion disc becomes warped due to the Pringle instability, or due to Bardeen-Petterson effect, before it fragments to stars. We show that this is plausible if the star formation efficiency $\epsilon_{SF} \lesssim 1$, and the viscosity parameter $\alpha \sim 0.1$. After fragmentation, we model the disc as a collection of concentric, circular, mutually tilted rings, and construct warped disc models for mass ratios and other parameters relevant to the GC environment, but also for more massive discs. We take into account the disc's self-gravity and the torques exerted by a surrounding star cluster. We show that a self-gravitating low-mass disc ($M_d / M_{bh} \sim 0.001$) precesses in integrity in the life-time of the stars, but precesses freely when the torques from a non-spherical cluster are included. An intermediate-mass disc ($M_d / M_{bh} \sim 0.01$) breaks into pieces which precess independently in the self-gravity-only case, and become disrupted in the presence of the star cluster torques. For a high-mass disc ($M_d / M_{bh} \sim 0.1$) the evolution is dominated by self-gravity and the disc is broken but not dissolved. The time-scale after which the disc breaks scales almost linearly with ($M_d / M_{bh}$) for self-gravitating models. Typical values are longer than the age of the stars for a low mass disc, and are in the range $\sim 8 \times 10^4-10^5$ yr for high and intermediate-mass discs respectively. None of these models explain the rotation properties of the two GC discs, but a comparison of them with the clockwise disc shows that the lowest mass model in a spherical star cluster matches the data best., Comment: 16 pages, 19 figures, abstract abridged to meet arXiv requirements. Accepted for publication in MNRAS

We have compared the halo kinematics traced by globular clusters (GCs) and planetary nebulae (PNe) for two elliptical galaxies in the Fornax and Virgo clusters NGC 1399 and NGC 4649, and for the merger remnant NGC 5128 (Centaurus A). We find differences in the rotational properties of the PN, red GC, and blue GC systems in all these three galaxies. NGC 1399 PNe and GCs show line of sight velocity distributions in specific regions that are significantly different, based on Kolmogorov-Smirnov tests. The PN system shows multi-spin components, with nearly opposite direction of rotation in the inner and the outer parts. The GCs velocity field is not point-symmetric in the outer regions of the galaxy, indicating that the system has not reached dynamical equilibrium yet. In NGC 4649 PNe, red and blue GCs have different rotation axes and rotational velocities. Finally, in NGC 5128 both PNe and GCs deviate from equilibrium in the outer regions of the galaxy, and in the inner regions the PN system is rotationally supported, whereas the GC system is dominated by velocity dispersion. The observed different kinematic properties, including deviations from point-symmetry, between PNe and GCs suggest that these systems are accreted at different times by the host galaxy, and the most recent accretion took place only few Gyr ago.We discuss two scenarios which may explain some of these differences: i) tidal stripping of loosely-bound GCs, and ii) multiple accretion of low luminosity and dwarf galaxies. Because these two mechanisms affect mostly the GC system, differences with the PNe kinematics can be expected., Comment: 14 pages, 13 figures, Accepted for publication in MNRAS. This new version contains an improved analysis, which includes the study of point-symmetry in the velocity fields and its implications for dynamical equilibrium

Adaptive optics observations of the flattened nuclear star cluster in the nearby edge-on spiral galaxy NGC 4244 using the Gemini Near-Infrared Integral Field Spectrograph (NIFS) have revealed clear rotation. Using these kinematics plus 2MASS photometry we construct a series of axisymmetric two-component particle dynamical models with our improved version of NMAGIC, a flexible Chi^2-made-to-measure code. The models consist of a nuclear cluster disc embedded within a spheroidal particle population. We find a mass for the nuclear star cluster of M=1.6^+0.5_-0.2 x 10^7 M_sun within ~42.4 pc (2"). We also explore the presence of an intermediate mass black hole and show that models with a black hole as massive as M_bh = 5.0 x 10^5 M_sun are consistent with the available data. Regardless of whether a black hole is present or not, the nuclear cluster is vertically anisotropic (beta_z < 0), as was found with earlier two-integral models. We then use the models as initial conditions for N-body simulations. These simulations show that the nuclear star cluster is stable against non-axisymmetric perturbations. We also explore the effect of the nuclear cluster accreting star clusters at various inclinations. Accretion of a star cluster with mass 13% that of the nuclear cluster is already enough to destroy the vertical anisotropy, regardless of orbital inclination., Comment: Replaced with the version accepted for publication in MNRAS. 13 pages, 10 figures, 3 tables

The accuracy and robustness of a simple method to estimate the total mass profile of a galaxy is tested using a sample of 65 cosmological zoom-simulations of individual galaxies. The method only requires information on the optical surface brightness and the projected velocity dispersion profiles and therefore can be applied even in case of poor observational data. In the simulated sample massive galaxies ($\sigma \simeq 200-400$ $\kms$) at redshift $z=0$ have almost isothermal rotation curves for broad range of radii (RMS $\simeq 5%$ for the circular speed deviations from a constant value over $0.5R_{\rm eff} < r < 3R_{\rm eff}$). For such galaxies the method recovers the unbiased value of the circular speed. The sample averaged deviation from the true circular speed is less than $\sim 1%$ with the scatter of $\simeq 5-8%$ (RMS) up to $R \simeq 5R_{\rm eff}$. Circular speed estimates of massive non-rotating simulated galaxies at higher redshifts ($z=1$ and $z=2$) are also almost unbiased and with the same scatter. For the least massive galaxies in the sample ($\sigma < 150$ $\kms$) at $z=0$ the RMS deviation is $\simeq 7-9%$ and the mean deviation is biased low by about $1-2%$. We also derive the circular velocity profile from the hydrostatic equilibrium (HE) equation for hot gas in the simulated galaxies. The accuracy of this estimate is about RMS $\simeq 4-5%$ for massive objects ($M > 6.5\times 10^{12} M_\odot$) and the HE estimate is biased low by $\simeq 3-4%$, which can be traced to the presence of gas motions. This implies that the simple mass estimate can be used to determine the mass of observed massive elliptical galaxies to an accuracy of $5-8 %$ and can be very useful for galaxy surveys., Comment: 15 pages, 14 figures, 1 table

Context. We investigate the stellar population and the origin of diffuse light around brightest cluster galaxies. Aims. We study the stellar population of the dynamically hot stellar halo of NGC 3311, the brightest galaxy in the Hydra I cluster, and that of photometric substructures in the diffuse light to constrain the origin of these components. Methods. We analyze absorption lines in medium-resolution, long-slit spectra in the wavelength range 4800-5800 angstrom obtained with FORS2 at the Very Large Telescope. We measure the equivalent width of Lick indices out to 20 kpc from the center of NGC 3311 and fit them with stellar population models that account for the [alpha/Fe] overabundance. Results. Stars in the dynamically hot halo of NGC 3311 are old (age > 13 Gyr), metal-poor ([Z/H] ~ -0.35), and alpha-enhanced ([alpha/Fe] ~ 0.48). Together with the high velocity dispersion, these measurements indicate that the stars in the halo were accreted from the outskirts of other early-type galaxies, with a possible contribution from dwarf galaxies. We identify a region in the halo of NGC 3311 associated with a photometric substructure where the stellar population is even more metal-poor ([Z/H] ~ -0.73). In this region, our measurements are consistent with a composite stellar population superposed along the line of sight, consisting of stars from the dynamically hot halo of NGC 3311 and stars stripped from dwarf galaxies. The latter component contributes < 28% to the local surface brightness. Conclusions. The build-up of diffuse light around NGC 3311 is on-going. Based on the observed stellar population properties, the dominant part of these stars may have come from the outskirts of bright early-type galaxies, while stars from stripped dwarf galaxies are presently being added., Comment: 8 pages, 4 figures. Accepted for publication in Astronomy & Astrophysics

Using a slitless spectroscopy method with (a) the 8.2 m Subaru telescope and its FOCAS Cassegrain spectrograph, and (b) the ESO Very Large Telescope (VLT) unit 1 (Antu) and its FORS2 Cassegrain spectrograph, we have detected 326 planetary nebulae (PNs) in the giant Virgo elliptical galaxy NGC 4649 (M 60), and we have measured their radial velocities. After rejecting some PNs more likely to belong to the companion galaxy NGC 4647, we have built a catalog with kinematic information for 298 PNs in M 60. Using these radial velocities we have concluded that they support the presence of a dark matter halo around M 60. Based on an isotropic, two-component Hernquist model, we estimate the dark matter halo mass within 3$R_{\rm e}$ to be 4$\times10^{11} M_{\odot}$, which is almost one half of the total mass of about $10^{12} M_{\odot}$ within 3$R_{\rm e}$. This total mass is similar to that estimated from globular cluster, XMM-Newton and Chandra observations. The dark matter becomes dominant outside. More detailed dynamical modeling of the PN data is being published in a companion paper. We have also measured the $m$(5007) magnitudes of many of these PNs, and built a statistically complete sample of 218 PNs. The resulting PN luminosity function (PNLF) was used to estimate a distance modulus of 30.7$\pm$0.2 mag, equivalent to 14$\pm$1 Mpc. This confirms an earlier PNLF distance measurement, based on a much smaller sample. The PNLF distance modulus remains smaller than the surface brightness fluctuation (SBF) distance modulus by 0.4 mag., Comment: Accepted for publication in ApJ; 14 figures

We present the third data release of the RAdial Velocity Experiment (RAVE) which is the first milestone of the RAVE project, releasing the full pilot survey. The catalog contains 83,072 radial velocity measurements for 77,461 stars in the southern celestial hemisphere, as well as stellar parameters for 39,833 stars. This paper describes the content of the new release, the new processing pipeline, as well as an updated calibration for the metallicity based upon the observation of additional standard stars. Spectra will be made available in a future release. The data release can be accessed via the RAVE webpage: http://www.rave-survey.org., Comment: AJ accepted. 54 pages, 20 figures. Figure 17 in low resolution mode

The study of the chemical abundances of HII regions in polar ring galaxies and their implications for the evolutionary scenario of these systems has been a step forward both in tracing the formation history of the galaxy and giving hints on the mechanisms at work during the building of disk by cold accretion process. It's now important to establish whether such results are typical for the class of polar disk galaxies as whole. The present work aims at checking the cold accretion of gas through a "cosmic filament" as a possible scenario for the formation of the polar structures in UGC7576 and UGC9796. If these form by cold accretion, we expect the HII regions abundances and metallicities to be lower than those of same-luminosity spiral disks, with values of the order of Z ~ 1/10 Zsun, as predicted by cosmological simulations. We have used deep long-slit spectra, obtained with DOLORES@TNG in the optical wavelengths, of the brightest HII regions associated with the polar structures to derive their chemical abundances and star formation rate. We used the "Empirical methods", based on the intensities of easily observable lines, to derive the oxygen abundance 12+log(O/H) of both galaxies. Such values are compared with those typical for different morphological galaxy types of comparable luminosity. The average metallicity values for UGC7576 and UGC9796 are Z = 0.4 Zsun and Z = 0.1 Zsun respectively. Both values are lower than those measured for ordinary spirals of similar luminosity and UGC7576 presents no metallicity gradient along the polar structure. These data, toghether with other observed features, available for the two PRGs in previous works, are compared with the predictions of simulations of tidal accretion, cold accretion and merging, to disentangle between these scenarios., Comment: 12 pages, 8 figures, accepted for publication on A&A

To investigate the origins of S0 galaxies, we present a new method of analyzing their stellar kinematics from discrete tracers such as planetary nebulae. This method involves binning the data in the radial direction so as to extract the most general possible non-parametric kinematic profiles, and using a maximum likelihood fit within each bin in order to make full use of the information in the discrete kinematic tracers. Both disk and spheroid kinematic components are fitted, with a two-dimensional decomposition of imaging data used to attribute to each tracer a probability of membership in the separate components. Likelihood clipping also allows us to identify objects whose properties are not consistent with the adopted model, rendering the technique robust against contaminants and able to identify additional kinematic features. The method is first tested on an N-body simulated galaxy to assess possible sources of systematic error associated with the structural and kinematic decomposition, which are found to be small. It is then applied to the S0 system NGC~1023, for which a planetary nebula catalogue has already been released and analyzed by (Noordermeer et al., 2008). The correct inclusion of the spheroidal component allows us to show that, contrary to previous claims, the stellar kinematics of this galaxy are indistinguishable from those of a normal spiral galaxy, indicating that it may have evolved directly from such a system via gas stripping or secular evolution. The method also successfully identifies a population of outliers whose kinematics are different from those of the main galaxy; these objects can be identified with a stellar stream associated with the companion galaxy NGC~1023A., Comment: 11 pages, 11 figures, accepted for publication in MNRAS