Detecting hot stars in the Galactic centre with combined near- and mid-infrared photometry

{The Galactic centre (GC) is a unique astrophysical laboratory to study the stellar population of galactic nuclei because it is the only galactic nucleus whose stars can be resolved down to milliparsec scales. However, the extreme and spatially highly variable interstellar extinction towards the GC poses a serious obstacle to photometric stellar classification.} {Our goal is to identify hot, massive stars in the nuclear stellar disc (NSD) region through combining near-infrared (NIR) and mid-infrared (MIR) photometry, and thus to demonstrate the feasibility of this technique, which may gain great importance with the arrival of the James Webb Space Telescope (JWST).} {We combined the GALACTICNUCLEUS NIR survey with the IRAC/Spitzer MIR survey of the GC. We applied the so-called Rayleigh-Jeans colour excess (RJCE) de-reddening method to our combined NIR-MIR data to identify potential hot stars in colour-magnitude diagrams (CMDs).} {Despite the very low angular resolution of IRAC we find 12 clear candidates for young massive stars among the $1\,065$ sources that meet our selection criteria. Seven out of these 12 stars are previously known hot, massive stars belonging to the Arches and Quintuplet clusters, as well as sources detected by the Hubble Space Telescope/NICMOS Paschen-$\alpha$ survey. Five of our massive star candidates have not been previously reported in the literature.} {We show that the RJCE method is a valuable tool to identify hot stars in the GC using photometry alone. Upcoming instruments with high angular resolution MIR imaging capabilities such as the JWST could surely make more substantial use of this de-reddening method and help establish a far more complete census of hot, young stars in the GC area than what is possible at the moment.}Comment: 6 pages, 6 figures, accepted for publication in A&

Compact Radio Sources within 30" of Sgr A*: Proper Motions, Stellar Winds and the Accretion Rate onto Sgr A*

Recent broad-band 34 and 44 GHz radio continuum observations of the Galactic center have revealed 41 massive stars identified with near-IR counterparts, as well as 44 proplyd candidates within 30" of Sgr A*. Radio observations obtained in 2011 and 2014 have been used to derive proper motions of eight young stars near Sgr A*. The accuracy of proper motion estimates based on near-IR observations by Lu et al. and Paumard et al. have been investigated by using their proper motions to predict the 2014 epoch positions of near-IR stars and comparing the predicted positions with those of radio counterparts in the 2014 radio observations. Predicted positions from Lu et al. show an rms scatter of 6 mas relative to the radio positions, while those from Paumard et al. show rms residuals of 20 mas, which is mainly due to uncertainties in the IR-based proper motions. Under the assumption of homogeneous ionized winds, we also determine the mass-loss rates of 11 radio stars, finding rates that are on average $\sim$2 times smaller than those determined from model atmosphere calculations and near-IR data. Clumpiness of ionized winds would reduce the mass loss rate of WR and O stars by additional factors of 3 and 10, respectively. One important implication of this is a reduction in the expected mass accretion rate onto Sgr A* from stellar winds by nearly an order of magnitude to a value of few$\times10^{-7}$ \msol\ yr$^{-1}$. Finally, we present the positions of 318 compact 34.5 GHz radio sources within 30\arcs\ of Sgr A*. At least 45 of these have stellar counterparts in the near-IR $K_s$ (2.18 $\mu$m) and $L'$ (3.8$\mu$m) bands.Comment: 30 pages, 4 figures, ApJ (in press

The distribution of stars around the Milky Way's central black hole II: Diffuse light from sub-giants and dwarfs

This is the second of three papers that search for the predicted stellar cusp around the Milky Way's central black hole, Sagittarius A*, with new data and methods. We aim to infer the distribution of the faintest stellar population currently accessible through observations around Sagittarius A*. We use adaptive optics assisted high angular resolution images obtained with the NACO instrument at the ESO VLT. Through optimised PSF fitting we remove the light from all detected stars above a given magnitude limit. Subsequently we analyse the remaining, diffuse light density. The analysed diffuse light arises from sub-giant and main-sequence stars with KS ~ 19 - 20 with masses of 1 - 2 Msol . These stars can be old enough to be dynamically relaxed. The observed power-law profile and its slope are consistent with the existence of a relaxed stellar cusp around the Milky Way's central black hole. We find that a Nuker law provides an adequate description of the nuclear cluster's intrinsic shape (assuming spherical symmetry). The 3D power-law slope near Sgr A* is \gamma = 1.23 +- 0.05. At a distance of 0.01 pc from the black hole, we estimate a stellar mass density of 2.3 +- 0.3 x 10^7 Msol pc^-3 and a total enclosed stellar mass of 180 +- 20 Msol. These estimates assume a constant mass-to-light ratio and do not take stellar remnants into account. The fact that no cusp is observed for bright (Ks 16) giant stars at projected distances of roughly 0.1-0.3 pc implies that some mechanism has altered their appearance or distribution.Comment: Accepted for publication A&

The distribution of old stars around the Milky Way's central black hole I: Star counts

(abridged) In this paper we revisit the problem of inferring the innermost structure of the Milky Way's nuclear star cluster via star counts, to clarify whether it displays a core or a cusp around the central black hole. Through image stacking and improved PSF fitting we push the completeness limit about one magnitude deeper than in previous, comparable work. Contrary to previous work, we analyse the stellar density in well-defined magnitude ranges in order to be able to constrain stellar masses and ages. The RC and brighter giant stars display a core-like surface density profile within a projected radius R<0.3 pc of the central black hole, in agreement with previous studies, but show a cusp-like surface density distribution at larger R. The surface density of the fainter stars can be described well by a single power-law at R<2 pc. The cusp-like profile of the faint stars persists even if we take into account the possible contamination of stars in this brightness range by young pre-main sequence stars. The data are inconsistent with a core-profile for the faint stars.Finally, we show that a 3D Nuker law provides a very good description of the cluster structure. We conclude that the observed stellar density at the Galactic Centre, as it can be inferred with current instruments, is consistent with the existence of a stellar cusp around the Milky Way's central black hole, Sgr A*. This cusp is well developed inside the influence radius of about 3 pc of Sgr A* and can be described by a single three-dimensional power-law with an exponent gamma=1.23+-0.05. The apparent lack of RC stars and brighter giants at projected distances of R < 0.3 pc (R<8") of the massive black hole may indicate that some mechanism has altered their distribution or intrinsic luminosity.Comment: Accepted for publication A&

KMOS view of the Galactic Centre - II. Metallicity distribution of late-type stars

Knowing the metallicity distribution of stars in the Galactic Centre has important implications for the formation history of the Milky Way nuclear star cluster. However, this distribution is not well known, and is currently based on a small sample of fewer than 100 stars. We obtained near-infrared K-band spectra of more than 700 late-type stars in the central 4 pc^2 of the Milky Way nuclear star cluster with the integral-field spectrograph KMOS (VLT). We analyse the medium-resolution spectra using a full-spectral fitting method employing the G\"ottingen Spectral library of synthetic PHOENIX spectra. The derived stellar metallicities range from metal-rich [M/H]>+0.3 dex to metal-poor [M/H]<-1.0 dex, with a fraction of 5.2(^{+6.0}+{-3.1}) per cent metal-poor ([M/H]<-0.5 dex) stars. The metal-poor stars are distributed over the entire observed field. The origin of metal-poor stars remains unclear. They could originate from infalling globular clusters. For the metal-rich stellar population ([M/H]>0 dex) a globular cluster origin can be ruled out. As there is only a very low fraction of metal-poor stars in the central 4 pc^2 of the Galactic Centre, we believe that our data can discard a scenario in which the Milky Way nuclear star cluster is purely formed from infalling globular clusters.Comment: 18 pages, 9 Figures, accepted for publication in MNRA

Near-Infrared Variability Study of the Central 2.3 arcmin x 2.3 arcmin of the Galactic Centre I. Catalog of Variable Sources

We used four-year baseline HST/WFC3 IR observations of the Galactic Centre in the F153M band (1.53 micron) to identify variable stars in the central ~2.3'x2.3' field. We classified 3845 long-term (periods from months to years) and 76 short-term (periods of a few days or less) variables among a total sample of 33070 stars. For 36 of the latter ones, we also derived their periods (<3 days). Our catalog not only confirms bright long period variables and massive eclipsing binaries identified in previous works, but also contains many newly recognized dim variable stars. For example, we found \delta Scuti and RR Lyrae stars towards the Galactic Centre for the first time, as well as one BL Her star (period < 1.3 d). We cross-correlated our catalog with previous spectroscopic studies and found that 319 variables have well-defined stellar types, such as Wolf-Rayet, OB main sequence, supergiants and asymptotic giant branch stars. We used colours and magnitudes to infer the probable variable types for those stars without accurately measured periods or spectroscopic information. We conclude that the majority of unclassified variables could potentially be eclipsing/ellipsoidal binaries and Type II Cepheids. Our source catalog will be valuable for future studies aimed at constraining the distance, star formation history and massive binary fraction of the Milky Way nuclear star cluster.Comment: has been accepted to be published in MNRAS, 64 pages, 26 figures. The complete lists of table 3, 4, 8 and 9 will be published onlin

Jeans modelling of the Milky Way's nuclear stellar disc

The nuclear stellar disc (NSD) is a flattened stellar structure that dominates the gravitational potential of the Milky Way at Galactocentric radii 30≲R≲300pc⁠. In this paper, we construct axisymmetric Jeans dynamical models of the NSD based on previous photometric studies and we fit them to line-of-sight kinematic data of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and silicon monoxide (SiO) maser stars. We find that (i) the NSD mass is lower but consistent with the mass independently determined from photometry by Launhardt et al. Our fiducial model has a mass contained within spherical radius r=100pc of M(r 1. Observations and theoretical models of the star-forming molecular gas in the central molecular zone suggest that large vertical oscillations may be already imprinted at stellar birth. However, the finding σz/σR > 1 depends on a drop in the velocity dispersion in the innermost few tens of parsecs, on our assumption that the NSD is axisymmetric, and that the available (extinction corrected) stellar samples broadly trace the underlying light and mass distributions, all of which need to be established by future observations and/or modelling. (iii) We provide the most accurate rotation curve to date for the innermost 500pc of our Galaxy