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.

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.