🔭 astrophysics

Chemical analysis of the Milky Way's Nuclear Star Cluster: Evidence for a metallicity gradient

By reanalyzing low-resolution KMOS spectra of 1,140 M giant stars with an improved line-list, this study establishes precise stellar parameters, identifies a dual metallicity population, and reveals a negative metallicity gradient in the Milky Way's Nuclear Star Cluster that supports an inside-out formation scenario.

Original authors: M. Schultheis, L. Serrano, B. Thorsbro, F. Nogueras-Lara, A. Feldmeier-Krause, G. Nandakumar, K. Fiteni, M. C. Sormani

Published 2026-01-28

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: M. Schultheis, L. Serrano, B. Thorsbro, F. Nogueras-Lara, A. Feldmeier-Krause, G. Nandakumar, K. Fiteni, M. C. Sormani

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the center of our Milky Way galaxy as a bustling, crowded city. At the very heart of this city lies a massive, dense neighborhood called the Nuclear Star Cluster (NSC). For a long time, astronomers have been trying to figure out how this neighborhood was built. Did it grow from the inside out, like a tree adding rings? Or did it form by massive star clusters crashing into the center and merging, like cars crashing into a pileup?

To answer this, a team of astronomers led by M. Schultheis went back to the drawing board with a specific set of "blueprints" (data) to get a clearer picture of the residents (stars) living there.

Here is what they found, explained simply:

1. The Problem: The "Fuzzy" Photos

The center of our galaxy is incredibly dusty and dark. To see the stars there, astronomers have to use infrared light (like night-vision goggles). However, the data they had from previous studies was a bit like looking at these stars through a foggy window. The "photos" (spectra) were low-resolution, making it hard to tell exactly how hot the stars were, how heavy they were, or what they were made of.

2. The Solution: Sharpening the Lens

The team took the existing data and applied a new, highly detailed recipe (a new "line-list" and model grid) specifically designed for the cool, giant stars that dominate this region.

The Analogy: Think of it like upgrading from a blurry, old map to a high-definition GPS. They also used a "cleaning" process to remove static and noise from the data before analyzing it.
The Check: They tested this new method against high-quality data from nearby stars (which they could see clearly) and confirmed that their new "GPS" was accurate. They could now measure the stars' properties with much higher precision.

3. The Discovery: A City with Two Distinct Neighborhoods

When they analyzed 1,140 stars in this cluster, they found the population wasn't just one uniform group. It was a mix of two distinct types of "residents":

The Rich Neighbors: The majority of the stars are "metal-rich." In astronomy, "metals" are elements heavier than hydrogen and helium (like iron, calcium, and oxygen). These stars are like people living in a wealthy district, full of heavy elements.
The Poor Neighbors: Surprisingly, they found a significant group (about 17-20%) of "metal-poor" stars. These are older, less enriched stars, similar to a poorer district on the outskirts.

4. The Big Clue: The "Slope" of the City

The most important finding is how these stars are arranged in space.

The Gradient: The team created a map showing that the metal content changes as you move away from the very center.
The Pattern: The stars closest to the center are the most metal-rich. As you move outward, the stars become progressively more metal-poor.
The Analogy: Imagine a coffee cup where the liquid at the very bottom is super strong and dark (high metallicity), and as you move up toward the rim, the coffee gets lighter and weaker (lower metallicity).

5. What This Means: The "Inside-Out" Story

This specific pattern—a gradient where the center is richer than the edges—tells a story about how the neighborhood was built.

The Theory: It supports an "inside-out" formation scenario.
The Story: Imagine gas flowing into the center of the galaxy like water down a drain. This gas formed stars right in the middle first. Over time, more gas flowed in, but by then, the gas had been "polluted" with heavy elements from the first generation of stars. This new, richer gas formed the next generation of stars slightly further out.
The Result: The cluster grew outward, layer by layer, creating the metal-rich center and the metal-poor edges. This suggests the Nuclear Star Cluster and the surrounding Nuclear Stellar Disc might be part of the same continuous structure, growing together like rings on a tree.

Summary

In short, the astronomers cleaned up old, blurry data to get a sharp new view of the Milky Way's core. They discovered that the stars there aren't all the same; there's a mix of rich and poor populations arranged in a specific way. This arrangement acts like a fingerprint, proving that the cluster likely grew from the inside out, rather than being a random pile-up of star clusters crashing together.

Technical Summary: Chemical Analysis of the Milky Way's Nuclear Star Cluster

Problem and Context
The formation and evolutionary history of the Milky Way's Nuclear Star Cluster (MWNSC) and its surrounding Nuclear Stellar Disc (MWNSD) remain subjects of debate. While photometric and kinematic studies suggest distinct formation histories or a smooth transition between the two components, chemical constraints have been limited by observational challenges. The extreme extinction toward the Galactic center ( $A_V \sim 30$ ) necessitates infrared observations, but previous spectroscopic surveys, such as those using KMOS (Feldmeier-Krause et al. 2017a, 2020), relied on low-resolution spectra ( $R \sim 3000-4000$ ). Analyzing cool M giant stars ( $T_{\text{eff}} < 4000$ K) at this resolution is challenging due to severe line blending with molecular features. Furthermore, previous analyses indicated a predominantly metal-rich population but suggested an anisotropic distribution of metal-poor stars in the Galactic North, a finding that required re-evaluation with improved modeling techniques.

Methodology
This study re-analyzes the KMOS dataset from Feldmeier-Krause et al. (2020) (FK1720) to derive improved stellar parameters ( $T_{\text{eff}}$ , $\log g$ , $[M/H]$ ) and $\alpha$ -element abundances. The methodology involves several key technical upgrades:

Improved Line Lists and Model Grids: The authors utilize a dedicated synthetic model grid generated with the Spectroscopy Made Easy (SME) code, optimized specifically for cool M giants. Unlike the PHOENIX library used in FK1720, this grid incorporates updated atomic line lists (e.g., from the solar center intensity atlas) and molecular line lists for CN, CO, SiO, H $_2$ O, and OH. Crucially, the synthesis accounts for Non-Local Thermodynamic Equilibrium (NLTE) effects using pre-computed departure coefficients for Mg, Si, Ca, and Fe.
Spectral Fitting and Bayesian Inference: The analysis employs the STARKIT code with the MultiNest nested sampling algorithm. The process includes a sixth-order polynomial continuum normalization with $2\sigma$ clipping prior to fitting. Uniform priors are applied for temperature ( $2800 < T_{\text{eff}} < 4600$ K) and surface gravity ( $-0.5 < \log g < 3$ ), avoiding the use of photometric constraints for $\log g$ due to large uncertainties.
Validation: The derived parameters are validated against high-resolution ( $R=45,000$ ) IGRINS spectra of cool M giants from the solar neighborhood (Nandakumar et al. 2024b) and specific MWNSC benchmark stars (Nandakumar et al. 2025). The high-resolution data were degraded to KMOS resolution to ensure a direct comparison.
Sample Selection: The study cross-matches the FK1720 sample with GALACTICNUCLEUS photometry to separate MWNSC members from foreground stars and the MWNSD using a color-magnitude diagram ( $H-K_s$ vs. $K_s$ ). Quality control filters exclude stars with large parameter variations between exposures, low signal-to-noise ratios ( $<10$ ), or parameters near the synthetic grid boundaries.

Key Results

Stellar Parameters and Uncertainties: The re-analysis yields a high-quality sample of 1,140 M giant stars. Validation against high-resolution spectra demonstrates that stellar parameters can be recovered with typical uncertainties of $\sim 150$ K in $T_{\text{eff}}$ , $\sim 0.4$ dex in $\log g$ , and $\sim 0.2$ dex in $[M/H]$ . The $\alpha$ -element abundances agree within $0.1$ dex with high-resolution data after applying a linear correction for systematic overestimation at high $\alpha$ -abundances.
Metallicity Distribution Function (MDF): The MDF reveals a bimodal distribution best described by a two-component Gaussian Mixture Model. The dominant component is metal-rich ( $[M/H] \simeq 0.26$ dex, $\sigma \simeq 0.10$ ), comprising the majority of the sample. A significant metal-poor population is identified centered at $[M/H] \simeq -0.77$ dex ( $\sigma \simeq 0.24$ ), accounting for approximately 17% of the sample. This fraction of metal-poor stars is notably higher than the $\sim 10\%$ reported in FK1720.
Spatial Distribution and Gradients: The metal-poor population is spatially concentrated in the southern outer parts of the MWNSC, contrasting with previous claims of anisotropy in the Galactic North. The study constructs a metallicity map and identifies a clear negative radial metallicity gradient of $-0.10 \pm 0.02$ dex/pc relative to Sgr A*. This gradient is statistically significant ( $p=0.0002$ , $R^2 \sim 0.84$ ) and confirmed via MCMC analysis.

Significance and Claims
The paper claims that the observed negative radial metallicity gradient provides strong evidence for an "inside-out" formation scenario for the MWNSC. In this scenario, the cluster forms from gas inflows that create circum-nuclear rings, with star formation progressing outward over time, leaving older, metal-poor stars at larger radii and younger, metal-rich stars closer to the center.

The authors assert that their improved synthetic grid and rigorous treatment of NLTE effects are responsible for the detection of a larger metal-poor fraction compared to previous studies. They conclude that the MWNSC and MWNSD may share a connected formation history, potentially as a single structure where the NSD represents the growing edge of the NSC. The study emphasizes that a fully consistent analysis of both the NSC and NSD using the same methodology is required to definitively resolve their interplay, a task anticipated for future instruments like MOONS at the VLT.