---
title: 'Westerlund 1: A Massive Galactic Star Cluster'
url: https://www.emergentmind.com/topics/westerlund-1-wd1
type: topic
---

# Westerlund 1: A Massive Galactic Star Cluster

Westerlund 1 (Wd1) is a young massive Galactic star cluster observed through extreme foreground extinction and containing one of the richest known concentrations of evolved massive stars in the Milky Way, including Wolf–Rayet stars, yellow and red hypergiants, a luminous blue variable, OB supergiants, massive eclipsing binaries, and a magnetar. Its stellar content spans unevolved and highly evolved phases, while its surrounding medium shows strong signatures of mechanical and high-energy feedback. For that reason, Wd1 is widely used as a benchmark for massive-star evolution, binarity, mass loss, cluster dynamics, and cluster–interstellar-medium coupling. At the same time, its fundamental parameters have remained unusually contentious, with the distance, extinction law, and especially the age and star-formation history repeatedly revised by Gaia astrometry, eclipsing-binary analysis, infrared spectroscopy, and multiwavelength studies [2204.09414][2210.04985][2602.24218].

## 1. Extinction and observational regime

The line of sight to Wd1 is heavily reddened. Damineli et al. derived an extinction law over $0.4$–$4.8~\mu$m characterized in the near-infrared by a single-index power law, $A_\lambda \propto \lambda^{-\alpha}$, with $\alpha = 2.126 \pm 0.080$, and found a total-to-selective ratio $R_V = 2.50 \pm 0.04$, lower than the canonical $R_V=3.1$ often assumed for the diffuse Galactic ISM [1607.04639]. For cluster members excluding dusty WC stars they obtained $\langle A_{K_s}\rangle = 0.736 \pm 0.079$, corresponding to $\langle A_V\rangle \simeq 11.4 \pm 1.2$ mag, with an interstellar foreground component $A_{K_s}({\rm ISM}) = 0.63 \pm 0.02$ or $A_V({\rm ISM}) = 9.7 \pm 0.3$ mag [1607.04639]. The same study emphasized that previous $A_{K_s}$ estimates were systematically larger by $\sim20$–$30\%$, implying systematic biases in older luminosity and distance estimates.

The obscuration is not spatially uniform. Across the $\approx3.5'$ cluster field, member extinctions span $A_{K_s}=0.55\ldots1.17$, equivalent to $A_V=8.5\ldots17.9$, so that $\Delta A_V \simeq 8.7$ mag, with a clear SW$\rightarrow$NE gradient that follows the warm-dust distribution mapped by Herschel and WISE [1607.04639]. This differential reddening is astrophysically significant because IMF reconstruction, mass segregation measurements, and HRD placement all become sensitive to the adopted extinction law and to local extinction corrections.

The same work also calibrated the $8620$ Å diffuse interstellar band as a high-extinction tracer. For the combined Wd1 and literature sample, the best-fit polynomial is
$$
A_{K_s} = 0.612\,EW - 0.191\,EW^2,
$$
with $EW$ in Å, extending the low-reddening Munari et al. relation into the non-linear regime [1607.04639]. A later infrared spectroscopic study reiterated that optical work had long been hampered by $A_V \approx 11$–$13$ mag, which had prevented direct detection of the unevolved main sequence before near-infrared spectroscopy with KMOS [2602.24218].

## 2. Distance scale

Wd1 has been the subject of unusually discrepant distance determinations. Near-infrared photometry and empirical fiducial-line fitting yielded $d \simeq 3.8$ kpc in Lim et al. [1211.5832]. A Bayesian analysis of Gaia DR2 parallaxes, explicitly modeling the cluster as a population embedded in a broader Galactic field-star distribution and including a parallax zero-point term, inferred a cluster parallax of $0.35^{+0.07}_{-0.06}$ mas, corresponding to $d = 2.6^{+0.6}_{-0.4}$ kpc, and argued that $5$ kpc was excluded at the $99\%$ confidence level [1901.06582]. That result would have reduced the turn-off mass from $\simeq40~M_\odot$ to about $22~M_\odot$ and lowered the total cluster mass by a factor of about four [1901.06582].

Subsequent Gaia EDR3 work moved the distance scale back upward. Using spectroscopically confirmed members and red-source zero-point corrections, Navarete et al. obtained $d_{\rm gedr3}=4.06^{+0.36}_{-0.34}$ kpc from 172 sources [2204.09414]. Independently, Rocha et al. modeled the eclipsing binary W36 as a double-lined O6.5 III + O9.5 IV system and derived $d_{\rm W36}=4.03 \pm 0.25$ kpc, consistent with the Gaia EDR3-based values [2210.04985]. The combined distance is
$$
d_{\rm Wd1} =
\frac{d_{\rm W36}/\sigma_{\rm W36}^2 + d_{\rm Gaia}/\sigma_{\rm Gaia}^2}
{1/\sigma_{\rm W36}^2 + 1/\sigma_{\rm Gaia}^2}
= 4.05 \pm 0.20~{\rm kpc},
$$
with corresponding distance modulus $(m-M)_0 = 13.04 \pm 0.11$ mag [2204.09414][2210.04985].

More recent work on the newly identified main sequence adopted $d = 4.23^{+0.23}_{-0.21}$ kpc and argued that its derived age is robust against reasonable changes in distance and extinction law [2602.24218]. Accordingly, the modern literature contains at least three distinct distance regimes: a DR2-based short distance near $2.6$ kpc, an intermediate near-infrared value near $3.8$ kpc, and an EDR3/eclipsing-binary scale near $4.0$–$4.2$ kpc [1901.06582][1211.5832][2204.09414]. This suggests that the treatment of crowding, red-source parallax zero points, and extinction remains decisive.

## 3. Age determinations and star-formation history

The age of Wd1 is the central controversy in its modern literature. Earlier cluster analyses based on the evolved population and rotating Geneva tracks placed the age near $5.0 \pm 1.0$ Myr [1211.5832], and Wd1 was long treated as a nearly coeval, $\sim4$–$5$ Myr system. That view was challenged when Beasor et al. integrated SEDs of the cool supergiants using Gaia EDR3, 2MASS, and new SOFIA/FORCAST photometry and found that the RSGs and YSGs are under-luminous by $\simeq0.4$ dex relative to $4$–$5$ Myr models, yielding a cool-supergiant age of $\tau_{\rm CSG}=10.4^{+1.3}_{-1.2}$ Myr; their re-analysis of the pre-main-sequence population gave $\tau_{\rm pre-MS}=7.2^{+1.1}_{-2.3}$ Myr, while the W13 eclipsing binary remained in tension at $\lesssim5$ Myr [2103.02609]. Navarete et al. independently derived $\tau_{\rm BPASS}=10.7 \pm 1.0$ Myr from the four red supergiants using BPASS v2.2 single-star tracks, in excellent agreement with the cool-supergiant result [2204.09414].

Rocha et al. then decomposed the upper HRD into multiple age components. Using Geneva rotating models with supersolar metallicity, they found that the secondary of W36 is matched by an isochrone of $t_{\rm W36B}=6.4 \pm 0.7$ Myr, the bulk of the OB giants and supergiants define a tight locus around $t_{\rm OB\,G/S}=7.1 \pm 0.5$ Myr, BA-type luminous evolved stars and yellow hypergiants scatter between $8$ and $11$ Myr, and the four bona-fide red supergiants lie closest to $t_{\rm RSG}=10.7 \pm 1.0$ Myr [2210.04985]. On that basis, they argued that Wd1 was not born in a single instantaneous burst but over at least three episodes spanning $\sim4$ Myr, with an earliest burst near $11$ Myr, a major burst near $7$ Myr, and a younger component at $<6.5$ Myr [2210.04985].

In direct contrast, the first spectroscopic identification of the main sequence with VLT/KMOS analyzed $\sim110$ members with CMFGEN atmospheres, identified 47 new O9–B1 III–V members, and derived a turn-off age of $5.5 \pm 1.0$ Myr at $4.23^{+0.23}_{-0.21}$ kpc [2602.24218]. That study reported a tight main-sequence band, found that evolved giants and supergiants fall along the same $5$–$6.3$ Myr isochrones, found no secondary “$10$ Myr” turn-off, and concluded that Wd1 formed in a single, moderately coeval burst $\approx5.5$ Myr ago [2602.24218].

| Diagnostic | Age result | Reported implication |
|---|---:|---|
| Cool supergiants | $10.4^{+1.3}_{-1.2}$ Myr | Single-age $4$–$5$ Myr models fail [2103.02609] |
| Red supergiants (BPASS) | $10.7 \pm 1.0$ Myr | Oldest population [2204.09414] |
| W36B | $6.4 \pm 0.7$ Myr | Younger component [2210.04985] |
| OB giants/supergiants | $7.1 \pm 0.5$ Myr | Major star-formation burst [2210.04985] |
| PMS re-analysis | $7.2^{+1.1}_{-2.3}$ Myr | Intermediate-age population [2103.02609] |
| NIR main-sequence turn-off | $5.5 \pm 1.0$ Myr | Moderate coevality, single burst [2602.24218] |

These incompatible chronometers have kept Wd1 at the center of discussions of extended star formation versus moderate coevality in young massive clusters. A plausible implication is that Wd1’s age is highly sensitive to the adopted tracer population and evolutionary prescription.

## 4. Stellar content, eclipsing binaries, and multiplicity

Wd1 contains an extraordinary evolved population. Recent summaries describe $\sim24$ Wolf–Rayet stars, numerous OB supergiants, yellow and red hypergiants, a luminous blue variable, and a magnetar, while earlier work also emphasized four massive eclipsing binaries: Wddeb, W13, W36, and WR77o [2602.24218][1009.4709]. This stellar inventory is the empirical basis for Wd1’s role as a laboratory for mass loss, binary interaction, and compact-remnant formation.

The eclipsing binary W36 has become central to both the distance scale and stellar-parameter calibration. Optical light curves and FLAMES/VLT radial velocities were modeled simultaneously with the Wilson–Devinney code coupled to an MCMC sampler, with IGRINS and SPAMMS constraints fixing the primary temperature at $T_1=37100$ K and the mass ratio at $q \approx 0.71$ [2210.04985]. The resulting solution identified W36 as O6.5 III + O9.5 IV with orbital period $P \simeq 3.1722$ d and yielded $M_1 = 22.9^{+2.7}_{-2.5}~M_\odot$, $R_1 = 12.7^{+0.6}_{-0.5}~R_\odot$, $M_2 = 16.3^{+2.2}_{-2.3}~M_\odot$, and $R_2 = 10.5^{+0.6}_{-0.7}~R_\odot$ [2210.04985]. Earlier spectroscopy of W13, W239, W43a, and W3003 had already shown that close and interacting binaries are common among Wd1’s evolved stars; that survey inferred a spectroscopic binary fraction of at least $\sim60$–$70\%$ among evolved massive members and highlighted binary mass transfer as a pathway into the Wolf–Rayet phase [1010.6270].

Multiplicity diagnostics from other wavelengths strengthen that view. A deep ATCA radio census detected 30 cluster stars, including 10 WR stars, 5 YHGs, 4 RSGs, 1 LBV, the sgB[e] star W9, and 7 O/B supergiants; negative spectral-index limits for W10, W15, W18, and W1031 were interpreted as evidence for non-thermal wind–wind shocks and new OB supergiant binary candidates [1909.04339]. The same survey measured mass-loss diagnostics across evolutionary phases, with inferred $\dot M$ values from $10^{-7}$ to $10^{-4}~M_\odot\,{\rm yr}^{-1}$ and clumping ratios for several WR winds that generally decline with radius [1909.04339].

X-rays now indicate that the WR population is even more binary-dominated than optical and radio work had implied. In the deepest Chandra study to date, 20 of the 24 known WR stars were detected, 17 showed Fe XXV emission near $6.7$ keV, nine showed clear variability over a year-long baseline, and the WR binary fraction was estimated as $\geq80\%$ from X-ray diagnostics alone and $\geq92\%$ when combined with other wavelength information [2408.11087]. The newly identified main-sequence and giant population also shows strong multiplicity signatures: multi-epoch KMOS spectroscopy found that $\sim65\%$ of the OB stars with multi-epoch coverage exhibit radial-velocity variability, with 63% of unevolved OB stars showing $\Delta RV > 20~{\rm km\,s^{-1}}$ at $\geq4\sigma$ significance [2602.24218]. Wd1 is therefore not only rich in massive stars; it is also rich in interacting and probable interacting systems.

## 5. Mass function, structure, and dynamics

The cluster’s stellar mass spectrum has been measured from the high-mass regime down into the subsolar domain. Lim et al. derived a high-mass IMF slope of $\Gamma = -0.8 \pm 0.1$, shallower than Salpeter, and obtained a total cluster mass in excess of $5.0 \times 10^4~M_\odot$, with stitched IMF integrations giving $M_{\rm tot}\simeq(7.3$–$7.8)\times10^4~M_\odot$ [1211.5832]. Andersen et al., using deep HST/WFC3 observations, extended the mass function to $0.15~M_\odot$ in the outer cluster and found a flattening consistent with that of the Galactic field and nearby embedded clusters. In the outer region they obtained a log-normal peak mass $m_c = 0.22 \pm 0.02~M_\odot$ and width $\sigma = 0.35 \pm 0.05~M_\odot$, with a total mass of $(4.4$–$5.7)\times10^4~M_\odot$ and the conclusion that Wd1 is sufficiently massive to remain bound and could potentially evolve into a low-mass globular cluster [1602.05918].

Mass segregation is present, but its strength depends on the stellar-mass range and method. Lim et al. reported a strong radial steepening of the high-mass IMF slope from the core outward, identifying clear mass segregation [1211.5832]. Andersen et al. likewise found strong segregation for supersolar stars but no significant segregation for $M \lesssim 1~M_\odot$ beyond $r\sim1.5$ pc [1602.05918]. A later HST proper-motion analysis measured kinematics for 10,346 stars, assigned kinematic memberships with a three-component Gaussian mixture model, and after extinction and completeness corrections modeled the surface-density profile with the Elson–Fall–Freeman form
$$
\Sigma(r)=\Sigma_0\,[1+(r/a)^2]^{-\gamma/2}.
$$
For stars with $1.00 \le M \le 12.14~M_\odot$ it found $\gamma = 2.62^{+0.32}_{-0.26}$, $r_c = 0.78 \pm 0.05$ pc, and only weak mass segregation, with $\Lambda_{\rm MSR}=1.11\pm0.11$ and a core radius that decreases modestly with mass [2501.16755]. That study interpreted the weak segregation in the $1$–$12~M_\odot$ range as dynamical rather than primordial [2501.16755].

The same HST kinematic work also showed that Wd1 is elongated, with axis ratio $\epsilon \simeq 0.70$ and eccentricity $e \simeq 0.71$, aligned at position angle $\simeq56^\circ$ East of North, approximately along the Galactic plane [2501.16755]. Its 1D velocity dispersion is $3.42 \pm 0.10~{\rm km\,s^{-1}}$, below the quoted virial expectation for $M_{\rm tot}\sim4$–$6\times10^4~M_\odot$ and $r_h \simeq 1$ pc, implying a subvirial state [2501.16755]. The derived crossing time is $0.30$ Myr and the relaxation time is $0.26$ Gyr [2501.16755].

The intermediate-mass regime has recently been filled in by EWOCS-VI. Using VVV and VVVX VIRAC2 astrometry and PSF photometry, HDBSCAN clustering in a six-dimensional parameter space identified 1,286 high-probability members with $12 < J < 18$ mag, corresponding to $\sim1.5$–$20~M_\odot$ for adopted PARSEC 5 and 6 Myr isochrones at $A_{K_s}=0.6$ mag and $d=4.23$ kpc [2606.27623]. Excluding the saturated central $\sim2'$, a King-profile fit yielded $r_c = 5.0 \pm 2.8'$ and $r_t = 9.44 \pm 2.23'$, while the spatial distribution showed an asymmetric halo extending to the north-east [2606.27623]. This bridges the observational gap between Gaia-selected high-mass stars and the HST/JWST low-mass population.

## 6. Feedback, cluster wind, and high-energy environment

Wd1 is not only a stellar aggregate but also an engine of cluster-scale feedback. XMM-Newton spectroscopy of the diffuse emission in the central $2'$ detected the He-like Fe line at $6.7$ keV and established that the hard diffuse component is predominantly thermal, not non-thermal [1106.2665]. A two-temperature thermal fit gave $N_H = (2.03^{+0.11}_{-0.15}) \times 10^{22}~{\rm cm^{-2}}$, $kT_1 = 0.68^{+0.12}_{-0.13}$ keV, $kT_2 = 3.07^{+0.60}_{-0.38}$ keV, and an unabsorbed $2$–$8$ keV luminosity of $(2.56 \pm 0.15)\times10^{33}~{\rm erg\,s^{-1}}$ [1106.2665]. Comparison with a thermalized cluster-wind model showed that the cluster wind likely contributes the majority of the hard emission, with an unresolved PMS population contributing $\lesssim15\%$ and supernova remnants unlikely to contribute significantly at the present epoch [1106.2665].

Radio imaging provides complementary evidence that this cluster wind shapes circumstellar ejecta and the intracluster medium. ATCA observations of cool hypergiants and red supergiants detected nine of ten targets, with eight spatially resolved; the nebulae of Wd1-4a, -12a, -265, -20, and -26 show cometary or tail-like morphologies directed away from the cluster core, which was interpreted as compelling evidence for ram-pressure sculpting by a cluster wind [1803.07008]. A broader radio census then revealed a pervasive asymmetric “sea” of free–free emission at $\sim0.6$–$2~{\rm mJy\,beam^{-1}}$ across the central $1.5'$, along with clumpy shells, filaments, and 53 ATCA-only knots lacking optical counterparts [1909.04339]. That morphology was attributed to a dense intracluster medium shaped by a collective cluster wind, photo-ionization, and non-thermal wind–wind shocks [1909.04339].

The circumcluster environment is also a high-energy source. H.E.S.S. observations of HESS J1646–458 show very extended TeV emission around Wd1, roughly circular with radius $\sim1^\circ$ and a complex multi-peaked surface-brightness distribution, with no pronounced peak at the cluster center and no significant energy-dependent morphological change up to $\geq10$ TeV [2108.03003]. A hadronic $\pi^0$-decay fit with naima yielded a proton index $\Gamma_p = 2.33 \pm 0.06$, cutoff energy $E_p^{\rm cut} = 400^{+250}_{-130}$ TeV, and total proton energy above $1$ TeV of $W_p \simeq 5\times10^{49}(n/10~{\rm cm^{-3}})(d/3.9~{\rm kpc})^2$ erg [2108.03003]. Earlier Fermi-LAT analysis of the GeV counterpart found a Gaussian extension $\sigma_S = (0.475 \pm 0.050)^\circ$, offset by $\sim1^\circ$ from the cluster, and concluded that a PWN interpretation can account for the GeV flux but leaves the TeV emission unexplained, whereas a hadronic scenario can reproduce both GeV and TeV data but requires slow diffusion and a high proton energy input of $\sim10^{51}$ erg [1306.5642].

Recent CO and H I mapping has tightened the case for a hadronic component. Molecular clouds at $V_{\rm LSR}\sim-32$ and $\sim-55~{\rm km\,s^{-1}}$ were found to be physically associated with Wd1, with the $-32~{\rm km\,s^{-1}}$ component showing a cavity-like structure of radius $\sim3.4$ pc and expansion speed $\simeq5~{\rm km\,s^{-1}}$, implying a dynamical age of $\simeq0.7$ Myr and a recently formed wind-blown bubble [2605.25467]. The $-55$ and $-32~{\rm km\,s^{-1}}$ clouds also show complementary morphology and an intermediate-velocity bridge at $V_{\rm LSR}\sim-44~{\rm km\,s^{-1}}$, consistent with a cloud–cloud collision scenario for triggering Wd1’s formation [2605.25467]. On larger scales, the total interstellar gas mass within the TeV shell is $\sim1.6\times10^6~M_\odot$, the total proton column density shows a positive correlation with TeV flux across 16 subregions with $r\approx0.52$, and under the hadronic assumption the required cosmic-ray proton energy is $\sim6\times10^{49}$ erg [2605.25467]. In that interpretation, Wd1 links massive-cluster formation, present-day feedback, and Galactic cosmic-ray acceleration within a single system.

Source: https://www.emergentmind.com/topics/westerlund-1-wd1