Open Access
Issue
A&A
Volume 682, February 2024
Article Number A49
Number of page(s) 22
Section Catalogs and data
DOI https://doi.org/10.1051/0004-6361/202347778
Published online 02 February 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

The star formation rate and the properties of the most common star-forming environments in galaxies vary over time. When considering cosmological timescales, the star formation rate is known to reach its peak at approximately 𝓏 ~ 2–3 and then gradually decline (e.g., Hopkins & Beacom 2006; Dunlop 2011; Madau & Dickinson 2014). In the local Universe, mergers play a dominant role in shaping the star formation process in galaxies (Rieke & Rujopakarn 2011) as they influence the overall properties of the interstellar medium. Such interactions happen frequently, and various studies have demonstrated that interacting galaxies undergo periods of intense star formation (e.g., Larson & Tinsley 1978; Smith & Struck 2010) because of the considerable impact interactions have on the stars formation process, for instance from the enhancement of the star cluster formation rate due to close encounters (e.g. in the Magellanic Clouds) and the ram pressure stripping enhancing star formation (e.g., jellyfish galaxies). Noticeable examples are the nearby interacting galaxies M51 and M82, where we observe extreme star formation taking place in very massive young clusters with masses reaching several times 105 M (known as super star clusters; de Grijs et al. 2001, 2003b,a). Generally, these highly massive star clusters constitute the dominant star-forming environments in starburst galaxies and are likely prevalent during the peak era of cosmic star formation (e.g., Figer 2008; Adamo et al. 2020).

In the Milky Way, current estimates of the star forming heavily rely on the methods employed. For example, Robitaille & Whitney (2010) derived a range of 0.68–1.45 M yr−1 based on the population of young stellar objects identified in the Spit𝓏er/IRAC survey of the Galactic plane GLIMPSE (Benjamin et al. 2003). On the other hand, Licquia & Newman (2015) applied a hierarchical Bayesian statistical method to previous analyses and determined a star formation rate of about 1.6 M yr−1. For comparison, recent estimates of the star formation rate in M51 range from 4.8 M yr−1 (from a 158 µm map of the galaxy; Pineda et al. 2018) to 2.7 M yr−1 (from combined

UV+optical spectral energy distribution fitting; Eufrasio et al. 2017), while in M82 star formation rates of 2–4 M yr−1 were observed (de Grijs et al. 2001). Nevertheless, all these studies indicate that our Galaxy does not currently have a high star formation rate. Consequently, it is not surprising that the Milky Way lacks a prominent population of super star clusters with masses exceeding 104 M. In order of distance from the Sun, the most massive clusters known are Westerlund 1 (2.6–5 kpc; Aghakhan-loo et al. 2020; Clark et al. 2005), Westerlund 2 (~4.2 kpc; Vargas Álvarez et al. 2013), NGC 3603 (7.6 kpc; Melena et al. 2008), the Arches and Quintuplet clusters (both at ~8.5 kpc; Figer et al. 2002, 1999), Mercer 81 (11 kpc; Davies et al. 2012), and Mercer 30 (12 kpc; de la Fuente et al. 2016). Similar regions in terms of mass, but with with a low stellar density, are the Cygnus OB2 association (1.4 kpc; Rygl et al. 2012) and the W3 complex (about 2 kpc; Hachisuka et al. 2006). Slightly older supermassive star clusters (10–20 Myr) are found in the Scutum-Crux arm (about 6 kpc from the Sun; Figer et al. 2006; Davies et al. 2007; Clark et al. 2009). Despite their limited number, these super star clusters hold significant importance as they enable the study of star and planet formation, as well as early stellar evolution, in a star-forming environment that was characteristic of epochs when the Milky Way had higher rates of star formation than today and most of the field stars in our Galaxy formed.

In this paper we present the Extended Westerlund 1 and 2 Open Clusters Survey (EWOCS) project, which is focused on studying star and planet formation and early stellar evolution in compact starbursts, using Westerlund 1 and 2 as first science cases. In particular, this paper focuses on the catalog of X-ray sources detected in the deep Chandra observations of Westerlund 1. The paper is organized as follows: we present Westerlund 1 and the EWOCS project in Sect. 2. The EWOCS observations are described in Sect. 3, the procedure for source detection is described in Sect. 4 and that of source validation and extraction in Sect. 5. The final catalog of the X-ray sources in Westerlund 1 is described in Sect. 6.

2 Westerlund 1 and the EWOCS project

Westerlund 1 is located at RAJ2000=16h47m04s and DecJ2000= −45°51′05″, corresponding to Galactic coordinates l=339.55º and b=−00.40°. The cluster was discovered by Westerlund (1961) through observations made with the 26-inch Uppsala-Schmidt telescope at Mt. Stromlo Observatory in Australia. From these initial observations, it became evident that Westerlund 1 is a very massive cluster. Today, it is considered to be the most massive young cluster known within the Milky Way, with mass estimates ranging from approximately 5 × 104 M to over 105 M (Clark et al. 2005; Brandner et al. 2008; Gennaro et al. 2011; Lim et al. 2013; Andersen et al. 2017).

Despite over 60 yr of studies and observations of Westerlund 1, the tension regarding the parameters of this distinctive cluster remains unresolved. This is primarily due to its compact nature and the significant extinction that has long hindered the ability to resolve its low-mass stars.

The distance to the cluster has been a subject of longstanding debate. The initial estimate by Westerlund (1961) was of 1.4 kpc. However, the same authors later presented a more distant estimate of 5 kpc based on photographic observations in the VRI bands (Westerlund 1968). The first study utilizing CCD imaging of the cluster (Piatti et al. 1998) reported a distance estimate of 1.0±0.4 kpc. However, this estimate was based on the incorrect assumption that all cluster members were on the main sequence.

Several authors have made distance estimates for Westerlund 1 based on the analysis of its rich population of massive stars. For example, Clark et al. (2005) based their estimate on six yellow hypergiants (YHGs), assuming that these stars have the standard luminosity for this class of objects (log(L/L)~5.7; Smith et al. 2004), and adopting an extinction of AV=11m, resulting in a distance range between 2 kpc and 5.5 kpc. A similar value was found by Crowther et al. (2006) through infrared analysis of WN and WC stars. Koumpia & Bonanos (2012) derived a distance of 3.7±0.6 kpc from the analysis of the dynamics and geometry of the eclipsing binary W13. By comparing the cluster locus in color-magnitude diagrams with suitable isochrones, Brandner et al. (2008) determined a distance of 3.55±0.17 kpc, while Gennaro et al. (2011) found a distance of 4.0±0.2 kpc, and Lim et al. (2013) reported a distance of about 3.8 kpc. An independent estimate (3.9+0.7 kpc) was provided by Kothes & Dougherty (2007) using the radial velocity of HI clouds in the direction of the cluster, assuming they were physically connected to Westerlund 1.

More recently, Gaia data have been extensively utilized to measure the distances of star clusters, providing precise values up to distances of about 1 kpc (Gaia Collaboration 2016). However, for more distant clusters, careful analysis and assumptions are required to obtain reliable distance measurements. Consequently, it is not surprising that different estimates of the distance to Westerlund 1 have emerged from authors who have analyzed Gaia data. Aghakhanloo et al. (2020) conducted a Bayesian analysis of Gaia data along the line of sight to Westerlund 1 and obtained a mean cluster parallax 0.35_{ - 0.06}^{ + 0.07}0.350.06+0.07$0.35_{ - 0.06}^{ + 0.07}$ mas, which corresponds to a distance of 2.60.4+0.6$2.6_{ - 0.4}^{ + 0.6}$ kpc and is in tension with the previous estimate of approximately 0.19 mas provided by Clark et al. (2020). Focusing on known members of Westerlund 1, Davies & Beasor (2019) found a distance of 3.90.64+1.0$3.9_{ - 0.64}^{ + 1.0}$ kpc. More recently, Negueruela et al. (2022) carried out a detailed determination of candidate members in Westerlund 1 using Gaia Early Third Data Release (EDR3; van Leeuwen et al. 2021) data and obtained a distance of 4.230.21+0.23$4.23_{ - 0.21}^{ + 0.23}$ kpc, suggesting that the cluster is located in the Norma arm. A similar estimate from the Gaia/EDR3 was obtained by Navarete et al. (2022).

Given the uncertainty surrounding the distance to Westerlund 1, it is not surprising that estimates of the cluster’s age provided by different authors also vary significantly. Age estimates in the range of 3.2 to 5 million years have been derived using isochrone fitting on the high-mass sequence and arguments based on the diverse population of massive stars, including Wolf-Rayet (WR) stars, YHGs, and red supergiants (RSGs; Clark et al. 2005; Crowther et al. 2006; Brandner et al. 2008; Ritchie et al. 2010; Gennaro et al. 2011; Koumpia & Bonanos 2012; Kudryavtseva et al. 2012; Mackey et al. 2015). These authors found a relatively narrow age spread, with an upper limit of 0.4 million yr indicating that Westerlund 1 likely formed in a single burst of star formation (Kudryavtseva et al. 2012). However, some of these estimates are based on arguments that strictly apply to single stars, whereas it is known that the binary fraction among the massive members of Westerlund 1 is very high (Crowther et al. 2006). More recent studies suggest a more complex star formation history and a slightly older age (Aghakhanloo et al. 2020; Beasor et al. 2021; Navarete et al. 2022; Negueruela et al. 2022). In particular, arguments based on spectral energy distribution fitting and the luminosity of individual RSGs support an age estimate exceeding 10 Myrs (Beasor et al. 2021; Navarete et al. 2022), although this estimate is in tension with other properties of the cluster.

There is a general consensus in the literature regarding other important properties of Westerlund 1, including its high extinction, large mass, and notably, its impressive population of massive stars. The significant extinction toward Westerlund 1 has been acknowledged since the initial publication on this cluster, where an approximate visual extinction of AV~12m has been found (Westerlund 1961). Subsequent estimates range from 10m to 13m of visual extinction (Negueruela et al. 2010; Lim et al. 2013; Damineli et al. 2016). There is some disagreement regarding the extinction law in the direction of Westerlund 1: while according to Negueruela et al. (2010) it follows the standard law in the VRI bands, Lim et al. (2013) and Damineli et al. (2016) suggested a steeper extinction law in the near-IR (RV=2.50±0.04).

The most remarkable characteristic of Westerlund 1 is its large population of massive stars (Clark et al. 2005, 2020; Ritchie et al. 2009), which includes 24 WR stars (Clark & Negueruela 2002; Negueruela & Clark 2005; Skinner et al. 2006; Groh et al. 2006; Crowther et al. 2006), the luminous blue variable (LBV) Wd1-243 (Clark & Negueruela 2004), ten YHGs with spectral classes ranging from A5Ia+ to F8Ia+ (Clark et al. 2005)1, two blue stragglers (Clark et al. 2019), four RSGs (Wright et al. 2014b), seven blue hypergiants (BHGs), and over 100 bright OB supergiants dominated by spectral classes O9-B1 (Negueruela et al. 2010). Most of these sources are concentrated in the inner region of the cluster, spanning approximately 1 arcmin, with only a few more isolated massive stars, such as WR77.

In particular, Westerlund 1 hosts examples of every known transitional evolutionary phase between H-rich OB supergiants and H-depleted WR stars. This makes the cluster a unique target for studying massive stars and, specifically, for understanding how binarity and mass loss impact the evolutionary paths of these stars and how the initial stellar masses are linked to the types of compact objects that form at the end of their evolution. Winds and mass loss in these stars have been extensively studied with radio and millimeter-continuum observations, which have detected individual bright sources, such as W9, surrounded by extended nebulae, providing evidence of intense mass loss in the past (up to several 10−4 M per year, Dougherty et al. 2010; Fenech et al. 2018; Andrews et al. 2019). These short-lived and episodic mass-loss events appear to be necessary to explain the diversity of evolved massive stars in Westerlund 1.

Westerlund 1 is rich in binary systems. A high binary fraction has been identified in massive stars through spectroscopic (Ritchie et al. 2022), radio (Dougherty et al. 2010), infrared (Crowther et al. 2006), and X-ray (Skinner et al. 2006; Clark et al. 2008, 2019) observations. For instance, the WR population of Westerlund 1 has an estimated binary fraction of at least 70% (Crowther et al. 2006; Clark et al. 2008). Isolated stars are primarily found among the mid-B to F hypergiants, with the exception of the LBV star W243, whose binarity is supported by interferometric (Clark et al. 2019), X-ray (Mahy et al. 2022), and spectroscopic (Ritchie et al. 2009) observations.

This harsh environment is expected to have effects on the star formation process, the evolution and dispersal of protoplanetary disks, and the formation and early evolution of planets and their atmospheres. While no studies to date have been able to identify the population of protoplanetary disks in Westerlund 1 and explore the feedback provided by the starburst environment on their evolution and dispersal, several authors have attempted to quantify the cluster’s initial mass function (IMF) to investigate possible deviations from the universal law. An IMF consistent with the Salpeter (1955) law has been found by Brandner et al. (2008) in the 3.4–27 M range, by Gennaro et al. (2011) extrapolated in the 0.5–120 M range, and Andersen et al. (2017) down to 0.15 M in the outer cluster, while a shallower IMF slope was found by Lim et al. (2013), integrated in the 0.08–85 M range.

2.1 Previous X-ray observations

X-ray observations of young clusters provide valuable diagnostics for selecting pre-main-sequence (PMS) stars independently of the presence of circumstellar disks (e.g., Montmerle 1996), down to low stellar masses (e.g., Getman et al. 2005; Barrado et al. 2011). Additionally, in a cluster rich in massive stars with a very compact configuration, X-ray observations can reveal a plethora of processes and physical mechanisms that play an important role in the evolution of massive stars (Seward et al. 1979; Berghoefer et al. 1996). It is also worth mentioning that only the Chandra X-Ray Observatory (Weisskopf et al. 2002) can currently provide the high spatial resolution required to resolve individual X-ray faint sources in a crowded cluster like Westerlund 1. Given the designs of future X-ray missions currently in development, such observations will likely be challenging for quite some time after the Chandra era.

Both Chandra and XMM have been used in the past to observe Westerlund 1. The initial observations performed with Chandra reached a depth of approximately 58 ksec (P.I. Skinner) and resolved numerous X-ray sources (Skinner et al. 2006; Muno et al. 2006; Clark et al. 2008). Skinner et al. (2006) focused on the WR stars and their spectral properties, detecting 12 out of 24 known stars and finding strong evidence for the existence of very hot plasma in the circumstellar environment in the two brightest objects (W72/A and WRB), strongly suggesting the presence of a colliding winds in these binary systems; Muno et al. (2006) studied the diffuse X-ray emission and its dominating hard spectral component, which was later confirmed by Kavanagh et al. (2011) from 48 ks XMM-Newton observations. These authors detected a strong Fe 6.7 keV line in the diffuse emission spectrum, indicating its thermal nature. Clark et al. (2008) found that 46 known high-mass members of Westerlund 1 were detected in X-rays, and they supposed that the remaining ~60 X-ray sources detected in these images are likely PMS stars with masses ≤1.5 M.

Since its discovery by Muno et al. (2006), the magnetar CXO J164710.2-455216 (CXOU J16) in Westerlund 1 — the brightest X-ray source in the cluster — has garnered significant attention. Dedicated observations using XMM-Newton and Chandra/ACIS-S have accumulated a total exposure of 273.14 ksec (P.I.s Israel, Muno and Schartel) and 94.65 ksec (P.I.s Israel and Rea), respectively. A typical property of this class of pulsars is their frequent bursts and recurrent outbursts. In fact, three distinct outbursts from CXOU J16 have been observed in the past 17 yr (Borghese et al. 2019). The first one occurred in September 2006 and was triggered by a short burst that released an energy of approximately 1039 erg in the 15–150 keV band (Krimm et al. 2006). It was followed by a second outburst in September 2011 (Israel et al. 2011), during which the pulsar exhibited a peculiar behavior: the pulse profile evolved from a single peak in the pre-outburst phase to an energy-dependent tri-peaked profile post-outburst. The overall spectrum evolved from a single blackbody to a more complex shape that was well modeled by including an additional hotter blackbody component. The most recent outburst was again triggered by a short burst detected by the Swift Burst Alert Telescope in May 2017 (D’Ai et al. 2017). During these intense outburst activities, the magnetic field strength was estimated to range from 7×1013 G (a value typical for low-field magnetars, Perna & Pons 2011) to ~l·014 G (An et al. 2013; Israel et al. 2007).

2.2 The EWOCS project

The pre-EWOCS Chandra observations of Westerlund 1 have been analyzed by Townsley et al. (2018) in the framework of the Second Installment of the Massive Star-forming Regions Omnibus X-ray Catalog (MOXC2), identifying 1721 X-ray sources. This work has confirmed that Westerlund 1 is rich in X-ray bright sources, even though its low-mass stellar content remained undetected in the pre-EWOCS observations. According to these authors, the X-ray luminosity limit in the broad band where half of the brighter population is detected was log(Lx)=30.69, with Lx in erg s−1, corresponding to a 1.5 M star2.

The need to unveil the low-mass population of Westerlund 1 has motivated the 1 Msec observation of Westerlund 1 (P.I. Guarcello) with the Chandra Advanced CCD Imaging Spectrometer (ACIS-I; Garmire et al. 2003), which, together with a 18.9 h Cycle 1 JWST/MIRI and NIRCam observation (program ID 1905, P.I. Guarcello) and a 48 ksec NICER observation (P.I. Borghese) of CXOU J16, constitutes the set of new observations of the EWOCS project3. The main objective of EWOCS is to use Westerlund 1 and 2 as a test cases for understanding how star and planet formation, early stellar evolution, and the production of compact objects occur in a starburst environment. Specifically, the project aims to achieve the following objectives:

  • Unveil the low-mass stellar population of Westerlund 1 and 2, both in their core and halo. X-ray observations are expected to be critically important for selecting cluster members in the halo, where contamination from background and foreground sources could affect membership determination based on photometric data.

  • Determine the actual stellar content of the clusters, down to the low-mass regime, mainly thanks to the JWST observations; calculate their IMF down to the brown dwarf regime, and understand whether the starburst environment impact the formation of low-mass and very-low mass stars.

  • Study the clusters properties, particularly age, age spread, morphology and dynamics. The project aims to understand whether the clusters formed in a single burst of star formation or through a process spanning several million years, as well as how and if they will disperse.

  • Identify the disk-bearing population of the clusters, mainly though the JWST observations. Combining this with the detection of disk-less stars from the Chandra/ACIS-I observations and modeling of disks dispersal, we will finally assess how disks evolve and how planet formation proceeds in a starburst environment.

  • If planets can form, understand how they evolve while immersed in such an environment characterized by high local fluxes of UV and X-ray radiation and relativistic particles.

  • Study how binarity and mass-loss affect the evolution of the massive stars in the clusters, and how their initial mass is mapped into the type of compact objects formed at the end of their evolution.

  • Determine whether binarity across stellar masses is different in a starburst environment.

  • Study for the first time the status of CXOU J16 far from bursts, which will allow us to estimate the intrinsic properties of the pulsar.

  • Search for the expected population of compact objects that have been suggested to exist in Westerlund 1, since, under specific assumptions, up to ~65 core-collapse supernovae could have already occurred in the cluster (Muno et al. 2006; Brandner et al. 2008). Besides, Westerlund 1 is one of the few known star clusters meeting the properties required for the formation of intermediate mass black holes from runaway coalescence (Portegies Zwart et al. 2004). As estimated by Clark et al. (2008), such objects, if present and if currently accreting mass, should be observable with a very deep Chandra observation.

  • Understanding how stellar winds from massive evolved stars can affect the ISM to produce diffuse X-ray emission, whether this hot gas could affect star formation throughout the region, and whether we can prove ongoing accumulation of polluted material in the cluster core.

Figure 1 shows the contours of the pre-EWOCS and EWOCS observations of Westerlund 1 and CXOU J16, plotted over the combined Chandra event file and a KS band image of the cluster and the surrounding area obtained with the FourStar infrared camera mounted on the Magellan 6.5 m telescopes.

3 Chandra observations and data reduction

3.1 Overview of the Chandra/ACIS observations

The EWOCS survey also includes eight pre-EWOCS observations performed with ACIS-S, two of which were pointed at Westerlund 1 and six at CXOU J16. These observations were conducted between June 2005 and February 2018 (Table 1). Additionally, 36 EWOCS observations were carried out with the ACIS-I detector from June 2020 to August 2021. The aim point of each ACIS-I observation was adjusted based on the nominal roll angle, as indicated in Table 2. This adjustment was crucial to for avoiding gaps that cover the cluster core, the pulsar, or some of the brightest massive members, while ensuring the cluster remained within the inner arcminute. By adopting this design, we maximized the benefits of the subarcsecond spatial resolution and sensitivity in the central part of the ACIS-I detector to observe the cluster core, which is highly compact and crowded.

The total exposure for the pre-EWOCS observations is 151.93 ksec, while for the EWOCS observations it is 967.80 ksec. The exposure times for the individual EWOCS observations range from 11.92 ksec (Obs.ID 23272) to 39.55 ksec (Obs.ID 22319), with a mean exposure of 26.33 ksec. The EWOCS observations span over 1 yr, providing a robust baseline for studying the X-ray variability of the brightest sources. All EWOCS observations were conducted using the ACIS-I detector in imaging mode, utilizing all four chips (I0–I3). The observations were performed in the VERY FAINT mode, which employs telemetry in 5×5 pixel event islands for improved background suppres-sion4. When combined with the pre-EWOCS observations, the total time baseline exceeds 16yr, which is particularly valuable for studying certain sources in the cluster, such as the magnetar. Figure 2 displays a composite RGB ACIS imageofWesterlund 1, where colors represent different photon energies (red: soft band, green: medium band, blue: hard band). The image shows both the entire EWOCS field and a central region of approximately ~3′. In the right panel, it is evident that the source density is high in the cluster core and it reveals that the majority of faint sources are predominantly hard, likely due to high absorption or since they have been observed during periods of intense magnetic activity such as flares. The list of Chandra datasets used in this paper, and obtained by the Chandra X-ray Observatory, are contained in the Chandra Data Collection (CDC) 1535.

thumbnail Fig. 1

Contours of the pre-EWOCS and EWOCS Chandra observations of Westerlund 1 overlaid on the combined ACIS event file (left panel) on an image in the Ks band obtained with the FourStar infrared camera mounted on the Magellan 6.5 m telescopes (right panel).

Table 1

Pre-EWOCS observations of Westerlund 1.

3.2 Data reduction

The Chandra observations were analyzed using the “pre-ACIS Extract workflow” procedure outlined in Townsley et al. (2003) and Broos et al. (2010). This procedure utilizes various tools integrated within the Chandra Interactive Analysis of Observations (CIAO) software (Fruscione et al. 2006). We employed versions 4.13 of CIAO along with the CALDB 4.9.5 calibration files.

The L1-to-L2 processing flow aims to generate calibrated Chandra event files from the L1 products provided by the Chandra X-Ray Center (CXC). It includes event energy calibration, refinement of event positions, and correction for contamination caused by bad pixels and cosmic-ray afterglow. This workflow utilizes a less aggressive bad pixel table compared to the one produced by CIAO and it incorporates the clean55 algorithm for background reduction. Additionally, cosmic-ray afterglows are removed, and the source point spread function (PSF) is improved by disabling the random ±0.25 pixel randomization. The standard grade filter is applied to events, retaining only AS CA grades 0, 2, 3, 4, and 6. However, events are not filtered for the standard status=0 requirement, which may result in the exclusion of a significant number of reliable events.

Afterglows, which are groups of events appearing at the same location in consecutive CCD frames, can often be mistaken for faint sources. To address this, the CIAO tool acis_detect_afterglow is typically employed to remove afterglows. This tool applies a relatively aggressive cleaning approach, eliminating several false positives. Another tool, acis_run_hotpix, is less aggressive but it may fail to detect afterglow series with fewer than ten counts. In this L1-to-L2 procedure, a bifurcated workflow is adopted, where we applied an aggressive cleaning to the files used for source detection and validation, and a less aggressive cleaning for the files used in spectral analysis.

The background light curves were examined to identify and exclude intervals with intense and fluctuating background. This correction was required only for Obs.ID 5411, as the background remained relatively stable throughout the other observations.

The astrometry of the event files was corrected in three steps. In the first step, we addressed the offset of each Obs.ID relative to Obs.ID 22319, which is the deepest observation. We utilized Wavdetect to identify the brightest sources in each observation and cross-matched their positions with those detected in the Obs.ID 22319 image. Subsequently, we employed the CIAO tools wcs_match and wcs_update to update the astrometry for each observation. In the second step, which was part of the L1-to-L2 workflow, we corrected the astrometry of each event file using the Gaia Third Data Release (DR3; Gaia Collaboration 2023) astrometric system. This process was repeated as the third step, but using the brightest sources from the final list of validated sources (Sect. 6).

Exposure maps were calculated using the standard CIAO tools implemented in the pre-ACIS Extract workflow for each observation in the broad (0.5–7.9 keV), soft (0.5–1.0 keV), medium (1.0–2.0 keV), hard (2.0–7.9 keV), and very hard (4.0–7.9 keV) bands, and subsequently combined. The resulting combined exposure map in the broad band is displayed in Fig. 3, revealing a deep and nearly uniform exposure in the central region. This region is sufficiently large to encompass both the core of Westerlund 1 and a portion of the expected halo of the cluster (as recently discovered, extended haloes are typically associated with stellar clusters; Meingast et al. 2021; Prisinzano et al. 2022).

Table 2

EWOCS observations.

4 Source detection

The strategy we employed for source detection aims to maximize the depth of the EWOCS catalog, even in the core of Westerlund 1. This region presents challenges due to source confusion and a bright, irregular background, making the detection of faint sources a complex task (see Fig. 4).

Source detection is implemented using four different methods:

  • The wavelet-based algorithm PWDetect (Damiani et al. 1997) is applied in the broad, soft, medium, hard, and very hard energy bands. The detection threshold we adopted roughly corresponds to 50 spurious sources. We excluded the outermost regions where the selection resulted in a very large number of false positives, resulting in 2306 detected sources.

  • The wavelet-based algorithm Wavdetect (Freeman et al. 2002) is applied to images in the broad, soft, medium, hard, and very hard energy bands. We set the sigmathreshold parameter equal to 10−4 and used only two small detection scales, resulting in 2509 detected sources.

  • The maximum likelihood reconstruction method developed by Townsley et al. (2006) is applied in the broad, soft, hard, and very hard energy bands. This algorithm operates over small tiles across the observed field, making it more sensitive to the spatial variation in the PSF and background and thus more capable of detecting faint sources in crowded fields (Broos et al. 2010). The reconstructed image is first calculated using the Lucy-Richardson algorithm (Lucy 1974), and then searched for peaks that identify the positions of point sources. This method resulted in 7585 detected sources.

  • A time-resolved deployment of PWDetect, described in more details in the following, is performed over segments of 10 ksec of the observations. This method is aimed at detecting faint and variable sources that may only be significant during specific short time segments in which they were detected. This method produced a list of 1147 detected sources.

Figure 5 shows a comparison of the spatial distribution of candidate sources detected using the four methods. In all cases, the cluster appears highly crowded, with the image reconstruction method being the only one capable of detecting a large number of sources in the central region of Westerlund 1, as expected.

thumbnail Fig. 2

RGB images of the whole ACIS-I field (left panel) and the central area (right panel) of the composite Chandra images. Soft band (0.5–1.0 keV) photons are marked in red, medium band (1.0–2.0 keV) photons in green, and hard band (2–7.9 keV) photons in blue. The brightest source in the southeast direction is CXOU J16. The two images were smoothed adopting a Gaussian kernel with a radius of 2 pixels.

thumbnail Fig. 3

Combined exposure map in the broad band.

4.1 Time-resolved PWDetect

We devised a simple time-resolved detection method, tailored to faint transient sources such as magnetically flaring low-mass stars. These may remain undetected in the full dataset because of the high background, but may be detected in a shorter time slice that includes the transient emission, thanks to the enhanced source-counts/background contrast.

We started by considering 10 ks time slices from each observation segment. Since exposure times are not multiple of 10 ks, the exposure time of the last frame was forced to range between 7 ks and 17 ks. Moreover, in order to fully capture transients that would otherwise be split between two frames, we also considered intervals shifted in time by half a frame (5 ks).

The 44 EWOCS and pre-EWOCS observations were thus split into 194 frames: 106 were 10 ks long, while the duration of the remaining ones are quite uniformly distributed between 5 and 16 ks. For each frame, event lists (and exposure maps) were then extracted in the following five energy bands: 0.5–7.0 keV, 4.0–7.0 keV, 0.5–1.2 keV, 1.2–2.0 keV, and 2.0–7.0 keV, resulting in a total of 970 event files.

We ran PWDetect twice on each of these 970 event files, once to estimate the background level and thus evaluate the significance thresholds to obtain the desired number of spurious sources, and once more for the final detection. For the first run we adopted a low significance threshold, 4.9σ, so to detect as many sources as possible, but also resulting in several spurious faint sources. The number of background photons was then estimated by subtracting the detected source photons from the total; we then chose the final detection threshold so to yields, on average, 0.1 spurious sources per frame. This was derived from the appropriate significance versus background curve provided by Damiani et al. (1997).

The final PWDetect runs produced 944 lists of sources6 for a total of 14178 sources, most of which are, obviously, repeated detections of the same source in multiple frames and/or bands. We started cleaning up this large sample by removing ~1600 extended sources, many of which were unresolved detections of multiple point sources (extent parameter, as given by PWDe-tect, larger than 2). We then screened for the remaining sources for possible cosmic rays afterglow events: for each source we extracted photons from a circle with radius twice the “detection scale” provided by PWDetect. In 284 cases the arrival times of all extracted photons were in subsequent 3.14s-long readout frames, and the detection, a likely afterglow artifact was discarded.

All detection lists were cross-identified and merged in a final source list using an iterative procedure: first we cross-identified and merged the first two catalogs. The resulting catalog was then merged with the third original catalog, and so on for all the 944 catalogs. Identifications were performed searching the close spatial coincidences with identification radii of each original detection taken as the 1σ uncertainties as estimated by PWDetect (rounded up to 0.5 arcsec if smaller). The coordinates and uncertainties or identification radii of merged sources were computed, at each step, as the uncertainty-weighted means of the coordinates/radii of cross-identified sources7. At the end of this process we are left with 1262 cross-matched sources.

Finally, we inspected all the final sources by eye, examining individual detections in the original event file, and the positions in the Hubble Space Telescope (HST) H-band image (when available, in the field center) and images from the Digitized Sky Survey (DSS) and the Two Micron All-Sky Survey (2MASS). Some cross-identifications were adjusted and a number of “sources,” which were not merged by the automatic process above, where merged as they clearly referred to the same star. The final list counts 1147 sources.

thumbnail Fig. 4

Inner region of Westerlund 1 observed with ACIS (left panel) in the broad band and with HST (right panel) using the F160W filter. In the ACIS image, source confusion and a high background intensity dominate the cluster core.

4.2 The merged list of candidate sources

At this stage, we generated a list of candidate sources that includes all the sources detected using the adopted methods. Contamination of this list by false positives is expected to be large, but we relied on the source-validation step, described in the next sections, to prune the catalog from these false and not significant sources.

The lists generated by the four detection methods were crossmatched by eye. In cases where it was difficult to confidently determine the presence of one or more nearby sources, we left those entries unmatched between the catalogs. Additionally, we included in the input list the following sources that were not detected by any of the aforementioned methods:

  • 446 faint sources from the catalog presented by Townsley et al. (2018), which are likely not detected in EWOCS observations because of the intrinsic variability of young stars;

  • 21 massive stars of Westerlund 1 from the list published by Clark et al. (2020);

  • 47 candidate sources added by eye corresponding to the positions of Gaia sources in or nearby the cluster center.

    The final list of candidate sources, which was used as input for the source validation process in ACIS-Extract (AE), consists of 9420 sources.

thumbnail Fig. 5

Spatial distribution of candidate sources detected with the four methods (from the top: Pwdetect, Wavdetect, image reconstruction, and time-resolved Pwdetect). The left panels show the whole ACIS field, those on the right the inner region. Different colors in the first and second rows mark sources detected at different energy bands.

5 Source extraction, validation, and photometry

Source validation and photometry were performed using the AE software in IDL (Broos et al. 2010)8, which has been successfully employed in previous X-ray surveys including the Chandra Carina Complex Project (Townsley et al. 2011), the Massive Young Star-forming Complex Study in Infrared and X-Rays (MYStIX) survey (Feigelson et al. 2013), the three Massive Star-forming Regions Omnibus X-ray Catalog (MOCX) data releases (Townsley et al. 2014, 2018, 2019), the Chandra Cygnus OB2 Legacy Survey (Wright et al. 2014a), and the Star Formation In Nearby Clouds (SFiNCs) project (Getman et al. 2017). AE enables the extraction and validation of sources across multiple observations, generating individual source spectra and light curves. It utilizes various data analysis software packages including CIAO, MARX (Davis et al. 2012), HEASoft9, and the IDL Astronomy User’s Library (Landsman 1993).

Following the guidelines provided by the authors and available on the AE website, we adopted a three-step procedure to compile the X-ray EWOCS catalog:

  • Initially, sources were extracted and validated using a parameter defined by AE, which helps distinguish between genuine and false sources. This step was repeated iteratively until no more false sources were identified and removed (Sect. 5.1).

  • Subsequently, source positions were updated, followed by another round of source validation process (Sect. 5.2).

  • Once the catalog reached a stable state, we performed the photometric procedure, to extract source events comprehensively and calculate the primary spectral and temporal properties for each source across multiple energy bands (Sect. 5.3).

5.1 Source validation

The AE procedure assesses the local PSF at the given position of each source and defines extraction regions based on the 1.5 keV local PSF, ensuring they do not overlap with neighboring sources. In the case of close pairs, the extraction region of the fainter source is progressively reduced to prevent overlap until it reaches 40% of its original size. Once this threshold is reached, if the two extraction regions still overlap, AE further reduces the size of the brighter source until the regions no longer overlap. If overlap persists even when both extraction regions are reduced to 40%, AE either discards the specific observation or automatically removes the fainter source.

The local background is determined within an optimized region surrounding the source. For isolated sources, this region is delimited by an inner radius, which is 1.1 times the radius encompassing 99% of the PSF, and an outer radius large enough to collect at least 100 background events not associated with nearby sources. AE adjusts the size of the background-extraction region to ensure that Poissonian noise contributes no more than 3% to the background uncertainty. However, in crowded regions, defining a region with 100 events may not be feasible. In such cases, AE employs a different calculation that incorporates the contribution from nearby bright sources and a model accounting for the spatial variation of the background.

Source validation relies on a parameter provided by AE called prob_no_source (PB ), which represents the probability that there is no real source at a given position. In our case, where multiple observations of a source are available, AE calculates PB based on the extractions with the highest source significance. To differentiate between valid and spurious sources, we applied a threshold of PB=0.01, consistent with previous studies. Since the removal of not valid sources could potentially impact the extraction region and background of valid sources, the procedure is iterated until the catalog reaches convergence and no further spurious sources are detected.

After the first iteration, we conducted a visual inspection of sources flagged by AE as potentially resulting from the hook-shaped feature of the PSF10. This feature can account for up to 5% of the source flux and its position is influenced by the rollangle, making it distinguishable from the actual source only in a few cases where the real source is both on-axis and sufficiently bright. Figure 6 illustrates an example of a source (MOXC2) that was flagged by AE as a potential PSF hook and subsequently removed after visual inspection.

AE also identifies sources that are expected to suffer from significant pileup11 (which is the loss of information due to different incident photons registered as a unique event by the detector). In our case, the only source affected by piled source is the magnetar.

thumbnail Fig. 6

Example of a source (label MOXC2) that was excluded as a potential product of a PSF hook near a brighter source (label c10200). The green contours outline the extraction regions of MOXC2 in all observations, while the red polygons indicate the locations where the PSF hook may appear in each observation, depending on the roll-angle.

5.2 Positions update and visual review

For each source, AE calculates three different position estimates. The first estimate is obtained by taking the mean value of the positions of the events associated with the source (mean-data position). However, this estimate may be inaccurate for large off-axis angles and in cases where there are significant offsets between the true source position and the extraction region (which can happen when the PSF is asymmetric). To obtain a more accurate estimate in these cases, AE correlates the source PSF with the spatial distribution of extracted events (PSF position). This calculation takes into account the combination of several Obs.IDs by using the PSF calculated in each observation. Both of these estimates can be influenced by nearby sources. In crowded fields, a third estimate is provided by AE using the reconstructed image of the source’s neighborhood. It identifies the position of the closest peak in the reconstructed image. According to AE’s recommendations, the mean-data position is used for on-axis sources, the PSF position is used for off-axis sources, and the image reconstruction position is used for sources in crowded regions. The repositioning of sources was performed twice, with each step followed by a new sequence of iterations for source validation, as described in the previous section.

Before conducting the visual review of validated sources, the astrometry of both the X-ray sources and the main products file was corrected using the Gaia/DR3 astrometric system. After this step, and when catalog stability was achieved again, we conducted a visual review of specific critical sources, including very faint sources that could affect the size of the extraction region of nearby bright sources and suspected afterglows. The decisions made during the visual review were guided also by the presence of high-probability optical and/or infrared counterparts. After the visual review, a new round of source validation was performed.

thumbnail Fig. 7

Extraction regions of the validated sources across the entire merged ACIS image (left panel), and extraction regions in the central area of approximately 3′ in size (right panel).

5.3 Spectral extraction

After 21 iterations of the source validation process, the catalog reached stability, with a total of 5963 validated X-ray sources. The final step involved the extraction of X-ray events and the estimation of X-ray properties in 17 energy bands, merging all available observations in a consistent manner. In addition, AE generates light curves and spectra for each source, although these will not be discussed in this paper. AE performs this calculation by excluding observations where the sources are observed off-axis to improve the overall signal-to-noise ratio. However, in our case, this correction was not necessary due to the design of our survey. The calculated quantities include source counts, net counts, photon flux in photon cm−2 s−1, and the quartiles of photons energy.

Figure 7 depicts the spatial distribution of the validated sources in the merged ACIS event files. In the left panel it is evident that there is a high concentration of validated sources toward the center of the cluster, as well as a significant number of sources surrounding the cluster core. This indicates that we have detected stars associated with the extended halo of Westerlund 1. This will be further investigated in upcoming papers of this series, which will focus on source classification and the identification of optical/infrared (OIR) counterparts. The right panel also highlights how this survey has pushed to the limits of Chandra in resolving individual stars within such a densely populated stellar cluster with a bright and irregular background. In fact, in the actual core of the cluster, where the background is both intense and variable, a few tens of sources that were initially included as input to AE were subsequently discarded during the validation process (see Fig. 8). The limited number of validated sources in the central region can be attributed to the intense background, and it is likely that many of these discarded sources are indeed genuine X-ray sources. Although we did not attempt to recover these stars, in future papers of this series their candidate OIR counterparts will be analyzed in order to estimate the fraction of real sources that we have excluded.

6 The final catalog

It is informative to analyze the number of sources detected using the various methods we employed and assess how many have survived the pruning process. Table 3 presents the total number of input sources for each detection method, as well as the fraction of these sources within 1″ and 3″ of a source in the final catalog (source positions changed during the pruning process and thus an exact position match was not possible). The image reconstruction method is the only one that experienced significant pruning of the input catalog, as it selects sources that are too faint according to the adopted PB threshold. According to Table 3, and considering also that the number of sources in the input PWDetect list of candidate sources more distant than 3″ from any source in the Image reconstruction input catalog is low, but not negligible (314, 688 for Wavdetect), it is evident that in complex fields like these, deploying different detection methods is crucial for optimizing the number of detected sources.

Given the design of the EWOCS survey and the compact nature of Westerlund 1, it is not surprising that the majority of sources are observed at low off-axis angles, as depicted in Fig. 9. Specifically, 63.7% of the sources (3485/5464) are located within 1 arcmin from the field center, and 87.7% are within 3 arcmin. Consequently, source positions are generally well-determined, with a median position error of 0.17″ and a 75% quantile position error of 0.27″. Position errors are estimated from the single-axis standard deviations of the PSF inside the extraction region and the number of counts extracted. This precision is crucial for the search of OIR counterparts and for dynamics studies.

As depicted in Fig. 10, the catalog is predominantly composed of faint sources. In the broad band, the median value of the net counts is 12.9 counts. There are 607 sources (10.2%) with fewer than 5 net counts and only 69 sources (1.2%) with fewer than 3 counts. It is well known that the sensitivity of ACIS-I decreases with the off-axis angle. This must be taken in consideration when comparing the spatial distribution of X-ray sources detected with ACIS-I to those detected with other instruments. Figure 11 illustrates the spatial distributions of EWOCS X-ray sources with fewer than 12.9 net counts and those with more net counts. The former sample exhibits a higher concentration in the center of the field, with only 37 sources having an off-axis angle larger than 7′. This region is considerably large compared to the size of Westerlund 1, so studies based on the spatial distribution of cluster members would not be significantly affected by the decline in sensitivity with the off-axis angle.

Given the design of the EWOCS observations and the intricate procedure we employed for source detection and validation, it is not currently feasible to provide a reliable estimate of catalog completeness without making strong and unverified assumptions about cluster properties, its morphology, and both mass and LX distributions. Instead, we prefer to discuss the achieved completeness in future papers of this series, once the identification of OIR counterparts and the determination of true cluster members have been accomplished. In Appendix C, however, we present a simplified analysis of completeness based on different assumptions regarding cluster morphology, along with simulations conducted using the MARX simulator.

The distribution of the median photon energy for the EWOCS X-ray sources is shown in Fig. 12. The median value of the distribution is 2.8 keV. In the case of young stellar populations in clusters with low extinction, the median photon energy serves as a reliable indicator of membership since the coronal plasma temperature in young stars is typically higher than in older stars. However, since interstellar absorption is significant in the direction of Westerlund 1, it becomes challenging to differentiate between the absorbed background population and the young stars within the cluster. However, the secondary peak observed at energies below 2 keV in the Emed distribution could potentially be attributed to a foreground population. Nevertheless, the spatial distribution of these soft sources does not differ significantly from that of the more energetic sources. Additionally, Fig. 12 includes a comparison between the photon median energy distribution of the EWOCS X-ray sources and the X-ray catalog published by Townsley et al. (2018) based on pre-EWOCS observations, which clearly exhibits a peak below 2 keV. The evident differences in the two distributions can be due to a combination of factors: the presence of a large population of stars associated with Westerlund 1 in the EWOCS catalog (which is a factor of ~5 deeper in X-ray photon flux compared to the catalog published by Townsley et al. 2018, considering the faintest sources in the two catalogs), as well as the decline in sensitivity in the soft band of the ACIS detector with the years, and the better sensitivity toward soft events of the ACIS-S detector compared with ACIS-I.

The catalog also includes a measure of source flux provided by AE: the photon flux (Fphotons), which is calculated as the ratio of the source net counts to the product of the mean effective area and nominal exposure time (thus expressed in units of photons cm−2 s−1). A model-independent estimate of the apparent source energy flux can be calculated as 1.602×10−9 × Emed × Fphotons. The coefficient is derived from the conversion between keV and erg, as determined by Getman et al. (2010).

In Appendix A, we show ten rows of the X-ray EWOCS catalog, which is available in full at the CDS. We have also made the output table produced by AE available on the EWOCS website in its original IDL format12.

thumbnail Fig. 8

Extraction regions of the validated sources and positions of the input candidate sources (crosses) within the central 1 arcmin region

Table 3

EWOCS sources and detection methods.

thumbnail Fig. 9

Distribution of the off-axis angles of the EWOCS X-ray sources.

thumbnail Fig. 10

Distributions of source net counts in the broad, soft, medium, hard, and very hard energy bands. The number of sources with more than 50 counts is indicated in the top-right corner of each panel.

thumbnail Fig. 11

Spatial distribution of EWOCS X-ray sources, sorted into two bins based on their net counts in the broad band.

thumbnail Fig. 12

Distribution of the median photon energy in the broad band for the validated EWOCS X-ray sources (black) and the catalog published by Townsley et al. (2018), in red.

7 Specific sources

Not surprisingly, CXOU J 16 is the brightest source in the EWOCS X-ray catalog, with a total of 71601 ± 268 net counts collected and a photon flux of 3.45 × 10−4 photons cm−2 s−1. The source, which is strongly piled-up, is quite isolated and produces a surrounding bright background, with the closest source being at about 8 arcsec. The pulsar will be analyzed in detail in future papers of this project.

As explained in Sect. 2, Westerlund 1 hosts a unique ensemble of massive stars caught in different evolutionary stages. Understanding the mechanisms responsible for the emission of X-rays and studying both binarity and the circumstellar environment in these stars is of primary importance for EWOCS. We visually inspected the X-ray counterparts of massive stars published by Clark et al. (2020) and found 126 coincidences out of the 166 listed massive stars. The results are listed in Table B.1. In the vast majority of cases, there was a clear one-to-one correspondence between the sources in the two catalogs. There are a few uncertain cases, which can be easily identified by repeated massive star IDs, EWOCS objects, or large separations.

The brightest X-ray massive star in the EWOCS catalog is the SgB[e] star W9 (Clark et al. 2014), with nearly 8000 net counts collected in the broad band. The intense X-ray brightness (LX ~ 3.6 × 1033 erg s−1) and the hardness of the spectrum, as previously reported by Clark et al. (2008), are consistent with the evidence of intense mass loss rate, estimated to be around 10−5 M yr−1 (Andrews et al. 2019), and strong indications of binarity (Ritchie et al. 2022). W9 is the brightest source in the cluster also at radio wavelengths (Andrews et al. 2019), millimeter band (Fenech et al. 2018), and it shows very bright mid-IR emission (Clark et al. 1998).

In terms of X-ray luminosity, W9 is followed by the post-binary blue straggler W30 (04-5Ia+; Clark et al. 2019, 2008), for which a putative orbital period of approximately 6.2 days has been identified from a radial velocity series analyzed by Ritchie et al. (2022). We have detected nearly 6000 net counts in the broad band for W30. After W9 and W30, the list of sources with net counts ranging between 190 and 5400 photons includes most of the known WR stars and OB supergiant binary systems.

The deep EWOCS observations provide, for the first time, candidate X-ray detections for some normal giant and subgiant stars, such as W50b and W1051. The lack of detection in the pre-EWOCS observations has been explained as a natural consequence of the lower intrinsic bolometric luminosity of these stars compared to more evolved massive stars in the cluster (Clark et al. 2019). The detection in the EWOCS observations supports this hypothesis. Faint X-ray counterparts have also been found for the two YHGs W4 and W8 (with 22 and 73 net counts, respectively), whose nature has been recently discussed by Beasor et al. (2023), who classified them as yellow super-giants. Additionally, a faint counterpart has been detected for the BHG W1049 (with 14.168.420.4$14.16_{8.4}^{20.4}$ net counts). For the first time, faint counterparts have been found for the four O9.5II SB1 stars W1022, W1050, W1056, and W1060, as well as for the B1.5II star W1048. We also confirm the relatively faint (75.7565.885.7$75.75_{65.8}^{85.7}$ net counts) and soft (median photon energy of 1.9 keV) X-ray emission from the SB2 star (B0.5I+OB) W10, as previously reported by Clark et al. (2008). Ritchie et al. (2022) attributed the X-ray properties of this star to the possibility that the pre-EWOCS observations were made at a phase where the wind collision zone was weak or obscured. However, given the length of the EWOCS observations, it is more likely that these properties are intrinsic to the star.

It is quite interesting that out of the 124 sources in our catalog with more than 100 net counts, 94 do not readily match any known massive stars in the cluster. This subset will be studied in detail in future works of this series to determine their nature. It is intriguing that this sample does not seem to follow the distribution of photon energy shown in Fig. 12. In fact, its photon median energy distribution exhibits three distinct peaks: one below 2 keV (which may be dominated by foreground stars), one between 2.5 keV and 3 keV (compatible with cluster stars), and one between 3.8 keV and 4.3 keV (which could be influenced by background sources or flaring low-mass cluster members).

We also provide a list of positions for the unvalidated candidate X-ray sources on the EWOCS website13. This list will be cross-matched with existing optical and infrared catalogs of Westerlund 1 to determine the fraction of rejected sources that could potentially be true counterparts of cluster members. Likewise, identifying optical and infrared counterparts will enable us to assess the level of contamination and the fraction of expected spurious sources in the EWOCS X-ray source catalog, as well as determine the completeness limit achieved by our survey.

8 Conclusions

In this paper, we present the EWOCS project and a new list of X-ray sources in the young supermassive star cluster Westerlund 1 and its surrounding area. The EWOCS project aims to investigate the impact of the starburst environment on the formation process of stars and planets, the dispersal of protoplanetary disks, and the evolutionary pathway of massive stars.

Here we present the 1 Msec Chandra/ACIS-I EWOCS observations of Westerlund 1, the workflow for data reduction, the procedure for source detection and validation, and the spectral extraction of the validated sources. Initially, we generated a preliminary list of 9420 candidate X-ray sources using the image reconstruction method, PWDetect, WAVDETECT, and a specific deployment of PWDetect focused on identifying flaring stars that exhibited a significant signal above the background for a brief duration. Additionally, a few sources were manually added or obtained from existing catalogs of Westerlund 1 sources. From these input sources, we compiled the EWOCS catalog of X-ray sources in Westerlund 1 of 5963 sources successfully validated using the IDL-based software AE.

The median value of net counts in the EWOCS X-ray catalog is approximately 13 counts, with about 10% of sources having fewer than 5 net counts detected in the broad energy band. The distribution of the median photon energy of the sources peaks at approximately 2.8 keV, with a contribution from unrelated (foreground and background) sources that is challenging to distinguish from the candidate cluster members. The brightest source in the catalog is the magnetar CXO J164710.2-455216, with over 70 000 net counts detected in the broad band. It is followed by several massive stars in Westerlund 1, including the SgB[e] star W9, the post-binary blue straggler W30, and some WR stars and supergiants in binary systems. Out of the 166 known very massive stars in Westerlund 1, we have identified a reliable X-ray counterpart for 126 of them. Additionally, we have made the first detection of an extended and rich halo surrounding the core of Westerlund 1, which will be crucial in assessing the cluster’s true mass content, formation, and evolution.

Acknowledgements

We acknowledge the referee for his/her careful reading of our paper and suggestions. M.G.G. is also indebted to Leisa K. Townsley and Patrick S. Broos for their valuable suggestions and advice on data reduction and ACIS extraction. M.G.G., C.A., R.B., E.F., G.L.I., L.P., and S.S. acknowledge the INAF grant 1.05.12.05.03. K.M. acknowledges support from the Fundação para a Ciência e a Tecnologia (FCT) through the CEEC-individual contract 2022.03809.CEECIND and research grants UIDB/04434/2020 and UIDP/04434/2020. Support for this work was also provided by the National Aeronautics and Space Administration through Chandra Proposal 21200267 issued by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060. I.N. is partially supported by the Spanish Government Ministerio de Ciencia e Innovación (MCIN) and Agencia Estatal de Investigación (MCIN/AEI/10.130 39/501 100 011 033/FEDER, UE) under grant PID2021-122397NB-C22, and also by MCIN with funding from the European Union NextGenerationEU and Generalitat Valenciana in the call Pro-grama de Planes Complementarios de I + D + i (PRTR 2022), project HIAMAS, reference ASFAE/2022/017. M. G. G. and R. B. acknowledge partial support from the Grant INAF 2022 YODA; R. B. also acknowledges partial support from the project PRIN-INAF 2019 “Spectroscopically Tracing the Disk Dispersal Evolution”. The scientific results reported in this article are based on observations made by the Chandra X-ray Observatory. M.M. acknowledges financial support from the European Research Council for the ERC Consolidator grant DEMOBLACK, under contract no. 770017, and from the German Excellence Strategy via the Heidelberg Cluster of Excellence (EXC 2181 – 390900948) STRUCTURES.

Appendix A Extract of the EWOCS X-ray Westerlund 1 sources catalog

Ten rows of the EWOCS catalog of the X-ray sources in Westerlund 1 are shown here. The full catalog is available at the CDS.

Table A.1

Ten rows extracted from the EWOCS catalog of the X-ray sources in Westerlund 1

Appendix B EWOCS X-ray counterparts of the massive stars in Westerlund 1

Table B.1 shows the EWOCS X-ray counterparts of the massive stars in Westerlund 1 listed in Clark et al. (2020).

Table B.1

Known massive stars in the EWOCS X-ray catalog

Appendix C Estimate of catalog completeness

The resulting completeness of our survey depends not only on the total exposure, but also on source crowding and the bright and irregular background. A full understanding of completeness will only be possible after the identification and classification of the OIR counterparts of the X-ray sources in order to distinguish between cluster members and sources in the foreground and background. However, we conducted some simple simulations using MARX, which, despite being based on strong assumptions, can provide some hints about completeness.

In order to simulate the cluster population, since the true shape of the IMF of Westerlund 1 is still a subject of debate, particularly in the low-mass regime, we made the assumption that the cluster IMF follows the law proposed by Kroupa (2001), which is applicable to most known young stellar clusters. We understand that the starburst environment can influence the distribution of stellar masses, leading to different mass functions. However, at this level of approximation, this is considered a secondary effect. To accommodate the compact morphology of the cluster, we assumed that cluster members are distributed according to a Gaussian function with a full width at half maximum of 4 arcminutes. Therefore, we did not account for the asymmetric morphology of Westerlund 1, as suggested by previous authors (e.g., Gennaro et al. 2011). Additionally, we assumed a total cluster mass of 45000 solar masses, encompassing stars with masses as low as 0.08 solar masses.

To convert the mass distribution into an LX distribution, we utilized the LX versus mass distribution derived from the Chandra Orion Ultradeep Project (COUP) conducted in the Orion Nebula Cluster (Preibisch & Feigelson 2005), accounting for its observed spread. We chose this distribution because the COUP survey provides the most complete X-ray observation of a young stellar cluster. However, it should be noted that this distribution may not accurately represent the population of Westerlund 1 due to differences in age and the presence of a distinct massive stellar population in this cluster. To account for this massive stellar population, we simply added the massive sources identified by Clark et al. (2005) with their corresponding measured LX values to the simulated cluster population. Additionally, we normalized the COUP LX versus mass distribution to account for the decline in stellar X-ray luminosity with age (Preibisch & Feigelson 2005), and we used the specified values for cluster distance and absorption to convert luminosity into flux.

We simulated a 1 Msec ACIS-I observation of this fake cluster, taking into account instrumental background14, and performed source detection using Wavdetect (thus not accounting for the source validation procedure we adopted with AE). By comparing the input and output lists of sources, we determined that the completeness in the 0.8–2 solar mass range is approximately 40% within the central 4 arcminute region, decreasing by approximately 10% in the inner 1 arcminute region. For more massive stars, the estimated completeness is around 85% regardless of the distance from the cluster center. It is important to note that this is a preliminary estimation of the completeness of the EWOCS X-ray catalog, which will be further validated through the identification of OIR counterparts and source classification.

References

  1. Adamo, A., Zeidler, P., Kruijssen, J. M. D., et al. 2020, Space Sci. Rev., 216, 69 [Google Scholar]
  2. Aghakhanloo, M., Murphy, J. W., Smith, N., et al. 2020, MNRAS, 492, 2497 [NASA ADS] [CrossRef] [Google Scholar]
  3. An, H., Kaspi, V. M., Archibald, R., & Cumming, A. 2013, ApJ, 763, 82 [NASA ADS] [CrossRef] [Google Scholar]
  4. Andersen, M., Gennaro, M., Brandner, W., et al. 2017, A&A, 602, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Andrews, H., Fenech, D., Prinja, R. K., Clark, J. S., & Hindson, L. 2019, A&A, 632, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Barrado, D., Stelzer, B., Morales-Calderón, M., et al. 2011, A&A, 526, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Beasor, E. R., Davies, B., Smith, N., Gehrz, R. D., & Figer, D. F. 2021, ApJ, 912, 16 [NASA ADS] [CrossRef] [Google Scholar]
  8. Beasor, E. R., Smith, N., & Andrews, J. E. 2023, ApJ, 952, 113 [NASA ADS] [CrossRef] [Google Scholar]
  9. Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953 [Google Scholar]
  10. Berghoefer, T. W., Schmitt, J. H. M. M., & Cassinelli, J. P. 1996, A&AS, 118, 481 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Borghese, A., Rea, N., Turolla, R., et al. 2019, MNRAS, 484, 2931 [NASA ADS] [CrossRef] [Google Scholar]
  12. Brandner, W., Clark, J. S., Stolte, A., et al. 2008, A&A, 478, 137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Broos, P. S., Townsley, L. K., Feigelson, E. D., et al. 2010, ApJ, 714, 1582 [NASA ADS] [CrossRef] [Google Scholar]
  14. Clark, J. S., & Negueruela, I. 2002, A&A, 396, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Clark, J. S., & Negueruela, I. 2004, A&A, 413, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Clark, J. S., Fender, R. P., Waters, L. B. F. M., et al. 1998, MNRAS, 299, L43 [NASA ADS] [CrossRef] [Google Scholar]
  17. Clark, J. S., Negueruela, I., Crowther, P. A., & Goodwin, S. P. 2005, A&A, 434, 949 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Clark, J. S., Muno, M. P., Negueruela, I., et al. 2008, A&A, 477, 147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Clark, J. S., Negueruela, I., Davies, B., et al. 2009, A&A, 498, 109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Clark, J. S., Negueruela, I., & González-Fernández, C. 2014, A&A, 561, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Clark, J. S., Najarro, F., Negueruela, I., et al. 2019, A&A, 623, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Clark, J. S., Ritchie, B. W., & Negueruela, I. 2020, A&A, 635, A187 [EDP Sciences] [Google Scholar]
  23. Crowther, P. A., Hadfield, L. J., Clark, J. S., Negueruela, I., & Vacca, W. D. 2006, MNRAS, 372, 1407 [Google Scholar]
  24. D’Ai, A., Evans, P. A., Krimm, H. A., et al. 2017, GRB Coordinates Network, 21095, 1 [Google Scholar]
  25. Damiani, F., Maggio, A., Micela, G., & Sciortino, S. 1997, ApJ, 483, 350 [NASA ADS] [CrossRef] [Google Scholar]
  26. Damineli, A., Almeida, L. A., Blum, R. D., et al. 2016, MNRAS, 463, 2653 [NASA ADS] [CrossRef] [Google Scholar]
  27. Davies, B., & Beasor, E. R. 2019, MNRAS, 486, L10 [NASA ADS] [CrossRef] [Google Scholar]
  28. Davies, B., Figer, D. F., Kudritzki, R.-P., et al. 2007, ApJ, 671, 781 [NASA ADS] [CrossRef] [Google Scholar]
  29. Davies, B., de La Fuente, D., Najarro, F., et al. 2012, MNRAS, 419, 1860 [NASA ADS] [CrossRef] [Google Scholar]
  30. Davis, J. E., Bautz, M. W., Dewey, D., et al. 2012, SPIE Conf. Ser., 8443, 84431A [NASA ADS] [Google Scholar]
  31. de Grijs, R., O’Connell, R. W., & Gallagher, John S., I. 2001, AJ, 121, 768 [NASA ADS] [CrossRef] [Google Scholar]
  32. de Grijs, R., Anders, P., Bastian, N., et al. 2003a, MNRAS, 343, 1285 [NASA ADS] [CrossRef] [Google Scholar]
  33. de Grijs, R., Fritze-v Alvensleben, U., Anders, P., et al. 2003b, MNRAS, 342, 259 [NASA ADS] [CrossRef] [Google Scholar]
  34. de la Fuente, D., Najarro, F., Borissova, J., et al. 2016, A&A, 589, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Dougherty, S. M., Clark, J. S., Negueruela, I., Johnson, T., & Chapman, J. M. 2010, A&A, 511, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Dunlop, J. S. 2011, Science, 333, 178 [NASA ADS] [CrossRef] [Google Scholar]
  37. Eufrasio, R. T., Lehmer, B. D., Zezas, A., et al. 2017, ApJ, 851, 10 [NASA ADS] [CrossRef] [Google Scholar]
  38. Feigelson, E. D., Townsley, L. K., Broos, P. S., et al. 2013, ApJS, 209, 26 [Google Scholar]
  39. Fenech, D. M., Clark, J. S., Prinja, R. K., et al. 2018, A&A, 617, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Figer, D. F. 2008, in Massive Stars as Cosmic Engines, 250, eds. F. Bresolin, P. A. Crowther, & J. Puls, 247 [NASA ADS] [Google Scholar]
  41. Figer, D. F., McLean, I. S., & Morris, M. 1999, ApJ, 514, 202 [Google Scholar]
  42. Figer, D. F., Najarro, F., Gilmore, D., et al. 2002, ApJ, 581, 258 [Google Scholar]
  43. Figer, D. F., MacKenty, J. W., Robberto, M., et al. 2006, ApJ, 643, 1166 [Google Scholar]
  44. Freeman, P. E., Kashyap, V., Rosner, R., & Lamb, D. Q. 2002, ApJS, 138, 185 [Google Scholar]
  45. Fruscione, A., McDowell, J. C., Allen, G. E., et al. 2006, Proc. SPIE, 6270, 62701V [Google Scholar]
  46. Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Garmire, G. P., Bautz, M. W., Ford, P. G., Nousek, J. A., & Ricker, Jr., G. R. 2003, SPIE Conf. Ser., 4851, 28 [NASA ADS] [Google Scholar]
  49. Gennaro, M., Brandner, W., Stolte, A., & Henning, T. 2011, MNRAS, 412, 2469 [NASA ADS] [CrossRef] [Google Scholar]
  50. Getman, K. V., Flaccomio, E., Broos, P. S., et al. 2005, ApJS, 160, 319 [NASA ADS] [CrossRef] [Google Scholar]
  51. Getman, K. V., Feigelson, E. D., Broos, P. S., Townsley, L. K., & Garmire, G. P. 2010, ApJ, 708, 1760 [NASA ADS] [CrossRef] [Google Scholar]
  52. Getman, K. V., Broos, P. S., Kuhn, M. A., et al. 2017, ApJS, 229, 28 [Google Scholar]
  53. Groh, J. H., Damineli, A., Teodoro, M., & Barbosa, C. L. 2006, A&A, 457, 591 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Hachisuka, K., Brunthaler, A., Menten, K. M., et al. 2006, ApJ, 645, 337 [NASA ADS] [CrossRef] [Google Scholar]
  55. Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 [NASA ADS] [CrossRef] [Google Scholar]
  56. Israel, G. L., Campana, S., Dall'Osso, S., et al. 2007, ApJ, 664, 448 [NASA ADS] [CrossRef] [Google Scholar]
  57. Israel, G. L., Esposito, P., & Rea, N. 2011, ATel, 3653, 1 [NASA ADS] [Google Scholar]
  58. Kavanagh, P. J., Norci, L., & Meurs, E. J. A. 2011, New A, 16, 461 [NASA ADS] [CrossRef] [Google Scholar]
  59. Kothes, R., & Dougherty, S. M. 2007, A&A, 468, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Koumpia, E., & Bonanos, A. Z. 2012, A&A, 547, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Krimm, H., Barthelmy, S., Campana, S., et al. 2006, GRB Coordinates Network, 5581, 1 [NASA ADS] [Google Scholar]
  62. Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
  63. Kudryavtseva, N., Brandner, W., Gennaro, M., et al. 2012, ApJ, 750, L44 [NASA ADS] [CrossRef] [Google Scholar]
  64. Landsman, W. B. 1993, in ASP Conf. Ser., 52, Astronomical Data Analysis Software and Systems II, eds. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes, 246 [NASA ADS] [Google Scholar]
  65. Larson, R. B., & Tinsley, B. M. 1978, ApJ, 219, 46 [Google Scholar]
  66. Licquia, T. C., & Newman, J. A. 2015, ApJ, 806, 96 [NASA ADS] [CrossRef] [Google Scholar]
  67. Lim, B., Chun, M.-Y., Sung, H., et al. 2013, AJ, 145, 46 [NASA ADS] [CrossRef] [Google Scholar]
  68. Lucy, L. B. 1974, AJ, 79, 745 [Google Scholar]
  69. Mackey, J., Castro, N., Fossati, L., & Langer, N. 2015, A&A, 582, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 [Google Scholar]
  71. Mahy, L., Lanthermann, C., Hutsemékers, D., et al. 2022, A&A, 657, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Meingast, S., Alves, J., & Rottensteiner, A. 2021, A&A, 645, A84 [EDP Sciences] [Google Scholar]
  73. Melena, N. W., Massey, P., Morrell, N. I., & Zangari, A. M. 2008, AJ, 135, 878 [NASA ADS] [CrossRef] [Google Scholar]
  74. Montmerle, T. 1996, in ASP Conf. Ser., 109, Cool Stars, Stellar Systems, and the Sun, eds. R. Pallavicini, & A. K. Dupree, 405 [NASA ADS] [Google Scholar]
  75. Muno, M. P., Clark, J. S., Crowther, P. A., et al. 2006, ApJ, 636, L41 [NASA ADS] [CrossRef] [Google Scholar]
  76. Navarete, F., Damineli, A., Ramirez, A. E., Rocha, D. F., & Almeida, L. A. 2022, MNRAS, 516, 1289 [NASA ADS] [CrossRef] [Google Scholar]
  77. Negueruela, I., & Clark, J. S. 2005, A&A, 436, 541 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Negueruela, I., Clark, J. S., & Ritchie, B. W. 2010, A&A, 516, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Negueruela, I., Alfaro, E. J., Dorda, R., et al. 2022, A&A, 664, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Perna, R., & Pons, J. A. 2011, ApJ, 727, L51 [NASA ADS] [CrossRef] [Google Scholar]
  81. Piatti, A. E., Bica, E., & Claria, J. J. 1998, A&AS, 127, 423 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Pineda, J. L., Fischer, C., Kapala, M., et al. 2018, ApJ, 869, L30 [NASA ADS] [CrossRef] [Google Scholar]
  83. Portegies Zwart, S. F., Baumgardt, H., Hut, P., Makino, J., & McMillan, S. L. W. 2004, Nature, 428, 724 [Google Scholar]
  84. Preibisch, T., & Feigelson, E. D. 2005, ApJS, 160, 390 [NASA ADS] [CrossRef] [Google Scholar]
  85. Preibisch, T., Kim, Y.-C., Favata, F., et al. 2005, ApJS, 160, 401 [NASA ADS] [CrossRef] [Google Scholar]
  86. Prisinzano, L., Damiani, F., Sciortino, S., et al. 2022, A&A, 664, A175 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Rieke, G. H., & Rujopakarn, W. 2011, in ASP Conf. Ser., 446, Galaxy Evolution: Infrared to Millimeter Wavelength Perspective, eds. W. Wang, J. Lu, Z. Luo, et al. 3 [NASA ADS] [Google Scholar]
  88. Ritchie, B. W., Clark, J. S., Negueruela, I., & Crowther, P. A. 2009, A&A, 507, 1585 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Ritchie, B. W., Clark, J. S., Negueruela, I., & Langer, N. 2010, A&A, 520, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Ritchie, B. W., Clark, J. S., Negueruela, I., & Najarro, F. 2022, A&A, 660, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Robitaille, T. P., & Whitney, B. A. 2010, ApJ, 710, L11 [NASA ADS] [CrossRef] [Google Scholar]
  92. Rygl, K. L. J., Brunthaler, A., Sanna, A., et al. 2012, A&A, 539, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
  94. Seward, F. D., Forman, W. R., Giacconi, R., et al. 1979, ApJ, 234, L55 [NASA ADS] [CrossRef] [Google Scholar]
  95. Skinner, S. L., Simmons, A. E., Zhekov, S. A., et al. 2006, ApJ, 639, L35 [NASA ADS] [CrossRef] [Google Scholar]
  96. Smith, B. J., & Struck, C. 2010, AJ, 140, 1975 [NASA ADS] [CrossRef] [Google Scholar]
  97. Smith, N., Vink, J. S., & de Koter, A. 2004, ApJ, 615, 475 [NASA ADS] [CrossRef] [Google Scholar]
  98. Townsley, L. K., Feigelson, E. D., Montmerle, T., et al. 2003, ApJ, 593, 874 [NASA ADS] [CrossRef] [Google Scholar]
  99. Townsley, L. K., Broos, P. S., Feigelson, E. D., et al. 2006, AJ, 131, 2140 [Google Scholar]
  100. Townsley, L. K., Broos, P. S., Corcoran, M. F., et al. 2011, ApJS, 194, 1 [Google Scholar]
  101. Townsley, L. K., Broos, P. S., Garmire, G. P., et al. 2014, ApJS, 213, 1 [NASA ADS] [CrossRef] [Google Scholar]
  102. Townsley, L. K., Broos, P. S., Garmire, G. P., et al. 2018, ApJS, 235, 43 [NASA ADS] [CrossRef] [Google Scholar]
  103. Townsley, L. K., Broos, P. S., Garmire, G. P., & Povich, M. S. 2019, ApJS, 244, 28 [NASA ADS] [CrossRef] [Google Scholar]
  104. van Leeuwen, F., de Bruijne, J., Babusiaux, C., et al. 2021, Gaia EDR3 documentation, European Space Agency; Gaia Data Processing and Analysis Consortium. Online at https://gea.esac.esa.int/archive/documentation/GEDR3/index.html [Google Scholar]
  105. Vargas Álvarez, C. A., Kobulnicky, H. A., Bradley, D. R., et al. 2013, AJ, 145, 125 [CrossRef] [Google Scholar]
  106. Weisskopf, M. C., Brinkman, B., Canizares, C., et al. 2002, PASP, 114, 1 [NASA ADS] [CrossRef] [Google Scholar]
  107. Westerlund, B. 1961, PASP, 73, 51 [NASA ADS] [CrossRef] [Google Scholar]
  108. Westerlund, B. E. 1968, ApJ, 154, L67 [NASA ADS] [CrossRef] [Google Scholar]
  109. Wright, N. J., Drake, J. J., Guarcello, M. G., et al. 2014a, ArXiv e-prints [arXiv:1408.6579] [Google Scholar]
  110. Wright, N. J., Wesson, R., Drew, J. E., et al. 2014b, MNRAS, 437, L1 [NASA ADS] [CrossRef] [Google Scholar]

1

A different classification for six YHGs has recently been presented by Beasor et al. (2023).

2

As stated by Townsley et al. (2018), the corresponding X-ray flux has been calculated using PIMMS6 assuming a limit of five-counts detection on-axis, for a source with an APEC thermal plasma with kT=2.7 keV and abundance 0.4 × Z, which are typical values for a PMS star (Preibisch et al. 2005).

6

For the observations in standard ACIS-I configuration, detection was performed only on the most on-axis CCD (CCD.ID=7). In 26 cases PWdetect crashed or no sources were found. We did not investigate these cases further.

7

Since most detections are not independent (they may share the same photons because of energy band or overlapping time frames), we computed weighted means only among values belonging to independent groups of positions. Within each group of dependent detections we chose the coordinates and radii of the source with the smallest positional uncertainty.

All Tables

Table 1

Pre-EWOCS observations of Westerlund 1.

Table 2

EWOCS observations.

Table 3

EWOCS sources and detection methods.

Table A.1

Ten rows extracted from the EWOCS catalog of the X-ray sources in Westerlund 1

Table B.1

Known massive stars in the EWOCS X-ray catalog

All Figures

thumbnail Fig. 1

Contours of the pre-EWOCS and EWOCS Chandra observations of Westerlund 1 overlaid on the combined ACIS event file (left panel) on an image in the Ks band obtained with the FourStar infrared camera mounted on the Magellan 6.5 m telescopes (right panel).

In the text
thumbnail Fig. 2

RGB images of the whole ACIS-I field (left panel) and the central area (right panel) of the composite Chandra images. Soft band (0.5–1.0 keV) photons are marked in red, medium band (1.0–2.0 keV) photons in green, and hard band (2–7.9 keV) photons in blue. The brightest source in the southeast direction is CXOU J16. The two images were smoothed adopting a Gaussian kernel with a radius of 2 pixels.

In the text
thumbnail Fig. 3

Combined exposure map in the broad band.

In the text
thumbnail Fig. 4

Inner region of Westerlund 1 observed with ACIS (left panel) in the broad band and with HST (right panel) using the F160W filter. In the ACIS image, source confusion and a high background intensity dominate the cluster core.

In the text
thumbnail Fig. 5

Spatial distribution of candidate sources detected with the four methods (from the top: Pwdetect, Wavdetect, image reconstruction, and time-resolved Pwdetect). The left panels show the whole ACIS field, those on the right the inner region. Different colors in the first and second rows mark sources detected at different energy bands.

In the text
thumbnail Fig. 6

Example of a source (label MOXC2) that was excluded as a potential product of a PSF hook near a brighter source (label c10200). The green contours outline the extraction regions of MOXC2 in all observations, while the red polygons indicate the locations where the PSF hook may appear in each observation, depending on the roll-angle.

In the text
thumbnail Fig. 7

Extraction regions of the validated sources across the entire merged ACIS image (left panel), and extraction regions in the central area of approximately 3′ in size (right panel).

In the text
thumbnail Fig. 8

Extraction regions of the validated sources and positions of the input candidate sources (crosses) within the central 1 arcmin region

In the text
thumbnail Fig. 9

Distribution of the off-axis angles of the EWOCS X-ray sources.

In the text
thumbnail Fig. 10

Distributions of source net counts in the broad, soft, medium, hard, and very hard energy bands. The number of sources with more than 50 counts is indicated in the top-right corner of each panel.

In the text
thumbnail Fig. 11

Spatial distribution of EWOCS X-ray sources, sorted into two bins based on their net counts in the broad band.

In the text
thumbnail Fig. 12

Distribution of the median photon energy in the broad band for the validated EWOCS X-ray sources (black) and the catalog published by Townsley et al. (2018), in red.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.