Unveiling very young O stars. Two Galactic O2V((f*))z in Westerlund 2
A. Roman-Lopes
 

    
   Two Galactic O2V((f*))z in Westerlund 2
   

   Department of Astronomy, Universidad de La Serena,
              Av. Raul Bitran #1302, La Serena, Chile

              aroman@userena.cl
             

   



 
  
   O-type stars are known to significantly contribute to both the dynamics and evolution of galaxies. Massive and luminous, they probably control and regulate the galaxies star formation rates. The stellar feedback generated by such cosmic beasts can strongly affect the local star formation rate, with effects in the current (and future) generations of low and intermediate mass stars, and possibly also in the disruption process of the galaxies' giant gas reservoirs.  
  
   For this work I performed a redetermination of the spectral types and effective temperatures of the Galactic O-type stars MSP182, MSP183, MSP199, VPHAS-01338, and VPHAS-01273.
  
   From a careful examination of the spectral features present in the blue optical spectral region, it was possible to identify several nitrogen lines usually only seen in the blue optical spectra of O2-O3 stars. From the nitrogen ionic equivalent width ratios measured in the spectra of MSP182, MSP183, MSP199, VPHAS-01338, and VPHAS-01273, and in those of standard stars of the O2-O4 spectral types, earlier spectral types and hotter effective temperature values were derived.
   Two O2V((f*))z, together with three new O3 V stars are now firmly identified in the Westerlund 2 region.  Besides RFS1 in NGC3603, the O2 stars found in Westerlund 2 are the only other exemplars known to date in the Milky Way. From the nitrogen equivalent width line ratios measured in the spectra of standard stars of the O2-O4 spectral types, linear relations  between the Nivλ4058/Niiiλ4640 ratio and the effective temperature in the 47000K-51000K range were derived. Based on my spectroscopic analysis of the science targets and the use of a HRD, a mean heliocentric distance of 5kpc to Westerlund 2 was computed, a result that is in line with the mean heliocentric distance of 5.3±1.5 kpc obtained from the associated Gaia DR3 parallaxes and distances.
  
  The Westerlund 2 massive stars studied in this work probably share a common evolutionary process that might be representative of the evolutionary ages of a large fraction of the cluster's O-type stellar population, which seems to be much younger than 1 Myrs.
  


   

   Unveiling very young O stars
    A. Roman-Lopes
          1
          
    Received ; accepted 
=========================================






§ INTRODUCTION




  
O-type stars are known for their elevated temperature, substantial size, and high luminosities, playing a pivotal role in shaping the evolution and dynamics of galaxies. Their lifespans, which are relatively brief on astronomical timescales (typically lasting a few million years), offer valuable insights into the temporal progression of galactic populations. Functioning as stellar powerhouses, O-type stars emit substantial energy in the form of intense ultraviolet radiation, stellar winds, and eventual supernova events. Consequently, they exert a pronounced influence on the interstellar medium (ISM), significantly contributing to the dynamics and evolutionary processes within galaxies. In this sense, they probably regulate the star formation rates within galaxies since the stellar feedback generated by such cosmic beasts can affect both the formation of subsequent star generations and the dispersion process of giant molecular clouds <cit.>.

The O2 is the earliest spectral type among O stars. It was introduced by <cit.>, designating early O stars with distinct characteristics in their spectral lines. Their spectra show very weak (or absent) Hei and Niii 4640 lines, which – combined with a higher excitation manifested through prominent Niv 4058Å  and Nv 4602-4621Å transition lines – sets O2 stars apart from their later type siblings. The unique spectral fingerprint among normal O-type stars suggests that the O2 type boasts the highest effective temperatures among dwarf, giant, and super-giant stars.
In the Galaxy, only a couple of O2 super-giant stars are known to date: HD93129A in Trumpler 14 and RFS3 in the periphery of NGC 3603, with both being classified as O2 If* <cit.> and <cit.>. Surprisingly, no O2 III stars are known in the Milky Way to date, and the only O2 dwarf is RFS1 <cit.> which is found at about 1 arcmin south of the center of the super-star cluster NGC3603.
As can be noticed, the observed paucity of exemplars of the O2 type underscores the rarity of such stars that formed from metal-rich molecular gas clouds, as normally found in the Galaxy.
On the other hand, in the lower metallicity environment of the Magellanic Clouds, there are at least ten known O2 stars: three O2 If* (Brey77, RMC136a5, and cl* NGC 2070 MH 59), one O2 IIIf* (VFTS16), and six dwarfs (Walborn 2, BI237, BI253, VFTS169, VFTS468, and VFTS512), with many of them spread around the 30 Doradus starburst region <cit.> and <cit.>. The higher number of O2 stars found there when compared to the Galaxy can be understood if we consider that, in the presence of a lower metallicity gas environment, the most massive stars produced in the associated molecular clouds tend to be hotter for a given spectral type <cit.>.
For this work I performed a reexamination of the spectral types and effective temperatures of the Galactic O-type stars MSP182, MSP183, MSP199, VPHAS-01338, and VPHAS-01273, communicating the identification of two O2 and three O3 main sequence stars in the Westerlund 2 cluster, which were previously misclassified as being of later spectral types.



§ NEW O2 AND O3 MAIN SEQUENCE STARS IN WESTERLUND 2


The Galactic O-type stars MSP182, MSP183, MSP199, VPHAS-01338, and VPHAS-01273 (hereafter the O-star sample) are found in the direction of the Westerlund 2 stellar cluster. They were recently studied by <cit.> using non-local thermodynamic equilibrium model-atmosphere fitting techniques applied to the hydrogen and helium line profiles present in the associated blue optical spectra. The authors concluded that the earliest source of the sample would be MSP183, a O3 V star that was previously classified as such by <cit.>. In the case of the remaining sources of the O-star sample, MSP182, MSP199, VPHAS-01338, and VPHAS-01273, the authors concluded that all are O4 V stars. However, from a careful examination of the spectral features present in the 4000Å-5000Å wavelength range of the blue optical spectra of the sample (which are not presented in their entirety in the original article), it is seen that there are (in all of them) absorption and emission nitrogen lines in a combination that is usually only present in the spectra of stars of the O2-O3 types. This is of the utmost importance as it is well known that the canonical method of using the Hei/Heii line ratio <cit.> and <cit.> might fail in the determination of the effective temperature (or spectral types) of the hottest O-type stars. Instead, <cit.> and <cit.> suggested complementing the spectral analyses using specific nitrogen ionic lines such as the Niii 4634-4640-4642, Niv 4058, Nv 4602, and/or Nv 4621 transitions present in the associated blue optical spectral range.



 §.§ Spectral types of the O-star sample based on blue optical spectra of very early-O standard stars


In order to reestimate the spectral types of the O-star sample sources, in this section their XSHOOTER blue optical spectra are compared with those of the standard stars HD269810 - O2 III((f*)), BI237 - O2V((f*))z, BI253 - O2 V((f*)), HD64568 - O3 V((f*))z, and HD46223 - O4 V((f)) <cit.>, taken from the ESO phase-III archive[http://archive.eso.org/cms.html]. Table 1 lists the instrument name and the ESO program ID from which the associated spectral data were obtained. All spectra used in the comparison were normalized using the SPLOT task of IRAF, with all of them matching the same spectral resolution R=6655 that corresponds to the smallest spectral resolution of the spectra (in the 4000Å-5000Å wavelength range) used in this work. Finally, in the case of spectral data for which the original resolution is superior to that of the O-star sample, the template spectrum was degraded to match the spectral resolution R=6655. 
The results are shown in Figure 1. The  spectra of the O-star sample are shown in black, while those from the O-type standards are seem in gold and red. Also, the most important spectral features are indicated by the vertical black labels. Next, I describe and comment on the results obtained for each source of the O-star sample in detail.

   




 



 §.§ VPHAS-01338: O3 V((f*))


The blue optical spectrum of VPHAS-01338 is shown in Figure 1(a) together with the spectra of the Galactic stars HD64568 (O3 V((f*))z) and HD46223 (O4 V((f))). It can be observed that its Hei 4471Å absorption line is much weaker than that of HD46223, a behavior that indicates that  VPHAS-01338 is hotter. Its blue optical spectrum resembles that of HD64568 well, the Galactic standard of the class, with the NIV 4058Å emission line showing a similar intensity to the N III 4640Å one <cit.>. Taking into account the observed similarity between their He I, N III, N V, and C IV lines (see the insets shown in Figure 1(a), the spectral type O3 V((f*)) is assigned for VPHAS-01338.



 §.§ VPHAS-01273: O2 V((f*))z

  
In Figure 1 the rectified blue optical  spectrum of the VPHAS-01273 source (black) is shown, together with those of BI237 (gold) and BI253 (red), which are both standards of the O2 V((f*)) class <cit.>. On the one hand, the former shows very intense Heii absorption lines at 4200Å, 4542Å, and 4686Å, and the absence of any of the Hei lines normally seen in the optical spectra of most O-type stars indicate that the star's spectral type is certainly earlier than O4 <cit.>. Its Heii absorption line at 4686Å is seen in pure absorption (the horizontal dashed line indicates the minimum of its observed intensity), so a dwarf class is assigned. On the other hand, the presence of the Niv 4058Å line in emission, together with absence of He I lines indicate that this object probably belongs to the hottest lineage of O-type stars. Last but not the least, the spectral feature seen at 4686Å corresponds to the most intense He absorption line in the 4000Å-5000Å wavelength range, which is a behavior characteristic of the Vz class <cit.>, so VPHAS-01273 is likely very young. From the comparison of the VPHAS-01273 spectrum with those of the BI237 and BI253, one can see that the three are very similar, which suggests that VPHAS-01273 is a Galactic exemplar of the O2 V((f*))z type.



 §.§ MSP183: O2 V((f*))z


The blue optical spectrum of MSP183 shown in Figure 1(c) presents very intense Heii 4200Å, 4542Å, and 4686Å absorption lines. The absence of any of the Hei lines that are usually seen in the optical spectra of most O stars indicates that MSP183's spectral type is possibly earlier than O3 <cit.>.
The comparison of its blue optical spectrum with that of the O3 V((f*))z star HD64568 (the Galactic standard of the type) shows that it presents a stronger Hei 4471Å absorption line, which, together with a C IV 4658 line much weaker than the one in the science spectrogram <cit.>, leaves no doubt that MSP183 is a hotter star. 
On the other hand, the comparison of its spectrum with that of the O2 star HD269810 shows a rather similar behavior with one of the few noticeable differences between them being related to their Heii 4686Å absorption lines. In the case of HD269810, its Heii line is seen partially (weakly) in emission (so luminosity class III); whereas, in the MSP183 optical spectrogram, it is observed in pure absorption (luminosity class  V). As in the case of VPHAS-01273, its Heii 4686Å absorption line is the most intense He transition in the 4000Å-5000Å wavelength range, and the Vz class is secured <cit.>. Finally, the NIV 4058Å emission line in the MSP183 data is much stronger than NIII 4640Å, so the O2 type is the best choice in that case <cit.>.


   




 §.§ MSP182 (O3.5 V(f+)) and MSP199 (O3 V((f*))


The two remaining sources of the O-star sample are MSP182 and MSP199. Their blue optical spectra are shown in Figure 1(d) and 1(e), together with those of HD64568 (O3 V((f*))z) and HD46223 (O4 V((f))). The Hei 4471Å absorption lines are seen much weaker in the MSP182 and MSP199 spectrograms; although, in the case of MSP182, they are stronger than that seen in the O3 standard spectrum. Also, and as in the case of VPHAS-01338, the Niii 4640Å emission lines appear much weaker than that seen in the HD46223 spectrogram. In the case of MSP182, its Niv 4058Å emission line appears weaker than that of Niii 4640Å, which together enables its spectral types to be estimated to be O3.5 <cit.>. On the other hand, in the case of MSP199, its Niv 4058Å emission line appears with a similar intensity to that of Niii 4640Å (see the horizontal dashed line positioned in the peak of the Niii 4640Å emission line), so O3 type is assigned in its case. Finally, their Heii 4686Å lines are both seen in pure absorption, and the dwarf classes are assigned to MSP182 and MSP199.
In what follows, I further explore the optical spectra of the O-star sample by using the information provided by the Niiiλ4640, Nivλ4058, Nivλ6380, and Nvλ4603 ionic lines.



§ SPECTRAL TYPES AND EFFECTIVE TEMPERATURES FROM NITROGEN IONIC EQUIVALENT WIDTH RATIOS


The canonical method of using the Hei/Heii line ratios <cit.> fails in the determination of the effective temperatures for the earliest O stars, with a degeneracy of possible solutions in the associated parameters' space. In order to brake this degeneracy, <cit.> proposed that the scheme in which ionic transitions of the nitrogen atom, namely the ratio Nvλ4058/Niiiλ4640, be used instead. Later, <cit.> provided theoretical support to it by modeling the observed nitrogen ion line ratios. They concluded that once the luminosity class is known, the effective temperatures of O2-O4 stars can be estimated well from the Nvλ4058/Nivλ4640 line ratios.



 §.§ Spectral types






It is possible to provide new estimates for the spectral types of the O-star sample sources by comparing their observed Nivλ4058/Niiiλ4640 and Nvλ4603/Nivλ6380] equivalent width (EW) ratios, with those from known spectroscopic standards of the O2, O3, and O4 types. In Table 1 the list of the ESO Phase 3[http://archive.eso.org/cms.html] dataset containing the optical spectra of O-type standards used in this work is shown.
All spectra used in my measurements were carefully continuum rectified and degraded to the same spectral resolution as previously explained in Section 2.



   

In Figure 2 the Log[Nvλ4603/Nivλ6380] × Log[Nivλ4058/Niiiλ4640] diagram is presented in which the empirical positions of O2 (violet crosses), O3 (blue star), O3.5 (gold triangles), and O4 (red triangles) spectroscopic standard stars are shown. They are labeled by the bold blue text. It can be seen that the distinct spectral types appear clearly separated from each other, mainly in the Log[Nivλ4058/Niiiλ4640] dimension of the diagram, with the O2 stars' group being clearly displaced from the later types by about 0.5 dex. The two new O2V((f*))z stars VPHAS-01273 and MSP183 are found just in this part of the Log[Nvλ4603/Nivλ6380] × Log[Nivλ4058/Niiiλ4640] diagram.
Also, accordingly to <cit.>, the O3 spectral type should correspond to EW[Nivλ4058]/EW[Niiiλ4640] line ratio = 1 (e.g., Log[Nivλ4058/Niiiλ4640] = 0). From Figure 2, it can be seen that besides HD64568 (the Galactic standard of the O3 class), sources VPHAS-01338, MSP182, and MSP199 all have Log[Nivλ4058/Niiiλ4640] ratio values close to zero, which suggests that they are very likely also O3 stars <cit.>. These results are in line with those obtained in Section 2, in which the O3 V spectral type was assigned for VPHAS-01338, MSP182, and MSP199.



 §.§ Effective temperatures


As mentioned above, <cit.> performed a theoretical study of the behavior of the nitrogen ion line ratios in the case of O2-O4 stars. Particularly
in their Figure 3, they present the theoretical relations between the EW of Niiiλ4640 and Nivλ4058 lines, as a function of the effective temperature of O-type stars in the Galaxy (Milky Way), and in the Small Magellanic Cloud (SMC), for models with gravity and nitrogen abundance values of log g = 4.0 and [N] = 7.78 dex, respectively.

From the work of <cit.>, I proceeded by tabulating the Niiiλ4640 and Nivλ4058 EW values for the Milky Way and SMC models considering the thin wind case for constant temperature values in the interval 47000K-51000K. In Figure 3, the Log[Nivλ4058/Niiiλ4640] versus effective temperature diagram obtained from the computed ratio values is shown. There, the theoretical relations are represented, respectively, by the blue (Milky Way) and red (SMC) symbols. It can be seen that in both cases, there is a linear relation between the effective temperature and the Log[Nivλ4058/Niiiλ4640]. The linear trend appears steeper in the case of the SMC models, so, as consequence, any O2 star there should be hotter than its sibling in the Galaxy. However, one can also see that in the mentioned temperature range, the maximum difference between the effective temperature values of SMC and Milky Way stars (occurring close to the value Log[Nivλ4058/Niiiλ4640] = 1.0) is only  1000K. On the other hand, it can be observed that for a O3 star (Log[Nivλ4058/Niiiλ4640]=0), both Milky Way and SMC stars should have similar effective temperatures. This occurs because the use of EW ratios minimizes the effects of metallicity differences between stars in the Milky Way and in the SMC. This is a very fortunate behavior as it occurs just in the temperature range in which the canonical method of using the Hei/Heii line ratios <cit.> starts to fail. Linear fittings to the theoretical points were performed in each case, and the results are shown in Figure 3 by the dashed red and blue lines, together with the corresponding linear relations derived from the linear fits. They can be used to estimate the effective temperatures of the earliest O-type stars, in the ranges -0.75<Log[Nivλ4058/Niiiλ4640]<1.5 in the case of the Galaxy and in the range -0.5<Log[Nivλ4058/Niiiλ4640]<0.75 in the case of SMC early O-type stars. As a complement, the linear relation in the intermediate case between the Milky Way and SMC is also provided. It is represented in Figure 3 by the dashed gold line which corresponds to the linear fit applied to the average of the Log[Nivλ4058/Niiiλ4640] values.
   
One can use the derived linear relations to estimate approximate values for the effective temperatures of the standard O-type stars shown in the Log[Nvλ4603/Nivλ6380] × Log[Nivλ4058/Niiiλ4640] diagram of Figure 2. I used the corresponding linear equations shown in Figure 3 to compute the average SMC and Milky Way effective temperature values for O2, O3, O3.5, and O4 stars. They are represented in the plot by the vertical light red lines. The different widths of the lines indicate the approximate uncertainty ranges for the temperatures 46500K (O4), 47500K (O3.5), 48500K (O3), and 49500K-50500K (O2). Interestingly, the three standards of the O2V class present Log[Nivλ4058/Niiiλ4640] values around 0.6 dex, while the O2 III class members present higher values around 0.9 dex, which is equivalent to an increase in the effective temperature of about 1000K. This behavior is probably related to the fact that for stars with stronger stellar winds, the Nivλ4058/Niiiλ4640 ratio could be less reliable as a tool for the determination of effective temperatures of giants and super-giants of the O2 spectral type since the derived values  would  probably be overestimated by at least  1000K.



 §.§ The Hertzsprung-Russel diagram of the O-star sample



   

The Hertzsprung-Russel diagram (HRD) of the O-star sample is shown in Figure 4, with the stars represented by the golden-filled circles. Their effective temperatures (shown in Table 2) correspond to the values calculated using the Galactic linear relation derived in Section 3, applied to the observed Nivλ4058 and Niiiλ4640 EW values presented in Table 3.
Their absolute magnitudes M_V were computed using the V-band magnitudes taken from the work of <cit.> (VPHAS-1273 and VPHAS-1338), and the V-band magnitudes of MSP183 and MSP199 were computed from the F555W and F814W magnitudes provided by the Hubble Space Telescope Legacy Archive[https://hla.stsci.edu/hlaview.html], using the transformation equations derived by <cit.>. I chose to use this procedure in the case of these two stars because, from an inspection of the HST PC F336W image of the central part of the Westerlund 2 cluster, I realized that MSP183 appears blended (see Figure 5) in the HST WFPC2 F814W lower spacial resolution image (and probably also in those taken from ground-based telescopes). This fact provides a further explanation of why MSP183 appears so luminous (and massive) in previous studies of Westerlund 2 <cit.>.
On the other hand, similar to MSP199, MSP182 is also seen relatively isolated in the HST archive images of Westerlund 2, and in its case, I chose to use the V-band magnitude value derived by <cit.>. The visual absorption values of the sources of the O-star sample were corrected using the A_V values computed by <cit.>. In the case of MSP183, its computed A_V value was probably affected by its fainter companions. Fortunately, as it was placed only at a few arc-seconds from MSP199 (see Figure 5), it is reasonable to assume that its visual absorption is probably similar, so I used A_V=6.4 mag instead. 
The positions of the 0.10 Myrs, 0.25 Myrs, 0.5 Myrs, and 1.0 Myrs PAdova and TRieste Stellar Evolution Code (PARSEC))[http://stev.oapd.inaf.it/cgi-bin/cmd] isochrones <cit.> (40M_⊙-120M_⊙ mass range) are indicated by the dashed gray, red, blue, and green lines. Regarding the uncertainties on the HRD parameters, in the case of the effective temperatures, it is shown in Section 3.2 that the uncertainties in the Teff for values in the range 48000K-50000K are probably on the order of less than  500K. However, in order to be conservative, 1500K was assumed to be the total uncertainty instead. On the other hand, in the case of absolute visual magnitudes, an uncertainty of 0.25 magnitudes was assumed. It can be seen that the five stars present (virtually) the same visual absolute magnitudes, with an observed dispersion on M_V 0.5 mag. Regarding the mean heliocentric distance used in the computations, it was set considering that the observed position of the stars in the HRD must be consistent with the most recent results dictated by the massive star evolution theory. This, in the case of dwarf stars of the O2 and O3 types, corresponds to the mass range 60M_⊙-70M_⊙ <cit.>, which in turn sets a mean heliocentric distance of 5kpc. 

Although the heliocentric distance of Westerlund 2 is an important problem, it is well beyond the main goal of this work. However, it is useful and necessary to provide some comments on it. The heliocentric distance of 5kpc used in the construction of the HRD is the natural consequence of the consistency between the overall results presented here, and the need to satisfy the range of masses dictated by the current theory of massive star formation, through the use of state-of-the-art isochrones. As result, the 5 kpc value is the most probable mean heliocentric distance to the five stars of the O-star sample. On the other hand, it is possible to check the consistency of my result, by comparing it with the heliocentric distances provided by Gaia DR3.
In order to do that, the heliocentric distances for the five stars of the O-star sample were taken from the Gaia EDR3 catalog of <cit.>, and the mean values for rgeo and rpgeo were computed. The results of the calculations are 5293±1523 pc and 5257±1555 pc, respectively, with the quoted errors corresponding to the standard deviation values of the mean.




§ DISCUSSION


The O2 stars VPHAS-01273 and MSP183, together with RFS1 <cit.> are the only other dwarfs of the type identified to date in the Milky Way. They are placed in Westerlund 2, which, with the NGC 3603 starburst cluster, are the only known sites harboring Galactic dwarfs of the O2 type. 
As very young massive stars, MSP183, together with its siblings MSP182 and MSP199, are found exactly at the place predicted by star formation theory: in the center of their parental cluster. On the other hand, the odd novelty comes from the fact that VPHAS-01273 and VPHAS-01338 are not found in a stellar cluster. They are indeed found isolated in the field, at 11.6and 21.2from the closest known massive star cluster, respectively. How 60M_⊙-80M_⊙ stars can be found so far from their likely parental cluster is unclear. A further study on this theme is beyond the context of this Letter; however, next, I briefly comment on some aspects (and consequences) of the subject.



 §.§ The earliest runway dwarf stars known in the Galaxy


<cit.> proposed an explanation for the current location of VPHAS-01273 and VPHAS-01338, based on their study of radial velocities and proper motions of candidate runaway stars placed around the Westerlund 2 cluster. The authors postulated that, in the case of VPHAS-01273 and VPHAS-01338, the stars traveled from the central parts of the nearest massive star cluster, which in their case is Westerlund 2. The authors concluded that in both cases, their observed space velocities meet the currently accepted standard runaway criterion, calculating that the stars were possibly ejected from Westerlund 2 at  500,000 and 800,000 years ago,  respectively. Massive star formation theory predicts that stars such as VPHAS-1273 and VPHAS-01338 should evolve very quickly, leaving the main sequence in no more than a few million years. Considering the estimated travel times, it is reasonable to assume that if the stars were expelled from their parental cluster, it is very likely that it occurred in the very beginning of their main sequence lifetimes. Under such conditions, some of them should present some sort of spectral signature evidencing their very early evolutionary stage. In fact, such spectral evidence does exist, at least in the case of VPHAS-1273. Its Heii 4686Å line is the most intense He absorption line in the 4000Å-5000Å wavelength range and, as previously mentioned, a behavior characteristic of stars of a very young age <cit.>. From an astrophysical perspective, it is desirable to check the consistency of a very small age for the stars, by comparing their positions in the HRD relative to those of the PARSEC theoretical isochrones presented in Figure 4. There they are seen as very young stars with evolutionary ages in the 100,000-500,000 year range.
From this reasoning, it is possible to conclude that VPHAS-1273 and VPHAS-01338 are probably the earliest runway O-dwarf stars known in the Galaxy.



 §.§ Consequences for the age of the massive stellar population of Westerlund 2


From an astrophysical perspective, it is useful to return to the HRD of Figure 4. The correct understanding of the physical processes that impact (or potentially govern) the temporal evolution of very massive Galactic clusters such as NGC3603 and Westerlund 2 can be partially deepened if one adequately considers the results obtained here for the O-star sample as a whole.
We shall start by looking at the positions of the stars in the HRD, not only relative to the evolutionary isochrones, but also how the ages of the stars compare to each other. On the one hand, VPHAS-01273 and MSP183 are quite young with ages inferred from the theoretical isochrones in the range of 100,000-500,000 years. On the other hand, it turns out that MSP182, MSP199, and VPHAS-01338 appear to not only show very similar absolute visual magnitudes, but very similar evolutionary ages, a youth condition that is in line with the Vz behavior observed in both O2 stars (see section 2). In fact, the entire O-star sample appears in the HRD as a very cohesive group of sources, which, from an astrophysical point of view, indicates that they probably share a common evolutionary process, perhaps representative of the entire massive stellar population of the cluster. This could suggest that the Westerlund 2 very massive stellar population is probably as young (or even younger) than the less massive members of the cluster, which have estimated mean ages in the range 0.5Myrs-2.0Myrs <cit.>, <cit.>.





§ CONCLUDING REMARKS


For this work I performed a redetermination of the spectral types and effective temperatures of the Galactic O-type stars MSP182, MSP183, MSP199, VPHAS-01338, and VPHAS-01273. The main results are summarized as follows:

   
      
  * Two O2V((f*))z (MSP183 and VPHAS-01273), together with three new O3-type stars (MSP182, MSP199, and VPHAS-01338) are spectroscopically confirmed in Westerlund 2.  Besides RFS1 in NGC3603, the two O2 stars found in Westerlund 2 are the only other exemplars known to date in the Milky Way. The O2 star MSP183 is found to be very close to several fainter sources in the HST PC F336W image of the central part of Westerlund 2, providing an explanation on why it was found so luminous (and massive) in previous studies of Westerlund 2 <cit.>, <cit.>.
      
      
  * From the nitrogen EW line ratios measured in the spectra of standard stars of the O2-O4 spectral types, linear relations  between the Nivλ4058/Niiiλ4640 ratio and the effective temperature in the 47000K-51000K range were derived. As expected, the linear trend appears steeper in the case of the low metallicity models. However, the maximum difference between the effective temperature of SMC and Milky Way O2 stars, which occurs close to the edge of the mentioned temperature range, corresponds to only  1000K. On the other hand, in the case of O3 stars, both Milky Way and SMC stars should have similar effective temperatures. This occurs because the use of EW ratios probably minimizes the effects of metallicity differences between stars in the Milky Way and in the SMC. This is a very fortunate behavior as it occurs just in the temperature range in which the canonical method of using the Hei/Heii line ratios  starts to fail.
      
  * Based on the results obtained from my spectroscopic analysis of the science targets and the use of a HRD, a mean heliocentric distance of 5kpc to Westerlund 2 was computed, a result that is in line with the mean heliocentric distance of 5.3±1.5 kpc obtained from the associated Gaia DR3 parallaxes and distances.
      
  * The Westerlund 2 massive stars studied in this work probably share a common evolutionary process that might be representative of the evolutionary ages of a large fraction of the cluster's O-type stellar population, which seems to be much younger than 1 Myrs.
   



      ARL thanks the anonymous referees for the constructive reports that contribute to improve the interpretation of the results obtained in this work.
      This work is based (in part) on data obtained from the ESO Science Archive Facility with DOI(s): https://doi.org/10.18727/archive/71, https://doi.org/10.18727/archive/50, https://doi.org/10.18727/archive/24. It is also based [in part] on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). ARL acknowledges the developers of the TOPCAT software (http://www.starlink.ac.uk/topcat/), as well as the PARSEC team (ttp://stev.oapd.inaf.it/cmd) for providing the state of the art isochones models used in this work.


  

  [Bailer-Jones et al.(2021)]bailer21 Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., Andrae, R. 2021, AJ, 161, 147
  [Bressan et al.(2012)]bressan12 Bressan, Alessandro, Marigo, Paola, Girardi, Léo., Salasnich., Bernardo, Dal Cero, Claudia, Rubele, Stefano, Nanni, Ambra 2012, MNRAS, 427, 127 
  [Carraro et al.(2013)]carr13 G. Carraro, G., Turner, D., Majaess, D., Baume, G. 2013, A&A, 550, 50
  [Conti & Alschuler(1971)]conti71 Conti P. S., Alschuler W. R., 1971, ApJ, 170, 325   
  [Dale & Bonnell(2011)]dale11 Dale, James E., Bonnell, Ian 2011, MNRAS, 414, 1, 321
  [Drew et al.(2018)]drew18 Drew, J. E., Herrero, A., Mohr-Smith, M., Monguió, M., Wright, N. J., Kupfer, T., Napiwotzki, R. 2018, MNRAS, 480, 2109
  [Ferrara & Tolstoy(2000)]ferr00 Ferrara, Andrea, Tolstoy, Eline 2000, MNRAS, 313, 2, 291
  [Harris(2018)]harris18 Harris, W. E. 2018, AJ, 156, 296 
  [Hopkins et al.(2012)]hopkins12 Hopkins, Philip F., Quataert, Eliot, Murray, Norman 2012, MNRAS, 421, 4, 3522
  [Lasker et al.(2008)]lasker08 Lasker, Barry M., Lattanzi, Mario G., McLean, Brian J., Bucciarelli, Beatrice, Drimmel, Ronald, Garcia, Jorge, Greene, Gretchen, Guglielmetti, Fabrizia, Hanley, Christopher, Hawkins, George, Laidler, Victoria G., Loomis, Charles, Meakes, Michael, Mignani, Roberto, Morbidelli, Roberto, Morrison, Jane, Pannunzio, Renato, Rosenberg, Amy, Sarasso, Maria, Smart, Richard L., Spagna, Alessandro, Sturch, Conrad R., Volpicelli, Antonio, White, Richard L., Wolfe, David, Zacchei, Andrea 2008, AJ, 136, 735  
  [Maiz-Apellaniz et al.(2014)]maiz14 2014, ApJS, 211, 10
  [Martins & Palacios(2017)]mart17 F. Martins, A. Palacios 2017, A&A, 598, 56
  [Massey et al.(2009)]mass09 Massey, Philip, Zangari, Amanda M., Morrell, Nidia I., Puls, Joachim, DeGioia-Eastwood, Kathleen, Bresolin, Fabio, Kudritzki, Rolf-Peter 2009, ApJ, 692, 618
  [Mathys(1989)]mat89 Mathys G., 1989, A&AS, 81, 237
  [Mokiem et al.(2007)]mok07 M. R., Mokiem, A. de Koter, C. J. Evans, J. Puls, S. J. Smartt, P. A. Crowther, A. Herrero, N. Langer, D. J. Lennon, F. Najarro, M. R. Villamariz, and J. S. Vink 2007, A&A, 465, 1003
  [Rauw et al.(2007)]rauw07 Rauw, G., Manfroid, J., Gosset, E., Nazé, Y., Sana, H., De Becker, M., Foellmi, C., Moffat, A. F. J. 2007, A&A, 463, 981
  [Rauw, Sana & Nazé(2011)]rauw11 Rauw, G., Sana, H., Nazé, Y. 2011, A&A, 535, 40
  [Rivero-Gonzalez et al.(2012)]river12 Rivero González, J. G., Puls, J., Najarro, F., Brott, I. 2012, A&A, 537, 79
  [Rivero-Gonzalez et al.(2012)]river12b Rivero González, J. G., Puls, J.,  P. Massey, Najarro, F. 2012, A&A, 543, 95
  [Roman-Lopes et al.(2011)]roman11 Roman-Lopes A., Barba, R. H., Morrell, N. I. 2011, MNRAS, 416, 501
  [Roman-Lopes(2013)]roman13 Roman-Lopes A. 2013, MNRAS, 435, 73
  [Roman-Lopes et al.(2016)]roman16 Roman-Lopes A., Franco, G.A.P., Sanmartin, D. 2016, ApJ, 823, 96  
  [Smith et al.(2010)]smith10 Smith, Nathan, Povich, Matthew S. , Whitney, Barbara A., Churchwell, Ed, Babler, Brian L., Meade, Marilyn R., Bally, John, Gehrz, Robert D., Robitaille, Thomas P., Stassun, Keivan G. 2010, MNRAS, 406, 2, 952
  [Sota et al.(2011)]sota11 Sota, A., Maíz Apellániz, J., Walborn, N. R., Alfaro, E. J., Barbá, R. H., Morrell, N. I., Gamen, R. C., Arias, J. I. 2011, ApJS, 193, 24
  [Walborn et al.(2002)]walb02 Walborn, Nolan R., Howarth, Ian D., Lennon, Daniel J., Massey, Philip, Oey, M. S., Moffat, Anthony F. J., Skalkowski, Gwen, Morrell, Nidia I., Drissen, Laurent, Parker, Joel Wm. 2002, AJ, 123, 2754
  [Walborn et al.(2014)]walb14 Walborn, Nolan R., Sana, H., Simón-Díaz, S., Maíz Apellániz, J., Taylor, W. D., Evans, C. J., Markova, N., Lennon, D. J., de Koter, A. 2014, A&A, 564, 40
  [Zacharias et al.(2013)]zach13 Zacharias, N., Finch, C. T., Girard, T. M., Henden, A., Bartlett, J. L., Monet, D. G., Zacharias, M. 2013, AJ, 145, 44 
  [Zeidler et al.(2015)]zeid15 Peter Zeidler, p., Sabbi, E., Nota, A., Grebel1, E. K., Tosi, M., Bonanos, A. Z., Pasquali1, A., Christian, C., de Mink, S. E., Ubeda, L. 2015, AJ, 150, 78
  [Zeidler et al.(2018)]zeid18 Zeidler, p., Sabbi, E., Nota, A., Pasquali, A., Grebel, E. K., McLeod, A. F., Kamann, S., Tosi, 6 M., Cignoni, M., Ramsay, S. 2018, AJ, 156, 211
  




