Context. The Galactic centre (GC) is of fundamental astrophysical interest, but existing near-infrared surveys fall short covering it adequately, either in terms of angular resolution, multi-wavelength coverage, or both. Here we introduce the GALACTICNUCLEUS survey, a JHKs imaging survey of the centre of the Milky Way with a 0.2″ angular resolution. Aims. The purpose of this paper is to present the observations of Field 1 of our survey, centred approximately on SgrA* with an approximate size of 7.95′ × 3.43′. We describe the observational set-up and data reduction pipeline and discuss the quality of the data. Finally, we present the analysis of the data. Methods. The data were acquired with the near-infrared camera High Acuity Wide field K-band Imager (HAWK-I) at the ESO Very Large Telescope (VLT). Short readout times in combination with the speckle holography algorithm allowed us to produce final images with a stable, Gaussian PSF (point spread function) of 0.2″ FWHM (full width at half maximum). Astrometric calibration is achieved via the VISTA Variables in the Via Lactea (VVV) survey and photometric calibration is based on the SIRIUS/Infrared Survey Facility telescope (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. Results. We reach 5σ 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 can 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 high- and low-extinction star groups at the GC. We use our data to analyse the near-infrared extinction curve and find some evidence for a possible difference between the extinction index between J − H and H − Ks. However, we conclude that it can be described very well by a power law with an index of αJHKs = 2.30 ± 0.08. 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 (interstellar medium) at the GC. Finally, we estimate that the majority of the stars have solar or super-solar metallicity by comparing our extinction-corrected colour-magnitude diagrams with isochrones with different metallicities and a synthetic stellar model with a constant star formation.

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öttingen 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.

No abstract

Context. 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. Aims. We aim to infer the distribution of the faintest stellar population currently accessible through observations around Sagittarius A*. Methods. We used adaptive optics assisted high angular resolution images obtained with the NACO instrument at the ESO VLT. Through optimised PSF fitting we removed the light from all detected stars above a given magnitude limit. Subsequently we analysed the remaining, diffuse light density. Systematic uncertainties were constrained by the use of data from different observing epochs and obtained with different filters. We show that it is necessary to correct for the diffuse emission from the mini-spiral, which would otherwise lead to a systematically biased light density profile. We used a Paschen α map obtained with the Hubble Space Telescope for this purpose. Results. The azimuthally averaged diffuse surface light density profile within a projected distance of R ≲ 0.5 pc from Sagittarius A* can be described consistently by a single power law with an exponent of Γ = 0.26 ± 0.02stat ± 0.05sys, similar to what has been found for the surface number density of faint stars in Paper I. Conclusions. The analysed diffuse light arises from sub-giant and main-sequence stars with Ks ≈ 19−22 with masses of 0.8−1.5 M⊙. 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 γ = 1.13 ± 0.03model ± 0.05sys. The stellar density decreases more steeply beyond a break radius of about 3 pc, which corresponds roughly to the radius of influence of the massive black hole. At a distance of 0.01 pc from the black hole, we estimate a stellar mass density of 2.6 ± 0.3 × 107 M⊙ pc-3 and a total enclosed stellar mass of 180 ± 30  M⊙. These estimates assume a constant mass-to-light ratio and do not take stellar remnants into account. The fact that a flat projected surface density is observed for old giants at projected distances R ≲ 0.3 pc implies that some mechanism may have altered their appearance or distribution.

Context. The existence of dynamically relaxed stellar density cusp in dense clusters around massive black holes is a long-standing prediction of stellar dynamics, but it has so far escaped unambiguous observational confirmation. Aims. In this paper we aim to 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. Methods. We used judiciously selected adaptive optics assisted high angular resolution images obtained with the NACO instrument at the ESO VLT. Through image stacking and improved PSF fitting we pushed the completeness limit about one magnitude deeper than in previous, comparable work. Crowding and extinction corrections were derived and applied to the surface density estimates. Known young, and therefore dynamically not relaxed stars, are excluded from the analysis. 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. Results. We focus on giant stars, with observed magnitudes K = 12.5 − 16, and on stars with observed magnitudes K ≈ 18, which may have similar mean ages and masses than the former. The giants 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 their 3D density distribution is not inconsistent with a shallow cusp if we take into account the extent of the entire cluster, beyond the radius of influence of the central black hole. 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 good description of the cluster structure. Conclusions. We conclude that the observed density of the faintest stars detectable with reasonable completeness at the Galactic Centre, is consistent with the existence of a stellar cusp around the Milky Way's central black hole, Sagittarius A*. This cusp is well developed inside the influence radius of Sagittarius A* and can be described by a single three-dimensional power-law with an exponent γ = 1.43 ± 0.02 ± 0.1 sys . This corroborates existing conclusions from Nbody simulations performed in a companion paper. An important caveat is that the faint stars analysed here may be contaminated significantly by dynamically unrelaxed stars that formed about 100 Myr ago. The apparent lack of giants at projected distances of R 0.3 pc (R 8") of the massive black hole may indicate that some mechanism may have altered their distribution or intrinsic luminosity. We roughly estimate the number of possibly missing giants to about 100.

Context. The environment of Sagittarius A* (Sgr A*), the central black hole of the Milky Way, is the only place in the Universe where we can currently study the interaction between a nuclear star cluster and a massive black hole and infer the properties of a nuclear cluster from observations of individual stars. Aims. This work aims to explore the star formation history of the nuclear cluster and the structure of the innermost stellar cusp around Sgr A*. Methods. We combined and analysed multi epoch high quality AO observations. For the region close to Sgr A* we apply the speckle holography technique to the AO data and obtain images that are ≥50% complete down to Ks ≈ 19 within a projected radius of 5″ around Sgr A*. We used H-band images to derive extinction maps. Results. We provide Ks photometry for roughly 39 000 stars and H-band photometry for ∼11 000 stars within a field of about 40″ × 40″, centred on Sgr A*. In addition, we provide Ks photometry of ∼3000 stars in a very deep central field of 10″ × 10″, centred on Sgr A*. We find that the Ks luminosity function (KLF) is rather homogeneous within the studied field and does not show any significant changes as a function of distance from the central black hole on scales of a few 0.1 pc. By fitting theoretical luminosity functions to the KLF, we derive the star formation history of the nuclear star cluster. We find that about 80% of the original star formation took place 10 Gyr ago or longer, followed by a largely quiescent phase that lasted for more than 5 Gyr. We clearly detect the presence of intermediate-age stars of about 3 Gyr in age. This event makes up about 15% of the originally formed stellar mass of the cluster. A few percent of the stellar mass formed in the past few 100 Myr. Our results appear to be inconsistent with a quasi-continuous star formation history. The mean metallicity of the stars is consistent with being slightly super solar. The stellar density increases exponentially towards Sgr A* at all magnitudes between Ks = 15−19. We also show that the precise properties of the stellar cusp around Sgr A* are hard to determine because the star formation history suggests that the star counts can be significantly contaminated, at all magnitudes, by stars that are too young to be dynamically relaxed. We find that the probability of observing any young (non-millisecond) pulsar in a tight orbit around Sgr A* and beamed towards Earth is very low. We argue that typical globular clusters, such as they are observed in and around the Milky Way today, have probably not contributed to the nuclear cluster’s mass in any significant way. The nuclear cluster may have formed following major merger events in the early history of the Milky Way.

Context. 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. Aims. Here, we present the GALACTICNUCLEUS catalogue, a near infrared JHKs high-angular-resolution (0.2″) survey of the nuclear bulge of the Milky Way. Methods. 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. Results. We obtained accurate JHKs photometry for ∼3.3 × 106 stars in the GC detecting around 20% in J, 65% in H, and 90% in Ks. The survey covers a total area of ∼0.3 deg2, which corresponds to ∼6000 pc2. The GALACTICNUCLEUS survey reaches 5σ detections for J ∼ 22 mag, H ∼ 21 mag, and Ks ∼ 21 mag. The uncertainties are below 0.05 mag at J ∼ 21 mag, H ∼ 19 mag, and Ks ∼ 18 mag. The zero point systematic uncertainty is ≲0.04 mag in all three bands. We present colour–magnitude diagrams for the different regions covered by the survey.

Context. 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.4pc x 21pc 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érsic 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 ± 0.10, an effective radius of R e = (5.1 ± 1.0) pc, and a Sérsic index of n = 2.2 ± 0.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 γ = 1.38 ± 0.06 sys ± 0.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.