The characterisation of the extinction curve in the near infrared (NIR) is fundamental to analyse the structure and stellar population of the Galactic centre (GC), whose analysis is hampered by the extreme interstellar extinction ($A_V\sim 30$ mag) that varies on arc-second scales. Recent studies indicate that the behaviour of the extinction curve might be more complex than previously assumed, pointing towards a variation of the extinction curve as a function of wavelength. We aim at analysing the variations of the extinction index, $\alpha$, with wavelength, line-of-sight, and absolute extinction, extending previous analysis to a larger area of the innermost regions of the Galaxy. We analysed the whole GALACTICNUCLEUS survey, a high-angular resolution ($\sim 0.2''$) $JHK_s$ NIR survey specially designed to observe the GC in unprecedented detail. It covers a region of $\sim 6000$\,pc$^2$, comprising fields in the nuclear stellar disc, the inner bulge, and the transition region between them. We applied two independent methods based on red clump (RC) stars to constrain the extinction curve and analysed its variation superseding previous studies. We used more than 165,000 RC stars and increased significantly the size of the regions analysed to confirm that the extinction curve varies with the wavelength. We estimated a difference $\Delta\alpha = 0.21\pm0.07$ between the obtained extinction indices, $\alpha_{JH}=2.44\pm0.05$ and $\alpha_{HK_s} = 2.23\pm0.05$. We also concluded that there is no significant variation of the extinction curve with wavelength, with the line-of-sight or the absolute extinction. Finally, we computed the ratios between extinctions, $A_J/A_H = 1.87\pm0.03$ and $A_{H}/A_{K_s} = 1.84\pm0.03$, consistent with all the regions of the GALACTICNUCLEUS catalogue.

We have started an ambitious program to determine if the full diversity of extinction laws is real or if some of it is due to calibration or methodological issues. Here we start by analyzing the information on NIR extinction in a 2MASS stellar sample with good quality photometry and very red colours. We calculate the extinction at 1 $\mu$m, $A_1$, and the power-law exponent, $\alpha$ ($A_\lambda = A_1 \lambda^{-\alpha}$), for the 2MASS stars located in the extinction trajectory in the $H-K$ vs. $J-H$ plane expected for red giants with $A_1 > 5$ mag. We test the validity of the assumption about the nature of those stars, whether a single or multiple values of $\alpha$ are needed, and the spatial variations of the results. Most ($\sim$83%) of those stars can indeed be explained by high-extinction red giants and the rest is composed of extinguished AGB stars (mostly O-rich), blended sources, and smaller numbers of other objects, a contaminant fraction that can be reduced with the help of Gaia} DR2 data. Galactic red giants experience a NIR extinction with $\alpha\sim 2.27$ and an uncertainty of a few hundredths of a magnitude. There is no significant spread in $\alpha$ even though our sample is widely distributed and has a broad range of extinctions. Differences with previous results are ascribed to the treatment of non-linear photometric effects and/or the contaminant correction. Future research should concentrate in finding the correct functional form for the NIR extinction law. In the appendix we detail the treatment of non-linear photometric effects in the 2MASS bands.

We present stellar metallicity measurements of more than 600 late-type stars in the central 10 pc of the Galactic centre. Together with our previously published KMOS data, this data set allows us to investigate, for the first time, spatial variations of the nuclear star cluster's metallicity distribution. Using the integral-field spectrograph KMOS (VLT) we observed almost half of the area enclosed by the nuclear star cluster's effective radius. We extract spectra at medium spectral resolution, and apply full spectral fitting utilising the PHOENIX library of synthetic stellar spectra. The stellar metallicities range from [M/H]=-1.25 dex to [M/H]> +0.3 dex, with most of the stars having super-solar metallicity. We are able to measure an anisotropy of the stellar metallicity distribution. In the Galactic North, the portion of sub-solar metallicity stars with [M/H]<0.0 dex is more than twice as high as in the Galactic South. One possible explanation for different fractions of sub-solar metallicity stars in different parts of the cluster is a recent merger event. We propose to test this hypothesis with high-resolution spectroscopy, and by combining the metallicity information with kinematic data.

The Milky Way nuclear star cluster (MWNSC) is a crucial laboratory for studying the galactic nuclei of other galaxies, but its properties have not been determined unambiguously until now. Aims. We aim to study the size and spatial structure of the MWNSC. Methods. This study uses data and methods that address potential shortcomings of previous studies on the topic. We use $0.2''$ angular resolution $K_{s}$ data to create a stellar density map in the central $86.4$pc x $21$pc at the Galactic center. We include data from selected adaptive-optics-assisted images obtained for the inner parsecs. In addition, we use IRAC/Spitzer mid-infrared (MIR) images. We model the Galactic bulge and the nuclear stellar disk in order to subtract them from the MWNSC. Finally, we fit a S\'ersic model to the MWNSC and investigate its symmetry. Results. Our results are consistent with previous work. The MWNSC is flattened with an axis ratio of $q=0.71 \pm0.10$, an effective radius of $R_{e}=(5.1 \pm1.0)$ pc, and a S\'ersic index of $n=2.2 \pm0.7$. Its major axis may be tilted out of the Galactic plane by up to $-10$ degrees. The distribution of the giants brighter than the Red Clump (RC) is found to be significantly flatter than the distribution of the faint stars. We investigate the 3D structure of the central stellar cusp using our results on the MWNSC structure on large scales to constrain the deprojection of the measured stellar surface number density, obtaining a value of the 3D inner power law of $\gamma=1.38\pm0.06_{sys}\pm0.01_{stat}$. Conclusions. The MWNSC shares its main properties with other extragalactic NSCs found in spiral galaxies. The differences in the structure between bright giants and RC stars might be related to the existence of not completely mixed populations of different ages. This may hint at recent growth of the MWNSC through star formation or cluster accretion.

Proper motion studies of stars in the centre of the Milky Way have been typically limited to the Arches and Quintuplet clusters and to the central parsec. Here, we present the first results of a large-scale proper motion study of stars within several tens of parsecs of Sagittarius A* based on our $0.2''$ angular resolution GALACTICNUCLEUS survey (epoch 2015) combined with NICMOS/HST data from the Paschen-$\alpha$ survey (epoch 2008). This study will be the first extensive proper motion study of the central $\sim 36' \times 16'$ of the Galaxy, which is not covered adequately by any of the existing astronomical surveys such as Gaia because of its extreme interstellar extinction ($A_{V} \gtrsim 30$ mag). Proper motions can help us to disentangle the different stellar populations along the line-of-sight and interpret their properties in combination with multi-wavelength photometry from GALACTICNUCLEUS and other sources. It also allows us to infer the dynamics and interrelationship between the different stellar components of the Galactic Centre (GC). In particular, we use proper motions to detect co-moving groups of stars which can trace low mass or partially dissolved young clusters in the GC that can hardly be discovered by any other means. Our pilot study in this work is on a field in the nuclear bulge associated by HII regions that show the presence of young stars. We detect the first group of co-moving stars coincident with an HII region. Using colour-magnitude diagrams, we infer that the co-moving stars are consistent with being the post-main sequence stars with ages of few Myrs. Simulations show that this group of stars is a real group that can indicate the existence of a dissolving or low to intermediate mass young cluster. A census of these undiscovered clusters will ultimately help us to constrain star formation at the GC in the past few ten Myrs.

The nuclear disc is a dense stellar structure at the centre of the Milky Way, with a radius of $\sim$150 pc. It has been a place of intense star formation in the past several tens of millions of years but its overall formation history has remained unknown up to now. Here we report the first detailed star formation history of this region. The bulk of its stars formed at least eight billion years ago. This initial activity was followed by a long period of quiescence that was ended by an outstanding event about 1 Gyr ago, during which roughly 5% of its mass formed in a time window $\sim$100 Myr, in what may arguably have been one of the most energetic events in the history of the Milky Way. Star formation continued subsequently on a lower level, creating a few percent of the stellar mass in the past $\sim$500 Myr, with an increased rate up to $\sim$30 Myr ago. Our findings contradict the previously accepted paradigm of quasi-continuous star formation at the centre of the Milky Way. The long quiescent phase agrees with the overall quiescent history of the Milky Way and suggests that our Galaxy's bar may not have existed until recently, or that gas transport through the bar was extremely inefficient during a long stretch of the Milky Way's life, and that the central black hole may have acquired most of its mass already in the early days of the Milky Way.

Due to the extreme extinction towards the Galactic centre ($A_{V} \sim 30$ mag), its stellar population is mainly studied in the near-infrared (NIR) regime. Therefore, a proper analysis of the NIR extinction curve is necessary to fully characterise the stellar structure and population of the inner part of the galaxy. We studied the dependence of the extinction index ($\alpha_\lambda$) in the NIR on the line of sight, wavelength, and extinction. We used the GALACTICNUCLEUS imaging survey, a high angular resolution catalogue ($0.2''$) for the inner part of the Galaxy in $JHK_s$, and studied the spatial variation in the extinction index. We also applied two independent methods based on red clump stars to compute the extinction index between different bands and its variation with wavelength. We did not detect any significant line-of-sight or extinction variation in $\alpha$ within the studied region in the nuclear stellar disc. The extinction index between $JH$ and $HK_s$ differs by $0.19 \pm 0.05$. We obtained mean values for the extinction indices $\alpha_{JH} = 2.43\pm0.03$ and $\alpha_{HK_s} = 2.23\pm0.03$. The dependence of the extinction index on the wavelength could explain the differences obtained for $\alpha_\lambda$ in the literature since it was assumed constant for the NIR regime.

The high extinction and extreme source crowding of the central regions of the Milky Way are serious obstacles to the study of the structure and stellar population of the Galactic centre (GC). Existing surveys that cover the GC region (2MASS, UKIDSS, VVV, SIRIUS) do not have the necessary high angular resolution. Therefore, a high-angular-resolution survey in the near infrared is crucial to improve the state of the art. Here, we present the GALACTICNUCLEUS catalogue, a near infrared $JHK_s$ high-angular-resolution ($0.2''$) survey of the nuclear bulge of the Milky Way. We explain in detail the data reduction, data analysis, calibration, and uncertainty estimation of the GALACTICNUCLEUS survey. We assess the data quality comparing our results with previous surveys. We obtained accurate $JHK_s$ photometry for $\sim 3.3\times10^6$ stars in the GC detecting around 20 \% in $J$, 65 \% in $H,$ and 90 \% in $K_s$. The survey covers a total area of $\sim0.3$ square degrees, which corresponds to $\sim 6,000$ pc$^2$. The GALACTICNUCLEUS survey reaches 5\,$\sigma$ detections for $J \sim 22$ mag, $H \sim 21$ mag, and $K_s \sim 21$ mag. The uncertainties are below 0.05 mag at $J \sim 21$ mag, $H \sim 19$ mag, and $K_s \sim 18$ mag. The zero point systematic uncertainty is $\lesssim0.04$ mag in all three bands. We present colour-magnitude diagrams for the different regions covered by the survey.

The study of the inner region of the Milky Way's bulge is hampered by high interstellar extinction and extreme source crowding. Sensitive high angular resolution near-infrared imaging is needed to study stellar populations in such a complex environment. We use the 0.2$''$ angular resolution $JHK_s$ data from the GALACTICNUCLEUS survey to study the stellar population within two $8.0'\times 3.4'$ fields, about 0.6$^\circ$ and 0.4$^\circ$ to the Galactic north of the Milky Way's centre and to compare it with one in the immediate surroundings of Sagittarius A*. We also characterise the extinction curve of the two fields. The average interstellar extinction to the outer and the inner field is $A_{K_s} \sim 1.20 \pm 0.08$ mag and $\sim 1.48 \pm 0.10$ mag, respectively. We present $K_{s}$ luminosity functions that are complete down to at least 2 mag below the red clump (RC). We detect a feature in the luminosity functions that is fainter than the RC by $0.80\pm0.03$ and $0.79\pm0.02$ mag, respectively, in the $K_s$-band. It runs parallel to the reddening vector. We identify the feature as the red giant branch bump. Fitting $\alpha$-enhanced BaSTI luminosity functions to our data, we find that a single old stellar population of $\sim12.8 \pm 0.6$ Gyr and $Z = 0.040 \pm 0.003$ provides the best fit. We obtain that the stellar population in the innermost bulge is old, similar to the one at larger distances from the Galactic plane, and that its metallicity increases down to distances as short as about 60 pc from the centre of the Milky Way. Comparing it with previous known values at larger latitudes ($|b|>2^\circ$), our results favour a flattening of the gradient at $|b|<2^\circ$. As a secondary result we obtain that the extinction index in the studied regions agrees the value of $\alpha = 2.30\pm0.08$, derived in Nogueras-Lara et al. 2018 for the very Galactic centre.

(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.

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.

The Galactic Centre is of fundamental astrophysical interest, but existing near-infrared surveys fall short to cover it adequately. Here we introduce the GALACTICNUCLEUS survey, a JHKs imaging survey of the centre of the Milky Way with a 0.2" angular resolution. We present the observations of Field 1 of our survey, centred approximately on SgrA* with an approximate size of 7.95 ' x 3.43 '. We describe the observational set-up and data reduction pipeline and discuss the quality of the data. Finally, we present some preliminary analysis of the data. The data were acquired with the near-infrared camera HAWK-I at the ESO VLT. Short readout times in combination with the speckle holography algorithm allowed us to produce final images with a stable, Gaussian PSF of 0.2" FWHM. Astrometric calibration is achieved via the VVV survey and photometric calibration is based on the SIRIUS/IRSF survey. The quality of the data is assessed by comparison between observations of the same field with different detectors of HAWK-I and at different times. We reach 5 sigma detection limits of approximately J = 22, H = 21, and Ks = 20. The photometric uncertainties are less than 0.05 at J < 20, H < 17 and Ks < 16. We distinguish five stellar populations in the colour-magnitude diagrams; three of them appear to belong to foreground spiral arms, and the other two correspond to a high- and a low-extinction star groups at the Galactic Centre. We use our data to analyse the near-infrared extinction curve and conclude that it can be described very well by a power-law with an index of JHKs = 2.31 +- 0.03. We do not find any evidence that this index depends on the position along the line-of-sight, or on the absolute value of the extinction. We produce extinction maps that show the clumpiness of the ISM at the Galactic Centre. Finally, we estimate that the majority of the stars have solar or super-solar metallicity.

Because of strong and spatially highly variable interstellar extinction and extreme source crowding, the faint (K>15) stellar population in the Milky Way's nuclear cluster is still poorly studied. RR Lyrae stars provide us with a tool to estimate the mass of the oldest, relative dim stellar population. Recently, we analyzed HST/WFC3/IR observations of the central 2.3'x2.3' of the Milky Way and found 21 variable stars with periods between 0.2 and 1d. Here, we present a further comprehensive analysis of these stars. The period-luminosity relationship of RR Lyrae is used to derive their extinctions and distances. Using multiple approaches, we classify our sample as four RRc, four RRab and three candidates, ten binaries. Especially, the four RRabs show sawtooth light curves and fall exactly onto the Oosterhoff I division in the Bailey diagram. Compared to the RRabs reported by Minniti et al, 2016, our new RRabs have higher extinction (A_K>1.8) and should be closer to the Galactic Centre. The extinction and distance of one RRab match those for the nuclear star cluster given in previous works. We perform simulations and find that after correcting for incompleteness, there could be no more than 40 RRabs within the nuclear star cluster and in our field-of-view. Through comparing with the known globular clusters of the Milky Way, we estimate that if there exists an old, metal-poor (-1.5<[Fe/H]<-1) stellar population in the Milky Way nuclear star cluster on a scale of 5x5pc, then it contributes at most 4.7x10^5 solar mass, i.e., ~18% of the stellar mass.

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.

We present new [Ne II] (12.8 micron) IRTF/TEXES observations of the Galactic Center HII regions H1 and H2, which are at a projected distance of ~11 pc from the center of the Galaxy. The new observations allow to map the radial velocity distributions of ionized gas. The high spectroscopic resolution (~4 km/s) helps us to disentangle different velocity components and enables us to resolve previous ambiguity regarding the nature of these sources. The spatial distributions of the intensity and radial velocity of the [Ne II] line are mapped. In H1, the intensity distributions of the Paschen-\alpha (1.87 micron) and [Ne II] lines are significantly different, which suggests a strong variation of extinction across the HII region of A_K~0.56. The radial velocity distributions across these HII regions are consistent with the predictions of a bow-shock model for H1 and the pressure-driven model for H2. Furthermore, we find a concentration of bright stars in H2. These stars have similar H-K_s colors and can be explained as part of a 2 Myr old stellar cluster. H2 also falls on the orbit of the molecular clouds, suggested to be around Sgr A*. Our new results confirm what we had previously suggested: the O supergiant P114 in H1 is a runaway star, moving towards us through the -30-0 {km/s} molecular cloud, whereas the O If star P35 in H2 formed in-situ, and may mark the position of a so-far unknown small star cluster formed within the central 30 pc of the Galaxy.

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.

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.