1. NGC 346 and metal-poor massive stars
The open cluster NGC 346 (Fig. 1) is the most massive star cluster in the Small Magellanic Cloud (SMC), with a metallicity 5-7 times lower than the Milky Way (Dufour et al. Reference Dufour, Shields and Talbot1982; Larsen et al. Reference Larsen, Clausen and Storm2000; Trundle et al. Reference Trundle, Dufton and Hunter2007). For this reason, it makes a perfect laboratory to study the properties of massive stars at low metallicity. Its O-star population has been previously studied by Rickard et al. (Reference Rickard, Hainich and Hamann2022), but the weakness of the wind diagnostics in the UV and optical resulted in only upper limits for the fainter targets. We aim to complement this study with unprecedented JWST/NIRSpec observations. In this work, we present the first infrared spectroscopic observations of some of its O-stars. We describe the observations and data reduction in Sect. 2, while in Sect. 3 we show the observations of the Br
$\alpha$
line of our sample and describe its nature as a mass-loss diagnostic. Afterward, we present the results of a first comparison with stellar atmosphere models (Sect. 4). Finally, we discuss what the future looks like regarding massive stars and infrared spectroscopy (Sect. 5).
NIRCam image of NGC 346 with the NIRSpec spectra observed for this work as described in Mang Román et al. (Reference Mang Román, Zeidler and Hamann2025).

2. Observations and data reduction
The data was obtained with the NIRSpec spectrograph (Böker et al. Reference Böker, Beck and Birkmann2023) using the S200A1 slit, the high-resolution grating G395H, and the F290LP filter. This configuration enables a wavelength coverage between 2.8 and 5.1
$\mu$
m and a resolution of
$R \sim 2700$
. The data reduction was performed to remove cosmic rays and detector artifacts. However, the critical part for these observations was the background subtraction. As seen in Fig. 1, the targets are embedded in a dense H ii region, resulting in a strong contamination around the lines of interest. As described in Mang Román et al. (Reference Mang Román, Zeidler and Hamann2025), for some of the targets, the background subtraction was performed manually by extracting the nebular signal from the 2D spectrum. The reduced spectra can be seen in Fig. 1.
3. Br
$\alpha$
: theory and observations
The Br
$\alpha$
line at 4.05
$\mu$
m was predicted to be in emission even for hydrostatic atmospheres due to a non-LTE effect (Auer & Mihalas Reference Auer and Mihalas1969). More recently, Najarro et al. (Reference Najarro, Hanson and Puls2011) showed, using the non-LTE stellar atmosphere code CMFGEN (Hillier & Miller Reference Hillier and Miller1998), that the emission core of the line is sensitive to the mass-loss rate,
$\dot{M}$
. They found that the emission core reacts inversely to the value of
$\dot{M}$
, making it a very powerful tool to determine weak winds. The authors were able to apply this line to measure
$\dot{M}$
in a sample of Galactic O-stars. However, until the arrival of JWST, this was impossible for targets in the Magellanic Clouds. The observed Br
$\alpha$
spectra of our targets are shown in Fig. 2. We see that the line is in emission for all of them. For the brightest targets, the profile has a broad component, while only a narrow emission core can be seen for the fainter ones.
Normalized observed Br
$\alpha$
spectra for 12 of the 13 targets.

4. First comparison with PoWR models
For a preliminary analysis, we compared the observed NIRSpec spectra with models computed with the non-LTE stellar atmosphere code (Gräfener et al. Reference Gräfener, Koesterke and Hamann2002; Hamann & Gräfener Reference Hamann and Gräfener2003; Sander et al. Reference Sander, Shenar and Hainich2015). In Fig. 3, we show the PoWR spectrum of Br
$\alpha$
for two stars with weak and strong winds. We see that for the luminous star, Br
$\alpha$
reacts proportionally to
$\dot{M}$
, while the opposite occurs for the fainter one. Additionally, we fitted the model SED to multi-wavelength observations, including UV spectra and optical photometry (Fig. 4). As extinction has a negligible effect on the IR part of the spectrum, the new NIRSpec observations can be used to strongly constrain luminosity. The revised luminosities are systematically lower than when only considering the UV and the optical ranges (Mang Román et al. Reference Mang Román, Zeidler and Hamann2025).
PoWR model Br
$\alpha$
spectra of two stars with different
$\dot{M}$
(left panels) and clumping density contrast, D (right panels). Top and bottom panels show, respectively, models for thick (luminous supergiant) and thin (faint dwarf) winds.

Synthetic SED (red solid line) adjusted to UV, optical and IR observations (black lines and boxes). The blue dashed line is the model SED without considering the IR.

5. Future of IR spectroscopy in massive stars
The results presented in this work show that IR spectroscopy can be a powerful tool for massive star astrophysics. However, data reduction reminds extremely critical, especially when observing targets embedded in ionized nebulae. The capabilities of the JWST will allow to create a bridge between stars at different metallicities, and reveal useful spectral signatures in reddened or obscured stars. The presence of Br
$\alpha$
in low-Z stars opens a new window into the study of stellar winds.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S174392132610595X

α




α
α
M˙