eROSITA studies of the Carina nebula

Context. During the first four all-sky surveys eRASS:4, which was carried out from December 2019 to 2021, the extended Roentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spektrum-Roentgen-Gamma (Spektr-RG, SRG) observed the Galactic H II region, the Carina nebula. Aims. We analysed the eRASS:4 data to study the distribution and spectral properties of the hot interstellar plasma and the bright stellar sources in the Carina nebula. Methods. The spectral extraction regions of the diffuse emission were defined based on the X-ray spectral morphology and multi-wavelength data. The spectra were fit with a combination of thermal and non-thermal emission models. The X-ray bright point sources in the Carina nebula are the colliding wind binary η Car, several O stars, and Wolf–Rayet (WR) stars. We extracted the spectra of the brightest stellar sources, which can be well fit with a multi-component thermal plasma model. Results. The spectra of the diffuse emission in the brighter parts of the Carina nebula are well reproduced by two thermal models, a lower-temperature component ( ∼ 0.2 keV) and a higher-temperature component (0.6–0.8 keV). An additional non-thermal component dominates the emission above ∼ 1 keV in the Central region around η Car and the other massive stars. Significant orbital variation in the X-ray flux was measured for η Car, WR 22, and WR 25. η Car requires an additional time-variable thermal component in the spectral model, which is associated with the wind-wind collision zone. Conclusions. Properties such as temperature, pressure, and luminosity of the X-ray emitting plasma in the Carina nebula derived from the eROSITA data are consistent with theoretical calculations of emission from superbubbles. This confirms that the X-ray emission is caused by the hot plasma inside the Carina nebula that has been shocked-heated by the stellar winds of the massive stars, in particular, of η Car.


Introduction
The Carina nebula complex (NGC 3372, hereafter CNC) is the brightest emission nebula in the southern Galactic plane at a distance of 2.3±0.1 kpc (Smith 2006a;Shull et al. 2021;Göppl & Preibisch 2022).It harbours hundreds of protostars and several young stellar clusters (Trumpler 14,15,and 16;Collinder 228 and 232;and Bochum 10 and 11) and is a giant H ii region, which has been created by the radiation and winds from more than 60 OB stars, three Wolf-Rayet (WR) stars, and the massive binary η Carinae (η Car hereafter) in Trumpler 16 (Smith & Brooks 2008).In addition, at about 1 • northwest of the CNC, an H ii region around the stellar cluster NGC 3324 called the Gum 31 region is located, with a similar distance as to the CNC (Kharchenko et al. 2005).Thanks to its proximity, the CNC is the ideal place to study star formation and massive star feedback in the Milky Way.
Inside the CNC, there are dark shells observed in the optical, in which stars are being formed (Smith et al. 2000(Smith et al. , 2010)).Images taken with the Hubble and Spitzer Space Telescopes have revealed large pillars of dust in the nebula around η Car.A large amount of dust has also been found in observations in the sub-mm-range (Preibisch et al. 2011b) and the far-infrared (FIR, Preibisch et al. 2012).Observations of molecular line emission have revealed some dynamics in the dust shells with velocities that are different from that of the Carina arm, in which they are located (Whiteoak & Otrupcek 1984;Cox & Bronfman 1995;Brooks et al. 1998;Yonekura et al. 2005).
In the CNC, there is also molecular gas which may develop into new star-forming regions in the future (Smith & Brooks 2008).Three regions with different star formation properties have been detected in the CO data: Near η Car, there are small globules, which are most likely surviving cores in the giant molecular cloud complex.Around Trumpler 14, there is the Northern Cloud, which shows little indication for star formation.In the south, there is the Southern Cloud, which has been affected by the massive stars in Trumpler 16 and in which a second generation of stars is forming.
Active star formation is in particular observed in the dust pillars.In optical and infrared (IR) observations, more than 900 young stellar objects (YSOs) have been detected in the dust pillars (Smith et al. 2010) including protostars with ages of only 10 5 yr.The star formation in these dust pillars was most likely triggered by winds from the massive stars.In the entire CNC, more than 60,000 young stars have been detected (Preibisch et al. 2011a).With such a high rate of star formation, the CNC contributes significantly to the current star formation in the Milky Way (Povich et al. 2011;Preibisch et al. 2011a).
Soft X-ray emission has also been observed from the CNC with the Einstein Observatory (Seward & Chlebowski 1982) and ROSAT (Corcoran et al. 1995).Observations with XMM-Newton and Suzaku have shown a thermal spectrum with enhanced element abundances particularly of iron (Hamaguchi et al. 2007b;Ezoe et al. 2009).With the Chandra X-ray Observatory, diffuse X-ray emission was confirmed and analysed in a 1 Ms survey of the CNC (Chandra Carina Complex Project, CCCP, Townsley et al. 2011a,c).Spectral analysis of the diffuse emission has shown that the emission is from thermal plasma, with indications for charge exchange, most likely between the hot plasma inside the nebula and the cold pillars, ridges, and clumps.The diffuse X-ray emission is most likely caused by the stellar winds of massive stars (Castor et al. 1975), however, observations have also shown, that the X-ray emitting plasma is not concentrated around the massive stars but is distributed in more broad areas.There is a bright 'hook' in X-rays, which does not coincide with any prominent emission at longer wavelengths.In the east, the emission of the 'X-ray hook' is aligned with the eastern arm of the striking V-shaped dust lane seen in the optical.The diffuse X-ray emission is brightest south of the western arm of the Vshaped dust lane.Some anti-correlation is seen between the diffuse X-ray emission and the 8 µm IR emission, suggesting that the X-rays shine through holes in the distribution of heated dust.Alternatively, the bulk of the X-ray plasma is too hot and too tenuous to be seen, such that it only appears where that hot plasma pushes through the holes and mixes with gas from the hole rim.The mixing of both reduces the temperature of the hot gas and increases its emission measure, and the gas appears in soft X-rays.This process is predicted from 3D hydrodynamic simulations to occur frequently whenever massive-star wind bubbles merge.
In this paper we report on the studies of the diffuse emission of the CNC using the extended ROentgen survey with an Imaging Telescope Array (eROSITA; Merloni et al. 2012;Predehl et al. 2021) on board the Spektrum-Roentgen-Gamma (Spektr-RG, SRG) spacecraft.eROSITA has performed all-sky surveys in the soft to medium X-ray energy band of 0.2 -10 keV, which is in particular suited for the study of the hot phase of the interstellar medium (ISM).The all-sky surveys allow us to study the emission from the entire nebula with a uniform exposure and resolution.In addition, as the emission in the surroundings is observed as well, we are able to constrain the local X-ray background well.We present the study of the morphology and the spectra of the diffuse emission in different regions.We compare the physical parameters of the ISM derived from the eROSITA data with theoretical calculations.We have also analysed the emission of the X-ray brightest massive stars in the CNC and study the time variability.

eROSITA
The Carina nebula complex was observed in the all-sky surveys of eROSITA on SRG, which was launched on July 13, 2019 from the Baikonur cosmodrome.There are seven telescope modules (TMs) on eROSITA, which are each equipped with a CCD detector.After the launch, it became obvious that two TMs suffer light leak (TM5 and 7) and therefore have a higher background.The effect is not well understood yet, therefore, data taken with TM5 and 7 are not included in the analysis presented here.eROSITA has a large field of view (FOV) with a diameter of ∼1 degree and a high sensitivity in the energy range up to 10 keV, especially below 2.3 keV.It is thus the ideal telescope to study the distribution and the spectral properties of the large extended X-ray emission from the ISM.eROSITA is planned to perform a total of eight surveys of the entire sky (eROSITA all-sky survey, eRASS1-8, the total survey called eRASS:8), out of which four surveys have been completed now (eRASS:4).We present the analysis of the eRASS:4 data of the CNC.
For data processing and analysis, we have used the eROSITA Science Analysis Software System (eSASS, Brunner et al. 2022).We have used the energy calibrated event files from the eSASS pre-processing pipeline c020 and eSASS users version 211214.
The Carina nebula complex has an extent of about 2 • × 3 • .In order to fully cover the CNC and also to have data for estimating the local X-ray background, we have used eRASS:4 data from a region with a size of ∼ 10 • × 5 • centered on the CNC.

HI
Archival data of Hi 4π (HI4PI) all-sky Hi survey (HI4PI Collaboration et al. 2016) were used in this work to compare the distribution of the hot interstellar plasma emitting X-rays with that of the cold diffuse gas.It is likely that the massive stars in the CNC have produced a cavity in the distribution of the cold atomic gas, which is filled with photo-and shock-ionised gas seen in the optical to X-rays.The HI4PI survey combines the Effelsberg-Bonn Hi Survey (EBHI; Kerp et al. 2011;Winkel et al. 2016) and Galactic All-Sky Survey (GASS; McClure-Griffiths et al. 2009;Kalberla et al. 2010;Kalberla & Haud 2015) obtained with the Parkes telescope.The angular resolution is 16.′ 2, corresponding to ∼ 11 pc at the distance of the Carina nebula (2.3 kpc; Smith 2006a).The brightness temperature noise level is ∼43 mK at a velocity resolution of 1.49 km s −1 .

Images
Using the energy calibrated data of eRASS:4, we created images in the energy bands of 0.2 -10.0 keV, 0.2 -0.5 keV, 0.5 -1.0 keV, 1.0 -2.0 keV, and 2.0 -10.0 keV.The data are stored in partly overlapping sky tiles, each of which has a size of 3.6 • × 3.6 • .First, we merged the event files of the sky tiles and recalculated the sky coordinates with the CNC in the centre (RA = 10:43:72.8,Dec = -59:49:48).Images were created for each energy band using a bin size of 400, i.e., one image pixel has a size of 20 ′′ × 20 ′′ .Based on the images, vignettingcorrected exposure maps and exposure-corrected count-rate images were created.The average exposure time is ∼ 620 s.In addition, we created event files and images without point sources (called 'point-source-free' data hereafter).We excluded all significantly detected point sources from the data.We use the standard eRASS point source catalogue of eRASS:4 produced by the consortium with a circular extraction region with a radius of 3 times the half-energy width (HEW) of 14.5 ′′ (Merloni et al. 2023, for the catalogue of the first all-sky survey eRASS1).In the ∼10 degr 2 field around the CNC, we excluded ∼ 3000 sources.For η Car, which is the brightest X-ray source in the CNC, we used a circular extraction region with a radius of 64 ′′ , i.e., 3 times larger than the HEW of eROSITA at higher energies (29 ′′ below 3 keV and 42 ′′ at 4.5 keV, Dennerl et al. 2020), to remove as much emission from this source as possible.Using these filtered data, we binned images again with a bin size of 400 and created exposure maps and count-rate images in the same energy bands as for the complete data.The count-rate image in the full band (0.2 -10.0 keV) is shown in Fig. 1, while a three-colour image of the point-source-free data in the bands 0.2 -0.5 keV, 0.5 -1.0 keV, and 1.0 -2.0 keV is shown in Fig. 2a with extraction regions for the spectral analysis (Sect.3.2).
As can be seen in the eROSITA images, there is significant diffuse emission in the central region, where the open clusters Trumpler 14, 16, and Collinder 232 as well as η Car and the Wolf-Rayet star WR25 are located.Around the massive stars, there is significant emission above 1 keV (appearing blue in the three-colour image in Fig. 2a in the Central region 1 (C1) around Trumpler14), indicating emission from the stellar winds.The emission at the position of the stellar cluster Trumpler 15 is also bright above 1 keV.
There is bright and soft extended X-ray emission south of Trumpler 14, which extends to the east, beyond Trumpler 16 and Collinder 232 (Central region 2, C2).The Carina nebula is very bright in this region in the optical (see Fig. 2b) and in the radio continuum (Rebolledo et al. 2021, see Sect. 5.1).The gas seems to be photoionised strongly from the radiation of the massive stars of the stellar clusters, η Car, and WR25.
The most pronounced structure in X-rays is the 'X-ray hook'.It appears green in the three-colour image, indicating that the emission is dominated around 1 keV, with more contribution of emission above 1 keV than in the adjacent C2 region.
Next to the X-ray hook (Hook region 1 and 2, H1 and H2, respectively), there is emission mainly between 0.5 keV and 1.0 keV, indicating hot thermal plasma (Outskirt regions 1 and 3, OS1 and OS3).In these region, shell-like structures are seen in the photoionised gas (Fig. 2b), traced also by dust (Fig. 2c).
Towards lower Galactic latitudes, the emission is soft (around 0.5 keV) and rather faint and diffuse (Outskirt region 2, OS2, and further out, being a little brighter in Outer region 2, O2).Similar emission is also observed to the north (Outer region 1, O1).As can be seen in the HI map (Fig. 2d), soft faint emis-sion seems to indicate hot gas escaping into cavities in the cold gas distribution.

Spectra
For the spectral analysis, we used the point-source-free data.Therefore, the emission from the massive stars, e.g., in η Car is not included in the diffuse emission but will be addressed separately (see Sect. 4).For the analysis of the diffuse emission, however, one has to take into account that there might be contamination by the emission of the very bright X-ray sources in areas around them.
We first defined extraction regions based on the photon statistics using the Voronoi tessellation algorithm (Cappellari & Copin 2003).However, the diffuse emission from the CNC is very faint.Even selecting a signal-to-noise ratio as low as 40, no reasonable spectra were obtained, which would have allowed the analysis of the emission of the Carina nebula above the background.Therefore, we decided to define the extraction regions manually in the next step, based on the comparison to multiwavelength data (see Fig. 2 and Sect.3.1).The criteria we applied are the following: contiguous emission in X-rays with sufficient photon statistics, similar colour in the three-colour X-ray image indicating a similar spectrum, corresponding to known astrophysical objects/structures (stellar clusters, molecular clouds, interstellar cavities), well correlated with distribution observed at other wavelengths (photoionised gas seen in the optical or radio, dust, molecular clouds, atomic gas).The extraction regions are shown in Fig. 2. The region from which the spectrum was extracted to estimate the local X-ray background is shown with a dashed circle in the upper left panel.
We analysed the spectra using XSPEC version 12.12.1.All spectra are contaminated with particle-induced, non-X-ray background as well as with local astrophysical X-ray background (see Freyberg et al. 2020;Yeung et al. 2023).An additional spectrum was extracted in a region near but well outside the CNC to estimate the local X-ray background.The background was first modelled and fit, then the source spectrum (for regions in the CNC) was simultaneously fit with the background spectrum, also including the background components.The particle background can be described with a multi-component powerlaw spectrum with emission lines arising in the instruments.In order to esti-mate the particle background we used the data up to 9.0 keV, even though no diffuse emission was observed above ∼3.0keV in most cases (above ∼5.0 keV in the C1 region).The astrophysical X-ray background consists of emission from the Local Hot Bubble, Galactic disc and halo, and the extragalactic X-ray background.The spectral model parameter values, which were obtained from the fit of the spectrum of the background region, were also used for the fit of the source spectrum, from which the background had not been subtracted.The background model components were scaled with a multiplicative constant parameter, which should mainly account for the different size of the extraction regions, but was free in the fit.For example, the area of the Central regions C1 and C2 is 0.019 deg 2 and 0.076 deg 2 , while the area of the local background region is 0.785 deg 2 , the theoretical scaling factor thus being 0.024 and 0.096 for regions C1 and C2, respectively.The best fit yielded a scaling factor of 0.02 (< 0.046) for C1 and 0.11 (0.06 -0.23) for C2 with the 90% confidence limits given in brackets.
X-ray observations of the ISM in the Milky Way and nearby galaxies have shown that the emission is well reproduced with at least two thermal plasma components with different temperatures (e.g.Kuntz & Snowden 2010;Kavanagh et al. 2020, and references therein).The lower temperature component has ∼0.2 keV and indicates emission from the hot phase of the ISM in equilibrium.The higher temperature component has temperatures > 0.5 keV and is particularly significant in regions with recent additional heating (H ii regions, superbubbles, and supernova remnants, SNRs).To model the plasma emission we used the collisional ionisation equilibrium (CIE) model vapec 1 and non-equilibrium ionisation model vnei (Borkowski et al. 2001, and references therein).
We have also fit the element abundances, in particular for O and Fe.In all regions, the fits of the eROSITA spectra suggest that the element abundances are consistent with solar values.At first glance this looks contradictory to the results of Townsley et al. (2011a) who found enhanced abundances of O, Ne, and Mg in a few smaller regions, and super-solar abundances of Si, S, and Fe in the CNC.However, since eRASS is a shallow survey, as our first tessellation approach has shown, the photon statistics are not sufficient to analyse small regions.Therefore, enhancements of O, Ne, or Mg seen in small regions in the deep Chandra data will be dominated by the emission from plasma with solar abundances in larger regions.Moreover, there is an Al-K instrumental line at 1.5 keV in the eROSITA data and the sensitivity drops significantly above 2 keV, which makes it impossible to constrain the abundances of Si, S, and Fe.
In the Central regions C1 and C2, and to some extent also in the Hook regions H1 and H2, the emission appears harder than in the other regions.Townsley et al. (2011a) have shown that this emission is well reproduced by one or two additional thermal components with higher temperatures, accounting for possible emission from, e.g., pre-main-sequence stars or η Car.We therefore fit the spectra in these regions with additional thermal components.An alternative scenario for the hard emission is that particles are accelerated in shocks of stellar winds and can give rise to a non-thermal component also in the X-ray regime.To verify the existence of such a component, we also included a powerlaw emission component (powerlaw) in some regions.Townsley et al. (2011a) also considered the existence of a non-thermal component in the diffuse emission of Carina nebula, however, did neither present nor discuss any spectral fit results including a powerlaw component.
To account for absorption, we used the model tbvarabs (Wilms et al. 2000).The spectra of the local background and the source regions with the best-fit models are shown in Fig. 3 and the spectral fit parameters are listed in Table 1.
The extraction region C1 covers the emission around Trumpler 14 and Collinder 232, which are seen as bright point-like sources and have been removed.The diffuse emission is hard, and the spectrum shows significant emission above 1 keV without emission line features.For the source spectrum, an additional emission component is necessary to reproduce the emission between 1 and 5 keV.We first added a third thermal plasma component for the source model.As an alternative, we also fit the spectrum including a powerlaw component instead of the third thermal component.Both the model with a third thermal 1 http://atomdb.org/component and with a non-thermal component yield similarly good fits.A model including two additional thermal components with kT = 2.5, 6.0 keV similar to the model of Townsley et al. (2011a) for the hard emission also yields comparable fit statistics.If we assume one additional thermal component, the temperature becomes very high (kT = 6.8 (4.7 -12) keV, with 90% confidence range in brackets), rather unphysical for shocked interstellar plasma.The spectrum is fit well with a combination of a two-component thermal plasma model assuming non-equilibrium ionisation (NEI) with kT = 0.24 keV and 0.82 keV and a powerlaw component with photon index Γ = 1.7 (1.3 -2.2), which can be explained by a population of non-thermal electrons.A cavity is seen in the HI map (Fig. 2d).Most likely, the cold material around the stellar clusters has been ionised and blown away by the radiation and the stellar winds of the massive stars in the clusters.
In the region C2, there is diffuse emission around Trumpler 16.The bright X-ray emission from the massive star binary η Car and from the star WR25 along with other detected point sources inside the region has been removed.The emission appears yellow-ish in the three-colour-image (Fig. 2a), with dominant emission in the medium band, but significant contributions from the soft and the hard bands.The best-fit model requires two thermal NEI plasma components (kT = 0.28 keV and 0.82 keV with n e t ≈ 10 11 cm −3 s for both) and an additional harder component.In this region, the fit with a non-thermal component is preferred to a fit with two additional thermal components as in Townsley et al. (2011a) with χ 2 = 513 and degrees of freedom (dof) = 468 compared to χ 2 = 561 and dof = 466 for a pure thermal model.If we again assume one additional thermal component, the fit yields a very high temperature of kT = 15 keV with a lower limit of 9 keV (unconstrained to higher values).In the best-fit model, the non-thermal component has a photon index of Γ = 1.8 (1.5 -2.3).This component again suggests the existence of non-thermal electrons.
Next to the region C2, we have defined the X-ray hook region H1, which is located south of the triangular shaped bright region in the optical and IR, which also in part coincides with our region C2 (Fig. 2b and c).Relative to the massive star clusters in the CNC, the X-ray hook is located outside the bright photoionised nebula, indicative of a shock of the expanding H ii region.The emission appears yellow-green, with dominant emission between 0.5 -1.0 keV.The spectrum is well reproduced with two thermal components with kT = 0.26 keV and 0.63 keV, with the first being consistent with NEI (n e t= 2.7 ×10 11 cm −3 s), while the hotter component seems to have reached collisional ionisation equilibrium (CIE).In this region, the contribution of the hotter component is as strong as that of the cooler component (norm 2 /norm 1 = 1), while in the other regions, it is less than half (norm 2 /norm 1 < 0.42) or almost negligible.Both the high ionisation timescale and the relatively high brightness (and hence high norm 2 ) for the hotter component suggest that the gas might be highly compressed.
The tip of the X-ray hook (region H2) appears greener than the region H1.While the region H1 traces the edge of a cavity seen in HI (Fig. 2d), the region H2 coincides with a larger HI structure seen in the north.The spectrum of the region H2 suggests a similar intrinsic spectrum as in the region H1 (two thermal components with two temperatures), while the hotter component also seems to be in NEI in the region H2.The foreground absorption is higher in this region (N H = 2.5 (2.0 -3.0) × 10 21 cm −2 , suggesting that the tip of the X-ray hook is located inside or behind the HI structure.In both Hook regions, there is no indication for a non-thermal component.The other regions, in which there is diffuse X-ray emission with sufficient surface brightness that allows us to perform spectral analysis reasonably well, the emission appears red to orange, with most of the emission below 1 keV (regions OS1, OS2, OS3 as well as O1 and O2).These regions coincide with lowerdensity regions in the HI map.Most likely it is emission from hot plasma, escaping from the central active regions.The spectra are well reproduced with two thermal components for plasma in CIE with a lower-temperature component with kT ≈ 0.2 keV and a higher-temperature component with kT > 0.8 keV.The region O2 extends north of and along a broad optical filament, likely photoionised gas delineating the end of the cavity.The hot gas streaming through the HI cavity will likely be adiabatically compressed (hence heated) because it is being pushed against the HI cavity wall.Mixing with HI gas then creates the more luminous, colder component.The fact that the region O2 is spatially anti-correlated with the optical region to the south suggests that the cavity elongates southwards approximately in the plane of the sky.The region O1 is also associated with optical emission, but in this case, X-rays and optical emission are more co-spatial, which suggests that the associated HI cavity is significantly inclined to the plane of the sky.

Stellar sources
The Carina nebula complex contains several X-ray bright point sources associated with massive stars and stellar systems, which have been studied extensively over the last decades, see, e.g., the Chandra based CCCP (Townsley et al. 2011c;Gagné et al. 2011;Nazé et al. 2011).While the CNC harbours thousands of young stars, only the X-ray brightest ones are detected as individual sources in eRASS (see Table A.1).In our analysis of the eRASS data we thus focus on the prominent X-ray sources in the region, i.e., η Car itself, several O-type stars and binaries, and the Wolf-Rayet stars WR 22, WR 24 and WR 25.An intensity scaled eRASS image with the studied stars and used source regions is shown in Fig. 4.
For the point source analysis, events from all TMs were used.Source and background photons were extracted from respectively adapted circular and annular regions.While some minor contributions from other nearby stars and diffuse emission are present, for the eRASS sources discussed here, the photon distribution is in most cases strongly dominated by the X-ray emission from a single star or a stellar system.With the exception of η Car we use the eRASS:4 combined spectra in our analysis and thus derive the average properties of the sources.The Xray spectra are fit with absorbed multi-temperature CIE plasma models (apec) with errors denoting the 90% confidence range.A foreground column density of 2 × 10 21 cm −2 is used to account for the interstellar foreground absorption and the adopted distance to the CNC is 2.3 kpc.
During the eRASS mapping, the region was observed four times with a long-term sampling of about 0.5 yr.The exposure times per individual eRASS, obtained in multiple scans over 2 -3 days, are around 0.5 ks and provide sufficient photons to study the X-ray variability for our targets.Example light curves of three selected X-ray sources are shown in Fig. 5.

η Carinae
η Carinae (HD 93308), the most prominent stellar object in the CNC, is a very massive stellar binary in a highly eccentric 5.5 yr orbit located in Trumpler 16.It consists of an eruptive primary in its luminous blue variable (LBV) phase and a secondary, potentially evolved high-mass star with mass estimates of M 1 ≳ 80−100 and M 2 ≳ 30−40 M ⊙ .A discussion of the η Car binary parameters can be found, e.g., in Corcoran et al. (2001); Kashi & Soker (2016).In this system, the collision of the fast wind from the secondary and the slower but more massive wind from η Car A leads to the complex and variable X-ray emission and absorption scenarios over the orbit, especially around periastron (see, e.g., Pittard & Corcoran 2002;Hamaguchi et al. 2007a;Corcoran et al. 2017).
X-rays from η Car include emission from the outer structures created by the ejecta of the last eruptions, e.g., the famous Homunculus nebula and a more absorbed and variable hard component from the inner wind shock regions (Corcoran et al. 2022).The source is not spatially resolved by eROSITA and the data provide a combined source spectrum.Due to the optical/IR-brightness of the star (G = 4.0 mag, J = 1.7 mag), a mild contamination from optical loading is present at energies of ≲ 0.3 keV, and is included as a background component in the modelling.The combined eRASS data are first used to derive a 'baseline model' that describes the non-variable spectral components in our η Car data.We then study the variability between the individual eRASS visits.Particularly in eRASS1, the data taken on January 02 -05, 2020, stand out, as this observation covers the X-ray bright state before the periastron passage.The obtained spectra are shown in Fig. 6.
The X-ray emission from the stable component, associated to the outer ejecta/nebula around η Car, is well described by a twotemperature plasma with components at kT 1,2 = 0.15 keV and 0.7 keV, with the ratio of the emission measures EM1/EM2 ≈ 1, a high N overabundance by a factor of 15 and an absorption column fixed at the ISM value of N H = 2 × 10 21 cm −2 .From the observed F X = 1.8 × 10 −12 erg cm −2 s −1 we derive an emitted 0.2 -10.0 keV X-ray luminosity of 3.7 × 10 33 erg s −1 for the emission from the surrounding region.The outer ejecta is well described by a constant emission component in our data, albeit a slow temporal evolution is expected for the expanding nebula.Comparing our eRASS results with older X-ray measurements, we find that the derived flux from the combined eRASS 2020/21 data matches the long-term trend seen in X-rays, which indicates a t −3 flux decline, as deduced from a fit to XMM-Newton EPICpn data from observations performed between 2003 and 2015 (Corcoran et al. 2022).For the 0.3 -1.0 keV energy band we get F X = 1.3 ± 0.1 × 10 −12 erg cm −2 s −1 from our spectral model, well matching the fitted X-ray trend over two decades and supporting their Great Eruption and blast wave evolution scenario.
Further, a stronger absorbed, variable plasma component is detected.It is thought to be associated with the wind-windcollision zone and one or two hot plasma components are required to model the eRASS spectra.Being most prominent during the bright state, Fe Kα 6.4 keV emission is detected, e.g., with a line flux of 4.7 × 10 −4 photons s −1 in eRASS1.To put our data into context, we use the orbital X-ray phasing, MJD_min = 50799.4+ 2023.4 × n, following Corcoran et al. (2017).
The object η Carinae has been monitored in X-rays over several decades.The flux variability and spectral changes observed during the eRASS are similar to those observed over previous orbits (Corcoran et al. 2017) and overall agree with findings from NICER observations around the 2020 periastron passage (Espinoza-Galeas et al. 2022).However, the eRASS data give about 30 % lower fluxes in the 2.0 -10.0 keV band when compared to the NICER values.Although the data were not taken strictly simultaneously, this indicates cross-calibration issues potentially caused by effective area uncertainties at energies above ≈ 2 keV for the eROSITA survey data.A comparison of the derived fluxes from the two instruments is shown in Fig. 7. Given the errors, the eRASS data are well described by a two-temperature model with kT 1 = 0.8 − 1.0 keV and kT 2 = 4.0 − 4.5 keV.Notably, the eRASS2 data at phase 0.068 are better described by a temperature of roughly 3 keV and similarly the NICER spectral models give 3.0 -3.6 keV for data taken within a few days.However, it is also during a phase of quite strong absorption, while kT and N H are interdependent parameters in spectral modelling.
The eROSITA spectrum of η Car is predominantly of thermal origin as shown in our model fits and evidenced by the strong Fe xxv 6.7 keV emission line.Surely, a fainter non-thermal contribution can always be 'hidden', especially in the spectra with a moderate signal-to-noise ratio, as its basic spectral shape for reasonable values of the photon index is similar to those of a hot thermal plasma, but it is not required to adequately describe our data nor favored in the fits.On the other hand, non-thermal emission from η Car is detected at hard X-rays and soft Gammaray energies (Hamaguchi et al. 2018).However, the expected contribution is only at a percent level in the eROSITA energy range, consistent with our results and conclusions from the spectral analysis.

Wolf-Rayet stars
The CNC harbours three additional X-ray bright Wolf-Rayet stars, WR 22, WR 24, and WR 25, which are also eRASS sources.We discuss the data and summarise the spectral modelling results in Tab. 3.
WR 25 (HD 93162): This star is a massive colliding wind binary of O2.5If*/WN6+O type.It has component masses of about 75+30 M ⊙ , an orbital period of 208 d period and has been studied in pointed X-ray observations, see e.g.Raassen et al. (2003); Gamen et al. (2006); Pollock & Corcoran (2006).In our survey data, WR 25 is the X-ray brightest among the Wolf-Rayet stars and after η Car the second brightest X-ray source in the studied Carina field.During eRASS, WR 25 was observed at orbital phases of 0.90, 0.77, 0.66, and 0.57, and its X-ray brightness varies by a factor of about two as shown in Fig. 5.We measure an observed F X = 3.3 − 6.9 × 10 −12 erg cm −2 s −1 , with the highest flux detected in the pre-periastron phase.However, changes in the spectral properties are minor and beside a re-normalisation, a uniform spectral model is suitable to describe the data.We find an average log L X = 33.6 and plasma components of comparable strength at temperatures of kT ≈ 0.8 and 2.0 keV.

WR 22 (HD 92740):
The Wolf-Rayet star WR 22 is a WN7h+09III-V system with 70+25 M ⊙ and P orb = 80 d.A discussion of the WR 22 system and a phase-resolved X-ray study with XMM-Newton is presented in Gosset et al. (2009).Over the eRASS, X-ray flux variations due to the orbital motion of about a factor of two are seen, see Fig. 5.The phase dependence (eRASS1-4 phases are 0.59, 0.86, 0.15, 0.50) of the observed rate is comparable to the one observed with XMM-Newton.The spectra of WR 22 show a dominant 0.6 -0.7 keV plasma component.In addition, we find indications for contributions from hotter plasma with kT ≳ 2.0 keV as in WR 25, but this component is poorly constrained by our data.
WR 24 (HD 93131): WR 24 is a Wolf-Rayet star of spectral type WN6ha-w and the X-ray faintest among the studied WR stars.Its plasma temperature is comparable to those of WR 22 and similarly indications for a hotter plasma component are present.These properties are quite typical for a colliding wind binary, however so far no companion has been identified.While strong intrinsic absorption is seen in all WR stars, it is highest in WR 24.

X-ray bright O stars
Other X-ray bright stellar sources are detected among the Otype stars, which are typically the more massive members with very strong stellar winds in a binary or a multiple system.An overview of the population of massive stars and a description of individual sources is given in Gagné et al. (2011), from where we adopt the stellar parameters if not denoted otherwise.We model and discuss the eRASS data from the most outstanding X-ray bright stellar sources in the following.
HD 93129 (Trumpler 14-1): The O-type star HD 93129 is among the earliest-type stars in the Milky Way and a bright X-ray source.It was classified as an binary of the type O2If*+O3.5V((f))5;however, HD 93129 A was discovered to be a binary itself consisting of an O2 supergiant and an O3.5 star, making HD 93129 at least a triple system.There is another massive binary HD 93128 (Trumpler 14-2) nearby, which is a spectroscopic binary of the spectral type O3.5V((fc)).Both stars are separated by 24 arcsec and not fully resolved in the eRASS data, with HD 93129 being the dominant X-ray source by a factor of about four.Strong X-ray emission from hot plasma at kT of a few keV is detected from this source.
HD 93160 and HD 93161: These are O-type multiples of the spectral type O7III((f)) SB and O8V+O9V+O6.5(f)located in Trumpler 14.They are separated by 15 arcsec and blended in eRASS, the flux ratio between the two sources is about one with HD 93161 being moderately brighter by 20%.
HD 93205 (V560 Car) is an O-type binary of the spectral type O3.5V((f))+O8V.Its X-ray spectrum is among the softest of the studied sources and we find that plasma at temperatures up to 1 keV dominates its X-ray emission.
HD 93250 is an O4III(fc) star with suspected binarity, but a secondary was not unambiguously identified so far.The eRASS data confirms it as a bright X-ray source at L X = 10 33 erg s −1 with a strong hot plasma component at kT of a few keV.
HD 93403 is an O5.5III(fc)+O7V eccentric binary.It is a bright X-ray source with L X = 1.3 × 10 33 erg s −1 and has a hard X-ray spectrum that shows strong emission from hot plasma at kT of a few keV.
HD 93206 (QZ Car) located in Collinder 228 is a O9.7I+O8III double binary, i.e., at least a quadruple system.Bright X-ray emission with a hot plasma component is observed from this massive O-type multiple system.
Coll 228-103 is an OB star in the southern cluster Collinder 228, which is part of an OB mini-cluster consisting of Coll 228 101-103 and HD 305453.
In the studied O-type stars the X-ray emission typically exhibits a strong and often dominant plasma components at temperatures of about 0.6 -0.8 keV.In addition, in several targets, significant contributions from hotter plasma at a few keV is present, most prominently in HD 93129, HD 93206, HD 93250, and HD 93403, which all have kT av ≳ 1 keV.These high temperatures are specifically found in colliding wind binaries or other multiple systems.In contrast, a few targets have predominantly cooler plasma in addition to the 0.6 -0.8 keV component (HD 93205, and exhibit kT av ≲ 0.5 keV in our models.Comparing the four eRASS visits for the O-type stars, we find only mild variability within about 20 % in the observed X-ray fluxes.All spectral modelling results of the studied sources are summarised in Tab. 3.

Discussion
The CNC-Gum 31 region is a very bright H ii region and a very active star-forming complex, but not as extreme as 30 Doradus (Dor) in the Large Magellanic Cloud (Massey & Hunter 1998) or the Arches cluster near the Galactic Centre (Najarro et al. 2004).However, the close location of the CNC-Gum 31 complex and its properties make it an excellent target for the study of high-mass star-forming regions and therefore also for the understanding of the star formation in the early Universe.

Multi-wavelength comparison
In a new H i study at the Australia Telescope Compact Array, complex filamentary structures were found in the distribution of the cold neutral gas over a wide range of velocities (Rebolledo et al. 2017).Several 'bubbles' were detected in the CNC, which have obviously been produced by the impact of the massive stars.
High-resolution observations of the 12 CO(1-0) and 13 CO(1-0) molecular lines in the Carina nebula and the Gum 31 region obtained with the 22 m Mopra telescope as part of The Mopra Southern Galactic Plane CO Survey was presented by Rebolledo et al. (2016).From the sample of the massive clumps observed with Mopra, two regions with very different properties have been observed at the Atacama Large Millimeter/Submillimeter Array (ALMA) by Rebolledo et al. (2020).One region is located at the centre of the nebula and has most likely experienced some feedback from the massive star clusters, while the other region was selected further south and is believed to be less disturbed by the massive stars.It was shown that the centre of the nebula is forming fewer cores than in the south, however, the cores are much more massive.These differences suggest that the stellar feedback has an effect on the formation of clumps which evolve into starforming cores.In addition, Fujita et al. (2021) have analysed data taken with the Mopra and NANTEN2 telescopes as well as with ALMA and have detected signatures of cloud-cloud collisions in the velocities and the morphologies at different velocities of the clouds.They suggest that the collisions between clouds have triggered massive star formation in CNC.
Moreover, the radio continuum map in the 1 -3 GHz range (Rebolledo et al. 2021) shows a complex distribution of emission across the nebula consisting of filaments, shells, and fronts on different size scales.Ionization fronts are found at distances of 80 pc from the stellar clusters in the centre of the CNC to the north and the south.It has been estimated that approximately 15% of the ionising flux has escaped from the nebula into the diffuse Galactic interstellar medium.
From the comparison of the radio continuum emission to the diffuse X-ray emission we can conclude that the cavities seen in the H i map in the north and the south have most likely been formed by the influence of massive stars in the star clusters at the southern end of the northern cloud, also confirmed by the bright diffuse X-ray emission in the Central region C2.As pointed out by Rebolledo et al. (2021), the radio continuum and IR maps show plumes at the edges of the H ii region, also around the northern tip of the X-ray hook.These plumes are indicative of cavities filled with over-pressurised hot gas, which broke out and is expanding out of the molecular cloud into the Galactic plane.The diffuse X-ray emission in the OS2 region apparently is due to the expanding hot plasma in the southern cavity.

Origin of the X-ray emission
The results of the spectral analysis indicate that the emission is mainly caused by thermal plasma in all regions.The element abundances are consistent with solar values, which indicates that the origin of the hot plasma is shocked ISM.Townsley et al. (2011a) found enhanced Fe abundance in their 'central arc' region, which corresponds to our Hook region H1.In young stellar clusters like those in the Carina nebula it is likely that there has already been a supernova explosion and that an SNR contributes to the X-ray emission.The combination of stellar winds and SNRs will create a superbubble around the stellar clusters.Therefore, in the context of the multi-wavelength study, we compare the observed emission with superbubble simulations (Krause et al. 2013(Krause et al. , 2014)).
The bright X-ray hook region is most likely caused by a shock in a mixing region at the edge of a superbubble driven mainly by η Car.The time variability of the stellar driving power of superbubbles leads to various instabilities of superbubble shells, which tend to produce more extended turbulent mixing regions compared to the gradual change in density and temperature from the hot bubble interior to the dense shell, which would be expected from thermal conduction alone.Dissipation in such mixing regions dominates the radiative losses in superbubbles (compare Fierlinger et al. 2016).Contrary to early 1D models with steady source power (e.g., Chu & Mac Low 1990) the losses can be strong enough, hence the internal energy density in the superbubble low enough, such that a significant increase in stellar source power can drive a shock wave through the entire superbubble.In the simulations of Krause et al. (2014), supernova shock waves encountered the optimal conditions in such mixing regions towards the inner edge of supershells to produce an order of magnitude enhancement in soft X-rays.In the present case, it seems likely that η Car, currently in the LBV phase, drives such a shock.An LBV phase is expected to last about 10 5 years Table 3. Spectral fit results for eRASS:4 data.F X is the observed 0.2 -5.0 keV flux, L X is the absorption corrected 0.2 -5.0 keV X-ray luminosity.
name Notes.N H is the total column density, with the minimum value set to the mean ISM value.A two-temperature fit was used when it improved the fit significantly.(Groh et al. 2014) with strong mass loss and powerful eruptive, possibly explosive winds (Gormaz-Matamala et al. 2022;Vink 2022).Massive star winds have typical velocities of the order of 1000 km s −1 .The distance of ∼10 pc between η Car and the X-ray hook are hence traversed in about 10 4 yr.It is therefore likely that the X-ray enhancement in the Hook region is caused by η Car's current LBV phase.
The Central regions around η Car are best explained by the classical inner shock of its wind (Weaver et al. 1977).Due to strong adiabatic expansion, the temperature of the gas before the shock is low, hence the shock's Mach number is high.It is therefore highly likely to find diffusive shock acceleration at such shocks (Vink & Yamazaki 2014).Consequently, we find emission from non-thermal electrons in the regions around the massive stars only, including η Car.
The radio continuum image (Rebolledo et al. 2021) suggests a superbubble around η Car, which is limited towards the southeast by the Hook region and extends significantly further towards the north-west.The latter is confirmed by the HI map in Fig. 2. Our Outer region O1 extends across a significant part of this superbubble, as defined by the radio continuum image, and is cospatial with the radio continuum as well as the R-band emission.Both, radio continuum and R-band image show photoionised gas, which is found towards the inner edge of the superbubble shells.
The situation is different for our Outer region O2.There, the X-ray-bright region as well as the photoionised gas region is much narrower in the direction of the major axis of the corresponding superbubble (i.e., the one south of the Hook region).The X-rays here are anti-correlated with the photoionised gas, with X-rays being further inside the superbubble.
Both our outer regions can be explained by hot gas being pushed against the shell of their respective superbubbles.The compression, and possibly some mixing, enhances the emission measure, and the regions are thus more prominent in X-rays.Southwards the superbubble extends close to the plane of the sky.It is the comparatively denser part of the shell that is photoionised.Hence, the X-rays appear "inside" the shell that is delineated by the photoionised gas, similar to the case of the Orion-Eridanus superbubble (Krause et al. 2014).The region O1 is seen almost "from the top" (or bottom), i.e., the superbubble extends at a significant inclination to the plane of the sky.Hence, we see the X-ray-bright and the radio continuum-bright regions on top of each other.
The fact that we see both Outer regions O1 and O2 at similar X-ray brightness at the same time is probably accidental.As shown in Krause et al. (2014), the temporal variability of the stellar sources inside superbubbles leads to gas sloshing, i.e., the hot gas moves around in superbubbles.The expected timescale of this variability is about 10 5 yr.Wherever the hot gas hits the dense shell, it becomes brighter due to the adiabatic compression.This interpretation is consistent with the many arcs in the radio continuum image that are not associated with X-ray-bright regions.It is also possible that the Outer region O1 is indicative of a temporal compression in the interior of the superbubble related to the gas sloshing.
The different bubbles around the individual massive star groups in the CNC are highly likely to have very different pressures.This follows from the strong dependency of the wind power on mass and evolutionary state of a star.The η Car bubble is clearly the highest-pressure system, which is not only evident from it having the most powerful stellar source, η Car, but also from the shock it drives, visible as the Hook region.
Due to the overpressure in the η Car bubble, a hot wind is driven through the molecular clouds in the centre of the system.The wind has swept away the more tenuous parts of the molecular cloud with which it is mixed southwards in our Outskirts regions.The mixing decreases the X-ray temperature and increases the emission measure, as observed.It is also interesting to compare the HI map (Fig. 2) to the molecular and radio continuum data (e.g., Rebolledo et al. 2021, their Fig.9).While the northwest superbubble and the southern superbubble seem connected in the HI map, molecular as well as ionised gas seems to be left in the interface between the two superbubbles.This could mean that the intermediate-density gas (HI) is blown away and mixed into the X-ray gas, whereas the denser, molecular parts remain.Some of the dense gas is photoionised and keeps being ablated.
The denser gas in the molecular clouds, which is not swept away, is likely to be surrounded by the hot X-ray gas and compressed by its high pressure (compare Krause et al. 2018).This may well be the mechanism that formed the massive cores detected by ALMA in the central part of the nebula (compare above).This process has been suggested to explain the sequential star formation in the Scorpius-Centaurus superbubble by Krause et al. (2018).Their simulations directly support the combination of sweeping away of more tenuous gas combined with the formation of dense gas cores.The plumes identified in the radio continuum image are also consistent with this picture, probably representing somewhat denser gas phases that are pushed out of the molecular region by the hot wind.

Thermal properties of the superbubble
If the high-mass stars in the central region blow strong stellar winds and encounter an obstacle in the X-ray hook, these regions should have a higher pressure than the surroundings.We use the results from the X-ray spectral analysis to estimate the density and pressure in these regions.
We assume that in the central region, the stars blow a bubble with a size similar to the Central regions C1 and C2.We approximate the bubble as a sphere with a diameter of 0.37 • .The Hook region (H1 and H2) is most likely a shell-like structure around the central region, which we see from the side.We assume that it can be modelled as a part of a spherical shell with an outer diameter of 1 • and an inner diameter of 0.65 • .We will also study the Outskirt region OS2, which is plasma flowing into the northern part of the southern cavity.We model the structure as a half ellipsoid with a major axis a of 0.64 • and a minor axis b of 0.5 • .For the axis parallel to our line of sight c, we assume the same value as for the minor axis b.
Based on these assumptions, we calculate the following volumes for the X-ray emission: for the Central regions (C1 and C2) V c = 4.2 ×10 58 cm 3 , for Hook regions (H1 and H2) V h = 2.9 ×10 59 cm 3 , and for the region OS2 V os = 5.8 ×10 59 cm 3 .
To calculate the density, we use the normalisation of the spectral models for thermal plasma: with D being the distance to the CNC and assuming homogeneous density in each of the regions and n e = 1.2nH .The index i = c, h, and os denotes the parameters for the Central regions (combined), Hook regions (combined), or the OS2 region.Therefore, for the hydrogen density and the pressure of the plasma, we get and Using Eq. 2 and 3 we get for the Central regions n H,c = 4.0 ×10 −2 cm −3 and P c /k = 9.1 ×10 5 cm −3 K, for the Hook regions n H,h = 5.4 ×10 −2 cm −3 and P h /k = 9.4 ×10 5 cm −3 K, and for the OS2 region n H,os = 2.3 ×10 −2 cm −3 and P os /k = 5.1 ×10 5 cm −3 K.The uncertainties of the temperature and the normalisation are about 5% and 10%, respectively, while our volume estimate is most likely uncertain by 20 -30 %.Therefore, the derived values have an uncertainty of ∼50%.The numbers still indicate that the density and the pressure are higher in the Central and Hook regions, while the lower pressure in the southern cavity indicates outflowing gas.The pressure in the CNC is one to two orders of magnitude higher than the average pressure in the ISM in the Milky Way (Jenkins & Tripp 2021, and references therein).

Cosmic ray pressure
In the shocks of the strong stellar winds of massive stars, electrons and protons are accelerated to relativistic energies and become cosmic rays.In young stellar clusters like the Carina Nebula, there also regions in which stellar winds collide and which are thus very efficient acceleration regions (see, e.g., Steinmassl et al. 2023, and references therein).The harder emission in the spectrum in the Central regions can be either explained by synchrotron radiation or inverse Compton scattering of the accelerated particles.
For the synchrotron hypothesis, we need to assume a magnetic field strength.As the thermal pressure is about two orders of magnitude higher than the Galactic average, we might expect a similarly enhanced magnetic pressure.Since the magnetic field in the Milky Way is generally of the order of µG (Sun & Reich 2010;Xu & Han 2019;Dickey et al. 2022), we may expect a magnetic field in the Carina superbubble system of several 10 µG.The synchrotron emissivity for a power-law distribution of electrons is given by (Longair 2011): Here, m e and e are, respectively, electron mass and charge, c is the speed of light and ϵ 0 the vacuum permittivity.The Gamma function is denoted by Γ f and the electron energy distribution is n(E) = κE −q .The electron distribution power-law index q is related to the photon index Γ by q = 2Γ − 1.For an assumed magnetic field strength B, the characteristic frequency for synchrotron emission is given by The non-thermal luminosity observed in the X-ray band is given by the integral: J(ν) dν , and we take ν 1 = 4.5 keV/h and ν 2 = 12 keV/h to be clearly dominated by the non-thermal part of the spectrum.This procedure allows us to infer κ for an assumed magnetic field by equating the predicted luminosity to the non-thermal luminosity we observe in regions Central 1 and 2, L obs,nt,4.5−12keV = 5.59 × 10 33 erg s −1 , which we use to infer the relativistic electron pressure by integrating over the energy distribution, P cr,e = (γ ad − 1) where the adiabatic index is γ ad = 4/3 for relativistic particles.We follow del Palacio et al. (2022) in adopting a leptonic contribution to the overall cosmic ray pressure of 2 %.The cosmic ray pressure as well as the ratio between the magnetic and the cosmic ray pressure are shown over the assumed magnetic field strength in Fig. 8.The thermal pressure is also shown for comparison.Under the synchrotron hypothesis, our measurements allow for two possible solutions: if the magnetic field strength is low, the region must be dominated by the cosmic ray pressure.If the magnetic field is high, the magnetic pressure will dominate.The measurements exclude the possibility that the thermal pressure dominates.A solution with all Fig. 8. Possible cosmic ray pressure in regions C1 and C2 under the assumption that the measured non-thermal X-ray radiation is synchrotron emission.The blue line shows the total cosmic ray pressure (hadrons + leptons, left vertical axis) and respectively the ratio between magnetic and total cosmic ray pressure (right vertical axis) over magnetic field strength.The magnetic field values used are not observationally known, but are reasonable values given other observables of the region.The orange line shows the thermal pressure of the same region for comparison (left axis).See Sect.5.4 for a discussion of this plot.three pressure components in approximate equipartition exists for B ≈ 60 µG.We note that for these assumptions, we predict a non-thermal radio flux at 2 GHz of 7 Jy, well consistent with the thermally dominated radio continuum emission of 2 kJy observed by (Rebolledo et al. 2021).There are generally very few detections of non-thermal radio components in superbubbles, and it is possible that special geometries or circumstances may be required to distinguish the non-thermal contributions (compare Heesen et al. 2015).
Strong magnetic field amplification is expected for shocks that produce the high-energy leptons (up to about 30 TeV) required in the synchrotron hypothesis (Schure et al. 2012;Vink 2012).The physics could be similar to the situation in the 30 Dor superbubble (Kavanagh et al. 2019), where TeV emission of the high-energy lepton population may be observed, though in the present case the accelerating shock would be the terminal wind shocks of the massive stars in the Central regions.Gammarays up to at least 400 GeV are also observed from the direction of η Car (Abdo et al. 2010;Martí-Devesa & Reimer 2021;H. E. S. S. Collaboration et al. 2020;Steinmassl et al. 2023).Given the extent of the point spread function in the HESS detection, some of this emission might be related to any inverse Compton component of the hypothetical TeV electrons, though the correlation of the GeV emission with η Car's orbital period observed with Fermi (Martí-Devesa & Reimer 2021) and the detection of the colliding wind binary in hard X-rays with NuS-TAR associate the bulk of this high-energy emission directly with η Car (Hamaguchi et al. 2018).The GeV γ-ray emission observed with Fermi (Abdo et al. 2010;Martí-Devesa & Reimer 2021) could also have a contribution from leptonic emission (H.E. S. S. Collaboration et al. 2020), but seems at least partly to be due to pion decay caused by the interaction of a hadronic cosmic ray population with the dense gas in the region (Ge et al. 2022).
The possible non-thermal X-ray spectrum we detect could be caused by inverse Compton upscattering of starlight or photons from the cosmic microwave background (CMB) by relativistic electrons with much lower Lorentz factors.The frequency shift is approximately given by (Longair 2011): Hence, upscattering from optical / UV (wavelength λ ≈ 400 nm) to X-rays requires electrons with Lorentz factors of γ ≈ 30.
From the direct measurements of the local cosmic ray electron population at MeV energies, we know that the energy spectrum becomes very flat, q ≈ 1.4 (after accounting for the solar modulation, Orlando 2018; Vittino et al. 2019;Krause & Hardcastle 2021).If the situation was similar in the CNC, one would expect a non-thermal photon index of around 1.2 rather than 1.8, which we observe.We estimate the cosmic ray pressure for this case in the following way: for a radiation energy density U rad the inverse Compton luminosity is given by (Longair 2011): with Thompson cross section σ T and the number of upscattering electrons N e .The latter can then be obtained by equating L IC to the observed luminosity.The total cosmic ray pressure for the same assumptions as above is then given by: P CR = (γ ad − 1) 51N e γm e c 2 V (10) Smith (2006b) gives a luminosity in ionising photons of Q = 10 51 s −1 , corresponding to a photon density of a few photons per cm 3 or U rad ≈ 10 −11 erg cm −3 .This would yield a cosmic ray pressure of P CR /k ≈ 10 9 cm −3 K, much higher than what we would expect in this region (compare above).
The CMB peaks at a frequency of ν = 160 GHz and has an energy density of U rad = 4 × 10 13 erg cm −3 .A similar analysis as before yields an estimate of the total cosmic ray pressure of P CR /k ≈ 10 9 cm −3 K, similar to the above estimate and again much higher than the thermal pressure, which we derive above.
In conclusion, if the non-thermal spectrum in the Central regions was due to inverse Compton scattering, the pressure in the region would be strongly dominated by cosmic rays, even if only the leptonic component was taken into account.These considerations seem to favour the synchrotron hypothesis.More sensitive TeV observations might be able to confirm this or otherwise.
As a further consistency check, we can work out the radius of the termination (inner) shock of the combined wind from the η Car binary.We adopt the values collected in Groh et al. (2012) (see Table 4).With ρ being the wind density and v the velocity of the wind, the wind ram pressure equals the ambient pressure p = ρv 2 at the termination shock.We take the thermal pressure in the central region as proxy for the ambient pressure.With this assumption, we find for the radius of the termination shock: Given the evidence for the additional cosmic ray pressure, we expect a termination shock somewhat closer to η Car than 15 pc, consistent with our assumption above that it is this shock that accelerates the non-thermal particles to TeV energies in the Central regions, and the synchrotron hypothesis above.Given this magnitude of the termination shock radius, it is possible that the shock has traversed the entire distance between η Car and the denser gas in the Hook region.5.5.Comparison of X-ray luminosity to wind power X-ray emission arises in the shock-heated intermediate-density gas and due to mixing between the hotter, shocked ejecta and the denser and colder gas phases.The X-ray luminosity can increase sharply after a strong energy increase; up to a factor of several tens is seen after a supernova in the superbubble, decaying over about 10 5 yr, as shown by Krause et al. (2014).On average, the simulations predict that about 2 − 3 × 10 −4 of the kinetic power injected by massive stars will be emitted in the 0.2 -12 keV band, and a fraction of about 10 −6 in the 2 -4.5 keV band.
The current kinetic power output by massive stars in the CNC was estimated by Smith (2006b) to be L kin = 3.4 × 10 38 erg s −1 .We collect partly updated estimates from the literature in Table 4 for the stars we detect in X-rays, which, by comparison to Smith (2006b), probably contribute the bulk of the kinetic power.Some values have increased in more recent analyses.For our sample only, we already find a combined wind power of L kin = 3.6 × 10 38 erg s −1 .
We derive a thermal X-ray luminosity of the entire CNC of 1.31 × 10 35 erg s −1 (0.2 -12 keV), or about ≲ 3 × 10 −4 of the wind power.This is for our adopted assumption for the foreground absorption (without this, the values would be about 40% higher).This compares very well to the predictions of Krause et al. (2014).In the 2 -4.5 keV band, we measure a thermal Xray luminosity of 2.35×10 33 erg s −1 , which corresponds to a fraction of ≲ 7 × 10 −6 of the wind power.This is again in reasonable agreement with the prediction from the 3D superbubble simu-lations, and indicates that the thermal structure in the emerging Carina superbubble is similar to the one in the simulations.
A comparison between Table 4 and Fig. 4 confirms that the strongest stellar winds are from the η Car binary and the surrounding stars, and that the gas dynamics observed in X-rays in the CNC is thus indeed driven by this region.

Summary
We have analysed the diffuse X-ray emission in the CNC observed in four all-sky surveys with eROSITA (eRASS:4).The complete coverage of the CNC and its surroundings with eROSITA has allowed us to constrain the local X-ray background emission in a more improved way than in the studies before and to obtain reliable fit results of the X-ray spectra.We created images in X-rays, which were compared to multiwavelength images, in particular with HI data to identify cavities in the cold ISM.These cavities have most likely been created by the winds of the massive stars and are filled with hot X-ray emitting, low-density plasma.In all regions the dominating emission is modelled best as emission from thermal plasma.Near the massive stars, in particular, η Car, additional non-thermal emission component is detected in the X-ray spectrum.
The most prominent feature in the CNC in X-rays is the Vshaped X-ray hook, which is anti-correlated with the bright triangular shape of photoionised gas seen in the optical.As the spectral analysis of the eROSITA data have shown, in the X-ray bright Hook region H1 located south of the massive stars, the emission of the hot plasma component is stronger than in other regions.In this part of the nebula, the ionisation timescale is also high, suggesting high density, which is also consistent with the high brightness.The X-ray hook is most likely the edge of a superbubble formed by the shocks of the stellar winds of massive stars, mainly by η Car.The fainter, more diffuse emission in the outer regions is caused by the hot gas which fills the superbubbles around the massive stars in the CNC.The pressure inside the superbubbles, as derived from the spectral fit parameters, is one to two orders of magnitude higher than in the ISM.We suggest that the Outer regions result from mixing of the hot gas driven by an overpressure in the Central regions through the cold filamentary gas in the Hook region, with the ablated cold gas.
The spectrum in the centre of the CNC is harder than in the outer regions, indicating that there is either an additional hot thermal emission as found by Townsley et al. (2011a) or a non-thermal component.The latter is favoured in the fit of the eRASS:4 spectrum in the region around η Car and can be explained as synchrotron emission from electrons accelerated in enhanced magnetic fields inside the CNC.From the photon index and the flux obtained from X-ray spectrum we derived the energy distribution and the pressure of the electron gas and a magnetic field strength of B ≈ 60µG assuming equipartition.The theoretical size of the termination shock around η Car is very much consistent with the distance of the X-ray hook from η Car.The luminosity of the diffuse X-ray emission in the entire Carina nebula is about 10 35 erg s −1 (0.2 -12 keV), which is consistent with the luminosity of the diffuse emission obtained with Chandra (Townsley et al. 2011b) and expectations from superbubble simulations (Krause et al. 2014).
We have also studied the X-ray emission from the massive stars and binaries in the Carina Nebula, in particular the LBV binary η Car and the Wolf-Rayet stars WR 22, WR 24, and WR 25, and bright O stars.Their X-ray spectra are well modelled with one or two components of thermal plasma emission, with relatively high temperatures (kT > 0.6 keV).The data of the four surveys show variation of the X-ray flux with orbital phase for η Car, WR 22 and WR 25. η Car was observed shortly before periastron in eRASS1, with brighter emission above 1 keV.This harder variable emission can be explained as emission from the wind-wind-collision zone in the binary.

Fig. 3 .
Fig.3.eROSITA TM1-4,6 spectra are shown in the colours black, red, green, blue, cyan, respectively.For each TM we show the spectral model components with the corresponding colours (the cyan lines are best visible) with dotted lines.The model also includes the particle background (seen as a straight line) and local X-ray background components.In addition, we mark the source components with thick, red, dash-dotted lines.The extraction regions are indicated in the upper left corner of each panel.

Fig. 5 .
Fig. 5. Light curves of stellar sources in the CNC from eRASS, here highlighting the X-ray brighter targets η Car, WR 25, and WR 22, and their variability.The legend denotes the used energy band in keV of the respective light curve (e.g., 05_10: 0.5 -10 keV).

Fig. 7 .
Fig. 7. X-ray lightcurve of η Car showing the 2.0 -10.0 keV flux as observed with NICER and eROSITA.Note that the eRASS flux was scaled up by 30 % to broadly match the NICER results taken from Espinoza-Galeas et al. (2022).

Table 1 .
Spectral fit parameters.90% confidence interval is given in brackets.
* Size of the extraction region

Table 4 .
Properties of some massive-star winds in the CNC